重力坝三维仿真分析.docVIP

重力坝三维仿真分析 363页 h=180 h1=100 h2=80 h3=100 gm=1000 fial=90-atan(0.75)*180/3.1415926 fia2=90-fial !定义单元类型及材料属性 /prep7 et,1,plane42 et,2,solid65 et,3,solid45 !1号材料 mp,ex,1,2.85e10 !100米以下混凝土的材料特性 mp,prxy,1,0.167 tb,conc,1,1,9, tbdata,,0.3,1,1.96e6,22e6 !张开剪切传递系数为0.3，抗拉强度fc=1.96e6pa，抗压强度fs=22e6pa mp,dens,1,2400 !2号材料 mp,ex,2,2.6e10 !100米以上混凝土的材料特性 mp,prxy,2,0.167 tb,conc,2,1,9, tbdata,,0.3,1,1.2e6,17.5e6 !张开剪切传递系数为0.3，抗拉强度fc=1.2e6pa，抗压强度fs=17.5e6pa mp,dens,2,2400 !3号材料 mp,ex,3,2.9e10 mp,prxy,3,0.3 mp,dens,3,2600 !建立断面模型。先建立5个小面和1个大面，再进行布尔操作得到7个面的几何模型 k,1 k,2,0.9*0.75*h+0.1*h k,3,0.1*h,0.9*h k,4,,0.9*h l,1,2 l,2,3 l,3,4 l,4,1 al,1,2,3,4 rectng,,0.1*h,0.9*h,h !画坝顶矩形 rectng,-1.5*h,0,-2*h,0 !坝基左面矩形 rectng,0,0.9*0.75*h+0.1*h,-2*h,0 !坝基正下方矩形 rectng,0.9*0.75*h+0.1*h,(2+0.9*0.75)*h,-2*h,0 !坝基右面矩形 rectng,-1.5*h,(2+0.9*0.75)*h,-2*h,h !总矩形 aplot !对断面进行布尔操作。使用布尔操作中的重叠技术，即将分离的同等级对象(线段、面积、体积等)变为数个对象的连!续体，其中对象所有重叠边界所围成的区域都自成一个对象。 aovlap,all !将面单元进行粘贴布尔操作 /pnum,line,1 !打开面，线的号码开关 /pnum,area,1 nummrg,all numcmp,all aplot !对面进行网格划分 !划分3号面的网格，对应于坝体的顶部 lsel,s,,,3,5,2 lesize,all,,,5 lsel,s,,,12,13,1 lesize,all,,,2 amesh,3 !划分1号面的网格，对应于坝体的下半部分 lsel,s,,,2,4,2 lesize,all,,,18 lsel,s,,,1 lesize,all,,,5 amesh,1 eplot !划分2号面的网格，对应于坝体的正下方基岩 lsel,s,,,11 lesize,all,,,5 lsel,s,,,9,10,1 lesize,all,,,8,4 amesh,2 !划分4号面的网格 lsel,s,,,14 lesize,all,,,5,4 lsel,s,,,7 lesize,all,,,8,4 lsel,s,,,6 lesize,all,,,5,0.25 amesh,4 eplot !划分5号面的网格 lsel,s,,,15,16,1 lesize,all,,,8,4 lsel,s,,,8 lesize,all,,,8,0.25 amesh,5 eplot !在对剩余部分的几何面划分网格之前，需要先将相关的部分线段组合在一起，以形成由四条边构成的四边形，然后再划分网格，完成后必须将!组合线段删除。 !连接4号线和13号线，同时划分7号面的网格 lsel,s,,,4,13,9 lccat,all lsel,s,,,19 lesize,all,,,5,4 lsel,s,,,20,21,1 lesize,all,,,20 amesh,7 eplot !连接2号线和12号线，并划分6号面的网格 lsel,s,,,2,12,10 lccat,all lsel,s,,,17 lesize,all,,,8,4 lsel,s,,,18,22,4 lesize,all,,,20 amesh,6 eplot !删除前面连接的线元素