高温下混凝土的本构模拟及破坏分析

针对已建立的高温下混凝土中化学-热-水力-力学耦合过程分析的分级数学模型,发展了混凝土的化学-热-水力-力学(CTHM)耦合本构模型。在已有的Willam-Warnke弹塑性屈服准则基础上发展了考虑脱水和脱盐引起的材料损伤及化学塑性软化、塑性应变硬化/软化和吸力硬化的广义Willam-Warnke本构模型,模拟高温下混凝土的材料非线性行为。为保证全局守恒方程的Newton迭代过程的二阶收敛率,导出了非线性化学-热-水力-力学(CTHM)耦合本构模型的一致性切线模量矩阵。数值结果显示了本文所发展的化学-热-水力-力学(CTHM)耦合本构模型在模拟高温下混凝土中复杂破坏过程的能力和有效性。     

A new thermodynamically consistent continuum model for hardening plasticity coupled with damage

A phenomenological model for hardening–softening elasto-plasticity coupled with damage is presented. Specific kinematic internal variables are used to describe the mechanical state of the system. These, in the hypothesis of infinitesimal changes of configuration, are partitioned in the sum of a reversible and an irreversible part. The constitutive equations, developed in the framework of the Generalised Standard Material Model, are derived for reversible processes from an internal energy functional, postulated as the sum of the deformation energy and of the hardening energy both coupled with damage, while for irreversible phenomena from a dissipation functional.Performing duality transformations, the conjugated potentials of the complementary elastic energy and of the complementary dissipation are obtained. From the latter a generalised elastic domain in the extended space of stresses and thermodynamic forces is derived. The model, which is completely formulated in the space of actual stresses, is compared with other formulations based on the concept of effective stresses in the case of isotropic damage. It is observed that such models are consistent only for particular choices of the damage coupling. Finally, the predictions of the proposed model for some simple processes are analysed.     

A generalized coupled viscoplastic-viscodamage-viscohealing theory for glassy polymers

A thermodynamic consistent, small-strain, non-unified model is developed to capture the irregular rate dependency included in the strain controlled inelastic responses of polymers at the glassy state. The model is considered as a generalized Frederick-Armstrong-Philips-Chaboche (FAPC) theory proposed by [Voyiadjis and Basuroychowdhury, 1998] and [Voyiadjis and Abu Al-Rub, 2003] which is based on a von Mises and Chaboche isotropic hardening type viscoplasticity formulation. Using the proposed model, different experimental results are simulated and the range of viscoplastic related material constants are obtained through a parametric study. The thermodynamic framework is used to incorporate the effect of coupling between viscodamage and viscohealing phenomena into the inelastic deformation of glassy polymers. This coupling effect is crucial for polymeric based self healing systems in which different damage mechanisms are active and the efficiency of the healing processes are highly dependent on the damage. The computational aspect for general coupled inelastic-damage-healing processes together with the required solution algorithms are elaborated and the inelastic-damage-healing response of a polymeric based self-healing system is simulated. The proposed viscoplasticity theory constitutes a physically consistent approach to model the irregular mechanical responses of glassy polymers and the viscodamage model provides an exquisite predicting tool to evaluate the ductile damage associated with the large inelastic deformation and low cycle fatigue in polymeric based material systems. In conclusion, a well structured viscohealing theory is formulated for polymeric based self healing systems.     

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

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

Thermodynamic-based cohesive zone healing model for self-healing materials

This work presents a thermodynamic-based cohesive zone framework to model healing in materials that tend to self-heal. The nominal, healing and effective configurations of continuum damage-healing mechanics are extended to represent cohesive zone configurations. To incorporate healing in a cohesive zone model, the principle of virtual power is used to derive the local static/dynamic macroforce balance and the boundary traction as well as the damage and healing microforce balances. A thermodynamic framework for constitutive modeling of damage and healing mechanisms of cracks is used to derive the evolution equations for the damage and healing internal state variables. The effects of temperature, resting time, crack closure, history of healing and damage, and level of damage on the healing behavior of the cohesive zone are incorporated. The proposed model promises solid basis for understanding the self-healing phenomena in self-healing materials.     

Shape memory behaviour: modelling within continuum thermomechanics

A phenomenological material model to represent the multiaxial material behaviour of shape memory alloys is proposed. The material model is able to represent the main effects of shape memory alloys: the one-way shape memory effect, the two-way shape memory effect due to external loads, the pseudoelastic and pseudoplastic behaviour as well as the transition range between pseudoelasticity and pseudoplasticity.The material model is based on a free energy function and evolution equations for internal variables. By means of the free energy function, the energy storage during thermomechanical processes is described. Evolution equations for internal variables, e.g. the inelastic strain tensor or the fraction of martensite are formulated to represent the dissipative material behaviour. In order to distinguish between different deformation mechanisms, case distinctions are introduced into the evolution equations. Thermomechanical consistency is ensured in the sense that the constitutive model satisfies the Clausius–Duhem inequality.Finally, some numerical solutions of the constitutive equations for isothermal and non-isothermal strain and stress processes demonstrate that the various phenomena of the material behaviour are well represented. This applies for uniaxial processes and for non-proportional loadings as well.     

温度梯度对热致型形状记忆聚合物板力学性能的影响

形状记忆聚合物作为一种新型的智能材料,由于质量轻、成本低以及变形回复率大等优势,已经在航空、医学等领域得到广泛应用。当前对于热致型形状记忆聚合物力学行为的研究,大都集中在整体变温的情况下,随着应用环境的越来越广泛,温度梯度对材料力学性能的影响效果越发重要。本文在均匀应力的假设下,给出材料在不同初始和传热条件下的温度梯度分布情况,结合传热学和热致型形状记忆聚合物相变理论本构模型,分别讨论了不同温度梯度对存储应变、弹性模量等力学性能的影响,通过理论分析和实验数据对比验证了模型正确性。本文研究可为不同导热情况下,对热致型形状记忆聚合物力学行为监测提供思路,也为形状记忆聚合物的进一步工程应用提供理论依据。     

On the constitutive modeling of coupled thermomechanical phase-change problems

This paper deals with a thermodynamically consistent numerical formulation for coupled thermoplastic problems including phase-change phenomena and frictional contact. The final goal is to get an accurate, efficient and robust numerical model, able for the numerical simulation of industrial solidification processes. Some of the current issues addressed in the paper are the following. A fractional step method arising from an operator split of the governing differential equations has been used to solve the nonlinear coupled system of equations, leading to a staggered product formula solution algorithm. Nonlinear stability issues are discussed and isentropic and isothermal operator splits are formulated. Within the isentropic split, a strong operator split design constraint is introduced, by requiring that the elastic and plastic entropy, as well as the phase-change induced elastic entropy due to the latent heat, remain fixed in the mechanical problem. The formulation of the model has been consistently derived within a thermodynamic framework. All the material properties have been considered to be temperature dependent. The constitutive behavior has been defined by a thermoviscous/elastoplastic free energy function, including a thermal multiphase change contribution. Plastic response has been modeled by a J2 temperature dependent model, including plastic hardening and thermal softening. The constitutive model proposed accounts for a continuous transition between the initial liquid state, the intermediate mushy state and the final solid state taking place in a solidification process. In particular, a pure viscous deviatoric model has been used at the initial fluid-like state. A thermomecanical contact model, including a frictional hardening and temperature dependent coupled potential, is derived within a fully consistent thermodinamical theory. The numerical model has been implemented into the computational finite element code COMET developed by the authors. Numerical simulations of solidification processes show the good performance of the computational model developed.     

A finite element model for shape memory alloys considering thermomechanical couplings at large strains

In the present work we propose a new thermomechanically coupled material model for shape memory alloys (SMA) which describes two important phenomena typical for the material behaviour of shape memory alloys: pseudoelasticity as well as the shape memory effect. The constitutive equations are derived in the framework of large strains since the martensitic phase transformation involves inelastic deformations up to 8%, or even up to 20% if the plastic deformation after the phase transformation is taken into account. Therefore, we apply a multiplicative split of the deformation gradient into elastic and inelastic parts, the latter concerning the martensitic phase transformation. An extended phase transformation function has been considered to include the tension–compression asymmetry particularly typical for textured SMA samples. In order to apply the concept in the simulation of complex structures, it is implemented into a finite element code. This implementation is based on an innovative integration scheme for the existing evolution equations and a monolithic solution algorithm for the coupled mechanical and thermal fields. The coupling effect is accurately investigated in several numerical examples including pseudoelasticity as well as the free and the suppressed shape memory effect. Finally, the model is used to simulate the shape memory effect in a medical foot staple which interacts with a bone segment.     

Thermomechanical behavior of shape memory alloy under complex loading conditions

The thermo-mechanical behavior of polycrystalline shape memory alloy (SMA) under multi-axial loading with varying temperature conditions has been studied by experiments. Recently the research has been extended theoretically and a mechanical model of polycrystalline SMA and the corresponding mesoscopic constitutive equations have been developed. The model presented in this paper is constructed on the basis of the crystal plasticity and the deformation mechanism of SMA. The variants in the crystal grains and the orientations of crystal grains in the polycrystal are considered in the proposed model; the constitutive equations are derived on the basis of the proposed model. The volume fraction of the martensite variants in the transformation process and the influence of the stress state on the transformation process are also considered. Some calculated results obtained by the constitutive equations are presented and compared with the experimental results. It is found that the deformation behavior of SMA under complex loading conditions can be well reproduced by the calculation of the constitutive equations.     

复杂加载下混凝土的弹塑性本构模型

混凝土材料在不同应力路径下或复杂加载条件下会表现出差异性显著的应力应变关系,在小幅循环加载条件下,其应力应变关系会表现出类似于弹性变形的滞回曲线.在不同应力水平下,混凝土的应力应变关系以及破坏特性都具有静水压力相关特点,即随着静水压力增大,各向异性强度特性弱化.此外,混凝土受压及受拉破坏机理不同,因而对应于混凝土硬化损伤亦有不同,即可分为受压硬化损伤,受拉硬化损伤及两者的混合硬化损伤类型.基于Hsieh模型,对该模型进行了三点改进.（1）针对小幅循环加载下混凝土无塑性变形的试验规律,而模型中在应力水平较低的循环加载条件下始终存在塑性变形的预测问题,采用在边界面模型框架下,设置了应力空间的弹性域,初始屈服面与后续临界状态屈服面几何相似的假定.（2）基于广义非线性强度准则将原模型采用变换应力方法将其推广为三维弹塑性本构模型,采用变换后模型可合理的考虑不同应力路径对于子午面以及偏平面上静水压力效应形成的影响,并避免了边界面应力点奇异问题.（3）分别对拉压两种加载损伤模式建议了相应的硬化参数表达式,可分别用于描述上述加载中产生的应变软化及强度退化行为.基于多种加载路径模拟表明:所建立的三维弹塑性本构模型可合理地用于描述混凝土的一般应力应变关系特性.      

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

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

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

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

A crystallographic model for the study of local deformation processes in polycrystals

Most engineering materials possess a polycrystalline structure. Under load the anisotropy of the constituent grains leads to strong inhomogeneities of stresses and strains on the grain level. In order to investigate the local deformation processes, a new crystallographic model for pure fcc metals in the low temperature range has been developed. It is based on the framework of crystal plasticity and uses the finite element method (FEM). The rate dependent constitutive equations consider isotropic as well as kinematic hardening, whereby the mutual interactions of dislocation processes on the different slip systems are taken into account. Comprehensive calculations show that the essential features of both single crystals—which serve as a test object for the constitutive equations—and polycrystals are reproduced correctly. Moreover the simulations allow a deeper understanding of the mechanisms that control the local deformation behaviour of metals, especially of the mutual interactions of slip system activity, local hardening and resulting local strain. Furthermore, the model may serve as a physically motivated base for a later inclusion of damage terms which allow investigations of damage and fatigue on the local scale.     

Constitutive modeling of shape memory polymer based self-healing syntactic foam

In a previous study, it was found that the shape memory functionality of a shape memory polymer based syntactic foam can be utilized to self-seal impact damage repeatedly, efficiently, and almost autonomously [Li G., John M., 2008. A self-healing smart syntactic foam under multiple impacts. Comp. Sci. Technol. 68(15–16), 3337–3343]. The purpose of this study is to develop a thermodynamics based constitutive model to predict the thermomechanical behavior of the smart foam. First, based on DMA tests and FTIR tests, the foam is perceived as a three-phase composite with interfacial transition zone (interphase) coated microballoons dispersed in the shape memory polymer (SMP) matrix; for simplicity, it is assumed to be an equivalent two-phase composite by dispersing elastic microballoons into an equivalent SMP matrix. Second, the equivalent SMP matrix is phenomenologically assumed to consist of an active (rubbery) phase and a frozen (glassy) phase following Liu et al. [Liu, Y., Gall, K., Dunn, M.L., Greenberg, A.R., Diani J., 2006. Thermomechanics of shape memory polymers: uniaxial experiments and constitutive modeling. Int. J. Plasticity 22, 279–313]. The phase transition between these two phases is through the change of the volume fraction of each phase and it captures the thermomechanical behavior of the foam. The time rate effect is also considered by using rheological models. With some parameters determined by additional experimental testing, the prediction by this model is in good agreement with the 1D test result found in the literature. Parametric studies are also conducted using the constitutive model, which provide guidance for future design of this novel self-healing syntactic foam and a class of light-weight composite sandwich structures.     

准脆性材料的弹塑性各向异性损伤耦合本构模型研究

针对准脆性材料的非线性特征：强度软化和刚度退化、单边效应、侧限强化和拉压软化、不可恢复变形、剪胀及非弹性体胀,在热动力学框架内,建立了准脆性材料的弹塑性与各向异性损伤耦合的本构关系。对准脆性材料的变形机理和损伤诱发的各向异性进行了诠释,并给出了损伤构形和有效构形中各物理量之间的关系。在有效应力空间内,建立了塑性屈服准则、拉压不同的塑性随动强化法则和各向同性强化法则。在损伤构形中,采用应变能释放率,建立了拉压损伤准则、拉压不同的损伤随动强化法则和各向同性强化法则。基于塑性屈服准则和损伤准则,构建了塑性势泛函和损伤势泛函,并由正交性法则,给出了塑性和损伤强化效应内变量的演化规律,同时,联立塑性屈服面和损伤加载面,给出了塑性流动和损伤演化内变量的演化法则。将损伤力学和塑性力学结合起来,建立了应变驱动的应力-应变增量本构关系,给出了本构数值积分的要点。以单轴加载-卸载往复试验识别和校准了本构材料常数,并对单轴单调试验、单轴加载-卸载往复试验、二轴受压、二轴拉压试验和三轴受压试验进行了预测,并与试验结果作了比较,结果表明,所建本构模型对准脆性材料的非线性材料性能有良好的预测能力。     

A thermodynamics based damage mechanics model for particulate composites

A micro-mechanical damage model is proposed to predict the overall viscoplastic behavior and damage evolution in a particle filled polymer matrix composite. Particulate composite consists of polymer matrix, particle fillers, and an interfacial transition interphase around the filler particles. Yet the composite is treated as a two distinct phase material, namely the matrix and the equivalent particle-interface assembly. The CTE mismatch between the matrix and the filler particles is introduced into the model. A damage evolution function based on irreversible thermodynamics is also introduced into the constitutive model to describe the degradation of the composite. The efficient general return-mapping algorithm is exploited to implement the proposed unified damage coupled viscoplastic model into finite element formulation. Furthermore, the model predictions for uniaxial loading conditions are compared with the experimental data.     

A constitutive theory for shape memory polymers. Part II: A linearized model for small deformations

A constitutive theory is developed for shape memory polymers. It is to describe the thermomechanical properties of such materials under large deformations. The theory is based on the idea, which is developed in the work of Liu et al. [2006. Thermomechanics of shape memory polymers: uniaxial experiments and constitutive modelling. Int. J. Plasticity 22, 279-313], that the coexisting active and frozen phases of the polymer and the transitions between them provide the underlying mechanisms for strain storage and recovery during a shape memory cycle. General constitutive functions for nonlinear thermoelastic materials are used for the active and frozen phases. Also used is an internal state variable which describes the volume fraction of the frozen phase. The material behavior of history dependence in the frozen phase is captured by using the concept of frozen reference configuration. The relation between the overall deformation and the stress is derived by integration of the constitutive equations of the coexisting phases. As a special case of the nonlinear constitutive model, a neo-Hookean type constitutive function for each phase is considered. The material behaviors in a shape memory cycle under uniaxial loading are examined. A linear constitutive model is derived from the nonlinear theory by considering small deformations. The predictions of this model are compared with experimental measurements.