如何在地震波时程分析过程中更改结构某一单元的弹性模量？
地震波时程分析中，时程为0.2s，（这里时程如此短只是为了试算验证结果）。时间间隔为0.01s。先输入前0.1s的地震波，求解计算，然后更改5号单元的弹性模量，以模拟结构的损伤。然后接下来利用重启动求解，输入后0.1s的地震波，可得到整个时程的位移响应。可是为什么得到的位移结果和没改变弹性模量而得到的位移结果是一样的呢？难道我用的更改弹性模量的命令是无效的么？求高手指点一下，我的输入地震波并进行求解的命令流如下：
/prep7
*dim,ace1,TABLE,11,1,1,time,accel
*dim,ace2,TABLE,10,1,1,time,accel
*create,ff
*tread,ace1,bo1,txt,,1
*tread,ace2,bo2,txt,,1
*end
/input,ff
finish
/config,nres,21
/solu
antype,trans
t1=0.01
t2=0.1
dt=0.01
*do,tm,t1,t2,dt
time,tm
alpha,0
beta,0.0039!阻尼比
acel,ace1(tm),0,0
solve
*enddo
esel,s,elem,,5
mpchg,2,all!把5号单元改成为2号材料属性，当然这个属性在建模时已经定义
save
finish
/filname,rtest1
/solu
antype,,rest
t1=0.11
t2=0.2
dt=0.01
*do,tm,t1,t2,dt
time,tm
alpha,0
beta,0.0039!阻尼比
acel,ace2(tm),0,0
solve
*enddo
nsel,s,loc,z,3
/POST26
nsol,2,77,U,x,ux_1
PRVAR,2, , , , , ,
plvar,2
finish
附件bo1和bo2是地震波数据。
求高手帮帮我，谢谢了！！
框架模型的命令流如下：
框架模型的命令流如下：finish
/clear
/filname,rtest1
/prep7
T=0.03
S=0.04
IZZ=0.0004/3
IYY=0.0004/3
B=0.2
H=0.2
EX=2.06E11
PRXY=0.3
DENS=7.85E3
EX1=2.06E11*0.95
EX2=2.06E11*0.85
EX3=10
ET,1,SHELL63
ET,2,BEAM4
R,1,T
R,2,S,IZZ,IYY,B,H
MP,EX,1,EX
MP,PRXY,1,PRXY
MP,DENS,1,DENS
MP,EX,2,EX1
MP,PRXY,2,PRXY
MP,DENS,2,DENS
/view,,-1.2,-1,1
K,1
K,2,,,3
K,3,,,6
K,4,,,9
K,5,4
K,6,4,,3
K,7,4,,6
K,8,4,,9
L,1,2
L,2,3
L,3,4
L,5,6
L,6,7
L,7,8
TYPE,2
MAT,1
REAL,2
ESYS,0
LESIZE,ALL,,,5,1
LMESH,ALL
L,2,6
L,3,7
L,4,8
TYPE,2
MAT,1
REAL,2
ESYS,0
LESIZE,ALL,,,2,1
LMESH,ALL
LGEN,2,ALL,,,,3,,,0 !复制
L,2,10
L,3,11
L,4,12
L,8,16
L,7,15
L,6,14
TYPE,2
MAT,1
REAL,2
ESYS,0
LESIZE,ALL,,,2,1
LMESH,ALL
A,2,6,14,10
TYPE,1
MAT,1
REAL,1
ESYS,0
ALLSEL
/PNUM,LINE,1
LSEL,S,,,13,16
LESIZE,ALL,,,2,1
AMESH,1
ALLSEL
AGEN,2,1,,,,,3
AGEN,2,ALL,,,,,3!复制
ALLSEL
NUMMRG,ALL
NUMCMP,ALL
NSEL,R,LOC,Z
D,ALL,,,,,,ALL
ALLSEL
/VIEW,,-0.2,-0.9,0.4
EPLOT
FINISH
页:
[1]