对转角场和剪应变场进行合理插值的厚板元

  • 格式:pdf
  • 大小:576.29 KB
  • 文档页数:14

第巧

卷第󰀁

󰀂 !

年!

月工程力学

∀#

协田∀%&∋#∃(

伍)∗+#,)−./0

󰀂2#

/󰀁

+34󰀂 !

转角

场和剪

应变场

进行合理插值的厚板

元’

岑松龙

志飞

5清华大学土

木系北京󰀂666!785中国矿大北研部

北京󰀂666!󰀁8

提要本文提出构造厚板三角形元

的一种合理方法先按9/:(

;<=

/

厚梁理论导出单

元各

边的转角

和剪

应变公式然后进行合理

插值

导出单

元的转角场曲率场和剪

应变

场当

板的厚度变小时厚梁理

论自动退化

为薄梁理论各边剪应变以及单

元剪应变插值函数自动

退化为零厚板单

元自动退化为薄

板单

元彻底消除了剪

切闭锁现象

此单

元对厚板和薄

都给

出高精度

的结果

关键

词有限

元厚板元三角形元厚梁理论

应用%;>:

:<;?≅>

><

理论构造厚板元

虽然具

有只需满足)6连续性构造方法简

便

的优点但常

出现剪切闭锁现象

而不能

用于分

析薄

板问题为了

克服闭锁

问题已

经提

不少方法例

如缩

减积分

法川选择性缩减积分法0Β&

代替剪应变

法∋:&

等这类方法的特

点是

采用从厚板

元过渡到厚一薄板元的步骤Χ

首先根

据厚板理论构造出适用

于厚板的单

元然后采

取措

施在退化

薄板元时设

法克服剪切闭锁

现象

另一类

方法的特点是采

用从薄

板元过渡到

薄一

厚板

元的步骤∋

材0Χ

首先根

据薄板理论构

出适用于薄

板的单元

然后

采取措施引人剪切应变的影响把薄

板元推广

成为对薄板

家自

然科

学基金

高校博

点基金及煤

炭系统

学回国人员科技基金

资助项目

本文收稿日

期Χ󰀂

Δ

年7月

Ε

程力学

和厚板

都通用的单

元这样构造

出来的单元自然不会

出现剪切闭锁

本文

提出构造厚板

角形元的

一种

合理方法此

单元记

为9≅9

元仃

血/:(

;<=

。≅>

0><

角形元8首先按9

;<=

。厚梁理论确定单

元各边的切向转角

和切向剪应变的分布函

数Φ

然后

根据各

边分布函

对单

元域内进行插值确定单元的

转角场和剪应变场进而

定单

元的曲率场和单

元刚度矩阵

当厚度变小时厚梁理论自动

退化为

薄梁理论单元

边剪应变以

及单

元剪应变插值函数自动退化为零厚板元自动退化为

薄板元因

此9≅

9

元不会

出现剪切闭锁现象

厚梁元

的挠度

转角和

剪应

变公式

对于图0

所示厚梁元0

’6&

可得出挠度Γ

转角笋

和剪应变Η

公式如

下Χ

卜Ιϑ

ΚΛ)󰀂0Μ0

ΜΜ

ΜΜ

ΓΝΓ,50

一ϑ8ΟΓΠϑΟ

二5

,一价Χ8

户Χ

二0ϑ0

一Β。8

户󰀁

艺‘

拼Ν

班050

一ϑ8Ο

班ΠϑΟ󰀁

50

一Ε占8

Ν

占?

其中

Ν

子‘

一ΟΘΠ

,一拼󰀂一梦Ε

卫‘

0Ο󰀂Ε兄5

五8

Ε

0

舟50一、8

·Ε5

今,’

250

刀80

’Ρ

Ν∀(

󰀂Ε

50

一召’ΣΝ

互∃(

Ν5Ε8

凡Νϑ50

一?8

凡Νϑ50

一ϑ850

一Π

?8

式508的推导过程如

下Χ

设沁哟Γ

分别为

常数二次式和三次式、

根据位移端点条件设

对转角场和剪

应变场进行合理

插值

的厚板元

Π

二?

/

班Ν

笋,50

一ϑ8Ο

俨ΠϑΟΚ

50

一ϑ8

ΓΝΓ

,50

一ϑ

8ΟΓΠϑΟ

刀。,

ϑ50

一ϑ8Ο

刀00

ϑ50

一ϑ850

一Π

ϑ85󰀁8

矜梦,

应满

足下式Χ

ΑΓ

下二,一班Ν

0ΑΤ578

将5󰀁8代

人578

夕。

将528

人5󰀁8

得Χ

Ν

喜5

梦0

一俨Ε8

刀,ΝΧ

。一

。。Ν一Υ

5

厂。一

粤?

8528

Π

乙乙

笋二?

󰀂

Μ

梦,50

一ϑ8Ο

尹Πϑ

一Υ

妙。一

育08ϑ

吸0

一?8

乙5Υ8

ΘΝΘ

,5

,一,8

一‘Ο

合5

班】一俨Ε8

‘5

,一‘8Ο

‘5Φ

。一

合?

8

‘5

‘一‘8“

一,‘8

其中

。为一

内部参数

可由应变能最小条件求

出如下Χ

曲率Χ‘

Ν一

些Ν

粤,价,一班ΕΟ

Υ5

。一

喜?85卜Ε

,8

、。Ο

粤5卜Ε

,8

0Αϑ0

乙0

ς

。ΝΩΧ

俨󰀂

一犷Ε

一,5

‘一,‘

,,

弯曲应变能Χ

Ξ

,Ρ,

ΨΡ0

ςΑϑΝ

—0

ΕΨΧ

、。ΕΟ

鱼、。‘󰀂

一Χ

,8Ο

0

”050

一Π

8

’&Αϑ

Ν3Χ

竿5?

。?

一厂

,

其中

Ξ

ΧΝ

警ΨΦς

。’

剪切应变能Χ

应变能

。Ν

。、。Ν

刃[

互卫上生Ο

奖5󰀂ΕΟ

粤8Χ

0艺0兀

爵一6

‘导Χ

Υ

。二

石石石石

0Ο0乙

Ν

占?

导毕∴

7

工程力学

三厚板元转角场

的插值

公式

考虑图Ε

所示厚板三角形元单元有

三个结点每个结点有

三个

工程自

由度单元的九

个自

度组成的向量为Χ

5ΓΠ

尹Π

班,Ε8

5

Γ0

叭Χ

竹Χ85

Γ󰀁

汽󰀁

梦,

,8

󰀁

户咬

乡了_

图Ε

叮ΩΝ∋Γ

,俨,班,

,ΓΠ

梦Ε

梦,ΠΓ

,梦,俨,󰀁&9

根据

各边的转角

公式进行插

值可得到

单元转角

场的插

值公式

,各边法向和切向

转角梦

树公式

各边的法向转角班

设为线性分布单元三个边5Ε󰀁󰀁󰀂󰀂Ε8的转角

式为Χ

汽。Ν5叭Ε󰀁8ΠαΧΟ5汽Ε󰀁8󰀁

俨󰀁󰀂Ν5

叭󰀁Χ8󰀁

几Ο5汽󰀁,8Φ

汽ΧΕΝ

5汽ΧΕ8ΧαΧΟ5汽󰀂Ε8ΕαΕ

各边的切向转角叭可由式

50Λ8

得出如下Χ

笋Ε,Ν5

俨Ε󰀁

人Ο5汽Ε󰀁8

󰀁

几Ο󰀁

50

一Ε

占Χ8矶α

Π

梦󰀁,Ν

5

笋󰀁Χ8󰀁α󰀁

Ο5

笋󰀁Χ8Φα,Ο󰀁50

一Ε占Ε8几几α

尹ΦΕΝ

5

梦ΧΕ

8

,

石Ο5

尹ΧΕ8ΕαΕΟ󰀁50

一Ε

占󰀁8几

αΧαΠ5Δ8

5:Κ8

5:Λ

8

其中(

叹了8

0[

2

Ο

丁5󰀂

一产8

Υ5>Ν0Ε󰀁

85

8

(

一护去

,乙

、6了∀

奇〔,‘[

一,

一‘俨

二班·,·

‘班刀·

梦诚

,,

Λ

‘ΝΗ

,

一Η

=

;Ν]

女一

],5>Ψ=Ν0Ε󰀁8

Ε

各边

转角梦班,

的表示式

用各

边5笋

汽8

和5俨俨,8之间的关

系可求得各边转角俨笋,

的表示式

如下Χ