用生死单元求解材料性质改变时的结构

大树357

finish
/clear
/prep7
l0=1000
b0=10
h0=20
e0=2.0e5
a0=b0*h0
i0=b0*h0*h0*h0/12
coef1=3
et,1,beam3
mp,ex,1,e0
mp,prxy,1,0.3
mp,ex,2,coef1*e0
mp,prxy,2,0.3
r,1,a0,i0,h0
k,1
k,2,l0
l,1,2
lesize,all,,,10
lmesh,all
egen,2,0,all,,,1
d,1,all
/solu
antype,0
nlgeom,on
outres,all,all
autots,off
nropt,full
time,1
nsubst,10
f,2,fy,2000
esel,s,mat,,2
ekill,all
allsel
solve
time,2
esel,s,mat,,2
ealive,all
allsel
f,7,fy,4440
nsubst,20
solve
finish
/post26
nsol,2,2,u,y
nsol,3,2,u,x
nsol,4,7,u,y
nsol,5,7,u,x
rforce,6,1,f,y
rforce,7,1,m,z
prvar,2,3,4,5,6,7
prod,8,6,,,,,,-1
xvar,2
plvar,8

大树357

这是王新敏书中的一个命令流，用来模拟在P1荷载作用时材料弹性模量为E0，之后材料弹性模量变为E1且又增加了荷载P2
在荷载步2时激活材料号为2的单元，再施加荷载。
但是这时候并没有杀死材料号1的荷载，整个结构是不是变成了在相同位置处，同时存在两种单元，材料1和材料2.这不是和实际情况不符吗？实际情况在time2时只有材料2，？？

大树357

我后来在time2时，分别激活单元材料2，杀死单元材料1，但是求出的结果和书中的结果却不同

大树357

而且按照上述命令流，当激活单元2进行求解时，相同位置处是存在两种单元的，ansys可以查看到，并不是只存在单元2。这个就奇怪了，不是和实际情况不符合了吗？

IRISXYF

亲呐，我最近也是发现这本书上的这个问题，本来想用王新敏老师的这个方法实现材料属性改变，但是后来感觉他这个有问题，正如楼主所说一样，材料2激活后计算的结果并非P2对应MAT2的结果，而是MAT1/2都活着了，我理解的是第二步变成材料属性的叠加了，毕竟单元共用所有的节点，也必然会对所有的节点都有刚度贡献。不知道楼主解决了没有？难道书上这个真的错了？？？？迫切求交流2623144801.