ansys之东山船闸闸室参数化结构应力分析
finish!本命令只适合正常运行工况下内水压力计算。
/clear
/filname,ship
/title,solve ship
/units,si
!参数
h=28.6-10+0.4!高度，底板底面算起
h1=2.5!底板厚度
h2=27.11-9.6!内水位高度，底板地面算起
h3=10!基础厚度
C1=10!基础外延宽度
gm=9.8*1e3!水的重度，N/m3
l1=20!顺水流方向长度
b1=2.5!边墩厚度
b2=1.5!外贴角厚度
b3=14!闸室内净宽
b4=1!内贴角尺寸
b5=1!外延宽度
H4=5.4!外贴角高度
H5=7.9!填土高度
H6=25-9.6!外检修水位高度
/prep7
et,1,solid45 !定义第一组元素，船闸闸室
mp,ex,1,2.55e10!C20砼杨氏模量2.8e4n/mm2
mp,dens,1,2.5e3
mp,nuxy,1,0.167
mp,ex,2,2.5e8!土的弹性模量
mp,nuxy,2,0.4
mp,dens,2,2000
tb,dp,2
tbdata,1,1,32,30!砂粘聚力，摩擦角，膨胀角
K,1,0,0,0
K,2,B5+B2+B1+B3/2,0,0
K,3,B5+B2+B1+B3/2,0,H1
K,4,B5+B2+B1+B4,0,H1
K,5,B5+B2+B1,0,H1+B4
K,6,B5+B2+B1,0,H
K,7,B5+B2,0,H
K,8,B5+B2,0,H1+H4
K,9,B5,0,H1
K,10,0,0,H1
A,1,2,3,4,9,10!A1
A,9,4,5,6,7,8!A2
!VEXT,1,,,,20!VEXT, NA1, NA2, NINC, DX, DY, DZ, RX, RY, RZ
!VEXT,2,,,,20
AADD,1,2!1+2=A3
local,11,0,B5+B2+B1+B3/2,0!LOCAL, KCN, KCS, XC, YC, ZC, THXY, THYZ, THZX, PAR1, PAR2
arsym,X,3,0,0,0!A1,VSYMM, Ncomp, NV1, NV2, NINC, KINC, NOELEM, IMOVE
VEXT,3,,,,L1!V1左边
VEXT,1,,,,L1!V2右边
VEXT,3,,,,L1!V3左边
VEXT,1,,,,L1!V4右边
!VADD,ALL!1+3=2
!csys,0
block,-C1,B5+B2+B1+B3/2,0,L1,-H3,H5!左边基础V5
block,B5+B2+B1+B3/2,(B5+B2+B1+B3/2)*2+C1,0,L1,-H3,H5!右边基础V6
VSBV,5,1!5-1=8+7
VSBV,6,2!6-2=1+5
VDELE,7
VDELE,1
VADD,3,4!3+4=1
VADD,8,5!8+5=2
!内侧
wpoffs,,,H2!切闸墩
asbw,27
asbw,45
!外侧
wpoffs,,,-H2+H6!检修水位25-9.6
asbw,29
asbw,43
VGLUE,ALL!1+2=3基础+1闸室
!基础网格化
vsel,S,,,3
vatt,2,,1,
mshape,1,3d!四边形0，三角形1
mshkey,0!自由网格0，对应网格1
!esize,1.5!控制线段长度
vmesh,3
!闸室网格化
vsel,all
vsel,s,,,1
vatt,1,,1
mshape,1,3d
mshkey,0
vmesh,1
!nummrg,node
finish
/solu
antype,0
!边界条件
csys,0
dsys,0
nsel,s,loc,x,-C1!x方向两个侧面施加x方向位移约束
nsel,a,loc,x,(B5+B2+B1+B3/2)*2+C1
nplot
d,all,ux
allsel
nsel,s,loc,y,0
nsel,a,loc,y,L1
!nplot
d,all,uy
allsel
nsel,s,loc,z,-H3
d,all,all
!施加外载
allsel
!gplot
acel,0,0,9.8
!正常使用工况
!底板水压力
sfa,35,1,pres,gm*(H2-H1)
sfa,16,1,pres,gm*(H2-H1)
!侧墙
sfgrad,pres,0,z,H1,-gm
sfa,11,1,pres,gm*(H2-H1)
sfa,26,1,pres,gm*(H2-H1)
sfa,38,1,pres,gm*(H2-H1)
sfa,46,1,pres,gm*(H2-H1)
sftran
LSWRITE!保存第一负载
!检修工况
!底板水压力
!sfa,33,1,pres,gm*(H6)
!sfa,14,1,pres,gm*(H6)
!侧墙左侧
!sfgrad,pres,0,z,0,-gm
!sfa,45,,pres,gm*(H6)
!sfa,30,,pres,gm*(H6)
!sfa,31,,pres,gm*(H6)
!sfa,32,,pres,gm*(H6)
!右侧
!sfa,29,,pres,gm*(H6)
!sfa,42,,pres,gm*(H6)
!sfa,41,,pres,gm*(H6)
!sfa,40,,pres,gm*(H6)
!sftran
!LSWRITE!保存第二负载
allsel
outres,all,all
LSSOLVE,1
!solve
save,ship_static_rst.db
!finish
!/post1
!只显示闸室下部
!esel,s,type,,1!选择1号solid45单元
!esel,r,mat,,1!从中选择1号砼材料
!wpof,,10
!wpro,,90.000000,
!/GRAPHICS,POWER!full(全模式)、图形增强模式
!/TYPE,1,7
!/CPLANE,1
!/SHADE,1,1
!/HBC,1,0
!/REPLOT
SZ应力
页:
[1]