三角形水坝的有限元分析

有限元程序

例题:

图示为一水坝示意图,水坝的上面1米处在水面之上,余下部分浸在水中.水的分布荷载gh

μ厚

=.E=20Gpa,泊松比17.0=

Pρ

t=1m,容重为0。

求：水坝的应力分布。

图一

一．单元划分：

1.单元划分：

图二

2.载荷分析：

3675

5.16

1

)21*8.1*(*31221===g g F ρρ

22050

5.1)21*8.1*(*13132222==⎪⎭

⎫

⎝⎛++=g g F ρρ

同理可得其他节点力的分布：

132300

8.16)2

1

*8.1*(*12110250

8.15)21

*8.1*(*10882008.14)21

*8.1*(*866150

8.13)21

*8.1*(*644100

8.12)21

*8.1*(*4227226225224223===============g g F g g F g g F g g F g g F ρρρρρρρρρρ

二．数据采集：

（1）基本参数：

1.单元数：NE＝36；

2.节点数：NJ＝28；

3.支承数：NZ＝14；

4.节点荷载数：NPJ＝7；

5.半带宽：DD＝（7+1）*2＝16；

6.节点位移数：NJ2＝NJ*2＝56；

（2）其它参数：

1.问题类型码：LMX＝0（平面问题）

2.弹性模量：E0＝20000000000；

3.泊松比：MU＝0.17；

4.容重：LOU＝24000；

5.板厚t: TE=1;

6.节点坐标数组：AJZ

{0.000000,0.000000,0.000000},//1节点坐标

{0.000000,1.000000,0.000000},//2节点坐标

{0.000000,2.000000,0.000000},//3节点坐标

{0.000000,3.000000,0.000000},//4节点坐标

{0.000000,4.000000,0.000000}, //5节点坐标

{0.000000,5.000000,0.000000},//6节点坐标

{0.000000,6.000000,0.000000},//7节点坐标

{0.000000,0.000000,1.500000},//8节点坐标{0.000000,1.000000,1.500000},//9节点坐标{0.000000,2.000000,1.500000},//10节点坐标{0.000000,3.000000,1.500000},//11节点坐标{0.000000,4.000000,1.500000},//12节点坐标{0.000000,5.000000,1.500000},//13节点坐标{0.000000,0.000000,3.000000},//14节点坐标{0.000000,1.000000,3.000000},//15节点坐标{0.000000,2.000000,3.000000},//16节点坐标{0.000000,3.000000,3.000000},//17节点坐标{0.000000,4.000000,3.000000},//18节点坐标{0.000000,0.000000,4.500000},//19节点坐标{0.000000,1.000000,4.500000},//20节点坐标{0.000000,2.000000,4.500000},//21节点坐标{0.000000,3.000000,4.500000},//22节点坐标{0.000000,0.000000,6.000000},//23节点坐标{0.000000,1.000000,6.000000},//24节点坐标{0.000000,2.000000,6.000000},//25节点坐标{0.000000,0.000000,7.500000},//26节点坐标{0.000000,1.000000,7.500000},//27节点坐标{0.000000,0.000000,9.000000},//28节点坐标7.节点码数组：JM

{0,2,9,8}, //单元2 {0,2,3,9}, //单元3 {0,3,10,9}, //单元4 {0,3,4,10}, //单元5 {0,4,11,10}, //单元6 {0,4,5,11}, //单元7 {0,5,12,11}, //单元8 {0,5,6,12}, //单元9 {0,6,13,12}, //单元10 {0,6,7,13}, //单元11 {0,12,13,18}, //单元12 {0,12,18,17}, //单元13 {0,11,12,17}, //单元14 {0,11,17,16}, //单元15 {0,10,11,16}, //单元16 {0,10,16,15}, //单元17 {0,9,10,15}, //单元18 {0,9,15,14}, //单元19 {0,8,9,14}, //单元20 {0,14,15,19}, //单元21 {0,15,20,19}, //单元22

{0,16,21,20}, //单元24

{0,16,17,21}, //单元25

{0,17,22,21}, //单元26

{0,17,18,22}, //单元27

{0,21,22,25}, //单元28

{0,21,25,24}, //单元29

{0,20,21,24}, //单元30

{0,20,24,23}, //单元31

{0,19,20,23}, //单元32

{0,23,24,26}, //单元33

{0,24,27,26}, //单元34

{0,24,25,27}, //单元35

{0,26,27,28}, //单元36

8.支承数组：NZC

NZC[NZ+1]={0,1,2,3,4,5,6,7,8,9,10,11,12,13,14};

9.节点载荷数组：PJ

PJ[NPJ+1][2+1]={{0,0,0},{0,3675,55},{0,22050,51},{0,44100,45},{0,66150,37 },{0,88200,27},{0,110250,15},{0,132300,1}};

3.主程序：

#include

#include

#define NE 36 //单元数

#define NJ 28 //节点数

#define NZ 14 //支承数

#define NPJ 7 //节点载荷数