马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
- 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应力 |