关于《弹性-蠕变体理论的广义变分原理》的几点讨论意见

邬瑞锋先生在《弹性-蠕变体理论的广义变分原理》中给出了四个广义变分原理的泛函: 小位移线性弹性-蠕变体理论广义变分原理: 大位移线性弹性-蠕变体理论广义变分原理: 文章中,作者列出了两个等式: 这等于说,I_1所对应的变分原理与I_2所对应的等价;而I_3所对应的变分原理与I_4所对应的等价,这是一个明显的错误。公式(1)中的两个等式并不是无条件成立的恒等式,而是在一定的条件下成立的,这个条件就是泛函I_1和I_3中的应力和应变函数满足弹性蠕变体的本构方程。     

混凝土单轴受拉的非局部本构模型

混凝土受拉本构行为存在很强的局部软化现象,使得单轴受拉试验无法给出应力-应变关系,而只能给出应力-位移关系。本文根据内变量理论和等效应变假设建立了基于真实应变的混凝土单轴受力本构方程,并根据Weibull分布可以描述混凝土等脆性材料断裂过程的试验现象,建立了关于弹性应变的损伤演化规律。然后,通过假设平均应变与真实弹性应变的函数关系,在应力-平均应变的本构关系中采用平均弹性应变以描述其非局部行为,而在材料的损伤演化规律中采用真实弹性应变以描述其局部行为,由此建立了单轴受拉荷载条件下的非局部本构模型。最后,对一个单调受拉试验和一个反复受拉试验的仿真结果表明所提出的非局部本构模型可以准确地模拟试验结果。     

求解非线性大变形的超级单元法

引言分析大位移大应变非线性问题,其基本关系存在着三个困难,第一,位移和应变关系是非线性函数,第二,应力和应变关系也是非线性函数,如橡皮这样的材料,应力应变关系是复杂的非线性函数,一般可用三个应变不变量I_1,I_2,I_3的应变能函数表示出来.第三,如橡皮这种材料,在变形过程中,体积没有明显的变化.某些金属材料,如软钢,在发生大塑性变形时,可以略去弹性应变不计,根据位错理论,塑性变形是由于沿晶格面的滑移而引起的,它是一种畸变,只引起形状的改变而不导致体积的变化,     

各向同性率无关材料本构关系的不变性表示

在内变量理论的框架下,针对各向同性率无关材料,使用张量函数表示理论建立了塑性应变全量及增量本构关系的最一般的张量不变性表示. 它们均由3个完备不可约的基张量组合构成,这3个基张量分别是应力的零次幂、一次幂和二次幂. 因此得出,塑性应变、塑性应变增量与应力三者共主轴. 通过对基张量的正交化,给出了本构关系式在主应力空间中的几何解释. 进一步,全量(或增量)本构关系中3个组合因子被表达为应力、塑性应变(或塑性应变增量)的不变量的函数. 当塑性应变(或塑性应变增量)的3个不变量之间满足一定关系时,所给出的本构关系将退化为经典的形变理论(或塑性势理论).最后,还讨论它与奇异屈服面理论的关系,当满足一定条件时,两者是一致的.      

高应变率下钛合金Ti-6Al-4V的热-粘塑性特性和绝热剪切变形

本文用SHPB 试验技术和金相观察相结合的方法研究了高应变率下钛合金 Ti6-6 Al-4 V的热-粘塑性本构关系和绝热剪切变形。发现用基于双曲型势垒热激活机制的热-粘塑性本构方程（1+（-s）0）-1=1-a0Tln（ep/es）能够很好地表述其应力、应变、应变率和温度间的关系。理论分析和试验还表明Ti-6Al-4V在室温下的绝热剪切临界条件是既与应变有关也与应变率有关的双变量准则,这和前人提出的各种单变量准则是有所区别的。     

有限应变下蠕变—弹塑性—损伤耦合本构模型的数值模拟

本文提出了一个模拟有限元应变下各向同性蠕变－弹塑性－损伤耦合本构行为的数值模型。模型基于不可逆过程的热动力学和内状态变量理论、以及应力就变本构关系的超弹性形式描述和变形梯度的弹性及非弹性部分的乘式分解。     

显式方法精确模拟形状记忆聚合物热力学行为

通过构建一个热耦合的多轴可压缩应变能函数,得到应力-应变、应力-温度和应变-温度之间的函数关系,建立形状记忆聚合物的本构方程.本文引入三个基于对数应变的不变量使得模型（i）可以模拟可压缩情况;（ii）适用于单轴拉伸和等双轴拉伸至少两个基准实验;（iii）多轴有效.通过显式方法（i）给出自由能和熵的具体表达,证明模型热力学定律;（ii）给出应变-应力,温度-应力以及,温度-应变的形函数具体表达.多轴模型在特定的情况下可以自动退化到各自的单轴情况. 通过调节形函数的参数,最终得到的模型结果和实验结果能够精确匹配.新方法建立的本构模型得到的结果能更加准确地指导形状记忆聚合物的工程设计。     

显式方法精确模拟形状记忆聚合物热力学行为

通过构建一个热耦合的多轴可压缩应变能函数，得到应力-应变、应力-温度和应变-温度之间的函数关系，建立形状记忆聚合物的本构方程.本文引入三个基于对数应变的不变量使得模型（i）可以模拟可压缩情况；（ii）适用于单轴拉伸和等双轴拉伸至少两个基准实验；（iii）多轴有效.通过显式方法（i）给出自由能和熵的具体表达，证明模型热力学定律；（ii）给出应变-应力，温度-应力以及，温度-应变的形函数具体表达.多轴模型在特定的情况下可以自动退化到各自的单轴情况. 通过调节形函数的参数，最终得到的模型结果和实验结果能够精确匹配.新方法建立的本构模型得到的结果能更加准确地指导形状记忆聚合物的工程设计。     

形状记忆合金伪弹性行为模拟

基于对伪弹性形状记忆合金(SMA)典型应力-应变曲线的特征分析,在原Graesser本构模型中增加简洁多项式来描述SMA应力诱发马氏体相变完成后在变形马氏体相下继续加载阶段的变形特征;并引入应变幅值与混相下SMA弹性模量的关系来改进不同应变幅值下卸载时SMA的应力-应变关系,从而提出了一种新的SMA一维本构关系模拟其伪弹性力学行为。该模型对直径为0.5mm的NiTi合金丝的拉伸加载、卸载试验曲线的模拟结果表明:改进本构模型与原Graesser模型相比,其能够准确地模拟SMA在不同应变幅值下加载和卸载应力-应变关系。此外,通过研究SMA本构模型的物理关系,推导出了控制SMA滞回曲线特征的关键参数fT与相变临界应力、弹性常数之间的明确关系,可利用该关系直接确定参数fT,摆脱了只靠试算获取该参数的传统做法,其准确性得到了试验验证。     

基于深度学习和细观力学的颗粒材料本构关系研究

颗粒材料的本构关系对岩土工程等众多领域至关重要. 不同于传统的唯象本构理论, 本文基于机器学习模型探索了一种细观力学理论指导下的数据驱动型颗粒材料本构关系预测方法. 根据Vogit均质化假设, 建立了小应变条件下颗粒材料应力?应变解析关系, 此关系唯一地确定了一组与颗粒材料本构行为相关的细观组构变量. 这些变量与反应颗粒材料宏观性质的主应力和主应变信息通过一系列离散元三轴压缩数值试验获得. 考虑到细观组构变量为内变量, 不能直接作为本构模型的输入. 本文基于有向图方法将颗粒材料微观结构信息隐式地包含在应力?应变的预测当中, 并采用门控循环单元(GRU)循环神经网络作为基础深度学习模型描述有向图中结点之间的映射关系. 通过将有向图从目标节点沿源节点展开, 整个应力?应变预测模型可由两个神经网络分别训练并组装而成. 将训练后的深度学习模型在全新的数据集上进行测试, 结果表明该训练策略能有效捕捉到颗粒材料在常规三轴任意加卸载, 等中主应力系数b的真三轴加载, 和等平均有效应力p的真三轴加卸载等复杂多轴加载工况下的应力?应变响应关系, 模型具有良好的内插和外推预测能力. 考虑到深度学习模型捕捉颗粒材料力学响应的能力及其开放式学习的特点, 充分结合数据驱动方法和理论本构模型可能是颗粒材料本构研究的一个重要方向.      

材料的弹粘塑性本构关系和应力波的演化

试验表明，大多数工程材料在冲击载荷作用之下的变形一般都同时包含有可恢复的瞬态性弹性变形和不可恢复的粘滞性塑性变形，即其本构关系可以用弹粘塑性模型来描述。本文从内变量理论出发，探讨了时率相关材料的弹粘塑性本构关系的一般特性，建立了增量型的弹粘塑性本构关系的一般理论框架和普适的表达式，并且对两种最常用的本构模型——Bodner-Partom模型和Johnson-Cook模型给出了在一维应变条件下的具体形式。通过计算和讨论一维应变粘塑性靶板中冲击波的衰减机制和应力波的演化规律，特别是考察各种粘塑性本构模型中的材料参数对冲击波的衰减和应力波的演化的影响，得出了一些可以直接应用或具有一定借鉴价值的结果，为研究应力波的其他衰减机制以及在人防工程中智能防护层设计时新材料的选取奠定了基础。     

微磁检测应力和塑性区的磁弹塑耦合理论

金属磁记忆微磁检测方法, 利用铁磁材料局部磁性状态的变化, 进行应力集中或塑性区域位置及程度的检测与评价. 面向微磁信号的定量理论分析可对其工程领域应用提供重要指导. 本文介绍铁磁材料微弱环境磁场下的磁弹塑性本构进展, 及其在微磁信号分析方面的应用. 力磁本构关系方面, 针对微磁检测弱磁化条件, 基于有效场理论构建了受弹塑性载荷铁磁材料的理想磁化本构的显式解析式, 并结合接近原理分析了恒定外加微弱磁场下应力-应变对材料磁化强度的影响. 检测信号分析方面, 基于弹性力学理论、静磁学理论和新建立的磁弹塑性本构关系, 建立并求解了微弱磁场下铁磁试件中弹性应力或塑性区诱导的表面磁信号的二维分析模型. 结合实验结果证实其在刻画弹塑性因素对微磁信号影响规律方面的能力, 并详细分析了微磁信号的特征量与局部弹性应力或塑性区的尺寸间的相关关系. 相比已有力磁本构关系, 本文建立的显式解析形式的理想磁化更加简洁, 有助于提升对力磁耦合效应的定量化理解和应用.      

一种孔隙介质力学模型及在水泥基材料的应用

基于细观力学和断裂力学的基本理论提出一个新的分析模型, 对孔隙介质的力学性能进行了分析. 依据孔隙介质内部孔隙的几何描述和状态参数,如孔隙率、形状、尺度及分布等,通过等效夹杂理论获得孔隙介质的等效本构方程,其最终变量为应力、应变和孔隙的形态参数. 根据断裂理论中材料承受载荷作用下破坏增长过程中的能量守恒,对孔隙介质变形过程中机械能、弹性应变能和载荷提供的势能进行分析, 根据能量守恒定律建立能量守恒方程,其最终变量也为应力、应变和孔隙的形态参数. 根据等效本构方程和能量守恒方程,获得孔隙介质承受载荷作用下的应力应变关系. 最后将该力学模型应用于水泥基材料,计算水泥基材料的力学性能并与文献中的结果进行对比分析,结果显示模型的计算结果准确有效.      

家兔肌腱的蠕变变形行为的实验研究

研究了家兔肌腱的蠕变规律，测定了稳定蠕变速率等力学参量，通过对试验结果的分析，建立了肌腱的蠕变本构关系，结果表明，蠕变应变与瞬时应力和时间有关，在一定的瞬时应力作用下，稳定蠕变速率近似为常数，稳定蠕变速率与瞬时应力近似成指数关系。     

聊聊动态塑性和黏塑性

固体力学研究者致力于具有应力-应变本构关系（以下简称为形变型本构关系）的变形体的力学响应研究，而流体力学研究者致力于具有应力-应变率本构关系（以下简称为流动型本构关系）的流动体的力学响应研究。当涉及结构和材料的动态塑性时，到底应该用“塑性变形”还是“塑性流动”来表示？本文从宏观塑性本构理论和微观位错动力学机理两个角度，分别讨论并指出塑性本构关系属于流动型黏塑性率相关本构关系，且同时适用于加载和卸载；因而不应该用应力-应变图来描述塑性加-卸载过程。弹塑性本构关系则是一种形变型和流动型本构关系的耦合。     

高应变率下航空透明聚氨酯的动态本构模型

采用低阻抗分离式霍普金森压杆对航空透明聚氨酯进行了高应变率下的动态力学性能测试,得到的应力应变曲线表现出了显著的非线性黏弹性特征。基于本构理论和实验数据,构建了航空透明聚氨酯的松弛时间应变相关的超黏弹性本构形式。该本构模型由2部分组成:一部分表征准静态下的超弹性行为,另一部分描述非线性应变率的相关特性。利用超黏弹性本构模型对不同应变率下航空透明聚氨酯的动态应力应变曲线进行拟合,拟合曲线与实验曲线一致性良好。     

有机玻璃在高应变率下的损伤型非线性粘弹性本构关系及破坏准则

对两种有机玻璃高应变率下的大变形和破坏行为进行了实验研究,通过改进文献[2]本构关系的非线性弹性项并引入损伤参量,建立了一个适用于更大变形范围、能描述应力平台及本构失稳的损伤型非线性粘弹性本构方程。相应地,从临界损伤量概念出发,提出以应变和应变率为控制变量的破坏准则。不论是本构关系还是破坏准则,理论计算均与试验结果吻合良好。     

基于塑性体积应变的岩石损伤变形特性实验研究

为分析岩石塑性变形与损伤的关系,在定义岩石的初始损伤和临界损伤,提出塑性体积应变分析方法,从而以塑性体积应变为损伤变量,采用归一化方法建立岩石的损伤本构模型。采用递增循环加载实验确定岩石损伤本构模型中的弹性卸载模量和弹性应变比例系数两个参数。通过实验和理论分析得出:当荷载较小时,普通单轴压缩状态下岩石损伤随荷载的增加具有减小趋势,荷载超过一定数值后,岩石损伤才开始增加;单轴递增循环压缩状态下当循环荷载大于约35%峰值强度后,卸载后岩石的损伤具有增加的趋势,小于该荷载之前具有减小的趋势。整个加载过程的理论应力-应变曲线能很好地与实验结果相吻合,在循环加载区间理论结果还能体现出岩石实验结果中的回滞环。     

超弹性镍钛形状记忆合金单轴相变棘轮行为的宏观唯象本构模型

超弹性镍钛形状记忆合金因其良好的力学性能以及独特的超弹性和形状记忆效应已广泛应用于土木工程、航空航天和生物医疗等多个领域,在实际服役环境中超弹性镍钛合金元件不可避免地会承受不同应力水平的循环载荷作用,亟待建立描述相变棘轮行为(即峰值应变和谷值应变随着正相变和逆相变循环的进行不断累积)的循环本构模型.为此,基于已有的超弹性镍钛形状记忆合金在不同峰值应力下的单轴相变棘轮行为实验研究结果,在广义黏塑性框架下,对Graesser等提出的通过背应力非线性演化方程反映超弹性镍钛形状记忆合金超弹性行为的一维宏观唯像本构模型进行了拓展,考虑了正相变和逆相变过程中特征变量的差异及其随循环的演化,以非弹性应变的累积量为内变量引入了正相变开始应力、逆相变开始应力、相变应变和残余应变的演化方程,同时通过峰值应力与正相变完成应力的比值来确定演化方程中的相关系数,建立了描述超弹性镍钛合金单轴相变棘轮行为的本构模型.将模拟结果与对应的实验结果进行对比发现,建立的宏观唯像本构模型能够合理地描述超弹性镍钛形状记忆合金的单轴相变棘轮行为及其峰值应力依赖性,模型的预测结果和实验结果吻合得很好.     

再论弹性力学中的广义变分原理——就等价定理问题和胡海昌先生商榷

本文就胡海昌先生提出的等价定理的论争,申述个人的观点和论证,与胡海昌先生商榷。 本文主要论证了下列三点: (1)通过待定的拉格朗日乘子法证明了胡海昌-鹫津久一郎原理(下文简称胡鹫原理)的三类变量之间并不独立,应力应变关系仍然是应力和应变之间应该首先满足的变分约束条件。这个变分原理只是在形式上有应力、应变、位移三类变量,在实际上,这些变量中只有两类是独立的。(2)通过高阶拉格朗日乘子法,我们求得了比胡海昌鹫津久一郎原理的泛函更一般形式的具有三类变量的变分泛函,而且证明有无穷个这样的变分泛函。利用唯一性定理,我们证明了这些泛函的变量中必须满足应力应变关系这个条件。同样也证明了胡鹫原理并不是三类变量都独立的和没有任何约束条件的完全的变分原理,而是一个以应力应变关系为变分约束条件的变分原理。(3)在应力应变关系的变分约束条件下,我们证明了Hellinger-Reissner原理和胡鹫原理的等价定理。 本文的结论是:等价定理是正确的,并非象胡海昌先生所指的那样是“误解”。郭仲衡、戴天民、陈至达、刘殿魁、张其洁、邬瑞铎、奚肖风等通过各自的努力,在各种变分问题上论证了等价定理,都是正确的,没有什么“误解”,更没有“误入迷途”。胡海昌先生认为大家都有“误解”的原因,似乎在于