收稿日期 :2007‐03‐03 ;修改稿收到日期 :2007‐06‐28畅
基金项目 :国家自然科学基金(10772035 ,10721062 ,
10472019) ;973 项目(2005CB321704) ;辽宁省
中青年学术带头人基金资助项目 .
作者简介 :王 凯(1982‐)男 ,硕士 ;
杨海天
倡
(1956‐)男 ,教授 ,博士生导师
(E‐mail :haitian@ dlut .edu .cn) .
第25卷第5期
2008 年 10 月
计 算 力 学 学 报
Chinese Journal of Computational Mechanics
Vol .25 ,No .5
October 2008
文章编号 :1007‐4708(2008)05‐0616‐05
高 斯 牛 顿 技 术 求 解 偶 应 力 反 问 题
王 凯
1
, 杨海天
倡
, 马莉英
2
(1 .大连理工大学 工程力学系 工业装备结构分析国家重点实验室 ,大连 116024 ;
2 .大连交通大学 土木与安全工程学院 ,大连 116028)
摘 要 :建立了便于敏度分析的偶应力反问题数值求解模型 ,给出了直接法和伴随法两种敏度计算格式 。 在反
演计算中采用了高斯牛顿技术对未知本构参数进行识别 ,探讨了测点数目 、初值选取和数据噪音对反演结果的
影响 ,数值算例给出了令人满意的结果 。
关键词 :偶应力 ;反问题 ;伴随法 ;高斯牛顿
中图分类号 :O302 文献标识码 :A
1 引 言
偶应 力 理 论 (Cosserat 理 论 )于 1909 年 由
Cosserat
[1]
兄弟提出 ,Topin
[2]
,Mindlin
[3]
和 Erin‐
g
en
[4]
等在 20 世纪 60 年代对该理论作了进一步发
展和完善 。 偶应力理论是研究具有微结构 、一定特
征尺度的介质的重要工具之一 ,如纤维材料 、颗粒
材料等力学行为 ,在实际问题中有许多应用 :如剪
切带分析 ,尺寸效应相关的硬化 、软化等分析 ,在岩
土工程中被用于描述考虑弯曲的层状岩体
[5]
、陡倾
角层状岩体中地下洞室围岩变形
[6]
等 。
本构参数的确定 ,特别是材料长度参数 l 的确
定(它是一个内变量 ,很难测量得到) ,是偶应力理
论应用的前提 。 除试验手段外 ,借助某些附加信
息 ,通过求解相应的偶应力连续介质场反问题以确
定本构参数 ,也是解决这类问题的方法之一 。 目
前 ,关于这方面直接的研究报道似不多见 。
基于以上考虑 ,本文借助有限元 、敏度分析及
高斯牛顿技术 ,考虑材料的非均质性 ,建立了一种
求解偶应力反问题的一般数值模式 ,并在高斯牛顿
求解中添加了一维搜索过程 ,以修正初值选取不当
对解的偏离 。 文中讨论了测点数目 、初值选取和数
据噪音对反演结果的影响 。 算例表明 ,所提算法具
有较高的精度和较好的抗噪性 。
2 控制方程
平面偶应力问题的平衡方程可写为
[7]
1
2
[
σ
i
j
+
σ
j
i
]
,
j
+
1
2
[
σ
i
j
-
σ
j
i
]
,
j
+
f
i
=
0
μ
i ,i
+
e
i
j
z
σ
i
j
+
m
=
0
(1)
式中 i ,
j
=
x ,
y
,
σ
i
j
为通常的应力部分 ,
μ
i
为偶应
力部分 ,
f
i
和 m 分别是体力和体力偶 ,e
i
j
z
为置换符
号 。
应变位移关系可写为
ε
i
j
=
u
i ,
j
+
e
i
j
z
φ
z
κ
i
=
φ
z ,i
(2)
式中 i ,
j
=
x ,
y
,
ε
i
j
为对应通常应力部分的变形张
量 ,
κ
i
是与偶应力相对应的曲率 ,u
i
是平动位移
φ
z ,i
是微观转角 。
本构方程为
σ
xx
σ
yy
σ
x
y
σ
y
x
μ
x
μ
y
=
A H 0 0 0 0
H A 0 0 0 0
0 0 B C 0 0
0 0 C B 0 0
0 0 0 0 Gl
2
0
0 0 0 0 0 Gl
2
·
ε
xx
ε
yy
ε
x
y
ε
y
x
κ
x
κ
y
(3)
式中 A
=
E/(1
-
υ
2
) ,H
=
E
υ
/(1
-
υ
2
) ,B
=
G
+
G
c
,
C
=
G
-
G
c
,G
=
E/2(1
+
υ
) ,E 为弹性模量 ,
υ
为泊
松比 ,l 为表征尺度效应的特征尺度 ,G
c
=
c G 为
Cosserat 剪切模量 ,根据文献[5] ,c 可取为 0畅5 。
边界条件为
u
=
u
-
, x
∈
Γ
u
(4)