计算力学4---加权残值法解矩形薄板的弹塑性问题

加权残值法解矩形薄板的弹塑性问题

4．加权残值法解矩形薄板的弹塑性问题

薄板理论的克希霍夫假设在板的弹塑性分析中仍可应用。采用增量形式表示，板的本构方程的矩阵形式为：

{ } ([De] [Dp]{ e}) （4-7-1）

式中

{ } [ x, y, xy]T

（a） T

{ e} [ x, y, xy]

分别为应力增量分量和应变分量增量。而弹性矩阵

1E

[De]

1 2

0

10

0

0 （b） 1 2

T

塑性矩阵

f f

[De] [De]

（4-7-2） [De] T f f H [De]

这里H d s/dep为硬化参数；f为屈服函数，对于密赛斯屈服条件

f s 0 （4-7-3）

式中

2221/2

( x y xy 3 xy) （c ）

式（4-7-1）中的β叫做塑性修正系数，在弹性区内β=0；在塑性区β=1；在弹

塑性过渡区，取

a s

a （d ） b

上标b,a分别表示加上载荷增量前后的值。板的应力偏量

11

Sx x ( x y)Sy y ( x y) （4-7-4）

33

将有关公式代入式（4-7-2）中，则得塑性矩阵：

(Sx Sy)2(Sx Sy)(Sy Sx)(1 )(Sx Sy) xy E 2

[Dp] (S S)(S S )(S S)(1 )(S S) xyyxyxyxxy

(1 )Q 22 (1 )(Sx Sy) xy (1 )(S S) (1 ) xyyxxy

（4-7-5）

其中

加权残值法解矩形薄板的弹塑性问题

Q S S 2 SxSy (1 ) 2x2y2xy

4H (1 2) 2 （e）

9E

分开应力增量{△σ}的弹性部分和塑性部分，沿板厚积分，即有

{ M} { Me} { MP} （f）

式中

{ M} [ Mx, My, Mxy]T (g)

弹性弯矩增量和挠度增量的关系

2 2

M D[2( w) 2( w)]

x y 22 e

My D[2( w) 2( w)] （4-7-6）

y x 2

e

Mxy D(1 )( w)

x y

e

x

Eh3

式中，D 为板的抗弯刚度。 2

12(1 )板的静力平衡条件有：

2 2 2e

( M) 2( M) ( M) q 0 （4-7-7） xxyy22

x y x y

将式（f）代入，同时考虑到方程（4-7-6），则得薄板弹塑性弯曲的控制方程：

4 4 4

D[4( w) 222( w) 4( w)

x x y y其中 qe为对应于弹性变形的外载荷增量，而

qe qp 0 （4-7-8）

2 2 2pp q 2( Mx) ( Mxy) 2( Myp) （4-7-9）

x y x y

p

称为塑性载荷项，它与当前的应力水平有关。

设有一正方形薄板，边长为a，取无量纲量

yxx

, ,hx ,hx aaaa （h ） p 22

V w(1 )apn , q , n 2

hEh Eh

取板弯曲的挠度函数

i 11j 1

A ( )

ij

i

N 1M 1

j

( ) [ ( )] [ ( )]{ A} （i）

加权残值法解矩形薄板的弹塑性问题

式中， 为克罗内斯克（Kronecker）乘积符号，而[ ( )],[ ( )]是由五次B样条函数构成的基函数。

由于位移试函数既不满足控制方程（4-7-8），又不满足边界条件，因此采用加权残值法中的混合法。将试函数（i）代入控制方程（4-7-8）及有关的边界条件，同时取样条结点为配点坐标，得残值方程：

{R} [H]{ A} { { p}} （j）

方差泛函

{I} {R}T{R} ({ A}T[H]T { }T { p}T)([H]{ A} { { p}) （k） 方差最小二乘极值条件

{I}

0 （I） { A}

求得

{ A} ([H]T[H]) 1[H]T({ p}) （4-7-10）

上式中{ p}是由塑性变形所引起的内力修正项，它随着载荷的增加不断变化，因此只能用数值逼近法求解。其基本思想是，对每一步载荷增量，{ p} {0}，则控制方程（4-7-8）变成：

4 4 4

D[4( w) 222( w) 4( w)] qe 0 （m） x x y y

上式为完全弹性方程。求解方程（m），计算应力，应变及位移，并进行屈服判断，当f( ij)＜0时，板属于弹变形；当f( ij)≥0，板发生塑性变形。对于应力超过屈服面的部分，求修正应力分量，进而求内力修正项{ p}。再将{ p}当作外载荷代入方程（m）中求解。如此循环反复，直到{ p}减小到允许范围，即可结束这一步加载。下一步的载荷增量类似处理，直到所有载荷加完。

作为具体算例，四边固定理想弹塑性矩形薄板， 承受均布载荷q作用，计算参数 0.3,E 3.0 104MPa,h 0.2m,a b 6.0m.取N M 8。图4-7-1给出了板中点处的载荷挠度曲线及屈服域扩展顺序图。

加权残值法解矩形薄板的弹塑性问题

最小二乘配点法分析圆形弹性地基板的弹塑性弯曲问题

用最小二乘配点法，对任意分布载荷作用下弹性地基上圆板的弹塑性弯曲问题进行了分析。方法简便可行，具有良好的精度。

已知文克勒（Winkler）地基上圆形薄板的静力平衡微分方程为：

2 1 1 21 (rM) 2[(M)] (M) (M ) kw(r, ) q 0r r 222 r rr r r r

（4-8-1）

圆板的Misses屈服条件为

r 3 2

r

2

2r

2

s （4-8-2）

相应于Misses屈服条件下的Levy-Missis弹塑性本构关系为

d r d rp d (2 r )

p

d d d (2 r)

（4-8-3） p

d r d r 6d r

这里忽略了弹性变形。为了求得应力分量，反解上式，并应用几何方程：

2w

r z2

r

1 w1 2w

z 2

2 r rr

（a）

1 w rr z r r

加权残值法解矩形薄板的弹塑性问题

则得

(2d r d )z 21 1 2

r 22(dw) (dw) 2(dw) …… 此处隐藏：7933字，全部文档内容请下载后查看。喜欢就下载吧 ……