类橡胶材料变形直到破坏的显式本构模型

本文通过直接、显式的方法提出一个多轴可压缩应变能弹性势来模拟类橡胶材料受载荷直到软化破坏的变形行为．首先，我们提出一个多轴可压缩应变能函数；其次，通过特定的不变量，该多轴应变能函数在单轴拉伸，平面应变和等双轴拉伸三个基准实验的情况下，可以退化为各自的单轴形函数形式；再次，我们显式给出带有软化破坏特性的形函数；最后，模型结果和试验数据可以精确匹配，同时可以预测材料临近破坏以后，接下来的变形行为．     

各向异性半平面问题的复应力函数基本解

应用复变函数Ｃａｕｃｈｙ积分的方法，对于各向异性半平面边界为自由或固定两种情形，分别导出了其在任意集中力或集中力偶作用下的复应力函数基本解，其特例与前人结果一致     

显式方法精确模拟形状记忆合金循环荷载下的变形行为

本文提出全新的有限弹塑性J2流方程,用来显式、精确地模拟SMAs（形状记忆合金）材料在循环加载-卸载条件下从塑性逐渐转变为伪弹性的变形行为．首先,改进流动法,使得本构方程耦合屈服中心的移动和屈服面的增大,并改进背应力演化方程,使模型可以产生强烈的包辛格效应,从理论上具备模拟SMAs独特变形行为的能力;其次,构造全过程下的统一硬化函数显式表达式,代入本构方程后能得到符合要求的形函数;再次,利用选定的数据点构造统一光滑的上屈服函数,再利用上下屈服应力之间的一种线性关系,推导得到下屈服阶段的形函数;最后,只需要给定一个参数就可以得到单个循环结果,利用拉格朗日插值方法构建参数随循环次数变化的函数,就可以模拟任意循环荷载下的变形行为．通过模型结果和实验数据对比证明新方法的有效性．     

雷诺应力各向异性涡黏模型的层析TRPIV测量

利用层析TRPIV测量水洞中平板湍流边界层3D-3C速度场的高分辨率时间序列数据库. 提出了空间局部平均多尺度速度结构函数的新概念, 描述湍流多尺度涡结构的空间拉伸、压缩、剪切变形和旋转. 用空间局部平均多尺度速度结构函数对湍流脉动速度进行了空间多尺度分解. 用空间流向局部平均多尺度速度结构函数, 根据湍流多尺度涡结构在流向的拉伸和压缩物理特征, 提出了新的湍流相干结构条件采样方法, 检测并提取了层析TRPIV数据中相干结构“喷射”和“扫掠”事件中的脉动速度、平均速度变形率、雷诺应力等物理量的空间拓扑形态. 通过研究平均速度变形率各分量与雷诺应力各分量之间的空间相位差异,肯定了壁湍流相干结构雷诺应力各向异性复涡黏模型的合理性.      

应力历史对颗粒材料强度和变形特征影响的离散元模拟

主要关注了颗粒材料前期所受的应力历史对其后期宏观力学响应的影响。该应力历史由一段等比例加载应力路径以及卸载垂直方向应力至与水平方向围压相同的卸载段描述。具体工作为:基于PFC2D双轴压缩数值实验,调查了应力历史对颗粒样本的强度、变形特征、细观参量如组构的影响,得出颗粒样本的偏应力-应变、体积应变曲线及其发生破坏时的名义应变云图。数值结果表明:高低围压下样本分别发生剪切破坏和弥散破坏,高围压下弹性阶段的刚度受应力历史影响较大,而低围压下样本在刚进入塑性至应力峰值点阶段的弹塑性刚度变化较为明显。随着应力历史中等比例加载系数的增大,剪胀加快,变形局部化范围有所不同;另一方面,颗粒形状的不规则性会增强颗粒材料的各向异性,导致样本强度更高。     

基于有限元模拟的SUS304奥氏体不锈钢摩擦耦合变形过程的应力分析

采用ANSYS/LS-DYNA软件对SUS 304奥氏体不锈钢薄板的摩擦耦合变形过程进行了数值模拟.采用隐式-显式序列求解法和分段线性塑性材料模型,分析了钢带摩擦耦合变形时的应力分布规律及载荷、下压量和滑动速度等因素对钢带剪应力、主应力及等效应力的影响.结果表明:摩擦耦合变形的试验参数显著影响钢带的应力分布,验证了钢带在低于其屈服强度的应力条件下发生塑性变形的摩擦诱发效应.对奥氏体不锈钢摩擦诱发马氏体转变行为的研究及其摩擦学性能的改善具有一定的指导意义.     

显式模拟类橡胶材料Mullins效应滞回圈

通过显式、直接的方法提出一个多轴可压缩应变能函数,用来模拟类橡胶材料在加载——卸载作用下,由于Mullins效应而产生的应力——应变滞回圈. 本文的创新点在于将表征能量耗散的变量引入到应变能函数.新的弹性势具有以下两个特点:第一,在加载情况下,新引入的变量不会对弹性势产生任何影响,因此,只要给出合适的形函数显式表达,3个基准实验,包括单轴拉伸和压缩,等双轴拉伸和压缩,以及平面应变,都可精确模拟;第二,新引入的变量在卸载情况下将被激活.在不同的卸载应力下,变量将发生改变,从而影响弹性势,使其最终产生不同的应力——应变关系卸载曲线,与对应的加载曲线共同构成应力——应变滞回圈.通过对Mullins效应实验数据进行分析和研究,得出了卸载形函数在不同卸载应力下变化的规律,并预测不同卸载应力下的应力——应变关系.最后,我们将得到精确匹配实验数据的数值模拟结果,从而证明本文方法不仅可以精确匹配至少3个基准实验,还可以模拟和预测类橡胶材料在加载——卸载作用下由于Mullins效应而产生的滞回圈.      

DAMAGE MODEL FOR HARDENING AND SOFTENING MATERIAL IN GENERAL STATE OF STRESS

The present paper is aimed to simulate progression of damage,hardeningand softening response in brittle materials such as concrete or rock in general state ofstress.Similar shape of surfaces for yield,failure and damage progressing areavailable,and softening strain is treated as plasticity.Then,the proposed model isapplied to solving several boundary value problems.     

THE RATE-INDEPENDENT CONSTITUTIVE MODELING FOR POROUS AND MULTI-PHASE NANOCRYSTALLINE MATERIALS

To determine the time-independent constitutive modeling for porous and multiphase nanocrystalline materials and understand the effects of grain size and porosity on their mechanical behavior, each phase was treated as a mixture of grain interior and grain boundary, and pores were taken as a single phase, then Budiansky's self-consistent method was used to calculate the Young's modulus of porous, possible multi-phase, nanocrystalline materials, the prediction being in good agreement with the results in the literature. Further, the established method is extended to simulate the constitutive relations of porous and possible multi-phase nanocrystalline materials with small plastic deformation in conjunction with the secant-moduli approach and iso-strain assumption. Comparisons between the experimental grain size and porosity dependent mechanical data and the corresponding predictions using the established model show that it appears to be capable of describing the time-independent mechanical behaviors for porous and multi-phase nanocrystalline materials in a small plastic strain range. Further discussion on the modification factor, the advantages and limitations of the model developed were present.     

土的各向同性化变换应力方法

在不同方向的力学参数、结构特性及应力应变关系的不同为材料的各向异性, 建立能够反映这种复杂特性的强度准则、本构模型, 对于材料本构关系的研究具有重要的理论意义. 但材料的各向异性一直是其力学特性描述的难点, 对此, 郑泉水院士提出了各向同性化定理, 为后续研究解决材料的各向异性问题提供了方向及思路. 作者等针对土的应力诱导各向异性提出了变换应力方法, 这种方法同样遵循对材料进行各向异性问题各向同性化处理的思路, 与郑泉水院士的各向同性化定理是一脉相承的, 也是对各向同性化定理的发展. 本文旨在通过分析各向同性化定理与变换应力方法明确两者间的内在联系, 并以土材料的应力诱导各向异性处理为例, 说明在具体材料的各向异性处理过程中面临的现实问题以及变换应力方法是如何解决这些问题的. 分析并给出了变换应力方法应用时的三个合理假设, 推导出了具体的变换应力数学公式, 阐明了在考虑土的应力诱导各向异性的具体函数已经给出的情况下, 在构造土的弹塑性本构模型中采用变换应力方法的必要性.      

弹塑性材料的平面应力非连续分岔

基于平面应力非连续分岔特性的一般描述,运用统一强度理论,得出了非相关流动情形的弹塑性材料平面应力非连续分岔的起始方位角以及相应的最大硬化模量的统一解析解,并且分析了材料拉压异性以及不同程度的中间应力对结果的影响,进而发现所得的结果一强度准则的选取有关,揭示了在分岔研究中正确选取符合材料特性的强度准则的重要性。最后,同特线理论比较发现平面应力剪切带型非连续分岔同平面应力特征线重合。     

循环软化45碳钢和循环硬化304不锈钢的棘轮行为实验研究

对循环软化45碳钢的单轴应力循环下的平均应力、应力幅值以及先前应变循环对棘轮效应的影响进行了实验研究;并对循环硬化的304不锈钢进行了多种非比例循环加载路径下路径形状、路径等效应力幅值、平均应变与平均应力对材料棘轮变形行为的影响实验.发现平均应力和应力幅值及其历史对于材料的棘轮行为都有很大的影响.     

基于深度强化学习算法的颗粒材料应力?应变关系数据驱动模拟研究

颗粒材料的宏观力学行为受颗粒组分等材料参数, 孔隙率、配位数等状态参数的影响, 同时又具备复杂的加载路径和加载历史相关性, 建立包含多个内变量以及各变量间相互关联的颗粒材料本构模型是一个重要的科学难题. 不同于传统的基于屈服面、流动法则和硬化函数框架下的唯象本构模型, 本文基于颗粒物质力学的研究基础, 以颗粒材料平均孔隙率、细观组构参数和弹性刚度参数作为内变量, 结合深度学习方法建立以有向图表征的数据本构模型. 有向图中以不同的链接网络表示不同的内变量信息流动方向, 各个内变量间的映射关系采用循环神经网络来建立, 将各个神经网络相互组合, 形成包含不同内变量且具有不同预测能力的本构模型. 该本构模型的建立过程等价于在众多可能的内变量链接关系空间中寻找最能描述实际材料宏观应力应变行为的优化问题. 因此, 可将有向图本构模型的建立过程看作“马尔可夫决策过程”, 采用深度强化学习算法构建有向图的内变量链接组合优化过程, 具体采用AlphaGo Zero算法自动寻找最优的颗粒材料数据驱动本构模型建模路径. 研究结果表明, 采用有向图和深度强化学习算法可建立起完全依靠“数据驱动”的颗粒材料应力?应变关系. 此外, 本方法提供了一种将不同理论模型从数据角度统一起来, 且基于人工智能算法发展更优模型的研究思路, 可为相似问题的研究提供借鉴.      

水库诱震机制新理论的探索──断层带弱化与岩体软化效应诱震理论

针对走滑断层诱发的水库地震，建立了尖点突变模型，提出了发震的充要力学条件判据及地震释放能量的简单表达式。以突变理论分析为依据，建立了水库诱震机制的新假说，简称为断层带弱化与岩体软化效应诱震理论，较合理地解释了水库地震与岩性、构造、水位的关系及水库地震序列ｂ值的特征。本文的讨论从整体上深化了对水库诱震机制的认识。     

层合材料的非傅里叶热传导及热应力

论文建立了一个双层材料层合板受瞬态加热情况下的非傅里叶热传导分析模型,用向后差分法得到了温度场的数值解,并对该差分格式的稳定性进行了讨论.给出了温度场随导热时间、热扩散率、空间与时间步长之比以及弛豫时间的变化趋势.同时,通过已经求得的温度场,求得了层合板内的应力场,给出了层合板内的热应力随时间的变化.     

大型液体火箭姿控与跷振大回路耦合动力学建模与分析

大型液体火箭结构模态的空间化分布特征导致结构振动、姿态运动和推进系统液路脉动存在相互耦合,进而影响传统姿控回路的稳定性. 针对大型液体火箭, 充分考虑姿态控制系统对箭体姿态动力学和弹性振动的影响, 以及箭体结构弹性振动与推进系统的耦合作用(跷振(POGO)), 建立了姿控与跷振大回路耦合模型. 该模型包含了推进系统、结构系统与姿控系统之间的耦合因素, 可进行姿控-结构-推进大回路耦合机理研究. 该模型具有非奇异的优点, 可以直接用于频域分析和时域仿真. 基于该模型研究了我国某型号液体捆绑火箭推进系统参数——泵增益和蓄压器能量值对姿态运动与结构振动稳定性的影响. 研究得出, 泵增益和蓄压器能量值的变化不仅导致了结构振动的不稳定, 而且也导致了姿态运动的发散. 因此, 对于大型液体捆绑火箭, 推进系统与姿控系统之间存在不可忽略的耦合作用, 在设计姿控系统时, 有必要考虑推进系统对姿控系统稳定性的影响.