这个我有一个小悬臂梁程序,你可以参考以下,希望对你有用!
/PREP7 !进入前处理模块PREP7
ET, 1, BEAM3 !定义第一类单元为平面梁单元BEAM3
ET, 2, MASS21, , ,4 !定义第二类单元为质量阻尼单元MASS21
R, 1, 0.003, 6.25e-7, 0.05 !定义单元的第一类实常数:Area,Inertia,Height
R, 2, 0.1 !定义单元的第二类实常数:集中质量
MP, EX, 1, 207e9 !定义第一类材料的弹性模量EX
N, 1, 0, 0 !定义各个结点
N, 2, 0.04, 0
N, 3, 0.08, 0
N, 4, 0.12, 0
TYPE, 1 !使用第一类单元
REAL, 1 !使用第一类实常数
MAT, 1 !使用第一类材料
E, 1, 2 !按上面设置定义单元
E, 2, 3
E, 3, 4
TYPE, 2 !使用第二类单元
REAL, 2 !使用第二类实常数
E, 4 !定义四号单元(集中质量)
FINISH !退出后模块
/SOLU !进入求解模块SOLUTION
ANTYPE, MODAL !申明求解类型是模态分析
MODOPT,LANB,5 !使用Block Lanczos方法求解前5阶振型和频率
D, 1, ALL, 0 !固定1号结点
M, 2, UY, 4, 1 !定义2号到4号结点的三个结点的Y方向为主自由度
SOLVE !开始求解
FINISH !退出后模块
/POST1 !进入后处理模块POST1
SET, 1, 1 !读入第一阶频率和振型
PLDISP ! 在图形窗口显示结构变形
ANMODE,10,0.05 !用10帧每隔0.05秒钟的动画显示振型
|