| 查看: 2453 | 回复: 1 | ||
[求助]
悬索桥主缆线形非线性分析找形
|
|
!!!主缆线形确定 !!!给变量和函数赋值 !!!设计变量 !!! finish /clear !!!定义桥梁设计基本信息 E_x=100e9 Are=0.12035 Rdens=7850 L_LEN=160 L_LEN2=70 f_LEN=26.666 TOWER=31.927 !!!定义单元类型、实常数个数 *dim,r_c,,L_LEN+L_LEN2 /PREP7 ET,1,LINK10,,1,0 !!!定义主缆材料 MP,EX,1,E_x MP,PRXY,1,0.3 MP,ALPX,1,0.000012 MP,DENS,1,Rdens !!!定义实常数，先按照数量级指定 *do,i,1,L_LEN+L_LEN2 r,i,Are,0.03 *enddo !!!定义节点坐标，按抛物线方程: Y=AA*X**2+BB*X+CC AA=f_LEN/((L_LEN/2)**2) BB=0 CC=TOWER-f_LEN *do,i,0,L_LEN Xd=i-L_LEN/2 n,i+1,Xd ,AA*Xd**2+CC *enddo !!!定义缆索单元 MAT,1 TYPE,1 *DO,I,1,L_LEN REAL,I EN,I,I,I+1 *ENDDO D,001,ALL D,L_LEN+1,ALL !!!定义吊杆力 NSEL,S,,,6,156,5 F,ALL,FY,-1958500 /SOLU ACEL,,10 ANTYPE,4 NLGEOM,1 SSTIF,ON TIME,1e6 ALLSEL SOLVE finish !!!提取单元内力 /post1 ETABLE,,LS,1 *VGET,r_c,ELEM,1,ETAB,LS1 finish u_y=1E-9 *dowhile,u_y /prep7 !!!按抛物线重新定义节点坐标 deltf=-uy(L_LEN/2+1) AA=-deltf/((L_LEN/2)**2) BB=0 CC=deltf *do,n,0,L_LEN Xd=n-L_LEN/2 delty=AA*xd**2+CC NMODIF,n+1,nx(n+1),ny(n+1)+uy(n+1)+delty *enddo !!!!重新定义实常数 *DO,n,1,L_LEN *GET,prrr,RCON,n,const,2 rrr=(prrr+r_c(n)/E_x)/2 r,n,Are,rrr *enddo finish /SOLU ACEL,,10 ANTYPE,4 TIME,1 ALLSEL SOLVE finish !!!!提取单元内力 /POST1 ETABLE,LS1,LS,1 *VGET,r_c,ELEM,1,ETAB,LS1 !!!!定义迭代终止条件 *do,i,1,L_LEN/2+1 u_yy=abs(uy(i)) *if,u_yy,ge,u_y,then u_y=u_yy *endif *enddo *if,u_y,le,0.00001,exit finish *enddo ANSYS dowhile循环语句只进行了第一步循环然后就一直停留在那一个位置不动了，控制变量u_y的数值也保持不变了，因此进入了死循环怎么破。 一样大神们能够帮忙解决 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
317一志愿华南理工电气工程求调剂
已经有5人回复
-
化工专硕348，一志愿985求调剂
已经有6人回复
-
292求调剂
已经有3人回复
-
304求调剂
已经有4人回复
-
290求调剂
已经有6人回复
-
295求调剂
已经有5人回复
-
26申博
已经有4人回复
-
材料学调剂
已经有5人回复
-
264求调剂
已经有3人回复
-
311求调剂
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2016-05-19 10:21:37
