简单重力坝ANSYS分析

简单混凝土坝的ANSYS分析计算

在上课和自习ANSYS之后，尝试用ANSYS建立一个简单的混凝土重力坝，并施加简单的约束和受力，最后尝试用ANSYS求解和后处理，查看大坝的变形，应力等情况，并做简要分析。

由于笔者本科初学ANSYS时老师教的是命令流方法建立模型也施加力，所以本篇全部采用命令流形式建立大坝和施力。

所建大坝喂混凝土重力坝，断面如下图，拔高180m,上游垂直，下游坝面洗漱0.75.坝基上游取270m，下游取360m，坝基深度取360m，坝顶长270m，坝顶宽1.8m，上游水位100m，下游水位80m。剖面图如下

一.首先，尝试对大坝的材料和材料类型形状等做出说明

/prep7

et,1,mesh200,6 ！划分用

et,2,solid65 ！混凝土单元

et,3,solid45 ！岩石单元

mp,ex,1,2.85e10 ！100m以下的材料弹性模量

mp,prxy,1,0.167 ！泊松比

mp,dens,1,2400 ！质量密度

mp,ex,2,2.6e10 ！100m以上的材料特性

mp,prxy,2,0.167

mp,dens,2,2400

mp,ex,3,2.9e10

mp,prxy,3,0.3 ！基岩特性mp,dens,3,2600

二：建立模型

k,1

k,2,139.5

k,3,18,162

k,4,0,162

a,1,2,3,4

rectng,0,18,162,180 rectng,-270,0,-360,0 rectng,0,139.5,-360,0 rectng,139.5,481.5,-360,0 rectng,-270,481.5,-360,180 aovlap,all

nummrg,all

numcmp,all

lsel,s,,,3,5,2

lesize,all,,,5

lsel,s,,,12,13,1

lesize,all,,,2

amesh,3

lsel,s,,,2,4,2

lesize,all,,,18

lsel,s,,,1

lesize,all,,,5

amesh,1

eplot

lsel,s,,,11

lesize,all,,,5

lsel,s,,,9,10,1

lesize,all,,,8,4

amesh,2

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

lsel,s,,,15,16,1 lesize,all,,,8,4 lsel,s,,,8 lesize,all,,,8,0.25 amesh,5

eplot

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

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 extopt,esize,8,0, type,2

mat,1

Vext,1,3,2,,,-135 extopt,esize,8,0, type,3

mat,3

Vext,2,,,,,-135 Vext,4,5,1,,,-135 allsel

extopt,esize,5,4, extopt,aclear,1 type,3

mat,3

Vext,1,7,1,,,180 local,11,0,,,-135 csys,11

/psymb,cs,1 dsys,11 vsymm,z,all nummrg,all nummcmp,all esel,s,mat,,1

eplot

nsle,s

nplot

nsel,r,loc,y,100,180 nplot

esln,s

eplot

mpchg,2,all

/solu

csys,0

dsys,0

nsel,s,loc,x,481.5 nsel,a,loc,x,-270 d,all,ux

nsel,s,loc,z,180

nsel,a,loc,z,-450

nplot

d,all,uz

nsel,s,loc,y,-360

d,all,uy

esel,s,type,,2

esel,s,type,,2

nsel,s,loc,x,0

nsel,r,loc,z,-270+0.1,-0.1 nsel,r,loc,y,0.1,100-5 esln,s

（模型图）至此，已经初步把模型建立并且划分了单元格，也施加了边界条件。再进行接下

来的步骤。

三．对建立好的模型施加荷载并求解

/psf,pres,norm,2,0,1 ！设定显示压力方向sfgrad,pres,0,y,0,-9.81 ！荷载梯度

sfe,all,2,pres,,9.81*100

allsel

esel,s,type,,2

nsle,s

fia2=36.9

local,12,0,139.5,,,fia2

csys,12

nsel,s,loc,y,0.1,80*1.25-25

nsel,r,loc,z,-270+0.1,-0.1

nsel,u,loc,x,-1000,-2

esln,s

/psf,pres,norm,2,0,1

sfgrad,pres,0,y,0,-9.81

sfe,all,4,pres,,9.81*80

csys,0

allsel

nsel,s,loc,y,0

nsel,r,loc,z,-135+0.1,0.1 ！选择

esln,s

esel,r,type,,2

sfcum,pres,add

sfe,all,5,pres,,9.81*80 ！下游浮托力

p0=9.81*100/139.5

sfgrad,pres,0,x,0,-p0 ！水压力荷载梯度

sfe,all,5,pres,,9.81*100

allsel

eplot

solve ！直接求解