24小时热门版块排行榜

返回列表

厚德载物le

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 62.5
帖子: 10
在线: 8.6小时
虫号: 4083036
注册: 2015-09-18
专业: 交通工程

[求助] 悬索桥主缆线形非线性分析找形

!!!主缆线形确定
!!!给变量和函数赋值
!!!设计变量
!!!
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客户端

回复此楼