ex3.13(1).inp

ansys有限元分析程序

! 3.13  塔机标准节内力分析(梁单元)
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/Clear
/FILNAME, EX3.13(1)	! 定义文件名称
/TITLE, EX3.13(1), STANDARD SECTION OF QTZ25 TOWER CRANE with BEAM4.
/PREP7
! （1）定义塔身标准节上的节点
N, 1, 650, 0, -650	$	N, 2, 650, 0, 650	$	N, 3, -650, 0, 650 
N, 4, -650, 0, -650	$	N, 5, 650, 100, -650	$	N, 6, 650, 100, 650 
N, 7, -650, 100, 650	$	N, 8, -650, 100, -650	$	N, 9, 650, 1200, -650 
N, 10, 650, 1200, 650	$	N, 11, -650, 1200, 650	$	N, 12, -650, 1200, -650 
N, 13, 650, 2300, -650	$	N, 14, 650, 2300, 650 
N, 15, -650, 2300, 650	$	N, 16, -650, 2300, -650 
N, 17, 650, 2400, -650	$	N, 18, 650, 2400, 650 
N, 19, -650, 2400, 650	$	N, 20, -650, 2400, -650 
! （2）定义单元类型和材料参数
ET, 1, BEAM4	! 定义第1类单元是三维梁单元BEAM4
MP, EX, 1, 210E3	! 定义第1类材料的弹性模量EX=210GPa=210E3 N/mm^2
MP, PRXY, 1, 0.3	! 定义第1类材料的泊松比PRXY=0.3
MP, DENS, 1, 7800e-12	! 定义第1类材料的密度DENS=7800 kg/m^3
Acel, 0, 9800, 0	! 定义沿着-Y方向的重力加速度9800 m/s^2
R, 1, 4536, 2104e4, 2104e4, 150, 150	! 塔身外立柱150@10
R, 2, 766, 910540, 910540, 80, 80	! 塔身水平杆，斜支撑 80@5
! （3）定义标准节的外框立柱杆件
Real, 1	! 塔身外立柱  150@10
E, 1, 5	$	E, 2, 6	$	E, 3, 7 
E, 4, 8	$	E, 5, 9	$	E, 6, 10 
E, 7, 11	$	E, 8, 12	$	E, 9, 13 
E, 10, 14	$	E, 11, 15	$	E, 12, 16 
E, 13, 17	$	E, 14, 18	$	E, 15, 19 
E, 16, 20 
! （4）定义标准节的水平杆和斜拉杆
Real, 2	! 塔身水平杆，斜支撑  80@5
E, 5, 6	$	E, 9, 10	$	E, 13, 14 
E, 6, 7	$	E, 10, 11	$	E, 14, 15 
E, 7, 8	$	E, 11, 12	$	E, 15, 16 
E, 5, 8	$	E, 9, 12	$	E, 13, 16 
E, 5, 10	$	E, 10, 13	$	E, 6, 11 
E, 11, 14	$	E, 7, 12	$	E, 12, 15 
E, 8, 9	$	E, 9, 16	$	E, 9, 11 
E, 13, 15	$	E, 5, 7
! （5）施加底部位移约束和集中力荷载
D, 1, UX, 0, , , , UY, UZ  
D, 2, UX, 0, , , , UY
D, 3, UY, 0  
D, 4, UY, 0, , , , UZ
F, 17, FY, 5E4	! 平衡臂一侧立柱的轴向拉力
F, 18, FY, 5E4	! 平衡臂一侧立柱的轴向拉力
F, 19, FY, -15E4	! 起重臂一侧立柱的轴向压力
F, 20, FY, -15E4	! 起重臂一侧立柱的轴向压力
FINISH	! 退出前处理模块
! （6）进入求解模块定义塔臂端部满载情况，并进入求解
/SOLU	! 进入求解模块
Solve	! 进行满载情况对称位置的静力分析求解
F, 17, FY, 3.07E4	! 平衡臂一侧立柱的轴向拉力
F, 18, FY, -4E4	! 中间立柱的轴向压力
F, 19, FY, -11.07E4	! 起重臂一侧立柱的轴向压力
F, 20, FY, -4E4	! 中间立柱的轴向压力
Solve 	! 进行满载情况45度方向的静力分析求解
F, 17, FZ, -2500/2.6	! 平衡臂一侧立柱的轴向拉力
F, 18, FX, 2500/2.6	! 中间立柱的轴向压力
F, 19, FZ, 2500/2.6	! 起重臂一侧立柱的轴向压力
F, 20, FX, -2500/2.6	! 中间立柱的轴向压力
F, 17, FX, -2500/2.6	! 平衡臂一侧立柱的轴向拉力
F, 18, FZ, -2500/2.6	! 中间立柱的轴向压力
F, 19, FX, 2500/2.6	! 起重臂一侧立柱的轴向压力
F, 20, FZ, -2500/2.6	! 中间立柱的轴向压力
Solve 	! 进行满载情况45度方向的静力分析求解
FINISH  	! 退出求解模块
! （7）进入后处理模块显示结果
/post1  	! 进入一般后处理模块
SET, 1, 1	! 读入第1个荷载步的计算结果
pldisp, 2	! 图形显示变形情况
Etable, LS1, LS, 1	! 起始端的轴向应力SDIR
Etable, LS6, LS, 6	! 终止端的轴向应力SDIR
Etable, LS2, LS, 2	! 起始端的＋Y一侧弯曲应力SBYT
Etable, LS3, LS, 3	! 起始端的－Y一侧弯曲应力SBYB
Etable, LS4, LS, 4	! 起始端的＋Z一侧弯曲应力SBZT
Etable, LS5, LS, 5	! 起始端的－Z一侧弯曲应力SBZB
Etable, V, Volu	! 用单元体积定义单元表V
Etable, smax, NMISC, 1	! 起始端的最大拉应力
Etable, smin, NMISC, 2	! 起始端的最大压应力
pretab, LS1, LS2, LS3, LS4, LS5, LS6, smax, smin	! 列表显示指定单元表的内容
/TITLE, Axial stress distribution in the first load case.
PLLS, LS1, LS6	! 图形显示单元轴向应力
! （8）显示45度情况下的结果
Set, 2, 1	! 读入第2个荷载步的计算结果
pldisp, 2	! 图形显示变形情况
Etable, SDIRI, LS, 1	! 起始端的轴向应力SDIR
Etable, SDIRJ, LS, 6	! 终止端的轴向应力SDIR
pretab, SDIRI, SDIRJ	! 列表显示所有单元表内容
/TITLE, Axial stress distribution in the second load case.
PLLS, SDIRI, SDIRJ	! 图形显示单元轴向应力
Set, 3, 1	! 读入第3个荷载步的计算结果
Etable, SDIRI3, LS, 1	! 起始端的轴向应力SDIR
Etable, SDIRJ3, LS, 6	! 终止端的轴向应力SDIR
/TITLE, Axial stress distribution in the third load case.
PLLS, SDIRI3, SDIRJ3	! 图形显示单元轴向应力
pretab, LS1, LS6, SDIRI, SDIRJ, SDIRI3, SDIRJ3	! 列表显示单元表内容
FINISH	! 退出后处理模块