文档介绍：静态分析 h=180 ! 大坝高度 fini_days=30 ! 完成大坝浇筑的总的天数 ht=h/fini_days ! 每天完成的大坝高度 h1=100 ! 上游水位高度 h2=80 ! 下游水位高度 h3=100 ! 坝体不同混凝土材料的分界线,为 100m gm=1000 !gm 为水的容重 fia1=90-atan()*180/ ! 计算下游斜面夹角 fia2=90-fia1 ! 计算坐标旋转角度/prep7 et,1,plane42 ! 辅助划分网格之用 et,2,solid65 ! 混凝土单元 et,3,solid45 ! 岩石单元 keyopt,2,5,2 ! 输出 solid65 单元的节点应力!1 号材料 mp,ex,1, !100m 以下混凝土的材料特性 MP,PRXY,1, TB,CONC,1,1,9, TBDATA,,,1,,22e6 ! 张开剪切传递系数为 , 抗拉强度 fc=, 抗压强度 fs=22e6pa, mp,dens,1,2400 !2 号材料 mp,ex,2, !100m 以上混凝土的材料特性 MP,PRXY,2, TB,CONC,2,1,9, TBDATA,,,1,, ! 张开剪切传递系数为 ,抗拉强度 fc=, 抗压强度 fs=, mp,dens,2,2400 !3 号材料 mp,ex,3, mp,PRXY,3, mp,dens,3,2600 ! 创建几何模型! 创建关键点 k,1 k,2,**h+*h k,3,*h,*h k,4,0,*h ! 创建线 l,1,2 l,2,3 l,3,4 l,4,1 ! 创建面 al,1,2,3,4 RECTNG,0,*h,*h,h ! 画坝顶矩形 RECTNG,-*h,0,-2*h,0 ! 坝基左面矩形 RECTNG,0,**h+*h,-2*h,0 ! 坝基正下方矩形 RECTNG,**h+*h,(2+*)*h,-2*h,0 ! 坝基右面矩形 RECTNG,-*h,(2+*)*h,-2*h,h ! 总矩形 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, amesh,4 eplot ! 划分 5 号面的网格 lsel,s,,,15,16,1 lesize,all,,,8,4 lsel,s,,,8 lesize,all,,,8, 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 ! 删除前面连接的线元素 allsel lsel,r,lcca ldele,all ! 拉伸成坝体单元, 采用 solid65 单元和 1 号材料, 随后改变 h3 以上坝体的材料属性 EXTOPT,ESIZE,8,0, ! 拉伸的份数 type,2 mat,1 VEXT,1,