水坝受力ansys模拟

合集下载

水坝受力ansys模拟

《大型结构分析软件》课程结业报告要求：1 课题来源和背景；2 模型基本情况和基本参数；3 分析主要步骤；4 后处理及结果分析；5 附录:本分析的精简的log 文件用小5号字；6报告正文以小四宋体，标题采用四号黑体加粗，有图片结果的需要保存为白底图。

分析题目ANSYS 重力坝受水压性能分析1、挡水建筑物用以拦截河流，形成水库等。

大型水利工程建设的首要目标是安全可靠，所以研究大坝的安全分析，评价等问题都需要解决，正确分析大坝性态急需解决。

用Ansys 软件可以很好的模拟大坝的受力等的性态对此解决大坝的安全问题。

2、（1）假设大坝的基础嵌入基岩中，且地基时刚性的。

（2）基本参数：弹性模量 E=35GPa 泊松比 v=0.2 容重 r=25KN/m3 重力加速度 g=9.8N/m (3) 计算分析大坝的水位为120米（4）水的质量密度为1000kg/m3 3、3.1对题目进行分析，画出简化模型，并对节点标号。

3.2设置解题类型启动ANSYS ，单击Preference,t 交互界面对话框，选中Structural ，来对后面的分析进行菜单及相应的图形界面过滤，单击OK 按钮。

3.3添加单元类型依次单击Main Menu ,Preprocessor ，Element Type ，Add/Edit/Delete ，弹出对话框，单击Add ，依次选择Soild ，Quad 4node4，单击OK 。

3.4设定单元选项：依次单击Main Menu ,Preprocesso ，Element Type Add/Edit/Delete 。

选中Type2PLANE42 ，单击Options.在Element behavior K3栏中选取Plane stress,其它为默认设置。

3.5添加材料属性依次单击Main Menu ,Preprocessor ，Material Props ，Material Models ，在对话框依次点击Structural ，Lingar,Elaastic,Isotropic ，在弹出对话框添加弹性模量3.5e10，泊松比0.2，单击OK 。

ansys重力坝建模

一、重力坝建模FINI/CLEAR,START/FILNAME,gravityDAM,1 !定义文件名/PLOPTS,DATE,0/TRIAD,LBOT/VIEW,1,1,1,1/prep7(1)定义单元材料及实常数et,1,plane42 !定义单元类型1 用于划分网格et,2,solid65 !定义单元类型2 钢筋混凝土，坝体et,3,solid45 !定义单元类型3 实体单元，基岩mp,ex,1,2.85e10 !定义材料1 弹性模量=2.85e10 mp,PRXY,1,0.167 !定义材料1 泊松比=0.167mp,dens,1,2400 !定义材料1 密度=2400mp,ex,2,2.9e10 !定义材料2 弹性模量=2.95e10 mp,PRXY,2,0.3 !定义材料2 泊松比=0.3mp,dens,2,2600 !定义材料2 密度=2600(2)画出坝体平面图形k,1,0,0 !设置关键点 1,2,3,4k,2,155 ,0k,3,20,180k,4,0,180l,1,2 !连接关键点l,2,3l,3,4l,4,1al,1,2,3,4 !由已知线生成面RECTNG,0,20,180,200 !画坝顶矩形RECTNG,-300,0,-400,0 !坝基左面矩形RECTNG,0,155,-400,0 !坝基正下方矩形RECTNG,155,535,-400,0 !坝基右面矩形RECTNG,-300,535,-400,200 !总矩形AOVLAP,all !将面单元进行粘贴布尔操作nummrg,all !合并重复元素numcmp,all(3)进行网格划分lsel,s,,,3,5,2 !选择线3,5lesize,all,,,5 !控制划分单元划分为5段lsel,s,,,12,13,1lesize,all,,,2amesh,3 !划分3号面的网格,坝顶矩形lsel,s,,,2,4,2 !同理划分1号面的网格，坝身矩形lesize,all,,,18lsel,s,,,1lesize,all,,,5amesh,1eplotlsel,s,,,11 !划分坝体正下方基岩网格lesize,all,,,5lsel,s,,,9,10,1lesize,all,,,8,4 !控制线9,10划分为8份，最后一份比第一份的比值为4 amesh,2eplotlsel,s,,,14 !同理划分4号面lesize,all,,,5,4lsel,s,,,7lesize,all,,,8,4lsel,s,,,6lesize,all,,,5,0.25amesh,4eplotlsel,s,,,15,16,1 !同理划分5号面lesize,all,,,8,4lsel,s,,,8lesize,all,,,8,0.25amesh,5eplotSAVElsel,s,,,4,13,9 !同理划分7号面LCCAT,all !合并线4,13lsel,s,,,19lesize,all,,,5,4lsel,s,,,20,21,1lesize,all,,,20amesh,7eplotlsel,s,,,2,12,10 !同理划分6号面lccat,alllsel,s,,,17lesize,all,,,8,4lsel,s,,,18,22,4lesize,all,,,20amesh,6EplotSAVEallsel ! 选择所有lsel,r,lcca !从中选择合并线ldele,all !删除选中的合并线(4)拉伸、镜像生成整个坝体模型EXTOPT,ESIZE,8,0, !控制拉伸的单元分数8type,2 !选择拉伸后单元类型mat,1 !选择拉伸后单元材料VEXT,1,3,2,,,-150 !将面1,3沿z轴负方向拉伸150EXTOPT,ESIZE,8,0, !同理拉伸2,4,5面生成坝基type,3mat,2VEXT,2,,,,,-150VEXT,4,5,1,,,-150allsel !选择所有面EXTOPT,ESIZE,5,4, !控制拉伸的份数为5，最后一份比第一份拉伸方向的尺寸比值为4EXTOPT,ACLEAR,1 !拉伸后删除源面上的网格type,3mat,2VEXT,1,7,1,,,200 !拉伸1~7面生成侧面基岩local,11,0,,,-150 !在z轴-150处建立标号为11的局部坐标系csys,11DSYS,11MainMenu>Preprocessor>Modeling>Reflect>Nodes >pick all !将全部节点沿x-y面镜像生成新节点，节点的标号增量为20000Menu>Preprocessor>Modeling>Reflect>Elements>User Numbered , pick all!在已生成的节点基础上镜像生成单元，单元标号增量为30000，节点增量为20000eplotnummrg,all !消除空号，使节点连续numcmp,all !合并重复项（5）施加约束及荷载csys,0dsys,0nsel,s,loc,x,535 !选择坝基x轴两侧面的节点施加约束nsel,a,loc,x,-300nplotMainmenu>Preprocessor>Loads>Define laods>Apply>Structural>Displacement>On nodes , pick all，选择UX， oknsel,s,loc,z,200 !同理施加z方向上的约束nsel,a,loc,z,-500nplotd,all,uznsel,s,loc,y,-400 !施加坝基底面的约束nplotd,all,uyallseleplotSAVEesel,s,type,,2 !选择单元类型为2的单元，即坝体eplotnsel,s,loc,x,0nsel,r,loc,z,-1.5*200+0.1,-0.1 !选择坝体垂直面上的节点nsel,r,loc,y,0.1,100-5!选择坝体受水压力高度的节点esln,s !选择已选节点上的单元nploteplot !显示单元/psf,pres,norm,2,0,1 !施加面荷载,压强, 垂直压强,打开轮廓线箭头线,打开载荷立即显示sfgrad,pres,0,y,0,-9810 !控制施加面荷载的梯度，沿y正方向的斜率为-9810sfe,all,2,pres,,9810*100 !在所选单元施加面荷载，9810*100allsel ! 选择施加下游面荷载ESEL,S,TYPE,,2 !选择坝体单元nsle,s !选择已选择单元上的节点eplot !显示单元local,12,0,155,,,90-53.1301024 !定义局部坐标系12，是一轴沿下游坝面向上csys,12/PSYMB,CS,1DSYS,12nsel,s,loc,y,0.1,80/sin(0.75)-25 ! 选择下游施加静水压力的节点、单元nsel,r,loc,z,-300+0.1,-0.1nsel,u,loc,x,-1000,-2esln,snploteplot/PSF,PRES,NORM,2,0,1 !施加下游静水压力sfgrad,pres,0,y,0,-9810sfe,all,4,pres,,9810*80/replotcsys,0 !选择坝体半边的底面单元dsys,0/psymb,cs,1allselnsel,s,loc,y,0nsel,r,loc,z,-150+0.1,0.1esln,sesel,r,type,,2eplotSFCUM,PRES,ADD !设置荷载为叠加荷载sfe,all,5,pres,,9810*80 !施加下游水位产生的浮托力p0=9810*100/(0.9*0.75*200+0.1*200) !施加上下游水位差产生的渗流压力sfgrad,pres,0,x,0,-p0sfe,all,5,pres,,9810*80/replotallsel !选择另外半边坝底单元施加扬压力nsel,s,loc,y,0nsel,r,loc,z,-300+0.1,-150-0.1esln,sesel,r,type,,2eplot*SET,p0,0sfgrad,pres,0,x,0,-p0sfe,all,3,pres,,9810*80*SET,p0,9810*80/(0.9*0.75*200+0.1*200) sfgrad,pres,0,x,0,-p0sfe,all,3,pres,,9810*100SFCUM,PRES,repl（6）求解及后处理AllselSolve。

拱坝ansys命令流

某砌石拱坝位于U型河谷中，坝高55.5m，为单曲等厚拱坝，顶宽5m，底宽16m，坝顶弧长115.65m，弧高比2.1。

首先输入基本参数，自底向上生成坝体。

!坐标系原点位于拱坝轴线与拱坝参考面的交点在坝基的投影，X轴指向右岸，Y轴指向下游，Z轴铅直向上!单位:长度-M,力-KN,加速度-M/S,密度-KN/M3!!FINI!退出当前处理器/CLEAR,START/FILNAME,ARCHDAM,1!指定新的工作文件名/PLOPTS,DATE,0/TRIAD,LBOT/VIEW,1,1,1,1/VUP,1,Z*AFUN,DEGZ_UP=53 !上游水位Z_DOWN=0 !下游水位Z_SAND=38 !泥沙淤积水位DENS_SAND=0.6 !淤积泥沙浮容重ANG_FRI=16 !淤积泥沙内摩擦角LAYER_NUM=8 !定义拱层数*DIM,ELEVATION,ARRAY,LAYER_NUM !定义拱圈高程*DIM,T_ARCH,ARRAY,LAYER_NUM !定义拱厚*DIM,ARCH_RAD,ARRAY,LAYER_NUM,2 !定义圆弧半径，上下游*DIM,RAD_CEN,ARRAY,LAYER_NUM !定义圆弧中心Y值*DIM,ARCH_ANGLE,ARRAY,LAYER_NUM,2 !定义圆弧左右中心角ELEVATION(1)=55.5,47,39,29,24,19,9,0 !拱圈高程H_DAM=ELEVATION(1)-ELEVATION(LAYER_NUM) !坝高T_ARCH(1)=5.00,6.68,8.27,10.25,11.24,12.23,14.22,16.00 !拱厚ARCH_RAD(1,1)=68.00,68.00,68.00,68.00,68.00,68.00,68.00,68.00 !外半径*VOPER,ARCH_RAD(1,2),ARCH_RAD(1,1),SUB,T_ARCH(1) !计算内半径RAD_CEN(1)=68.00,68.00,68.00,68.00,68.00,68.00,68.00,68.00 !圆心位置,Y坐标ARCH_ANGLE(1,1)=-47.72,-46.36,-45.06,-43.36,-41.88,-40.38,-37.32,-34.48 !左中心角ARCH_ANGLE(1,2)=49.73,47.13,44.68,41.57,40.44,39.29,36.92,35.72 !右中心角W_DAM1=-ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,1)) !坝顶左弦长W_DAM2=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,2)) !坝顶右弦长LOCAL,11,1,0,RAD_CEN(1),ELEVATION(1),-90 !局部柱面坐标,原点位于顶拱圆心,X轴为拱坝中心线,指向上游,Y轴为中心角!为查值计算坝体温度，重新定义一套数组*DIM,TEMPUP,TABLE,LAYER_NUM !上游面温度*DIM,TEMPDOWN,TABLE,LAYER_NUM !下游面温度*DIM,T2_ARCH,TABLE,LAYER_NUM !拱厚*DIM,Y_RAD,TABLE,LAYER_NUM !圆心到坝轴线距离*DIM,RADUP,TABLE,LAYER_NUM !外半径TEMPUP(1)=-8.84,-7.30,-5.66,-4.94,-4.99,-5.17,-5.73,-6.25TEMPUP(1,0)=0,6,14,24,29,34,44,53TEMPDOWN(1)=-8.84,-8.21,-7.36,-6.32,-5.86,-5.45,-4.72,-6.71 TEMPDOWN(1,0)=0,6,14,24,29,34,44,53T2_ARCH(1)=5.00,6.68,8.27,10.25,11.24,12.23,14.22,16.00T2_ARCH(1,0)=0,6,14,24,29,34,44,53Y_RAD(1)=68.00,68.00,68.00,68.00,68.00,68.00,68.00,68.00Y_RAD(1,0)=0,6,14,24,29,34,44,53RADUP(1)=68.00,68.00,68.00,68.00,68.00,68.00,68.00,68.00RADUP(1,0)=0,6,14,24,29,34,44,53CSYS,0!/PNUM,KP,1!/PNUM,LINE,1/PREP7KNN=0*DO,II,1,LAYER_NUM !生成拱坝控制关键点,每层6个点,上下游个三个点PX1=ARCH_RAD(II,1)*SIN(ARCH_ANGLE(II,1))PY1=RAD_CEN(II)-ARCH_RAD(II,1)*COS(ARCH_ANGLE(II,1))PX2=ARCH_RAD(II,1)*SIN(ARCH_ANGLE(II,2))PY2=RAD_CEN(II)-ARCH_RAD(II,1)*COS(ARCH_ANGLE(II,2))K,KNN+1,PX1, PY1,ELEVATION(II)K,KNN+2,0,RAD_CEN(II)-ARCH_RAD(II,1),ELEVATION(II)K,KNN+3,PX2, PY2,ELEVATION(II)PX3=ARCH_RAD(II,2)*SIN(ARCH_ANGLE(II,1))PY3=RAD_CEN(II)-ARCH_RAD(II,2)*COS(ARCH_ANGLE(II,1))PX4=ARCH_RAD(II,2)*SIN(ARCH_ANGLE(II,2))PY4=RAD_CEN(II)-ARCH_RAD(II,2)*COS(ARCH_ANGLE(II,2))K,KNN+4,PX3, PY3,ELEVATION(II)K,KNN+5,0,RAD_CEN(II)-ARCH_RAD(II,2),ELEVATION(II)K,KNN+6,PX4, PY4,ELEVATION(II)KNN=6*II*ENDDOLARC,1,3,2 !生成拱圈线*REPEAT,LAYER_NUM*2,3,3,3L,1,4*REPEAT,LAYER_NUM,6,6L,3,6*REPEAT,LAYER_NUM,6,6ASKIN,1,3,5,7,9,11,13,15 !生成上游面ASKIN,2,4,6,8,10,12,14,16 !生成下游面ASKIN,17,18,19,20,21,22,23,24 !生成侧面ASKIN,25,26,27,28,29,30,31,32AL,1,17,2,25AL,15,24,16,32VA,6,1,3,2,4,5==============================拉伸坝体侧面、底面，并进行切割，生成坝基!拉伸生成坝基ASEL,S,AREA,,4 !选择右侧面VEXT,ALL,,,4*H_DAM !拉伸成基岩ASEL,S,AREA,,3 !选择左侧面VEXT,ALL,,,-4*H_DAM !拉伸成基岩ALLSELASEL,S,LOC,Z,ELEVATION(LAYER_NUM)-1,ELEVATION(LAYER_NUM)+1 !选择底面VEXT,ALL,,,,,-NINT(H_DAM /10)*10ASEL,S,AREA,,10,20,5ASEL,A,AREA,,23,27,4VEXT,ALL,,,,4*H_DAMASEL,S,AREA,,8,18,5ASEL,A,AREA,,25,29,4VEXT,ALL,,,,-4*H_DAMALLSEL,ALLSAVEWPROTA,,90WPOFFS,,,NINT(H_DAM*1.5/10)*10VSBW,ALL,,DELETEVSEL,S,LOC,Y,-200,-H_DAM*1.4VDEL,ALL,,,1VSEL,ALLWPOFFS,,,-NINT(H_DAM*3.5/10)*10VSBW,ALL,,DELETEVSEL,S,LOC,Y,NINT(H_DAM*2/10)*10,NINT(H_DAM*2/10)*100VDEL,ALL,,,1VSEL,ALLWPCSYS,,0WPROTA,,,90WPOFFS,,,NINT((W_DAM2+H_DAM)/10)*10VSBW,ALL,,DELETEVSEL,S,LOC,X,NINT((W_DAM2+H_DAM)/10)*10,NINT((W_DAM2+H_DAM)/10)*50 VDEL,ALL,,,1VSEL,ALLWPOFFS,,,-NINT((W_DAM2+W_DAM1+2*H_DAM)/10)*10VSBW,ALL,,DELETEVSEL,S,LOC,X,-NINT((W_DAM1+H_DAM)/10)*50,-NINT((W_DAM1+H_DAM)/10)*10 VDEL,ALL,,,1WPCSYS,,0VSEL,ALLSAVE==============================用各控制高程切割整个模型，便于划分单元及计算成果和拱梁法比较*DO,II,LAYER_NUM,3,-1WPOFFS,,,ELEVATION(II-1)-ELEVATION(II)VSBW,ALL,,DELETE*ENDDOWPCSYS,,0ALLSELWPROTA,,,90VSBW,ALL,,DELETEWPCSYS,,0NUMMRG,ALLNUMCMP,ALLPX1=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,1))PX2=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,2))VSEL,S,LOC,X,PX1,PX2VSEL,R,LOC,Z,ELEVATION(LAYER_NUM),ELEVATION(1)CM,DAM,VOLUVSEL,INVECM,BEDROCK,VOLUALLSEL==============================生成溢流堰CMSEL,S,DAMVSEL,R,LOC,Z,ELEVATION(1)-0.1,ELEVATION(2)+0.1/VIEW,1,1,1,1VPLOT !生成溢流堰,为便于建模和施加荷载,与实际不太相符,根据圣维南原理,对大坝的总体应力不会有太大影响W_WEIR=8 !溢流堰单孔宽(坝轴线处)NUM_WEIR=5 !溢流堰孔数W_WALL=2.5 !溢流堰闸墩宽(坝轴线处)WPOFFS,,,ELEVATION(1)-ELEVATION(LAYER_NUM)WPOFFS,,RAD_CEN(1)WPROTA,,,90*AFUN,RADPI=2*ACOS(0)SITTA1=(W_WEIR*NUM_WEIR+W_WALL*(NUM_WEIR+1))/68*180/PI !溢流堰总中心角SITTA2=W_WEIR/ARCH_RAD(1,1)*180/PI !溢流堰单孔中心角SITTA3=W_WALL/ARCH_RAD(1,1)*180/PI !溢流堰单个闸墩中心角WPROTA,,-SITTA1/2VSBW,ALL,,DELETE*DO,II,1,NUM_WEIRWPROTA,,SITTA3VSBW,ALL,,DELETEWPROTA,,SITTA2VSBW,ALL,,DELETE*ENDDOWPROTA,,SITTA3VSBW,ALL,,DELETESAVECMSEL,S,BEDROCKVSEL,INVEVSEL,R,LOC,Z,ELEVATION(1),ELEVATION(2)CSYS,11VSEL,R,LOC,Y,-SITTA1/2,SITTA1/2CM,YLYY,VOLU*DO,II,1,NUM_WEIRCMSEL,S,YLYYVSEL,R,LOC,Y,-SITTA1/2+II*SITTA3+(II-1)*SITTA2,-SITTA1/2+II*SITTA3+II*SITTA2 VDEL,ALL,,,1,,,1*ENDDOCSYS,0CMSEL,S,BEDROCKVSEL,INVECM,DAM,VOLUCMSEL,S,DAMVSEL,R,LOC,Z,ELEVATION(1),ELEVATION(2)CSYS,11VSEL,R,LOC,Y,-SITTA1/2,SITTA1/2,CM,ZHADUN,VOLUCSYS,0WPCSYS,,0VSEL,ALLNUMMRG,ALLNUMCMP,ALLSAVE==============================划分单元!定义单元性质和材料ET,1,SOLID95MP,EX,1,1E7MP,NUXY,1,0.25MP,DENS,1,23 !坝体容重MP,ALPX,1,0.7E-5MP,REFT,1,0MP,EX,2,0.8E7MP,NUXY,2,0.21MP,ALPX,2,0.0MP,REFT,2,0!坝体单元划分CMSEL,S,DAMCMSEL,U,ZHADUNASLV,SLSLA,S*AFUN,DEGPX1=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,1))PX2=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,2))PX3=ARCH_RAD(LAYER_NUM,2)*SIN(ARCH_ANGLE(LAYER_NUM,1)) PX4=ARCH_RAD(LAYER_NUM,2)*SIN(ARCH_ANGLE(LAYER_NUM,2)) LSEL,U,LOC,X,PX3,PX1LSEL,U,LOC,X,PX4,PX2LSEL,U,LOC,X,0LSEL,U,LOC,Z,ELEVATION(1),ELEVATION(2)LESIZE,ALL, , ,20,CM,LTEMP1,LINECSYS,0CMSEL,S,DAMASLV,SLSLA,SCMSEL,U,LTEMP1CM,LTEMP2,LINELSEL,U,LENGTH,,T_ARCH(LAYER_NUM)+0.1,1000LSEL,R,TAN1,ZCSYS,11LSEL,U,LOC,X,ARCH_RAD(1,1)LSEL,U,LOC,X,ARCH_RAD(1,1)-T_ARCH(1)LSEL,U,LOC,X,ARCH_RAD(1,1)-T_ARCH(2)CSYS,0CM,LTEMP3,LINELESIZE,ALL, , ,5,ALLSELCMSEL,S,LTEMP2CMSEL,U,LTEMP3LSEL,U,TAN1,ZCM,LTEMP4,LINELESIZE,ALL, , ,3,CMSEL,S,DAMCMSEL,U,ZHADUNASLV,SLSLA,SLSEL,R,LOC,Z,ELEVATION(1),ELEVATION(2)CSYS,11LSEL,U,LOC,Y,-SITTA1/2-1,SITTA1/2+1LSEL,U,LENGTH,,0,T_ARCH(1)+0.1CSYS,0CMSEL,U,LTEMP3CMSEL,U,LTEMP4LESIZE,ALL, , ,10,CM,LTEMP5,LINECMSEL,S,DAMCMSEL,U,ZHADUNASLV,SLSLA,SLSEL,R,LOC,Z,ELEVATION(1),ELEVATION(2)CMSEL,U,LTEMP3CMSEL,U,LTEMP4CMSEL,U,LTEMP5LSEL,U,LENGTH,,0,3CM,LTEMP6,LINELESIZE,ALL, , ,4,CSYS,0ALLSELMSHAPE,0CMSEL,S,DAMVSWEEP,ALLSAVE!坝基单元划分CMSEL,S,BEDROCKASLV,SLSLA,SLSEL,U,TAN1,XLSEL,U,LOC,X,PX2,PX1CM,LTEMP7,LINELESIZE,ALL, , ,5,, , , ,1CMSEL,S,BEDROCKASLV,SLSLA,SLSEL,U,TAN1,YPY4=RAD_CEN(LAYER_NUM)-ARCH_RAD(LAYER_NUM,2)*COS(ARCH_ANGLE(LAYER_NUM,2))+20 LSEL,U,LOC,Y,-10,PY4CM,LTEMP8,LINELESIZE,ALL, , ,6,4, , , ,1CMSEL,S,BEDROCKASLV,SLSLA,SLSEL,U,TAN1,ZLSEL,R,LOC,Z,ELEVATION(LAYER_NUM),ELEVATION(LAYER_NUM)-200LESIZE,ALL, , ,5,5, , , ,1CMSEL,S,BEDROCKASLV,SLSLA,SLSEL,U,TAN1,ZLSEL,R,LOC,Z,ELEVATION(1),ELEVATION(LAYER_NUM)LESIZE,ALL, , ,3,CMSEL,S,BEDROCKASLV,SLSLA,SLSEL,R,LOC,Z,-50,-1000LSEL,R,LOC,Y,-10,PY4CSYS,11LSEL,U,LOC,X,ARCH_RAD(1,1)LSEL,U,LOC,X,ARCH_RAD(1,1)-T_ARCH(LAYER_NUM)CSYS,0LSEL,R,LENGTH,,T_ARCH(LAYER_NUM)-5,T_ARCH(LAYER_NUM)+3,LESIZE,ALL, , ,5,CSYS,0CMSEL,S,BEDROCKASLV,SLSLA,SLSEL,R,LOC,Z,ELEVATION(LAYER_NUM)-50,ELEVATION(LAYER_NUM)-1000LSEL,R,LOC,X,PX2,PX1LESIZE,ALL, , ,20,CMSEL,S,BEDROCKMAT,2VMESH,ALL!单元划分完毕ALLSELNUMMRG,ALLNUMCMP,ALLSAVE==============================定义边界条件，顺河向约束Y向位移，横河向约束X向位移，底面约束竖向Z位移FINI/SOLU!边界位移条件*GET,NXMAX,NODE,,MXLOC,X*GET,NXMIN,NODE,,MNLOC,X*GET,NYMAX,NODE,,MXLOC,Y*GET,NYMIN,NODE,,MNLOC,Y*GET,NZMIN,NODE,,MNLOC,ZNSEL,S,LOC,X,NXMAXNSEL,A,LOC,X,NXMIND,ALL,UX,0 !约束上下游边界X向位移NSEL,S,LOC,Y,NYMAXNSEL,A,LOC,Y,NYMIND,ALL,UY,0 !约束左右边界Y向位移NSEL,S,LOC,Z,NZMIND,ALL,,ALL,0 !约束底面边界Z向位移NSEL,ALL==============================施加水荷载，淤沙压力，坝体自重!大坝荷载ALLSELCMSEL,S,DAMCMSEL,U,ZHADUNASLV,SCSYS,11ASEL,R,LOC,X,ARCH_RAD(1,1)-0.1,ARCH_RAD(1,1)+0.1 NSLA,S,1CM,N_DAMUP,NODE !定义大坝上游面节点CSYS,0ALLSELCMSEL,S,DAMASLV,SASEL,R,EXTLSLA,SLSEL,R,LOC,X,0LSEL,U,LOC,Y,-0.1,0.1LSEL,U,LOC,Z,ELEVATION(2)LSEL,U,LOC,Z,ELEVATION(LAYER_NUM)ASLL,SASEL,U,LOC,X,0LSLA,SNSLA,S,1CM,N_DAMDOWN,NODE !定义大坝下游面节点ALLSELCMSEL,S,YLYYASLV,SCSYS,11ASEL,R,LOC,X,ARCH_RAD(1,1)-0.1,ARCH_RAD(1,1)+0.1 ASEL,U,,,291,326,326-291NSLA,S,1CM,N_ZHADUN1,NODE !定义溢流坝中墩节点ASEL,S,,,291,326,326-291NSLA,S,1CM,N_ZHADUN2,NODE !定义溢流坝边墩节点CSYS,0*IF,Z_UP,NE,ELEVATION(LAYER_NUM),THENCMSEL,S,N_DAMUPNSEL,R,LOC,Z,0,Z_UP/PSF,PRES,NORM,2,0.1SFGRAD, PRES, 0, Z, 0,-9.81SF,ALL,PRESS,Z_UP*9.81 !施加大坝上游面水压力*ENDIF*IF,Z_UP,GT,ELEVATION(2),THENCMSEL,S,N_ZHADUN1NSEL,R,LOC,Z,ELEVATION(2),Z_UPSFGRAD, PRES, 0, Z, 0,-(W_WEIR+W_WALL)/W_WALL*9.81SF,ALL,PRESS,Z_UP*(W_WEIR+W_WALL)/W_WALL*9.81 !施加中墩水压力,包含闸门传递的水压力CMSEL,S,N_ZHADUN2NSEL,R,LOC,Z,ELEVATION(2),Z_UPSFGRAD, PRES, 0, Z, 0,-(W_WEIR/2+W_WALL)/W_WALL*9.81SF,ALL,PRESS,Z_UP*(W_WEIR/2+W_WALL)/W_WALL*9.81 !施加边墩水压力,包含闸门传递的水压力*ENDIF*IF,Z_SAND,GT,ELEVATION(LAYER_NUM),THENSFCUM,PRES,ADD !荷载为累加方式CMSEL,S,N_DAMUPNSEL,R,LOC,Z,0,Z_SAND*AFUN,DEGSAND_GRADS=DENS_SAND*(TAN(45-ANG_FRI/2))**2*9.81 !计算淤沙压力梯度SFGRAD, PRES, 0, Z, 0,-SAND_GRADSSF,ALL,PRESS,Z_SAND*SAND_GRADS !施加淤沙压力SFCUM,PRES,REPL !荷载为替代方式*ENDIFALLSEL*IF,Z_DOWN,NE,ELEVATION(LAYER_NUM),THENCMSEL,S,N_DAMDOWNNSEL,R,LOC,Z,0,Z_DOWNSFGRAD, PRES, 0, Z, 0,-9.81SF,ALL,PRESS,Z_DOWN*9.81 !施加大坝下游面水压力*ENDIFALLSELACEL,0,0,1 !施加大坝自重荷载SAVE==============================施加温度体荷载!施加大坝温度荷载CMSEL,S,DAMNSLV,S,1*GET,NCOUNT_DAM, NODE,, COUNT*GET,NMIN_DAM, NODE,,NUM,MINTUNIF,0,TREF,0,NTT=NMIN_DAM*DO,II,1,NCOUNT_DAMZ_NTT=NZ(NTT) !ZZ_TEMP=ELEVATION(1)-Z_NTT !到坝基高度TL_ARCH=T2_ARCH(Z_TEMP) !拱厚CSYS,11X_NTT=NX(NTT) !X:距顶拱圆心距离CSYS,0RAD0=Y_RAD(Z_TEMP) !圆心到坝轴线距离RAD1=RADUP(Z_TEMP) !外半径LXUP=RAD_CEN(1)-RAD0+RAD1-X_NTTTTT1=TEMPUP(Z_TEMP)TTT2=TEMPDOWN(Z_TEMP)TTT=TTT1+(TTT2-TTT1)/TL_ARCH*LXUPBF,NTT,TEMP,TTT*IF,II,LT,NCOUNT_DAM,THENNTT=NDNEXT(NTT)*ENDIF*ENDDO==============================计算求解，采用PCG求解器ALLSELSAVEEQSLV,PCG,1E-8SOLVESAVEFINI最大变形2cm，位于拱冠梁顶，与拱梁分载法一致。

利用ANSYS 模拟大坝的运行期应力变形

利用ANSYS模拟大坝的运行期应力变形1.主要步骤：
S1：在Auto CAD里面按照给出尺寸绘图并导入到ANSYS软件中。

S2：通过命令流或者GUI操作来设置材料并给其赋值。

S3：通过APDL命令或者建模的“Operate”里面的“Extrude”沿Z轴将面拉成体，注意距离为10。

S4：按照题目要求，通过APDL或者“Meshing”里面的“Mesh Attributes”给体赋值材料号。

S5：划分网格，先将线分段再划分网格，注意在给侧面划分网格时划分份数为2，并将“KYNDIV SIZE,NDIV can be changed”勾选去掉，否则划分网格时划分份数会有变化。

S6：通过APDL或者“Solution”里面的“Define Loads”添加边界条件和自重。

S7：模拟初始状态，“杀死”坝体单元，求解。

S8：模拟完建状态，“激活”坝体状态，求解。

S9：运行期，施加水荷载，计算扬压力，求解运行期。

S10：创建荷载工况，查看结果。

2.导出结果：2.1初始状态模拟
2.2完建状态模拟
2.3运行期模拟
2.4运行期的变形图
2.5运行期的第一应力图
2.6运行期的第一应力等值线图
2.7运行期总位移变形图
2.8运行期总位移变形等值线图。

ANSYS水坝分析实例

如图为某水坝的截面及尺寸图，坝体为混凝土浇筑，材料的弹性10Pa,泊松比为0.25，密度为2500kg/3m,坝体挡水面承模量为2.14×10受静水压力作用，假设最危险的水平面刚好平坝底，计算重力和水压力作用下的坝体的承载状态，求解整个坝体的位移图及应力云图。

1、指定工作文件名“Utility Menu>File>Change Jobname”2、定义标题名“Utility Menu>File>Change Title”3定义单元“Preprocessor>Element Type >Add/Edit/Delete”4、定义弹性模量和泊松比“Preprocessor>Material Props >Material Models>Structural>Linear>Elastic>Isotropic”5、定义密度“Preprocessor>Material Props >Material Models>Structural> Density”6、创建关键点K1, “Preprocessor>Modeling>Create>Keypoints>In Active CS”分别定义K2、K,3、K,4、K,5、K67、依次连接各关键点“Preprocessor>Modeling>Create >Lines >Lines>StraightLine”8、用线合成面积“Preprocessor>Modeling>Create >Areas>Arbitrary>By Lines”9、定义关键点K7、K8，并分别与K3、K6连接，得到现L7、L810、用线L7、L8将A1切割“Preprocessor>Modeling>Operate>Blooleans>Divide>Area By Line”打开网格划分“Preprocessor>Meshing>MeshTool”11、将L9平均分割为8份12、同理将L7平均分割为12份,将L8平均分割为34份,将L2平均分割为26份13、采用映射划分网格划分A1 、A214得到A1 、A2的网格划分图15、采用线连接“Preprocessor>Meshing>Concatenate>Lines”16、将L3 、L7 连接成一条直线17、将L5 、L8 连接成一条直线18、将L7平均分割为12份、将L3平均分割为28份、将L8平均分割为34份，将L5平均分割为6份、将L4平均分割为6份19、采用映射划分网格划分A420、在底线L9、L11、L2上施加约束“Preprecessor>Loads>DefineLoads>Apply>Structuaral>Displacemengt>On Lines”Loads>Apply>Structural>Pressure>On Lines”Loads>Apply>Structural>Inertia>Gravity”23、进行求解“Solution>Solve>Current LS”Shape”24、绘制云图“General Postproc>Plot Results>Contour Plot>Nodal Solu”选择显示方式X方向节点云图X方向应力云图APDL命令：/PREP7 !进入前处理器ET,1,PLANE42 !定义单元体MP,EX,1,2.14e10 !定义弹性模量MP,NUXY,1,0.25 !定义泊松比MP,DENS,1,2500 !定义密度K,1,-20,0,0K,2,80,0,0K,3,0,30,0K,4,0,100,0K,5,15,100,0K,6,15,85,0 !创建关键点LSTR,1,3LSTR,1,2LSTR,3,4LSTR,4,5LSTR,5,6LSTR,6,2 !连点成线AL,1,2,3,4,5,6 !连线成面K,7,0,0,0K,8,15,0,0 !创建关键点LSTR,3,7LSTR,8,6 !连点成线ASBL,1,7ASBL,3,8 !切割图形ALLSEL,ALL !选择全部线LSEL,S,,,9 !选择L9 LESIZE,9,,,8 !将L9平均分成8份LSEL,S,,,7LESIZE,7,,,12ASEL,S,AREA,,2 !选择A2AATT,2,,1MSHKEY,1 !选择映射划分方式MSHAPE,0,2D !选择2D四边形划分形式AMESH,2 !划分网格ALLSEL,ALLLSEL,S,,,2LESIZE,2,,,26LSEL,S,,,8LESIZE,8,,,34ASEL,S,AREA,,1AATT,1,,1MSHKEY,1MSHAPE,0,2DAMESH,1ALLSEL,ALLLSEL,S,,,3,7LCCAT,ALLALLSEL,ALLLSEL,S,,,5,8LCCAT,ALL !连接ALLSEL,ALLLSEL,S,,,4LESIZE,2,,,6LSEL,S,,,3LESIZE,8,,,28LSEL,S,,,7LESIZE,8,,,12ASEL,S,AREA,,4AATT,4,,1MSHKEY,1MSHAPE,0,2DAMESH,4ALLSEL,ALLSFL,1, ,PRES,10000SFL,3, ,PRES,9000 !对两边施加分布荷载ACEL,,-9.8 !施加重力荷载DL,2,,ALL,0DL,9,,ALL,0DL,11,,ALL,0 !对底线进行约束ALLSEL,ALL/SOLU !进入求解SOLVE !求解当前荷载FINISH/POST1 ! 进入后处理器PLDISP,2 !绘制变形图PLNSOL,U,X !绘制X方向节点云位移图PLNSOL,U,Y !绘制Y方向节点云位移图PLNSOL,S,X !绘制X方向应力云图PLNSOL,S,Y !绘制Y方向应力云图。

基于Ansys对于坝体的三种工况研究讲解

基于Ansys对于坝体的研究分析报告坝体及相关建筑在使用过程中，会承受如重力、净水压力、淤泥荷载、浪压力、扬压力等各种作用，而我们在设计、建造这个建筑之前，要分析其产生的应力、应变进而选取材料和校核材料的安全性。

为分析所需，基于Ansys软件建立相应的模型，并施加荷载和作用，在三种工况下校核结构的安全性。

一：分析对象1：坝体的几何参数：2：基岩的几何参数：二：作用及荷载（1）约束基岩左右两端受x 方向的位移约束，基岩下端受x 、y 两个方向的位移约束。

（2）静水压力正常蓄水位高程91.75m ，防洪高水位97m ，校核洪水位101m 。

对应下游水位分别为15m ，20m 和25m 。

（3）泥沙荷载坝前泥沙淤积高程： 25m 。

坝前泥沙浮容重：6.0kN/m 3，淤沙内摩擦角：12°。

坝面上单位宽度上的泥沙压力为221(45)22s sk sb s p h tg ϕγ=︒- 式中： sk p ——淤沙压力标准值（KN/m ）；sb γ——泥沙的浮容重，取6kN/m 3；s h ——泥沙淤积深度（m ）；s ϕ——淤沙的内摩擦角，12°。

（4）浪压力50年一遇计算风速21m/s ，多年平均最大风速14m/s ，有效吹程1km 。

（5）扬压力取渗透压力强度系数α＝0.25，帷幕中心线坐标X=10m 。

三：选用单元及划分网格1) 单元选择：Solid –Quad 4node 422) 材料参数：坝体和基岩分别设置，见上图。

3) 划分网格：坝体部分-外围线按1m 每格划分，整体按自由网格划分。

基岩部分-靠近坝体网格密集，坝基面水平线上基岩外围线按20份、4的比率划分；垂直地表的按20份、0.25的比率划分；基底均分。

整体基岩自由网格划分。

四：三种工况的具体Ansys设置1)正常蓄水位（其中括号内为承载能力极限状态时的分项系数）上游水位高程为91.75m，下游水位高程为15m。

（1）上下游静水压力（分项系数为1.0）gradient 斜率为-9810 沿y轴方向，分别取91.75m和15m在各自位置。

ansys大坝受力分析

ansys水坝坝体受力分析图3-1 水坝截面υ操作步骤1. 清除内存，准备分析1) 清除内存，开始一个新的分析，选择菜单路径Utility Menu>File>Clear & Start New 弹出Clears database and Start New对话框，采用默认状态，单击ok按钮弹出Verify确认对话框，单击Yes按钮。

2) 指定工作文件名。

选取菜单路径Utility Menu>File>Change Jobname 弹出『Change Jobname』对话框，在enter New Jobname项中输入“PlaneStrain”，然后单击ok按钮。

3) 定义标题，选择菜单路径Utility Menu>File>Change Title，输入文字“The PlaneStrain Pressure Gradient Example：Dam Section”，单击OK按钮。

2. 创建有限元模型1) 进入前处理器，并定义单元类型1，选取菜单路径MainMenu>Preprocessor>Element Type>Add/Edit/Delete弹出『Element Types』对话框，单击Add弹出如图3-2所示『Library of Element Types』对话框，设置下列选项：图3-2 『Library of Element Types』对话框λ左边列表框中选择Solid。

λ右边列表框中选择Quad4node42。

λ Element type reference number：单元编号，输入1。

单击ok按钮，返回Element对话框，单击Option按钮，弹出如图3-3所示『PLANE42 element type option』设置对话框，将K3设置成Plane Strain，单击OK按钮返回Element Type对话框，单击Close按钮可关闭该对话框。

ansys大坝受力分析

ansys水坝坝体受力分析图3-1水坝截面: 操作步骤1. 清除内存，准备分析1) 清除内存，开始一个新的分析，选择菜单路径Utility Me nu>File>Clear & Start New 弹出Clears database and Start New 对话框，采用默认状态，单击ok按钮弹出Verify确认对话框，单击Yes按钮。

2) 指定工作文件名。

选取菜单路径Utility Men u>File>Cha nge Job name 弹出『Cha nge Job name』对话框，在enter New Job name 项中输入 "Pla neStrai n，"然后单击ok按钮。

3) 定义标题，选择菜单路径Utility Men u>File>Cha nge Title ，输入文字“ The Pla ne Strain Pressure Gradient Example : Dam Section ”,单击OK 按钮。

2. 创建有限元模型1) 进入前处理器，并定义单元类型1，选取菜单路径MainMenu >Preprocessor>Eleme nt Type>Add/Edit/Delete弹出『Element Types』对话框，单击Add弹出如图3-2所示『Library of ElementTypes』对话框，设置下列选项：图3-2 『Library of Element Types 』对话框■■■■■左边列表框中选择Solid。

- 右边列表框中选择Quad4node42 。

- Eleme nt type reference nu mber :单元编号，输入1。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

分析题目ANSYS 重力坝受水压性能分析1、挡水建筑物用以拦截河流，形成水库等。

大型水利工程建设的首要目标是安全可靠，所以研究大坝的安全分析，评价等问题都需要解决，正确分析大坝性态急需解决。

用Ansys 软件可以很好的模拟大坝的受力等的性态对此解决大坝的安全问题。

2、（1）假设大坝的基础嵌入基岩中，且地基时刚性的。

3.2设置解题类型启动ANSYS ，单击Preference,t 交互界面对话框，选中Structural ，来对后面的分析进行菜单及相应的图形界面过滤，单击OK 按钮。

3.3添加单元类型依次单击Main Menu ,Preprocessor ，Element Type ，Add/Edit/Delete ，弹出对话框，单击Add ，依次选择Soild ，Quad 4node4，单击OK 。

3.4设定单元选项：依次单击Main Menu ,Preprocesso ，Element Type Add/Edit/Delete 。

选中Type2PLANE42 ，单击Options.在Element behavior K3栏中选取Plane stress,其它为默认设置。

再选中Density ，在Density for Material Number1对话框，在DENS 中输入边坡土体材料的密度2500，单击OK 。

3.6建模和划分网格3.6.1创建节点依次单击Main Menu ,Preprocessor ，Modeling ，Create ，Keypoints ，In Active CS ，弹出对话框，依次输入节点编号，节点坐标（0，0,0），单击Apply ，输入节点编号2,3,4,4，节点坐标（76,0,0），(15.6,104.1,0),(15.6,120,0),(5.6,120,0)重复以上过程，输入其它节点的编号和坐标单击OK 。

3.6.2创建坝体线模型依次单击Main Menu ，Preprocessor ，Modeling ，Create ，Lines ，Straight line ，弹出对话框连,连接关键点1、2单击Apply ，然后依次连接23，3、4，4、5，5、1。

单击OK ，得到坝体线模型。

3.7创建坝体面模型： 3.7.1关键点标号显示：依次单击Utility Menu ,PiotCtrls,Numbering.选中Keypoint Numbers,将off 变为on 。

3.7.2创建坝体面模型：依次单击Main Menu ,Preprocessor ，Modeling ，Create ，Areas, Arbitrary, Through KPs. 选中关键点1,2,3,4,5单击Apply 即可。

3.8 划分网格3.8.1设置网格份数：Main Menu ,Preprocessor ，Meshing,Size Ctrls,ManualSize, Layers, Picked Lines.选取线L1单击OK ，在No.of line divivions 栏后面输入20单击OK 。

相同设置L2、L3、L4、L5的分割份数为32,6,4，40。

3.8.2划分单元网格：Main Menu ,Preprocessor ，Meshing ，Mesh ，Area ，Free ，拾取图形。

单击OK 。

得到如图3.9 施加载荷及边界条件 3.9.1施加重力加速度依次单击Main Menu ，Preprocessor ，Loads ，Define Loads ，Apply ，Structural ，Inertia ,Gravity ，弹出对话框。

在Global Cartesian Y-comp 栏中输入9.8。

3.9.2施加水压力：依次单击Main Menu ，Solution,Define Loads,Apply,Structural,Pressure,on Lines.选中L5，单击OK ，在对话框中输入0，1101370单击OK3.9.3坝体底部施加位移约束：依次单击Main Menu ，Preprocessor ，Loads ，Define Loads ，Apply ，Structural ，Displacement,on Nodes.在DOFs to be constrained 栏中远ALL DOF 。

在Ｄisplacement value 中输入0单击OK 。

3.10求解及后处理3.10.1指定求解类型3.10.2设置载荷步3.10.3设置线性搜索 3.10.4静力求解：依次单击Main Menu, Solution,Slove, Current LS 单击OK.等待至Solution is done 求解结束。

3.10.5后处理院系：理学院班级：理论与应用力学07—1 姓名：汤仁君学号： 2007135461．绘制坝体变形图：依次单击Main Menu, General Postproc, Plot Results, Defomed Shap.选中Def+undeformed出现下图2坝体X方向位移图：依次单击Main Menu, General Postproc, Piot Results, Contour Plot, Nodal solu.再依次点击Nodal Solution, DOF Soiution , X-Component of displacement.如下图3坝体X方向应力云图：依次单击Main Menu, General Postproc, Piot Results, Contour Plot, Nodal solu. 依次单击Nodal Solution,，Stress，X-Component of Stress如下图4：附log 文件/COM,ANSYS RELEASE 12.1 UP20091102 17:23:00 11/04/2010 /input,menust,tmp,'',,,,,,,,,,,,,,,,1/GRA,POWER/GST,ON/PLO,INFO,3/GRO,CURL,ON/CPLANE,1/REPLOT,RESIZEWPSTYLE,,,,,,,,0/CWD,'F:\听力'/FILNAME,ds,0/TITLE,ds/PREP7ET,1,PLANE42ET,2,PLANE42KEYOPT,2,1,0KEYOPT,2,2,0KEYOPT,2,3,0KEYOPT,2,5,0KEYOPT,2,6,0MPTEMP,,,,,,,, MPTEMP,1,0MPDATA,EX,1,,3.5e10 MPDATA,PRXY,1,,0.2 MPTEMP,,,,,,,, MPTEMP,1,0MPDATA,DENS,1,,2500K,1,0,0,0,K,2,76,0,0,K,3,15.6,104.1,0,K,4,15.6,120,0,K,5,5.6,120,0,LSTR, 1, 2 LSTR, 2, 3 LSTR, 3, 4 LSTR, 4, 5 LSTR, 5, 1 /PNUM,KP,1 /PNUM,LINE,0 /PNUM,AREA,0 /PNUM,VOLU,0 /PNUM,NODE,0 /PNUM,TABN,0 /PNUM,SVAL,0 /NUMBER,0/PNUM,ELEM,0 /REPLOTFLST,2,5,3 FITEM,2,1 FITEM,2,2 FITEM,2,3 FITEM,2,4 FITEM,2,5A,P51XFLST,5,1,4,ORDE,1 FITEM,5,1CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,S,_YLESIZE,_Y1,0, ,20,0,4,0,0 CMDELE,_YCMDELE,_Y1FLST,5,1,4,ORDE,1 FITEM,5,2CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,S,_YLESIZE,_Y1,0, ,32,0,4,0,0 CMDELE,_YCMDELE,_Y1FLST,5,1,4,ORDE,1 FITEM,5,3CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,S,_Y LESIZE,_Y1,0, ,6,0,4,0,0CMDELE,_YCMDELE,_Y1FLST,5,1,4,ORDE,1FITEM,5,4CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,S,_YLESIZE,_Y1,0, ,4,0,4,0,0CMDELE,_YCMDELE,_Y1FLST,5,1,4,ORDE,1FITEM,5,5CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,S,_YLESIZE,_Y1,0, ,40,0,4,0,0CMDELE,_YCMDELE,_Y1MSHKEY,0CM,_Y,AREAASEL, , , , 1CM,_Y1,AREACHKMSH,'AREA'CMSEL,S,_YAMESH,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2SAVEFINISH/SOLFLST,2,80,1,ORDE,6FITEM,2,1FITEM,2,-2FITEM,2,22FITEM,2,-60FITEM,2,64FITEM,2,-102/GOD,P51X, ,0, , , ,ALL, , , , ,ACEL,0,9.8,0,FLST,2,1,4,ORDE,1FITEM,2,5/GOSFL,P51X,PRES,0,1101370ANTYPE,0NSUBST,5,100,1CUTCONTROL,CRPLIMITexp,0.1,0CUTCONTROL,DSPLIMIT,10000000CUTCONTROL,PLSLIMIT,0.15LNSRCH,1/STATUS,SOLUSOLVEFINISH/POST1SET,LISTPLDISP,1/EFACET,1PLNSOL, U,X, 0,1.0/EFACET,1PLNSOL, U,Y, 0,1.0/EFACET,1PLNSOL, S,X, 0,1.0/EFACET,1PLNSOL, S,Y, 0,1.0FINISH! /EXIT,ALL。