下面是我的工作,我做了2个算例,其中一个是用初始应力来实现应力的叠加,用upgeom命令来实现节点坐标的更新;另一个算例是只用pgeom命令来实现节点坐标的更新,看看能否实现应力的叠加。
计算的结果发现:只用pgeom命令是不能实现应力的叠加的,也就是前步计算的应力不能作为初始应力来叠加到后续的计算。
命令流如下:
算例1:用初始应力把上一步的求解应力结果作为初始应力加到后续计算上。用upgeom命令更新节点坐标
/PREP7
ET,1,plane42
MP,EX,1,1e10
MP,PRXY,1,0.25
MP,DENS,1,2500
TB,BISO,1,1,2,
TBTEMP,0
TBDATA,,200e6,3.3e9,,,,
rectng,0,5,0,2
esize,0.1
amesh,all
FINISH
/SOL
iswrite,1
nlgeom,on
acel,0,9.8,0
lsel,,line,,4
dl,all,,all
allsel,all
lsel,,line,,1
dl,all,,uy
allsel,all
lsel,,line,,2
dl,all,,ux,-0.5
allsel,all
SOLVE
FINISH
/prep7
upgeom,1,,,test,rst
finish
/SOL
ISFILE,READ,test,ist, ,1
nlgeom,on
lsel,,line,,4
dl,all,,all
allsel,all
lsel,,line,,1
dl,all,,uy
allsel,all
lsel,,line,,2
dl,all,,ux,-0.5
allsel,all
SOLVE
FINISH
/prep7
upgeom,1,,,test,rst
finish
算例2:只用upgeom命令更新节点坐标,看看能否实现应力的叠加
/PREP7
ET,1,plane42
MP,EX,1,1e10
MP,PRXY,1,0.25
MP,DENS,1,2500
TB,BISO,1,1,2,
TBTEMP,0
TBDATA,,200e6,3.3e9,,,,
rectng,0,5,0,2
esize,0.1
amesh,all
FINISH
/SOL
nlgeom,on
acel,0,9.8,0
lsel,,line,,4
dl,all,,all
allsel,all
lsel,,line,,1
dl,all,,uy
allsel,all
lsel,,line,,2
dl,all,,ux,-0.5
allsel,all
SOLVE
FINISH
/prep7
upgeom,1,,,test,rst
finish
/SOL
nlgeom,on
lsel,,line,,4
dl,all,,all
allsel,all
lsel,,line,,1
dl,all,,uy
allsel,all
lsel,,line,,2
dl,all,,ux,-0.5
allsel,all
SOLVE
FINISH
/prep7
upgeom,1,,,test,rst
finish |