ansys初始应力文件编写
浏览:129531 回答:1
*cfopen,initial stress,IST
*do,i,1,6390
*vwrite,i
('!Stress for element',f5.0)
*vwrite
('! Sx Sy Sz Sxy Syz Sxz')
*vwrite,i
('eis',i,f5.0)
x=CENTRX(i)
y=CENTRX(i)
xx=SQRT(x*x+y*y)
yy=2*x*x/(x*x+y*y)-1
yy1=2*x*y/(x*x+y*y)
q=1e7
*if(xx.GT.6.AND.xx.LT.7) then
a=3/4*q-1/4*q*yy
b=3/4*q+1/4*q*yy
c=1/4*q*yy1
*else then
a=0
b=0
c=0
*endif
*vwrite,a,b,0,c,0,0
(6f15.3,/)
*enddo
*cfclos
这是我编写的一个圆孔隧洞圆周边界上单元的初始应力文件,但是在输入后出现了以下警告:'!Stress for element' is not a recognized prep7 command和do loop index same as outer index I.请问大家这是什么原因?
*do,i,1,6390
*vwrite,i
('!Stress for element',f5.0)
*vwrite
('! Sx Sy Sz Sxy Syz Sxz')
*vwrite,i
('eis',i,f5.0)
x=CENTRX(i)
y=CENTRX(i)
xx=SQRT(x*x+y*y)
yy=2*x*x/(x*x+y*y)-1
yy1=2*x*y/(x*x+y*y)
q=1e7
*if(xx.GT.6.AND.xx.LT.7) then
a=3/4*q-1/4*q*yy
b=3/4*q+1/4*q*yy
c=1/4*q*yy1
*else then
a=0
b=0
c=0
*endif
*vwrite,a,b,0,c,0,0
(6f15.3,/)
*enddo
*cfclos
这是我编写的一个圆孔隧洞圆周边界上单元的初始应力文件,但是在输入后出现了以下警告:'!Stress for element' is not a recognized prep7 command和do loop index same as outer index I.请问大家这是什么原因?