对转角场和剪应变场进行合理插值的厚板元
- 格式: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
一Π
ϑ858
矜梦,
应满
足下式Χ
ΑΓ
下二,一班Ν
厂
0ΑΤ578
将58代
人578
夕。
将528
代
人58
求
得Χ
Ν
喜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
俨Ε
8Ε
人Ο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之间的关
系可求得各边转角俨笋,
的表示式
如下Χ