| 查看: 2828 | 回复: 1 | ||||
[交流]
【求助】利用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 |
回复此楼» 收录本帖的淘帖专辑推荐
 ANSYS命令流 |
» 猜你喜欢
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
b口会评
已经有6人回复
-
信息学部会评时间
已经有3人回复
-
祈祷青基必中
已经有9人回复
-
在职考研
已经有7人回复
-
MSER送审了还被拒稿
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 利用3-氨基噻吩的草酸盐能够提纯出3-氨基噻吩吗？？？
已经有6人回复
- 有用ANSYS研究玻璃材料的么 交流下！
已经有16人回复
- 我想写一个像QQ魔法表情那样的动画，怎么实现啊？
已经有7人回复
- 利用电化学方法制备石墨烯/四氧化三钴
已经有3人回复
- ansys中正弦曲线如何创建？且最好是可以进行变化成非标准形式的...
已经有11人回复
- ansys 14.0安装失败：license 无法加载
已经有21人回复
- 寻求ANSYS电磁方面的书
已经有4人回复
- 滑模控制抖振问题
已经有7人回复
- 【求助】ANSYS或abaqus 等CAE软件的迭代求解的初值问题
已经有7人回复
- 【求助】关于ANSYS在体模型的表面施加面载荷的问题
已经有13人回复
- 【资源】ANSYS操作命令与参数化编程
已经有18人回复
2楼2013-11-16 11:25:33
50