某大坝的应力分析（附命令流）

chenfanjian

某大坝的应力分析（附命令流)
已知条件：
1、  大坝截面呈直角三角形,其中迎水面为垂直面，高100m，底边宽80 m。
2、  此题仅考虑自身的重力、上游水压力及扬压力。
3、  上游水位为100 m，下游无水。
4、  大坝的材料密度为2400kg/m3,杨氏模量为2.1e11，泊松比为0.3.
5、  扬压力的折减系数为0.4.
  
其下为本人的命令流，若有不对之处，望指正。

/PREP7   
ET,1,PLANE42
KEYOPT,1,1,0
KEYOPT,1,2,0
KEYOPT,1,3,2
KEYOPT,1,5,0
KEYOPT,1,6,0
MPTEMP,,,,,,,,   
MPTEMP,1,0   
MPDATA,EX,1,,2.1e11  
MPDATA,PRXY,1,,0.3   
MPDATA,DENS,1,,24000
K,1,0,0,,   
K,2,80,0,,   
K,3,0,100,,  
FLST,2,3,3   
FITEM,2,3   
FITEM,2,1   
FITEM,2,2   
FLST,2,3,3   
FITEM,2,1   
FITEM,2,2   
FITEM,2,3   
A,P51X   
FINISH   
/POST1   
FINISH   
/SOL
FLST,2,1,4,ORDE,1   
FITEM,2,1   
/GO  
DL,P51X, ,ALL,   
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*SET,_FNCNAME,'aa'   
! /INPUT,1.func  
*DIM,%_FNCNAME%,TABLE,6,7,1  
!   
! Begin of equation: (100-{Y})*9.81*1000
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0  
*SET,%_FNCNAME%(3,0,1), 0.0  
*SET,%_FNCNAME%(4,0,1), 0.0  
*SET,%_FNCNAME%(5,0,1), 0.0  
*SET,%_FNCNAME%(6,0,1), 0.0  
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 100, 0, 0, 3
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, -1, 2, 3  
*SET,%_FNCNAME%(0,3,1),   0, -1, 0, 9.81, 0, 0, -2   
*SET,%_FNCNAME%(0,4,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,5,1), 0.0, -1, 0, 1000, 0, 0, -3   
*SET,%_FNCNAME%(0,6,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,7,1), 0.0, 99, 0, 1, -2, 0, 0  
! End of equation: (100-{Y})*9.81*1000   
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*SET,_FNCNAME,'bb'   
! /INPUT,2.func  
*DIM,%_FNCNAME%,TABLE,6,13,1
! Begin of equation: (1-{X}/80)*9.81*1000*100*0.4   
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0  
*SET,%_FNCNAME%(3,0,1), 0.0  
*SET,%_FNCNAME%(4,0,1), 0.0  
*SET,%_FNCNAME%(5,0,1), 0.0  
*SET,%_FNCNAME%(6,0,1), 0.0  
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 80, 0, 0, 2  
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 2, 4, -1  
*SET,%_FNCNAME%(0,3,1),   0, -1, 0, 1, 0, 0, -2  
*SET,%_FNCNAME%(0,4,1), 0.0, -3, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,5,1), 0.0, -1, 0, 9.81, 0, 0, -3   
*SET,%_FNCNAME%(0,6,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 1000, 0, 0, -2   
*SET,%_FNCNAME%(0,8,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 100, 0, 0, -3   
*SET,%_FNCNAME%(0,10,1), 0.0, -2, 0, 1, -3, 3, -1   
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 0.4, 0, 0, -2   
*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, -2, 3, -1   
*SET,%_FNCNAME%(0,13,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: (1-{X}/80)*9.81*1000*100*0.4  
!-->
FLST,2,1,4,ORDE,1   
FITEM,2,3   
/GO  
SFL,P51X,PRES, %AA%  
FLST,2,1,4,ORDE,1   
FITEM,2,1   
/GO  
SFL,P51X,PRES, %BB%  
ACEL,0,1,0,  
FINISH   
/PREP7   
FLST,2,1,5,ORDE,1   
FITEM,2,1   
AESIZE,P51X,2,   
MSHAPE,0,2D  
MSHKEY,1
AMAP,1,3,1,2,   
FINISH   
/SOLU   
/STATUS,SOLU
SOLVE   
FINISH

chenfanjian

还得加一个条件是：认为大坝是固结于地面的

无奈和孤独的人

大坝的问题用二维建模说服力不够啊