请问谁用过这个命令:UPGEOM, FACTOR, LSTEP, SBSTEP, Fname, Ext, -- —命令?
我做的是一个三层框架(梁板柱型的)的地震时程分析。首先用ansys建立了一个三层框架模型,模型数据为:
框架总体尺寸是3m×4m×9m。层高3m。框架柱和梁均采用方钢,柱截面尺寸为200mm×200mm,梁截面尺寸为200mm×200mm,各层板厚30mm。所使用材料的的弹性模量E=2.06*10^11Pa,泊松比为0.3 ,密度为7.85*10^3kg/m^3。
向该框架结构体系输入时程为20s的地震波,时间间隔为0.01s;得到其位移响应。我做的是损伤识别。在输入地震波的时候,在中间10s的时候设定结构的25号单元弹性模量降低为原来的95%,以后的时间里该单元弹性模量就一直是原来的95%了,以此来模拟损伤。10s后继续输入地震波,这时候怎么样才能把10s前算的位移响应作为现在输入地震波时的位移初值?谁能帮帮俺?
我建框架模型的命令流是:
/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=2.06E11*0.75
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 !此处为了以后更改弹性模量设定的材料2.
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
下面是输入地震波并求解(地震波是一个名为jili.txt的文件;里面的时间列名为time,加速度列名为accel;
/prep7
*dim,acez,TABLE,2001,1,1,time,accel
*create,ff !
*tread,acez,jili,txt,,1
*end
/input,ff
finish
/config,nres,2001
/solu
antype,trans
t1=0.01
t2=20
dt=0.01
*do,tm,t1,t2,dt
time,tm
alpha,0
beta,0.0039
acel,acez(tm),0,0
solve
*enddo
算节点77的位移的命令流:
nsel,s,loc,z,3
/POST26
nsol,2,77,U,x,ux_1
PRVAR,2, , , , , ,
plvar,2
finish
这样算完后得到的是无损伤时结构的第77测点的位移响应。
但我要解决的是在10s时结构第5单元弹性模量降低5%后算出结构的位移响应。
显然需要用到一个命令将10s前分析所得的位移加到框架模型上,并生成一个已变形的几何形状。即要用到UPGEOM, FACTOR, LSTEP, SBSTEP, Fname, Ext, -- —命令,但不知道应该具体怎么使,哪位高手能帮我看看应该怎么做好么?多谢了!!
页:
[1]