反应谱分析.txt

是用apdl语言编的

!注意：在进行反应普分析时，若想要得到单元的应力状态，必须将solid65单元的非线性属性去掉


/prep7
resume,dam_model,db
keyopt,2,5,2			!输出solid65单元的节点应力
/solu
!施加边界条件
csys,0
dsys,0
nsel,s,loc,x,(2+0.9*0.75)*h
nsel,a,loc,x,-1.5*h
nplot
d,all,ux
allsel

nsel,s,loc,z,h
nsel,a,loc,z,-(1+1.5)*h
nplot
d,all,uz
allsel

nsel,s,loc,y,-2*h
nplot
d,all,uy
allsel
gplot

!模态分析
antype,modal
modopt,subsp,10
solve
save,dam_dynamic_rst11,db
finish

!谱分析
/solu
antype,spectr
spopt,sprs,10,yes
svtyp,2
sed,1,1,
freq,0.3444,0.3502,0.3546,0.3739,0.4175,0.4613,0.5114,0.5121,0.5811,0.5837
sv,1.226,1.2079,1.1945,1.1389,1.0312,0.9427,0.859,0.858,0.7657,0.7627
solve
save,dam_dynamic_rst2,db
finish

!模态扩展
/solu
antype,modal
expass,on
mxpand,10,,,yes,0.005
solve
save,dam_dynamic_rst3,db
finish

!合并模态
/solu
antype,spectr
srss,0.15,disp
solve
save,dam_dynamic_rst4,db
finish

!查看结果
/post1
/input,,mcom
esel,s,type,,2
eplot
SET,FIRST
PLNSOL,U,SUM,1,1
set,next
PLNSOL,U,SUM,1,1
set,next
PLNSOL,U,SUM,1,1

SET,FIRST
PLNSOL,S,1,0,1
PLNSOL,EPTO,1,0,1
SET,next
PLNSOL,S,1,0,1
PLNSOL,EPTO,1,0,1
finish

/post26
csys,0
nsel,s,loc,y,h
nsel,r,loc,z,-0.75*h
nsel,r,loc,x,0
nplot
NSOL,2,1466,U,X,nux
NSOL,3,1466,U,Y,nuy
NSOL,4,1466,U,z,nuz
XVAR,1  
PLVAR,2,3,4

ESOL,2,1,19,S,1,s1_19   	!变量2代表着：上游坝底侧面，1号单元上19号节点的第一主应力
ESOL,6,33,20,S,1,s1_20 		!变量3代表着：下游坝底侧面，33号单元上20号节点的第一主应力
ESOL,7,761,9,S,1,s1_9		!变量4代表着：上游坝顶侧面，761号单元上9号节点的第一主应力
ESOL,8,746,7,S,1,s1_7		!变量5代表着：下游坝顶侧面，746号单元上7号节点的第一主应力
/axlab,x,time
/axlab,y,stress1[pa]
XVAR,1  
PLVAR,5,6,7,8