海冰动力学中本构模型研究的若干进展

从20世纪70年代, 人们在不同尺度下建立了一系列的海冰动力学本构模型用于海冰数值模拟和预测. 将目前应用的海冰动力学本构模型分为弹塑性、黏塑性、各向异性和颗粒流体动力学中的黏弹塑性模型4类,并分别讨论了各类模型的特点和适用范围. 尤其对在极区及副极区大、中尺度下广泛应用的黏塑性及其改进的本构模型进行了深入地分析. 最后得出: 在大、中尺度下建立海冰动力学的黏弹塑性本构模型是提高海冰动力学计算精度的有效途径; 将小尺度下采用黏弹塑性本构模型的海冰颗粒流体动力学的计算时效进一步改进后, 可在大、中尺度下对海冰的重叠、堆积特性进行模拟; 进一步开展海冰动力学的尺度效应研究, 进行海冰本构模型的实验验证, 并建立不同尺度模型间的相互联系也应是海冰动力学本构模型研究的重要内容.      

心肌细胞大变形黏弹性模型及在实验中的应用

在衡量单个细胞力学行为的研究中,越来越多地采用结合实验的数值模拟方法. 在连续介质力学框架下,发展了一种新的心肌细胞本构模型,并与微管吮吸实验结合,探讨了心肌细胞的力学特性. 本构模型是对普遍使用的仅能用于小变形分析的标准线性固体模型的一种扩展,它将超弹性性能引入到黏弹性模型中,用以描述细胞的大变形黏弹性效应. 基于改进的本构模型,对心肌细胞微管吮吸实验过程进行了有限元模拟,并将计算结果与实验结果以及经典理论解进行了对比. 结果显示发展的本构模型适合细胞大变形问题的有限元数值模拟.      

沥青混合料的粘弹性疲劳损伤模型研究

从沥青混合料的黏弹性本构模型出发，按应变等效假设 将损伤引进本构模型，推导周期应力条件下含损伤的耗散能表达式，进而得到黏弹性损伤演 化方程和疲劳方程；该模型可以综合反映应力、加载频率和温度的影响，适合用来分析沥青 路面的疲劳损伤问题；利用MTS开展了沥青混合料的疲劳损伤试验，确定了相应的黏弹性疲 劳损伤模型参数，据此开展了模型验证，结果表明模型与实测值之间相关性很好.     

黏弹夹芯板的有限元建模及实验研究

从有限元分析和数值模拟及实验验证的角度研究了黏弹夹芯板的频率依赖振动特性。夹芯板中间层为黏弹性材料，其刚度和阻尼的频率依赖性行为直接影响系统的模态频率和阻尼，并导致振动模式求解的复杂化。采用三阶七参数Biot模型描述黏弹性材料频率相关的黏弹性行为。开发了三层四节点28自由度的夹芯板单元，基于经典板理论和哈密顿原理建立了黏弹夹芯板的有限元动力学方程。通过引入辅助耗散坐标，将Biot模型和黏弹夹芯板的有限元动力学模型结合起来，并将其转化为常规二阶线性系统形式，极大简化了求解非线性振动特性的过程。对一边固定、另三边自由的黏弹夹芯板进行了前三阶固有频率和损耗因子的预测，并与实验结果对比。数值模拟结果和实验结果吻合良好，说明所提有限元方法是正确有效的。     

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

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

基于实验的心肌被动力本构模型

在狗心脏外科手术和心肌力学性能试验的基础上,对心肌力学性能实验数据进行处理和拟合,并应用有限元方法对心肌的单拉实验、松弛实验进行数值模拟和计算. 得到不同心肌位点的本构模型的非线性弹性和黏性参数,从而发展了心肌被动力的非线性黏弹性本构模型,体现了心脏力学行为的非均匀性和黏弹性. 在理论模型和试验数据的基础上,开发了有限元软件ABAQUS的VUMAT材料子程序,应用于数值心脏的仿真研究.      

基于CDM的黏弹性薄板压屈失效分析

研究黏弹性薄板的压屈失效问题.基于连续介质损伤力学与Schapery对应原理,得到耦合各向异性蠕变损伤的黏弹性本构模型,构造板的分层模型以考虑损伤演化的空间特征,建立具初挠度板的非线性压屈控制方程与边界条件,并进行空间与时间域上的数值离散和迭代求解.计算结果与实验曲线的对比表明本文的分析方法是有效的,黏弹性材料的蠕变损伤将导致结构刚度不断削弱,在一定应力水平下,板发生变形快速增加的压屈失效.     

基于统一强度理论的岩石弹塑性本构模型及其数值实现

基于弹塑性力学理论，以统一强度准则为屈服准则，建立了考虑硬化/软化行为和应变率效应的岩石弹塑性本构模型；采用Fortran语言通过LS-DYNA的用户自定义材料接口(Umat)对该弹塑性本构模型进行编程，并把该程序生成求解器以达到对该模型进行应用的目的；通过岩石的单轴压缩实验和SHPB实验对所建的弹塑性本构模型进行验证，结果表明，该弹塑性本构模型能够反映岩石在准静态和动态下的力学行为。     

广义Maxwell黏弹性流体在两平板间的非定常流动

将分数阶微积分运算引入Maxwell黏弹性流体的本构方程，研究了黏弹性流体在两平板问的非定常流动．对于广义Maxwell黏弹性流体的分数阶导数模型，导出了对时间具有分数阶导数的特殊运动方程，利用分数阶微积分的Laplace变换理论，得到了流动的解析解．     

黏弹-Perzyna黏塑性有限元法应力更新隐式算法

黏弹-黏塑性耦合模型的黏弹性部分由弹簧、黏壶和Kelvin链串联而成,黏塑性部分为双曲线型DruckerPrager屈服函数、各向同性硬化和Perzyna黏塑性流动模型。基于黏弹性蠕变柔度,通过定义与弹性问题相对应的与时间增量相关的黏弹性剪切模量和体积模量,导出增量递推形式的本构方程。为保证算法的收敛和稳定性,把Perzyna黏塑性流动方程转化为与弹塑性相似的一致性条件,建立黏塑性增量因子单侧逼近其收敛值的N-R迭代算法。最后,给出应力更新完全隐式算法和最终计算公式。分别采用黏弹性、黏弹-塑性和黏弹-黏塑性本构关系对一地基蠕变模型进行三维有限元分析和比较,结果表明,本文算法具有较高的计算效率和稳定性。     

粘弹塑性统一本构模型理论

文章在已有的统一本构模型的基础上,将粘弹性变形引入到统一本构模型之中,成功地改善了材料过渡段的变形模拟情况。通过Hastelloy-X的变形模拟及与其它统一本构模型的变形模拟比较,证明了粘弹塑性统一本构模型的合理性。     

高应力高水压下砂岩三轴蠕变特性试验研究

对砂岩进行高围压高水压条件下的三轴压缩蠕变试验。试验表明,在整个加载过程中,孔隙水压力主要起到增强轴向变形和横向变形的作用;但在加载的初始阶段,孔隙水压力在一定程度上抑制了轴向变形。当应力水平高于屈服应力时,横向蠕变速率明显大于轴向蠕变速率,且横向蠕变率先进入加速蠕变阶段。本文提出一个新的非线性黏性元件,并引入一个能判断是否进入加速蠕变阶段的计时器元件,组建一个非线性黏塑性加速蠕变启动模型,将该黏塑性模型与Burgers模型串联,构建一个新的非线性黏弹塑性蠕变模型,推导了该模型在常规三轴应力状态下的本构方程。基于试验结果,通过对非线性优化算法(BFGS)的Matlab编程,实现对本文提出蠕变模型的参数识别,识别效果比较理想。对比试验曲线与拟合曲线,二者相当吻合,验证了新提出的非线性黏弹塑性蠕变模型的正确性。     

混凝土材料冲击特性的研究

基于混凝土材料强冲击加载下的试验研究，提出了两种损伤型动态本构模型： 损伤型黏弹性本构模型和损伤与塑性耦合的本构模型. 通过模型计算结果与冲击试验结果的 比较可发现，随着冲击速度的提高，混凝土材料内部产生了显著的塑性变形，由此 损伤型黏弹性本构模型的应用就存在一些不足. 而损伤与塑性耦合的本构模型由于考虑了裂 纹扩展引起的材料强度和刚度的弱化，以及微空洞缺陷塌陷引起的塑性变形，因而能更好地 用于模拟强冲击载荷作用下混凝土材料的冲击响应特性.     

高应变率大变形下的聚丙烯/尼龙共混高聚物损伤型本构特性

对两种采用不同相容剂的聚丙烯（PP）和尼龙（PA）共混高聚物材料在大变形下的粘弹性力学行为进行研究,着重考察应变率效应和损伤的演化,从而分析不同的界面分子设计对共混体系材料宏观性能的作用。在准静态及冲击实验研究的基础上,基于ZWT非线性粘弹性模型,并结合了遗传算法,分别得到了能有效描述两种共混高聚物大变形阶段计及损伤的非线性粘弹性本构关系。两种材料在不同加载条件下表现出明显不一致的性能,原因在于其损伤演化的率相关性,且两种材料的大变形机制存在一定的差别,能用ZWT方程进行描述的范围也不一样。     

模拟微可压粘弹性流体的WCCBS_SU方法

针对微可压缩粘弹性流动问题,发展了微可压缩流的WCCBS方法,详细推导了基于Oldroyd-B本构模型的WCCBS_SU方法的求解过程。在流场微可压的条件下,分别对平面Poiseuille流和4：1粘弹性收缩流进行了数值模拟。Poiseuille流在不同We数下数值结果与解析解的比较,验证了本文方法具有较高的精度和较好...     

索网天线的参变量变分及非线性有限元方法

针对大型周边桁架式索网天线由拉索拉压模量不同引起的本构非线性和结构大变形引起的几何非线性问题,给出了基于参变量变分原理的几何非线性有限元方法. 首先针对含预应力索单元拉压模量不同分段描述的本构关系,通过引入参变量,导出了基于参变量及其互补方程的统一描述形式,避免了传统算法需要根据当前变形对索单元张紧/松弛状态的预测,提高了算法收敛性. 然后利用拉格朗日应变描述索网天线结构大变形问题,结合几何非线性有限元法,建立了基于参变量的非线性平衡方程和线性互补方程;并给出了牛顿-拉斐逊迭代法与莱姆算法相结合的求解算法. 数值算例验证了本文提出的算法比传统算法具有更稳定的收敛性和更高的求解精度,特别适合于大型索网天线结构的高精度变形分析和预测.      

THE PARAMETRIC VARATIONAL PRINCIPLE AND NON-LINEAR FINITE ELEMENT METHOD FOR ANALYSIS OF ASTROMESH ANTENNA STRUCTURES

针对大型周边桁架式索网天线由拉索拉压模量不同引起的本构非线性和结构大变形引起的几何非线性问题,给出了基于参变量变分原理的几何非线性有限元方法. 首先针对含预应力索单元拉压模量不同分段描述的本构关系,通过引入参变量,导出了基于参变量及其互补方程的统一描述形式,避免了传统算法需要根据当前变形对索单元张紧/松弛状态的预测,提高了算法收敛性. 然后利用拉格朗日应变描述索网天线结构大变形问题,结合几何非线性有限元法,建立了基于参变量的非线性平衡方程和线性互补方程;并给出了牛顿-拉斐逊迭代法与莱姆算法相结合的求解算法. 数值算例验证了本文提出的算法比传统算法具有更稳定的收敛性和更高的求解精度,特别适合于大型索网天线结构的高精度变形分析和预测.     

A fractional non-linear creep model for coal considering damage effect and experimental validation

Accurate prediction of coal׳s creep behavior is of great significance to coalbed methane extraction. In this study, taking into account the visco-elastic–plastic characteristics and the damage effect, a fractional non-linear model is proposed to describe the creep behavior of coal. The constitutive and creep equations of the proposed fractional non-linear model are derived via the Boltzmann superposition principle and discrete inverse Laplace transform. Furthermore, uniaxial creep tests under different axial stress conditions were carried out to validate the proposed model. It is found that the present model can describe the experimental data from creep tests with better accuracy than classical models. Particularly, the present model can predict the accelerating creep deformation of coal which classical models fail to reproduce. Finally, the parametric sensitivity analysis is performed to investigate the effects of model parameters on the creep strain. It is verified that the introduction of fractional parameters and damage factor in the present model is essential to accurate prediction of the full creep stage of coal.