常量力移动载荷分析.txt

《ANSYS工程分析进阶实例》-王呼佳-源代码 看文件吧.不多说了.

!在移动载荷作用下，桥梁将发生振动，产生的变形和应力都比静载荷作用大。移动载荷的这种动力效应是不可忽视的，若在载荷处于最不利的静力作用位置同时满足共振条件，那么将会发生较大的动态响应，从而导致桥梁破坏。本例将移动载荷作为匀速移动的常量力来模拟悬索桥在移动载荷下的动态响应。为了提高分析效率，我们在前节初始位置确定分析中已经创建了宏“MODEL.MAC”，在本节只需要直接调用该宏，并将缆索单元的初应变值改为0.0043，即可得到移动载荷作用分析的有限元模型。
!参数定义
FINI								!退出
/CLE
/FILNAME,CABLE_BRIDGE,1
/TITLE,YIDONG HEZAI OF CABLE_BRIDGE
!建立模型
!"MODEL"宏文件应放在工作目录下。
MODEL								!运行宏生成模型
!加载及求解
!本例中的移动载荷视为匀速常量力载荷行分析。其中载荷MG＝2000N，简化成四个车轮上的相等集中载荷，车轮间距为8M，正好为两个桥面单元的长度，则简谐力F=500COS(10T)，载荷移动速度V＝120公里/小时。
!参数定义
!定义参数
DELTL=4							!桥面单元长度
V=120								!载荷移动速度
DELTT=DELTL/V*3.6				!荷载经过一个单元所用时间
F=1000								!常量力大小
W=10								!简谐力的圆频率
LF=8								!前后车轮间距
NUM=54							!桥面单元数
MM=NUM+1+LF/DELTL			!循环最大数目
!设置分析选项
FINI
/SOLU
ANTYPE,4							!指定分析类型为瞬态动力学分析
TRNOPT,FULL						!瞬态动力学分析采用FULL法
NLGEOM,ON
NROPT,FULL  
!计算重力的影响
TIMINT,OFF
TIME,1E-8
KBC,1
NSUBST,5
SSTIF,ON
ACEL,,9.8
ALLS
SOLVE                                 	!求解
!施加常量力移动载荷并求解
TIMINT,ON                             !打开时间积分开关
KBC,0
*DO,I,1,MM,1
OUTRES,ALL,ALL					!输出每一步的结果
TIME,I*DELTT
NSUBST,5

!***********常量力情况*******************
*IF,I,LT,LF/DELTL,THEN
	FDELE,ALL,ALL					!删除所有集中力载荷
	NSEL,S,LOC,X,I*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	NSEL,S,LOC,X,I*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,BW-1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	SOLV
*ELSEIF,I,LT,(NUM+1)
	FDELE,ALL,ALL
	NSEL,S,LOC,X,I*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	NSEL,S,LOC,X,I*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,BW-1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	NSEL,S,LOC,X,(I-2)*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	NSEL,S,LOC,X,(I-2)*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,BW-1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	SOLV
*ELSE
	FDELE,ALL,ALL
	NSEL,S,LOC,X,(I-2)*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	NSEL,S,LOC,X,(I-2)*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,BW-1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	SOLV
*ENDIF
*ENDDO
!************************
!**************简谐力情况******************
*IF,I,LT,LF/DELTL,THEN
	FDELE,ALL,ALL					!删除所有集中力载荷
	NSEL,S,LOC,X,I*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	NSEL,S,LOC,X,I*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,BW-1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	SOLV
*ELSEIF,I,LT,(NUM+1)
	FDELE,ALL,ALL
	NSEL,S,LOC,X,I*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	NSEL,S,LOC,X,I*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,BW-1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	NSEL,S,LOC,X,(I-2)*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	NSEL,S,LOC,X,(I-2)*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,BW-1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	SOLV
*ELSE
	FDELE,ALL,ALL
	NSEL,S,LOC,X,(I-2)*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	NSEL,S,LOC,X,(I-2)*DELTL
	NSEL,R,LOC,Y
	NSEL,R,LOC,Z,BW-1.89
	F,ALL,FY,-F					!施加常量力载荷
	ALLS
	SOLV
*ENDIF
*ENDDO
!***************************************
FINI
/POST1						!进入通用后处理
SET,30,2					!查看30步第2子步即载荷移动到跨中附近的结果
PLNSOL, U,Y				!显示y方向位移云图
PLNSOL, S,EQV				!显示等效应力云图
/POST26					!进入时间历程后处理
NUMVAR,20				!指定允许的变量数
NSOL,2,677,U,Y, UY_CENT  !取出中间节点的Y向位移数据
DERIV,3,2,1,,VY_CENT		!对中间节点的Y向位移微分运算，得到Y向速度
PLVAR,2					!显示中间节点的Y向位移随时间变化曲线
PLVAR,3					!显示中间节点的Y向速度随时间变化曲线