【讨论】某水库大坝应力计算分析！

西北郎

水压方向不正确！

西北郎

原帖由 fxlt513 于 2006-8-12 10:42 发表
我也做过大坝的应力分析,右侧为上游,水压力向左,那么最大的拉应力应该是在大坝的坝踵处,即右边角点处,尤其是拱坝更明显,出现楼主的结果是不是梯度荷载的方向弄反了呢,以上不太成熟的建议,有错的地方,还请原谅

赞成这个说法，肯定是这样的！

sls

原帖由 hero110 于 2006-5-9 11:10 发表


nod

应该没有错。楼主画的是等效应力的等值线，如画y方向正应力就是右下角受拉左下角受压了。

renfeifan2008

这是我作的,不过我没有加重力，大家看看对不对？

why79303633

大坝很长时采用平面应变模型: 选用单元时候,选择beam,还是选用PLANE42单元（选择PLANE STRAIN）? 水工建筑物上介绍因力分析,对重力坝视为悬臂梁的.(水工建筑物-第三版,天津大学主编p59叶,说的)
重力(是不是等效一个法向荷栽)
水压力 按照变化的均布力加上的(press,)
传一个例子:ANSYS进行悬索结构的找形和计算
/prep7
ET,1,LINK10
!找形时采用很小的弹性模量
R,1,0.001468,0.9,
MP,EX,1,0.60551e9
MP,NUXY,1,0.3
!定义节点
*do,i,1,5
*do,j,1,2*i-1
k=(i-1)*(i-1)+j
n,k,(j-i)*9.15,-36.6+(i-1)*9.15,0
*enddo
*enddo
*do,i,1,4
*do,j,1,2*i-1
k=(i-1)*(i-1)+j+25
n,k,(j-i)*9.15,36.6-(i-1)*9.15,0
*enddo
*enddo
!定义单元
*do,i,1,4
*do,j,1,2*i
k=i*i+j
e,k,k+1
*enddo
*enddo
*do,i,1,3
*do,j,1,2*i
k=i*i+j+25
e,k,k+1
*enddo
*enddo
*do,i,1,4
*do,j,1,2*i-1
k=(i-1)*(i-1)+j
e,k,k+2*i
*enddo
*enddo
*do,i,1,3
*do,j,1,2*i-1
k=(i-1)*(i-1)+j+25
e,k,k+2*i
*enddo
*enddo
*do,i,18,24
e,i,i+17
*enddo
!施加位移约束
d,1,ux,0
d,1,uy,0
d,1,uz,3.66
d,2,ux,0
d,2,uy,0
d,2,uz,1.83
d,4,ux,0
d,4,uy,0
d,4,uz,1.83
d,5,ux,0
d,5,uy,0
d,5,uz,0
d,9,ux,0
d,9,uy,0
d,9,uz,0
d,10,ux,0
d,10,uy,0
d,10,uz,-1.83
d,16,ux,0
d,16,uy,0
d,16,uz,-1.83
d,17,ux,0
d,17,uy,0
d,17,uz,-3.66
d,25,ux,0
d,25,uy,0
d,25,uz,-3.66
d,35,ux,0
d,35,uy,0
d,35,uz,-1.83
d,41,ux,0
d,41,uy,0
d,41,uz,-1.83
d,30,ux,0
d,30,uy,0
d,30,uz,0
d,34,ux,0
d,34,uy,0
d,34,uz,0
d,27,ux,0
d,27,uy,0
d,27,uz,1.83
d,29,ux,0
d,29,uy,0
d,29,uz,1.83
d,26,ux,0
d,26,uy,0
d,26,uz,3.66
!设定时间步
TIME,1
AUTOTS,0
NSUBST,20, , ,1
KBC,0
!输出每个子步的结果
OUTPR,BASIC,ALL,
OUTRES,ALL,ALL,
!设定大变形，应力刚化
NLGEOM,1
SSTIF,ON
FINISH
!求截
/SOLU
/STATUS,SOLU
SOLVE
FINISH
!后处理
/POST1
PRNSOL,DOF,
PRESOL,SMISC,1
PLESOL,SMISC,1
PLNSOL,U,X,0,1
finish
/PREP7
!更新几何形状
!a.rst为计算结果文件名，最后一个为目录
!这两个参数应根据你的计算情况定
UPGEOM,1,LAST,LAST,a,rst,E:\JZD\1\
!弹性模量恢复为真值
R,1,0.001468,0.0027248,
MP,EX,1,2.0e11
!重新施加位移约束
d,1,uz,0
d,2,uz,0
d,4,uz,0
d,5,uz,0
d,9,uz,0
d,10,uz,0
d,16,uz,0
d,17,uz,0
d,25,uz,0
d,35,uz,0
d,41,uz,0
d,30,uz,0
d,34,uz,0
d,27,uz,0
d,29,uz,0
d,26,uz,0
!求截
/SOLU
/STATUS,SOLU
SOLVE
FINISH
/POST1
!观察施加荷载之前结构的位移
!可见位移基本为零,预张力基本不变
PLNSOL,U,Z,0,1
PLESOL,SMISC,1
FINISH
!写荷载工况文件01
/PREP7
LSWRITE,01,
!施加节点荷载
*do,i,1,41
f,i,fz,-167445
*enddo
!设置第2荷载步
TIME,2
AUTOTS,0
NSUBST,20, , ,1
KBC,0
LSWRITE,02,
!求解
FINISH
/SOLU
LSSOLVE,1,2,1,
后处理
/POST26
!可以看到，该结果和书中结果一致
NSOL,2,21,U,Z,
PLVAR,2, , , , , , , , , ,
ESOL,3,33,1,F,Y,
PLVAR,3, , , , , , , , , ,

yangxuan1984

原帖由 西王赏功 于 2003-6-2 16:02 发表
chenhu兄，你的模型中是否没有考虑 重力？还有，忘了说的是，在坝顶有一小段共2米高是没有加水压力的（水不能蔓淹到顶的）：） 我后来在没有计算重力作用时，并且把未加水压力的那块面清除掉了——得到的应力分 ...

老兄，国际单位下 ，c30的ex应该等于3e10 n/m*m吧？

showniya

呵呵，我对ANSYS还是一窍不通的水平，但看了lz的结果分析，我个人认为是合理的，在对坝基岩土稳定性分析里面有两个检验：抗滑移和抗倾覆。其中抗倾覆验算就是以坝趾（临空面底端）为支点，检验水压力在这一点的力矩和坝体自重在这一点的力矩的平衡关系。

xugangseu

搞了半天，楼主的命令也没贴上来

wing

重力坝计算不加重力，算出的结果肯定是不能看得

yuanfang910

这个例子做的好像还是初步的？

mumu_ld

找人代做论文。
结构硕士论文。模型已建好，属性，荷载等都已完成，可直接计算。由于时间关系，需找人帮我挑出结论并分析撰写部分论文。
mumu_ld1981@163.com

dgupc

最大应力才1.1MPa，会不会不太符合实际情况啊！

jiangyu5291169

1.1Mpa 很符合实际情况啊
大了估计混凝土抗压不够了
而拉应力是根本不给出现的
除非是在溢流坝的挑流出口，可以配筋抗拉

令狐易

都是高手啊，小弟刚上路，路过，顶一下

blackmore

分布肯定是一样的，结果肯定不一样！

zuo-0000

学习，谢谢楼主！

忠诚憨厚

谢谢了，很感激不禁

hran.zh

大坝是属于平面应变问题