当前位置：文档之家 › 基于偶应力理论的高阶弹塑性本构模型的无网格法

基于偶应力理论的高阶弹塑性本构模型的无网格法

  • 格式：pdf
  • 大小：1.37 MB
  • 文档页数：5

下载文档原格式

  / 5

合集下载

弹塑性力学-弹塑性本构关系ppt课件

弹塑性力学-弹塑性本构关系ppt课件

为非负，即有 0
功，即 0
（应变硬化和理想塑性材料）
（应变软化材料）
工程弹塑性力学·塑性位势理论
(2) 德鲁克塑性公设的表述
德鲁克公设可陈述为：对于处在某一状态下的稳定材 料的质点(试件)，借助于一个外部作用在其原有应力状态 之上，缓慢地施加并卸除一组附加压力，在附加应力的施 加和卸除循环内，外部作用所作之功是非负的。
Ñ W
0 ij
ij
0 ij
d ij 0
Ñ 由于弹性应变εije在应力循环
中是可逆的，因而
( ij
0 ij
)
d
e
ij
0
0 ij
于是有：
Ñ WD WDp
( ij
0 ij
)d
p
ij
0
0 ij
工程弹塑性力学·塑性位势理论
(3) 德鲁克塑性公设的重要推论
Ñ WD WDp
( ij
0 ij
)d
势理论。他假设经过应力空间的任何一点M，必有一
塑性位势等势面存在，其数学表达式称为塑性位势函
数，记为：
g I1, J2, J3, H 0
或
g ij , H 0
式中， H 为硬化参数。
塑性应变增量可以用塑性位势函数对应力微分的表达
式来表示，即：
d
p ij
d
g
ij
工程弹塑性力学·塑性位势理论
不小于零，即附加应力的塑性功不出现负值， 则这种材料就是稳定的，这就是德鲁克公设。
工程弹塑性力学·塑性位势理论
在应力循环中，外载所作的 功为：
Ñ W
0 ij
ij
d ij
0
不论材料是不是稳定，上述 总功不可能是负的，不然， 我们可通过应力循环不断从 材料中吸取能量，这是不可 能的。要判断材料稳定必须 依据德鲁克公设，即附加应 力所作的塑性功不小零得出

第5章 弹塑性本构模型理论

第5章 弹塑性本构模型理论
?r3321??????232221??????op表示p点应力矢量的大小p为应力空间上一点代表某一应力状态过p点作与等倾线相垂直的面即为平面???32131???????oq??321321313?????????????rr由mioqr?331????平面上法向应力即为令????m3?????23212322212231??????????????oqoppq??????2222323122132231qj??????????????平面上的剪应力为令???qj3222??应力洛德角1?2?3?p?qr洛德参数313122???????????毕肖甫常数3132???????b洛德角??312332tan31312???????b????????洛德角与偏应力不变量之间的关系23232333sinjj????应变与应变增量ji??333231232221131211??????????????????zzyzxyzyyxxzxyx??????????????????212121212121??????????321??????应变状态体积应变增量321???????????v偏差应变增量????????3vijijijee?????应变张量不变量3211???????i偏差应变不变量3231212??????????i3211?????i01??j??232221121eeej??????????1233123213222271?????????????????j体积应变321???????v广义剪应变????????2123223122132?????????????应变洛德角??3131232tan???????????2增量弹塑性理论?弹性增量理论?以弹性模型与泊桑比表达????????????????????????????????????????????????????????????????????????????????????????zxyzxyzyxzxyzxyzyxvvvvvvvvvvvvvvvvvvvvve????????????????????????1221000000122100000012210000001110001110001112111?以剪切模型与体积模量表达?????????????????????????????????????????????????????????

本构模型

本构模型

2 应力-应变曲线
应力－应变反应与变形率无关的材料称为率无关；否则， 称为率相关。名义应变率定义为
x L0
因为
L 和
x x
L0 L L0 x
即名义应变率等于伸长率，例如 可以看出，对于 率无关材料的应力－ 应变曲线是应变率独 立的，而对于率相关 材料的应力－应变曲 线，当应变率提高时 是上升的；而当温度 升高时是下降的。
4 非线性弹性
对于一个由三个彼此正交的对称平面组成的正交材料(如木 材或纤维增强的复合材料)，仅有9个独立弹性常数，Kirchhoff 应力－应变关系为材料对称坐标平面，为正交各向异性体
对于各向同性材料，仅有3个常数
C11 C22 C33 C1 C21 C23 C31 C2 C44 C55 C66 C3
4 非线性弹性
小应变和大转动
许多工程应用包括小应变和大转动。在这些问题中，大变形 的效果主要来自于大转动，如直升机旋翼、船上升降器或者钓鱼 杆的弯曲。由线弹性定律的简单扩展即可以模拟材料的反应，但 要以PK2应力代替其中的应力和以Green应变代替线性应变，这称 为Saint-Venant- Kirchhoff材料，或者简称为Kirchhoff材料。 最一般的Kirchhoff模型为
3 一维弹性
2 2 应变能一般是应变的凸函数，例如， (w( 1 ) w( x ))( 1 x ) 0 x x
当
2 1 x x
公式的等号成立。
凸应变能函数的一个例子如图所示。在这种情况下，函数 是单调递增的，如果w 是非凸函数，则 s 先增后减，材料应变 软化，这是非稳定的材料反应， ds d x 0 如右下图。

非线性专题-本构模型

非线性专题-本构模型

非线性关系
4 非线性有限元的分类 非线性有限元的分类
非线性分析包含下列步骤： 非线性分析包含下列步骤： • 建立模型－ 建立模型－Pre-process • 基本方程的公式 • 离散方程 • 求解方法 • 表述结果－ 表述结果－Post-process
TSINGHUA UNIVERSITY
4 非线性有限元的分类 非线性有限元的分类
TSINGHUA UNIVERSITY
4 非线性有限元的分类 非线性有限元的分类
边界非线性
如果边界条件在分析过程中发生变化， 如果边界条件在分析过程中发生变化 ， 就会产生 边界非线性问题。 边界非线性问题。悬臂梁随着施加的载荷产生挠曲。 悬臂梁随着施加的载荷产生挠曲。 梁端点在接触到障碍物以前， 梁端点在接触到障碍物以前 ， 其竖向挠度与载荷 成线性关系（ 成线性关系 （ 如果挠度是小量）。 如果挠度是小量 ）。当碰到障碍物时梁 ）。 当碰到障碍物时梁 端点的边界条件发生了突然的变化， 端点的边界条件发生了突然的变化 ， 阻止了任何进一 步的竖向挠度， 步的竖向挠度 ， 因此梁的响应将不再是线性的。 因此梁的响应将不再是线性的 。 边界 非线性是极度的不连续； 接触 时 ， 结 非线性是极度的不连续 ； 当在模拟中发生接触 当在模拟中发生 接触时 构中的响应在瞬时会发生很大的变化。 构中的响应在瞬时会发生很大的变化。 另一个边界非线性的例子是将板材材料冲压入模 具的过程。 具的过程 。 在与模具接触前， 在与模具接触前 ， 板材在压力下比较容易 发生伸展变形。 发生伸展变形 。 在与模具接触后， 在与模具接触后 ， 由于边界条件的改 变，必须增加压力才能使板材继续成型。 必须增加压力才能使板材继续成型。
TSINGHUA UNIVERSITY

基于虚拟纤维的各向异性超弹性材料本构模型设计

基于虚拟纤维的各向异性超弹性材料本构模型设计

基于虚拟纤维的各向异性超弹性材料本构模型设计
卢子璇;何浩;吴笛;刘学慧
【期刊名称】《计算机辅助设计与图形学学报》
【年(卷),期】2024(36)4
【摘要】为了拓展图形学弹性体模拟中的各向异性超弹性虚拟材料种类,建立了基于虚拟纤维的本构模型.首先从能量可加性出发,将超弹性体应变能量密度函数分解为轴向、剪切、体积分量的纤维加和形式,然后建立单根纤维的轴向应变、剪切应变、体积应变的应变度量,最后推导出各分量的应力表示.仿真实验使用基于四面体的非线性有限元法(finite element method,FEM),半隐式时间积分进行解算,并采用CPU串行算法,测试了不同场景下非线性能量函数以及纤维权重组合对虚拟纤维材料刚度、泊松效应、轴向特性的影响.结果表明,虚拟纤维本构模型具有大变形稳定性,材料参数设置可良好地展现上述物理特性,相比现有的横观各向同性模型具有更丰富的可调节能力.
【总页数】10页(P533-542)
【作者】卢子璇;何浩;吴笛;刘学慧
【作者单位】中国科学院软件研究所计算机科学国家重点实验室;中国科学院大学计算机科学与技术学院;中国科学院软件研究所人机交互技术与智能信息处理实验室;中国科学院力学研究所微重力重点实验室
【正文语种】中文
【中图分类】TP391.41
【相关文献】
1.人体椎间盘纤维环各向异性超弹性材料本构模型的数值验证
2.帘线/橡胶复合材料各向异性黏-超弹性本构模型
3.基于Voigt模型和Reuss模型的三方晶粒各向异性集合的弹性本构关系
4.考虑纤维弯曲刚度的橡胶-帘线复合材料各向异性超弹性本构模型
5.涤纶增强橡胶基复合材料各向异性超弹性本构研究
因版权原因，仅展示原文概要，查看原文内容请购买。

非线性有限元9弹塑性本构关系ppt课件

非线性有限元9弹塑性本构关系ppt课件
单轴试验下材料的弹塑性性态 （1/3）
对塑性变形基本规律的认识来自于实验： • 从实验中找出在应力超出弹性极限后材料的特性; • 将这些特性进行归纳并提出合理的假设和简化模型，
确定应力超过弹性极限后材料的本构关系; • 建立塑性力学的基本方程; 1) 求解这些方程，得到不同塑性状态下物体内的应力和
应变。
• 塑性阶段：继续加载，材料可承受 更大应力，称为材料强化，并伴随 出现塑性应变。至A点以前卸载， 路径接近直线，即处于弹性卸载状 态，其斜率等于加载斜率E。
1) 破坏点：继续加载至可承受的最大 极限应力，试件出现颈缩而破坏，
称为强度极限。
在日常生活中，随处都可以看到浪费 粮食的 现象。 也许你 并未意 识到自 己在浪 费，也 许你认 为浪费 这一点 点算不 了什么
1913年：泰勒（Taylor）的实验证明，LevyMises本构关系是真实情况的一阶近似。
1924年：提出塑性全量理论，伊柳辛（Ilyushin） 等苏联学者用来解决大量实际问题。
1930年：罗伊斯（Reuss）在普朗特（Prandtle） 的启示下，提出包括弹性应变部分的三维塑性应力 -应变关系。至此，塑性增量理论初步建立。
（屈服点），描写多维问题的屈服条件就需要应力或应变空间的一个临界曲面，该
曲面称为屈服面。
考虑到塑性变形与静
水压力无关的特点
f1,2,3C
FJ2,J3C
至今已出现许多屈服理论。俞茂宏教授在这方面做出了重要贡献。 屈服函数：
是描写屈服条件的函数。不同屈服条件，其屈服函数不尽相同。
在日常生活中，随处都可以看到浪费 粮食的 现象。 也许你 并未意 识到自 己在浪 费，也 许你认 为浪费 这一点 点算不 了什么
基本实验有两个： • 简单拉伸实验：实验表明，塑性力学研究的应力与应变

应变梯度理论的新进展_一_偶应力理论和SG理论

第21卷第2期机 械 强 度V o l.21N o.2 1999年6月JOU RNAL　O F　M ECHAN I CAL　STR EN GTH June1999应变梯度理论的新进展(一)Ξ——偶应力理论和SG理论RECENT AD VANCES IN STRA IN GRAD IENT PLAST I C IT Y-——Couple stress theory and SG theory黄克智ΞΞ　邱信明　姜汉卿(清华大学工程力学系,北京100084)Hw a ng Ke hchih　Q iu X inm ing　J ia ng Ha nq ing(D ep a rt m en t of E ng ineering M echan ics,T sing hua U n iversity,B eij ing100084,Ch ina) 摘要　介绍两种应变梯度塑性本构模型:CS应变梯度塑性理论——偶应力理论、SG应变梯度塑性理论。

并对它们在断裂力学中的应用进行了评述。

给出一种考虑可压缩性的方法,并根据这种模型用薄梁弯曲的例子给出了可压缩性的影响。

本文的讨论虽限制在形变理论范围内,但按照相应的方法也可以得到流动理论的形式。

关键词　应变梯度　塑性　偶应力　高阶应力　断裂中图分类号　O344Abstract　In the paper tw o k inds of fram ew o rk of strain gradien t p lasticity recen tly developed and their app licati on s are review ed:strain gradien t p lasticity fo r CS so lid——the coup le stress2theo ry,strain gradien t p lasticity fo r SG so lid.T he app licati on s are m ain ly focu ssed on the fractu re p rob lem s.O ne w ay of accoun ting fo r m aterial comp ressib ility is suggested.T he review is confined to the defo rm ati on theo ry versi on,though the flow theo ry versi on can be parallelly con structed.Key words　stra i n grad ien t,pla stic ity,couple stress,h igher-order stress,fracture1　引言新近的试验表明,当非均匀塑性变形特征长度在微米量级时,材料具有很强的尺度效应。

金属材料在复杂应力状态下的塑性流动特性及本构模型

金属材料在复杂应力状态下的塑性流动特性及本构模型
秦彩芳;许泽建;窦旺;杜雨田;黄风雷
【期刊名称】《爆炸与冲击》
【年(卷),期】2022(42)9
【摘要】工程应用中,金属材料和结构往往处于复杂应力状态。

材料的塑性行为会受到应力状态的影响,要精确描述材料在复杂应力状态下的塑性流动行为,必须在本构模型中考虑应力状态效应的影响。

然而,由于在动态加载下材料的应变率效应和应力状态效应相互耦合、难以分离,给应力状态效应的研究和模型的建立造成很大困难。

通过对Ti-6Al-4V钛合金材料开展不同加载条件下的力学性能测试,提出了一个包含应力三轴度和罗德角参数影响的新型本构模型,并通过VUMAT用户子程序嵌入ABAQUS/Explicit软件。

分别采用新提出的塑性模型和Johnson-Cook模型对压剪复合试样的动态实验进行了数值模拟。

结果表明,新模型不仅在对材料本构曲线的拟合方面具有较强的优势,而且由该模型所得到的透射脉冲和载荷-位移曲线均更加准确。

因此,该模型能够更精确地描述和预测金属材料在复杂应力状态下的塑性流变行为。

【总页数】11页(P77-87)
【作者】秦彩芳;许泽建;窦旺;杜雨田;黄风雷
【作者单位】北京理工大学爆炸科学与技术国家重点实验室
【正文语种】中文
【中图分类】O347.3
【相关文献】
1.复杂应力状态下超细晶材料的塑性
2.金属材料在双拉应力状态下塑性变形过程的声发射检测
3.复杂应力状态下高强混凝土的损伤本构理论及应用
4.复杂应力状态下金属材料屈服模型及有限元验证
因版权原因，仅展示原文概要，查看原文内容请购买。

偶应力介质结构的最优强度拓扑优化

偶应力介质结构的最优强度拓扑优化苏文政【摘要】在等应变能密度分布的意义下给出了偶应力介质结构最优强度的拓扑优化设计方法.优化模型的设计变量为单元的密度,约束函数为许用材料的体积用量,目标函数通过等应变能密度准则隐式地给出.该方法的优点在于将强度约束/目标的显式应力形式进行转化,避免了应力函数优化问题中的奇异性、强非线性及约束的局部性等困难.数值算例表明,基于优化模型,偶应力介质的最优结果依赖于结构的最小局部尺寸与偶应力介质特征长度的比值.当宏观结构的尺寸远大于特征长度时,偶应力介质的结果趋于经典连续介质的相应结果.【期刊名称】《科学技术与工程》【年(卷),期】2014(014)035【总页数】3页(P157-159)【关键词】偶应力;强度;拓扑优化;准则法【作者】苏文政【作者单位】大连交通大学土木与安全工程学院,大连116028【正文语种】中文【中图分类】O342偶应力连续介质理论[1]由于包含了反映材料自身微结构尺度特征的材料特征长度参数，能够解释由于材料的微结构和宏观结构的尺寸相近而产生的尺寸效应，近年来广泛应用于超轻多孔固体以及多相复合材料的等效性能表征[2, 3]。

基于偶应力等效连续介质模型的结构分析和设计，需要新的实现方法。

同经典连续介质相比，偶应力连续介质具有高阶特点，如偶应力介质认为材料质点既具有平动自由度也具有转动自由度，质点之间既能够传递力也能够传递力偶。

这种高阶特点使得偶应力介质结构的分析和设计方法更为复杂。

目前，对于经典连续介质，结构拓扑优化已经成为结构设计领域的重要方法之一。

在结构的初始设计阶段，拓扑优化能够提供具有革命性的概念设计产品。

20世纪80年代，程耿东、Kikuchi、Bendsøe等人通过引入微结构的概念，最终创立了当代拓扑优化的主流方法——微结构法[4—7]。

此后，该领域也相继出现了一些其他的方法，如ESO方法[8]、ICM方法[9]以及水平集方法[10]等。

弹塑性本构模型理论

加工硬化
当材料中的应力状态处于某一个屈服面上时，如果因加荷 使它发生超越这个屈服面的应力变化，就会在材料中同时 引起新的弹性与塑性变形，形成新的屈服面。加荷使屈服 面膨胀、移动或改变形式，这些改变取决于材料的应力历 史与应力水平，这种现象称为加工硬化（软化）
等向硬化：屈服面大小不同
运动硬化：屈服面位置发生移动

剑桥模型
物态边界面
正常固结的饱和重塑黏土的孔隙比e和它所受的 力p与q之间存在一种固定关系，这一关系反映在 e-p-q空间中就形成了物态边界面
原始各向等压固结线AC（VICL）
在p
1
2


条件下的
3
e

p曲线
VICL表达式： e ea0 ln p
VICL回弹曲线：
Mises破坏条件
f
*

J2

k2 f
Mohr-Coulomb 破坏条件

cn
tan
f Drucker-Prager *
破坏条件
I1
J2 kf
屈服面：
定义：
特征
理想简单塑性材料：材料进入屈服状态，就可以认为材料 破坏了，屈服面与破坏面重合
加工硬化材料：屈服应力随荷载的提高与变形的增大而提 高，因此屈服面不同于破坏面，不是一种固定的面

1 2

3
应力不变量
3 I12 I2 I3 0 I1 1 2 3 I2 1 2 2 3 31
I3 1 2 3
偏差应力
sij ij ij (I1 / 3)
偏差应力不变量
E-V弹性模型 K-G弹性模型 南京水科所模型 剑桥模型 KW模型 LD模型 罗威剪胀模型
  1. 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

文 章 编 号 ：１６７１－６９０６（２０１９）０３－０００６－０５
基于偶应力理论的高阶弹塑性本模型的无网格法
王瑞昌，孙玉周
（中原工学院 建筑工程学院，河南 郑州 ４５０００７）
摘 要： 高阶连续理论包含材料内禀特征长度，可以反映材料尺度效应或微尺度特征对宏观性能的影 响，但 是，高 阶 连 续理论需要考虑位移的高阶导数，给数值模拟带来很大困 难。利 用 无 网 格 法 能 够 方 便 构 造 具 有 高 阶 连 续 特 征 形 函 数 的 优点，建立基于偶应力理论的高阶弹塑性本构模型的无网格 法，对 二 维 悬 臂 梁 弯 曲 变 形 进 行 数 值 模 拟，并 分 析 了 结 构 变 形中的尺度效应。结果表明：无网格法计算结果与有限元软件 ＡＮＳＹＳ计算结果相 吻 合，且 尺 度 因 子 对 结 构 变 形 有 一 定 的影响。 关 键 词 ： 高 阶 理 论 ；弹 塑 性 ；尺 度 效 应 ；无 网 格 法 ；数 值 模 拟 中 图 分 类 号 ： Ｏ３０２ 文 献 标 志 码 ： Ａ ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．１６７１－６９０６．２０１９．０３．００２
第３期
王 瑞 昌 ，等 ：基 于 偶 应 力 理 论 的 高 阶 弹 塑 性 本 构 模 型 的 无 网 格 法
·７·
１ 基 于 偶 应 力 理 论 的 高 阶 弹 塑 性 本 构 模 型
１．１ 本 构 模 型
偶应力理论的几何方程为
ε犻犼
＝
１ ２
（狌犻，犼
＋狌犼，犻），χ犻犼
＝ ω犼，犻
（１）
室 开 挖 过 程 进 行 了 研 究 。 ［８］ 整 体 上 说 ，由 于 数 值 离 散 时要 求位 移 函 数 满 足 Ｃ１ 连 续 性 ，高 阶 连 续 理 论 、特 别是高阶弹塑性理论的数值模拟仍然是工程数值仿 真领域的研究热点。
无网格法是一种 新 型 数 值 计 算 方 法，与 传 统 有 限 元 法 相 比 ，具 有 形 函 数 满 足 高 阶 连 续 特 征 、不 受 网 格 限 制等特点，特别适合 大 变 形、位 移 不 连 续、裂 纹 扩 展 和 高速碰撞等问题的数值模拟 。 ［９］ 本文尝试用无网 格 法 建立基于偶应力理论的高阶弹塑性本构模型的数值计 算框架 ，用 Ｆｏｒｔｒａｎ语言编写计算程序，对悬臂梁进行 数值模拟，以研究无 网 格 法 在 高 阶 弹 塑 性 问 题 数 值 模 拟中的应用。
ＷＡＮＧ Ｒｕｉｃｈａｎｇ，ＳＵＮＹｕｚｈｏｕ．Ｍｅｓｈｆｒｅｅｍｅｔｈｏｄｆｏｒｔｈｅｈｉｇｈｏｒｄｅｒｅｌａｓｔｉｃｐｌａｓｔｉｃｃｏｎｓｔｉｔｕｔｉｖｅｍｏｄｅｌｂａｓｅｄｏｎｔｈｅｃｏｕｐｌｅｓｔｒｅｓｓｔｈｅｏｒｙ［Ｊ］． ＪｏｕｒｎａｌｏｆＺｈｏｎｇｙｕａｎ ＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ，２０１９，３０（３）：６－１０（ｉｎＣｈｉｎｅｓｅ）．
收 稿 日 期 ：２０１９－０３－１０ 基 金 项 目 ：国 家 自 然 科 学 基 金 项 目 （１１４７２３１６）；河 南 省 高 校 科 技 创 新 人 才 支 持 计 划 项 目 （１６４２００５１００２０） 引 文 格 式 ：王 瑞 昌 ，孙 玉 周 ．基 于 偶 应 力 理 论 的 高 阶 弹 塑 性 本 构 模 型 的 无 网 格 法 ［Ｊ］．中 原 工 学 院 学 报 ，２０１９，３０（３）：６－１０．
近年来 ，各 种 特 殊 材 料 以 及 微／纳 米 量 级 新 型 构 件被 越 来 越 多 地 应 用 到 工 程 实 际 中 ，偶 应 力 理 论 和 ［１］ 应变 梯 度 理 论 等 ［２］ 高 阶 连 续 理 论 被 广 泛 用 在 微／纳 米 结构或者需 要 考 虑微尺 度 影响 的 宏观 结 构的 研 究 中［３］。ＢＯＲＳＴ 等 率 先 提 出 了 基 于 Ｖｏｎ－ Ｍｉｓｅｓ屈 服 准 则 的 Ｃｏｓｓｅｒａｔ弹 塑 性 理 论［４］；ＦＬＥＣＫ 等 分 别 从 统 计 存 储 位 错 和 几 何 必 需 位 错 两 个 角 度 出 发，改 进 了 Ｃｏｓｓｅｒａｔ弹 塑 性 模 型［５］；李 锡 夔 等 基 于 Ｄｒｕｃｋｅｒ－ Ｐｒａｇｅｒ屈服 准 则 给 出 了 Ｃｏｓｓｅｒａｔ弹 塑 性 理 论 的 塑 性 一 致 性 算 法［６］；冀 宾 等 基 于 偶 应 力 弹 塑 性 理 论 ，对 软 化行 为 的 剪 切 带 问 题 进 行 了 有 限 元 数 值 模 拟［７］；张 建 成 等 利 用 非线性 Ｃｏｓｓｅｒａｔ扩展模型对层状 岩体的洞
式中：ε犻犼为常规应 变，狌犻、狌犼 为 位 移，下 标 “，”表 示 对 其 后坐标求偏导，χ犻犼为微曲率，ω犼 为微转角。
第 ３０ 卷 第 ３ 期 ２０１９ 年 ６ 月
中原工学院学报 ＪＯＵＲＮＡＬ ＯＦＺＨＯＮＧＹＵＡＮ ＵＮＩＶＥＲＳＩＴＹ ＯＦＴＥＣＨＮＯＬＯＧＹ
Ｖｏｌ．３０ Ｎｏ．３ Ｊｕｎ．２０１９
孙 玉 周 ：男 ，汉 族 ，１９７４ 年 生 ，本 科 、硕 士 毕 业 于 兰 州 大 学 力 学 系 ，博 士 毕 业 于 香 港 城 市 大 学 建 筑 系 。 现 任 中 原 工 学 院 建 筑 工 程 学 院 教 授 ，研 究 生 处 处 长 。 主 要 研 究 方 向 为 微 纳 米 力 学 和 计 算 力 学 。 主 持 国 家 自 然 科 学 基 金 项 目 ２ 项 ，省 （部 ）级 科 研 项 目 ５ 项 ，发 表 学 术 论 文 ７０ 多 篇 ，其 中 ＳＣＩ收 录 ３０ 多 篇 ，出 版 学 术 专 著 １ 部 。 先 后 被 评 为 全 国 优 秀 力 学 教 师 、河 南 省 教 育 系 统 优 秀 教 师 、河 南 省 高 校 青 年 骨 干 教 师 、河 南 省 教 育 厅 学 术 技 术 带 头 人 ，入 选 河 南 省 高 校 科 技 创 新 人 才 、河 南 省 科 技 创 新 杰 出 人 才 和 教 育 部 新 世 纪 优 秀 人才。