结点矩形薄板弯曲元静力线性分析程序设计.docVIP

4结点矩形薄板弯曲元静力线性分析程序设计

------------------------------------------------------------------------------------------------

4结点矩形薄板弯曲元静力线性分析程序设计

有限单元法作业

(一):选题

(二):基本理论及公式

采用Kirchhoff理论

一、矩形单元的典型形式

将图示矩形薄板沿x,y方向划分成若干小矩形

从中取出一小矩形(单元)，共有四个结点，每个结点有三个位移分量:挠度w

绕x轴转角

?挠度W???绕x轴转角?x?上节为沿y方向倾角

?绕y轴转角??上节为沿x方向倾角?y?

即结点i的位移

???Wi??W??????di????xi????y??i?1,?4?

??????yi?????x?i?

同理,相应的结点力

?fi?竖向力??Fi???M?xi?绕x轴力偶(上节中的My)?M?绕y轴力偶(上节中的M)x?yi?

符号重新定义是为了有限元表示的方便，由此得单元结位移向量

?d?e?

?d1

——————————————————————————————————————

------------------------------------------------------------------------------------------------

?d4??W1?x1?y1?W4?x4?y4T??T1

节点力

?F?e??FT

T

1

?F4???f1Mx1My1?f4Mx4My4?

二、位移模式(函数)

,、位移模式的选取

插值多项式取为:

W?x,y???1??2x??3y??4x2??5xy??26y??7x3??8x2y??9xy2??10y3??11x3y

??12xy3

三、形函数和形函数矩阵

{W}??N??d?e

?Nx1Ny1?N4Nx4Ny4

?形函数矩阵:?N???

N1??N1/?y

??????N?

y4/?y???

?N1/?x?

?

??

——————————————————————————————————————

------------------------------------------------------------------------------------------------

?

?Ny4/?x??式中形函数:

Ni?18

?1??i???1??i???2??i???i

???2??2?

N1

xi??8

b?i?1??i???1??i???1??2?

Nyi?18

a?i?1??i???1??i

???1??2?

(i=1234)

局部坐标与整体坐标的关系为:??1?x?x?

a

??1b

?y?y0

?

??z

四、单元的几何矩阵[B]和内力矩阵[S]

1(几何矩阵[B]

2

???2N1?2N?2Ny1?2Ny4——————————————————————————————————————

------------------------------------------------------------------------------------------------

?x2x1?x2?x2??

???2?2?B????2N

?1?y2????N?

y4

?y2?

???2

?2?N1x?y???2?2Ny4?

x?y??

或以子块形式表示:[B]=[B1B2B3B4]。

式中:

2.内力矩阵[S]

?F???D??B??d?e??S??d?e

3

五、单元刚度矩阵

由一般公式得:?K????ab

?a?bt?B??D??B?dxdy。将几何矩阵[B]和弹性矩阵[D]的表达式代T

入，积分可得薄板弯曲问题矩形单元的单元刚度矩阵的显示:

4

六、荷载等效变换

由荷载等效变换的一般公式可得

?R???ab

——————————————————————————————————————

------------------------------------------------------------------------------------------------

?a??b[N(x,y)]Tq?x,y?dxdy

1(法向均布荷载q

代入上述公式得:

5

?1??b??3?N?1??a??N?????x1??3??Ny1??1????b??N2??3??N??a??x2??N?