ex8.2(2).inp

ansys有限元分析程序

! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!．ANSYS三维分析程序
FINISH
/CLEAR, NOSTART	! 清除数据，不加载初始化文件
! （1）定义工程选项、分析类型、单元类型、材料参数和实常数
/FILNAME, EX8.2(2)	! 定义工程文件名称为EX8.2(2)
/PREP7  	! 加载前处理模块
/TITLE, EX8.2(2), NONAXISYM. VIBR. OF A STRETCHED CIRCULAR MEMBRANE.
ANTYPE, STATIC	! 定义分析类型为静力分析
ET, 1, SHELL41, , 1	! 定义第1类单元薄膜单元SHELL41
MP, EX, 1, 207E9	! 定义第1类材料的弹性模量EX=207 GPa
MP, NUXY, 1, 0	! 定义第1类材料的泊松比NUXY=0
MP, ALPX, 1, 5.5556E-6	! 定义第1类材料的热膨胀系数ALPX=5.5556E-6
MP, DENS, 1, 7801.6	! 定义第1类材料的密度DENS = 7801.6 kg/m^3
R, 1, 1.27e-6	! 定义板厚t=1.27e-6 m
! （2）定义节点和单元
CSYS, 1	! 将当前坐标系切换为柱面坐标系
N, 1, .0001	! 在柱面坐标系中的（0.0001m, 0, 0）位置定义第1个节点
N, 10, 76.2e-3	! 在a=76.2mm位置定义节点10
FILL 	! 在1号节点和10号节点之间填充生成其余的8个节点
NGEN, 4, 10, 1, 10, 1, , 10	! 循环4次, 节点号增量10, 按照1到10号节点的格式, 所有Y坐标增加10度。
NROTAT, ALL	! 旋转所有节点的局部坐标系到当前的主坐标轴方向
E, 1, 2, 12, 11	! 过节点1, 2, 12和11定义第1个单元
EGEN, 3, 10, -1	! 循环3次, 每次单元的节点号增量为10, 按照前1个单元的模式生成单元
EGEN, 9, 1, -3	! 循环9次, 每次单元的节点号增量为1, 按照前3个单元的模式生成单元
! （3）选择节点, 施加位移边界条件
NSEL, S, LOC, Y, 0	! 选择所有Y＝0的直线边上的所有节点, 施加对称边界条件
! 定义用于循环对称分析的边界部件
CM, RIGHT, NODE	! 将这些选中的节点组合成一个部件, 叫作RIGHT
NSEL, S, LOC, Y, 30	! 选择另一条Y＝30度的直线边界上的节点, 施加对称边界条件
CM, LEFT, NODE	! 将这些选中的节点组合成一个部件, 叫作LEFT
NSEL, ALL	! 选择所有节点
/NOPR 	! 关闭输入数据的反馈和解释
CYCGEN	! 生成循环对称分析中的一个对称区
/GOPR	! 打开输入数据的反馈和解释
D, ALL, ALL, 0	! 固定所有节点, 为了便于施加预应力
TREF, 0	! 设定参考温度为0度
BFUNIF, TEMP, -12	! 将温度降低12度, 以生成预应力
FINISH 	! 退出前处理模块PREP7
! （4）进入求解模块, 打开预应力选项并开始求解静力分析
/SOLU	! 加载求解模块
PSTRES, ON	! 打开预应力选项
SOLVE	! 开始进行静力分析求解
FINISH	! 退出求解模块SOLU
! （5）再次进入前处理模块, 进行模态分析
/PREP7 	! 加载前处理模块
/TITLE, EX8-2(2), NONAXISYM. VIBR. OF A STRETCHED CIRCULAR MEMBRANE (MODAL)
ANTYPE, MODAL	! 指定分析类型为模态分析
MODOPT, LANB, 4, 1, 1000	! 使用LANCZOS循环求解方法, 展开前4阶模态, 频率范围1到1000 HZ
MXPAND, 4	! 设置展开模态个数为4
PSTRES, ON	! 打开预应力选项
! （6）修改位移边界条件
DDEL, ALL, ALL	! 放松所有节点的所有约束条件
NSEL, S, LOC, Y, 0	! 选择Y＝0的直线边上的所有节点
NSEL, A, LOC, Y, 30	! 添加选择Y＝30度的直线边上的所有节点
D, ALL, UX, , , , , UY	! 对以上选择的所有节点施加UX＝0和UY＝0的边界条件
NSEL, S, LOC, X, 76.2e-3	! 选择X＝76.2 mm的最外侧圆弧边上的所有节点
D, ALL, ALL, 0	! 固定被选择的最外边的这些节点的所有自由度
NSEL, ALL	! 选择所有节点
/REPLOT	! 重新绘制模型
FINISH	! 退出前处理模块PREP7
! （7）在求解模块, 进行轴对称问题的求解
/SOLU	! 加载求解模块
/OUTPUT, SCRATCH	! 定义输出文件的名称为SCRATCH
CYCSOL, 0, 1, 12, 'RIGHT'	! 波节直径上节点号从0到1, 将360度剖分为12份, 角度较小的直线边上节点是部件'RIGHT'
/OUTPUT	! 关闭文件输出
FINISH	! 退出求解模块SOLUTION
! （8）在一般后处理模块, 处理频率计算结果
/POST1	! 进入一般后处理模块
SET, 1, 1	! 读入第1个荷载步第1子步的计算结果
*GET, F01_1, ACTIVE, , SET, FREQ	! 提取0个节圆的第1阶频率
SET, 1, 2	! 读入第1个荷载步第2子步的计算结果
*GET, F02_1, ACTIVE, , SET, FREQ	! 提取0个节圆的第2阶频率
SET, 2, 1	! 读入第2个荷载步第1子步的计算结果
*GET, F11_1, ACTIVE, , SET, FREQ	! 提取1个节圆的第1阶频率
SET, 2, 2	! 读入第2个荷载步第2子步的计算结果
*GET, F11_2, ACTIVE, , SET, FREQ	! 提取1个节圆的第2阶频率
SET, 2, 3	! 读入第2个荷载步第3子步的计算结果
*GET, F12_1, ACTIVE, , SET, FREQ	! 提取1个节圆的第2阶频率
SET, 2, 4	! 读入第2个荷载步第4子步的计算结果
*GET, F12_2, ACTIVE, , SET, FREQ	! 提取1个节圆的第2阶重复的频率
*STATUS, PARM  	! 显示内存变量的值
FINISH	! 退出后处理模块
! （9）处理和显示模态
/NOPR	! 关闭输入数据反馈和解释
/POST1	! 进入一般后处理模块
EXPAND, 12	! 展开整个12个扇区模型的结果
/VIEW, , 1, 1, 1	! 从（1, 1, 1）到（0, 0, 0）的方向观察模型
/VUP, 1, Z	! 设置Z轴方向向上
/GLINE, , -1	! 不显示单元外轮廓线
/TRIAD, OFF	! 关闭坐标架显示
SET, 1, 2	! 读入第1个荷载步第2子步的计算结果
/TITLE, EX8_2(2), STRETCHED CIRCULAR MEMBRANE - NODAL DIAM 0, MODE 2
PLNSOL, U, Z	! 显示节点位移值UZ
SET, 2, 2	! 读入第1个荷载步第2子步的计算结果
/TITLE, EX8_2(2), STRETCHED CIRCULAR MEMBRANE - NODAL DIAM 1, MODE 2
PLNSOL, U, Z	! 显示节点位移值UZ
FINISH	! 退出后处理模块