ex8.2(1).inp

ansys有限元分析程序

! 8.2  薄膜二维和三维非轴对称自由振动
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!．ANSYS二维分析程序
FINISH
/CLEAR, NOSTART	! 清除数据, 不加载初始化文件
! （1）定义工程选项、分析类型、单元类型、材料参数和实常数
/FILNAME, EX8.2(1)	!  定义工程文件名称为EX8.2(1)
/PREP7	! 加载前处理模块
/TITLE, EX8.2(1), NONAXISYMMTRIC VIBRATION OF A STRETCHED CIRCULAR MEMBRANE.
ET, 1, SHELL61	! 定义第1类单元轴对称协调壳单元SHELL61
R, 1, 1.27e-6	! 定义第1类实常数, 板厚1.27e-3 mm
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
! （2）通过关键点的定义, 建立模型
K, 1	! 在坐标原点（0, 0, 0）定义第1个关键点
K, 2, 76.2e-3	! 在最外侧（a=76.2 mm, 0, 0）位置建立关键点2
L, 1, 2	! 过关键点1, 2定义直线
LESIZE, 1, , , 9	! 指定剖分线单元的格数为9
LMESH, 1	! 对1号线执行单元剖分
! （3）选择边界节点, 施加位移约束
NSEL, S, LOC, X, 0	! 选择X＝0位置处的所有节点
D, ALL, UX	! 给上面选中的所有节点施加UX＝0 的位移约束
NSEL, S, LOC, X, 76.2e-3	! 重新选择X＝76.2mm位置处的所有节点
D, ALL, UX, , , , , UY	! 给选中的所有节点施加UX＝0和UY＝0的位移约束
NSEL, ALL	! 选择所有节点
M, ALL, UY	! 指定所有节点的UY自由度为主自由度
D, ALL, UZ	! 位移约束所有节点的沿着Z方向的位移自由度UZ
TREF, 0	! 指定参考温度（无应力状态的温度）为0度
BFUNIF, TEMP, -12	! 将薄膜的温度降低12度, 以生成薄膜应力
FINISH 	! 退出前处理模块
! （4）进入求解模块, 开始求解
/SOLU  	! 加载求解模块
ANTYPE, STATIC	! 定义分析类型为静力分析
PSTRES, ON	! 打开预应力分析选项
OUTPR, BASIC, ALL	! 指定输出所有节点的基本数据
SOLVE  	! 开始进行静力分析的求解
FINISH 	! 退出求解模块
! （5）重新进入求解模块, 进行模态分析
/SOLU 	! 加载求解模块
ANTYPE, MODAL	! 指定分析类型是模态分析
PSTRES, ON	! 打开预应力选项, 进行预应力模态分析
MXPAND, 3	! 指定最大展开的模态数为3
MODOPT, REDUC	! 指定利用降阶（REDUCED）方法求解矩阵特征值
MODE, 0	! 指定求解第0阶谐波模态
SOLVE 	! 开始进行模态分析求解
FINISH	! 退出求解模块
! （6）在一般后处理模块, 图线显示模态
/POST1	! 进入一般后处理模块
/WINDOW, 1, LTOP	! 在窗口的上面左侧定义第1号子图形窗口
SET, 1, 1 	! 读入第1个荷载步第1个子步的计算结果（第1阶模态）, 在窗口1显示
PLDISP, 1	! 带原来模型图, 显示振动模态
/NOERASE	! 保持屏幕不擦除状态
*GET, F11, MODE, 0, FREQ	! 提取第0阶模态的频率, 存入变量F11
/WINDOW, 1, OFF	! 关闭前面定义的1号窗口
/WINDOW, 4, RBOT	! 在窗口的下面右侧定义第4号子图形窗口
SET, 1, 2	! 读入第1个荷载步第2个子步的计算结果（第2阶模态）, 在窗口4显示
PLDISP, 1	! 带原来模型图, 显示变形后的模型
*GET, F12, MODE, 2, FREQ	! 提取第2阶模态的频率, 存入变量F12
/WINDOW, 4, OFF	! 关闭前面定义的4号窗口
FINISH	! 退出后处理模块
! （7）重新进入求解模态, 求解第1阶谐波模态
/SOLU	! 加载求解模块
MODE, 1	! 指定求解第1阶谐波模态
SOLVE	! 开始求解第1阶模态
FINISH	! 退出求解模块SOLUTION
/POST1	! 进入一般后处理模块
/WINDOW, 2, RTOP	! 在窗口的上面右侧定义第4号子图形窗口
SET, 1, 1	! 读入第1个荷载步的第1个子步的计算结果
PLDISP, 1	! 带原来模型图, 显示振动模态
*GET, F21, MODE, 0, FREQ	! 提取第0阶模态的频率, 存入变量F21
/WINDOW, 2, OFF	! 关闭前面定义的2号窗口
FINISH    	! 退出前处理模块PREP7
! （8）重新进入求解模态, 求解第1阶谐波模态
/SOLU	! 加载求解模块
MODE, 2	! 指定求解第2阶谐波模态
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLUTION
! （9）在一般后处理模块, 图线显示谐波模态
/POST1	! 进入一般后处理模块
/WINDOW, 3, LBOT	! 在窗口的下面左侧定义第3号子图形窗口
SET, 1, 1 	! 读入第1个荷载步的第1个子步的计算结果
PLDISP, 1 	! 带原来模型图, 显示振动模态
*GET, F31, MODE, 0, FREQ	! 提取第0阶模态的频率, 存入变量F31
*STATUS, PARM  	! 列出所有变量的值
FINISH	! 退出后处理模块