24小时热门版块排行榜

返回列表

hugang0828

铁虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 13
在线: 14.1小时
虫号: 1087908
注册: 2010-09-02
专业: 结构工程

[交流] 【求助】利用ansys中apdl语言进行桥梁抖振时域分析的问题已有1人参与

本人利用ansys进行桥梁抖振的时域分析，需要利用APDL语言编程将脉动风荷载施加到桥梁结构上，但是编出的程序在运行过程中遇到了错误，下面是我加载的命令流，请各位高手帮小弟检查下，先谢了。
!***读入风速时程数据****
*DIM,FSX,ARRAY,6000,11,0, , ,
*DIM,FSY,ARRAY,6000,11,0, , ,
*CREAT,FSSC
*VREAD,FSX(1,1),'wind speed time history','txt','D:\Users\hp\desktop\BZ Yellow River',KIJ,11,6000
(11e12.9)
*VREAD,FSY(1,1),'wind speed time history','txt','D:\Users\hp\desktop\BZ Yellow River',KIJ,11,6000
(11e12.9)
*END
/INPUT,FSSC
!加载与求解
!*****************求解****************
/SOLU !进入求解器
ANTYPE,4 !指定分析类型为瞬态动力学分析
TRNOPT,FULL !瞬态动力学分析采用FULL法
TIMINT,OFF !关闭时间积分效应
OUTRES,BASIC,ALL !输出基本项，每一步都输出
KBC,1 !指定载荷
TIME,1E-10 !指定载荷步结束时间
NSUB,4,8 !指定载荷子步数为4，最大子步数为8
SSTIF,ON !打开应力刚化效应
ACEL,0,9.8,0 !施加重力加速度
ALLS !选中所有元素
SOLVE !求解第一个载荷步

!***在自重作用下，计算结构在风激励下的响应***
TIMINT,ON
OUTRES,BASIC,ALL
KBC,0
AUTOTS,ON
ALPHAD,ALPHA
BETAD,BETA
*DIM,Al,ARRAY,11
*DIM,Ad,ARRAY,11
*DIM,Lb,ARRAY,11
*DIM,Db,ARRAY,11
*DIM,Mb,ARRAY,11
*DIM,MMx,ARRAY,11
*DIM,Nod,ARRAY,11
*DIM,FFy,ARRAY,11
*DIM,FFz,ARRAY,11
Al(1)=304.371,279.00675,372.009,380.46375,380.46375,219.8235,380.46375,380.46375,372.009,279.00675,304.371
Ad(1)=162,148.5,198,202.5,202.5,117,202.5,202.5,198,148.5,162
Nod(1)=10,20,30,40,50,59,68,78,88,98,108
*SET,Ro,1.225
*SET,U,10
*SET,B,32.8
*SET,Cl,0.0222
*SET,Cd,0.8187
*SET,Cm,0.0093
*SET,Cll,4.776
*SET,Cdl,0
*SET,Cml,0.908
*DO,T,1,6000,1
*DO,n,1,11,1
*SET,v,FSX(T,n)
*SET,w,FSY(T,n)
TIME,0.1*T
Lb1=Ro*U*v*Al(n)*Cl
Lb2=1/2*Ro*U*w*(Ad(n)*Cd+Al(n)*Cll)
Lb(n)=Lb1+Lb2
Db1=Ro*U*v*Ad(n)*Cd
Db2=1/2*Ro*U*w*(Ad(n)*Cdl-Al(n)*Cl)
Db(n)=Db1+Db2
Mb1=Ro*B*U*v*Al(n)*Cm
Mb2=1/2*Ro*B*U*w*Al(n)*Cml
Mb(n)=Mb1+Mb2
FFy(n)=Lb(n)*cos(5/180*3.14)+Db(n)*sin(5/180*3.14)
FFz(n)=Db(n)*cos(5/180*3.14)-Lb(n)*sin(5/180*3.14)
MMx(n)=Mb(n)
!***********选择节点并施加风载荷*********
NSEL,S,NODE,,10,50,10    !通过点号选择节点
NSEL,A,NODE,,59,68,9
NSEL,A,NODE,,78,108,10
F,Nod(n),FY,,FFy(n)       !施加所选节点的风载荷
F,Nod(n),FX,,FFz(n)
NSUB,1                   !指定载荷子步数
SOLVE                !求解
*ENDDO
*ENDDO
SAVE
FINI

回复此楼