ex7.10.inp

ansys有限元分析程序

! 7.10  矩形截面薄壁悬臂梁的自由振动
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除数据, 不加载初始化文件
! （1）定义工程选项、单元类型、材料参数和实常数
/FILNAME, EX7.10	! 定义工程文件名称
/PREP7	! 加载前处理模块
/TITLE, EX7.10, TRANSVERSE VIBRATIONS OF A SHEAR BEAM.
ANTYPE, MODAL	! 指定分析类型为模态分析
ET, 1, SHELL28	! 使用考虑剪切和扭曲效应的壳单元SHELL28
R, 1, 0.1  	! 定义壳体厚度为0.1m
MP, EX, 1, 200E9	! 定义第1类材料的弹性模量EX=200GPa
MP, NUXY, 1, 0.27	! 定义第1类材料的泊松比NUXY=0.27
MP, DENS, 1, 7860.0	! 定义第1类材料的密度7860kg/m^3
! （2）定义节点位置和单元连接
/VIEW, 1, 1, 1, 1	! 改变窗口1的视点为等轴测位置
N, 1  	! 在(0, 0, 0)定义结点1
N, 2, 10 	! 在(10m, 0, 0)定义结点2
N, 3, 10, 10	! 在(10m, 10m, 0)定义结点3
N, 4, , 10	! 在(0, 10m, 0)定义结点4
NGEN, 9, 4, 1, 4, 1, , , 3.75	! 循环9次, 每次节点号增量4, 用前面定义的1号到4号节点模式, 增量为1, Z坐标增量3.75m, 自动生成其余节点
E, 1, 2, 6, 5	! 过节点1, 2, 6, 5定义单元
EGEN, 3, 1, -1	! 循环3次, 按照前面单元模式, 节点号增加1, 生成3个单元（包括前面一个单元）
E, 4, 1, 5, 8	! 过节点4, 1, 5, 8定义单元
EGEN, 8, 4, -4	! 按循环8次, 每次节点号增加4, 按照前面定义的4个单元模式, 生成其余单元
! （3）设置位移约束和耦合自由度
D, ALL, UZ 	! 限制所有节点的沿着Z方向的自由度UZ
D, 1, UX, , , 4, 1, UY	! 限制1号到4号节点的沿着X和Y方向的平移自由度UX和UY
NSEL, S, LOC, Z, 3.75	! 选择位于Z=3.75米位置上的所有节点
CP, 1, UX, ALL	! 将选择集中的节点在X方向的自由度UX耦合为一个自由度, 定义为集合1
CP, 2, UY, ALL	! 将选择集中的节点在Y方向的自由度UY耦合为一个自由度, 定义为集合2
NSEL, S, LOC, Z, 7.5	! 选择位于Z=7.5米位置上的所有节点
CP, 3, UX, ALL	! 将选择集中的节点在X方向的自由度UX耦合为一个自由度, 定义为集合3
CP, 4, UY, ALL	! 将选择集中的节点在Y方向的自由度UY耦合为一个自由度, 定义为集合4
NSEL, S, LOC, Z, 11.25	! 选择位于Z=11.25米位置上的所有节点
CP, 5, UX, ALL	! 将选择集中的节点在X方向的自由度UX耦合为一个自由度, 定义为集合5
CP, 6, UY, ALL	! 将选择集中的节点在Y方向的自由度UY耦合为一个自由度, 定义为集合6
NSEL, S, LOC, Z, 15	! 选择位于Z=15米位置上的所有节点
CP, 7, UX, ALL	! 将选择集中的节点在X方向的自由度UX耦合为一个自由度, 定义为集合7
CP, 8, UY, ALL	! 将选择集中的节点在Y方向的自由度UY耦合为一个自由度, 定义为集合8
NSEL, S, LOC, Z, 18.75	! 选择位于Z=18.75米位置上的所有节点
CP, 9, UX, ALL	! 将选择集中的节点在X方向的自由度UX耦合为一个自由度, 定义为集合9
CP, 10, UY, ALL	! 将选择集中的节点在Y方向的自由度UY耦合为一个自由度, 定义为集合10
NSEL, S, LOC, Z, 22.5	! 选择位于Z=22.5米位置上的所有节点
CP, 11, UX, ALL	! 将选择集中的节点在X方向的自由度UX耦合为一个自由度, 定义为集合11
CP, 12, UY, ALL	! 将选择集中的节点在Y方向的自由度UY耦合为一个自由度, 定义为集合12
NSEL, S, LOC, Z, 26.25	! 选择位于Z=26.25米位置上的所有节点
CP, 13, UX, ALL	! 将选择集中的节点在X方向的自由度UX耦合为一个自由度, 定义为集合13
CP, 14, UY, ALL	! 将选择集中的节点在Y方向的自由度UY耦合为一个自由度, 定义为集合14
NSEL, S, LOC, Z, 30	! 选择位于Z=30米位置上的所有节点
CP, 15, UX, ALL	! 将选择集中的节点在X方向的自由度UX耦合为一个自由度, 定义为集合15
CP, 16, UY, ALL	! 将选择集中的节点在Y方向的自由度UY耦合为一个自由度, 定义为集合16
NSEL, ALL	! 选择所有结点
FINISH	! 退出前处理模块PREP7
! （4）指定求解选项, 并求解
/SOLU	! 加载求解模块
MXPAND, 4	! 展开前面两阶重复的4个模态
MODOPT, REDUC	! 用降阶方法计算特征值
TOTAL, 8	! 设置自动主自由度生成到8
SOLVE	! 开始求解
! （5）提取频率计算结果
*GET, F1, MODE, 1, FREQ	! 获取第1阶频率, 存入变量F1
*GET, F2, MODE, 3, FREQ	! 获取第3阶频率, 存入变量F2
*STATUS, PARM	! 列表显示所有内存变量的值
FINISH    	! 退出求解模块SOLU
! （6）在一般后处理模块中, 显示模态形状, 列表显示内存变量的值
/POST1	! 进入一般后处理模块
SET, 1, 1	! 读入第1个荷载步的计算结果(第1阶振型)
/VUP, , Z	! 设置Z坐标轴向上
/VIEW, , , 1	! 从（0, 0, 1）位置向（0, 0, 0）方向观察模型
PLDISP, 1	! 图形显示第1阶模态
SET, 1, 3	! 读入第3个荷载步的计算结果(第3阶振型)
PLDISP, 1	! 图形显示第2阶模态
FINISH	! 退出后处理模块