化学-力学耦合理论与数值方法

该文研究了化学场中的质量扩散与力学耦合问题,构造了化-力耦合情况下的力学本构关系与质量扩散的本构关系,并由这些本构关系和化学场、力学场的控制方程,得到化-力场耦合的有限元方程.通过数值算例,详细分析了由应力场引起的质量重分布和由化学场引起的结构变形.研究表明,力学与化学之间存在明显的相互作用,并采用有限元数值方法进行了分析.     

Modeling of plasticity and aging as coupled phenomena

This paper concerns the study of the interaction between aging and plasticity in a 2024 aluminum-copper alloy. The framework of thermodynamics of irreversible processes is used to introduce the effect of aging in the constitutive equations by means of an internal variable. These constitutive equations are checked by comparing the prediction of uniaxial tests with the experimental responses. They are also introduced into a finite element code, and the evolution of notched structures shows that the coupling between aging and plasticity can strongly modify the distribution of the stresses.     

Formulation of thermo-hydro-mechanical coupling behavior of unsaturated soils based on hybrid mixture theory

Thermo-Hydro-Mechanical （THM） coupling pro- cesses in unsaturated soils are very important in both theoretical researches and engineering applications. A coupled formulation based on hybrid mixture theory is derived to model the THM coupling behavior of unsaturated soils. The free-energy and dissipative functions for different phases are derived from Taylor＇s series expansions. Constitutive relations for THM coupled behaviors of unsaturated soils, which include deformation, entropy change, fluid flow, heat conduction, and dynamic compatibility conditions on the interfaces, are then established. The number of field equations is shown to be equal to the number of unknown variables; thus, a closure of this coupling problem is established. In addition to modifications of the physical conservation equations with coupling effect terms, the constitutive equations, which consider the coupling between elastoplastic deformation of the soil skeleton, fluid flow, and heat transfer, are also derived.     

Restrictions on nonlinear constitutive equations for elastic shells

Constitutive equations for the resultant forces and moments applied to a shell-like body necessarily couple the influences of the shell geometry and the constitutive nature of the three-dimensional material from which the shell is constructed. Consequently, even when the nonlinear constitutive equation of the three-dimensional material is known, the complicated influence of the shell geometry on the constitutive response of the shell is not known. The main objective of this paper is to develop restrictions on the constitutive equations of nonlinear elastic shells which ensure that exact solutions of the shell equations are consistent with exact nonlinear solutions of the three-dimensional equations for homogeneous deformations. Since these restrictions are nonlinear in nature they provide valuable general theoretical guidance for specific constitutive assumptions about the coupling of material and geometric properties of shells. Examples of the linear theories of a plate and a spherical shell are considered.     

Restrictions on nonlinear constitutive equations for elastic rods

Constitutive equations for the resultant forces and moments applied to a rod-like body necessarily couple the influences of the rod geometry and the constitutive nature of the three-dimensional material from which the rod was constructed. Consequently, even when the nonlinear constitutive equation of the three-dimensional material is known, the influence of the rod geometry on the constitutive response of the rod is not known. The main objective of this paper is to develop restrictions on the constitutive equations of nonlinear elastic rods which ensure that exact solutions of the rod equations are consistent with exact nonlinear solutions of the three-dimensional equations for all homogeneous deformations. Since these restrictions are nonlinear in nature they provide valuable general theoretical guidance for specific constitutive assumptions about the coupling of material and geometric properties of rods. Also, an example of a straight beam clamped at one end and subjected to a shear force at the other end is used to examine the validity of the proposed value for the transverse shear deformation coefficient.     

一种化学-力学耦合连续介质理论与数值格式

对化学驱动的连续介质化学-力学耦合系统进行研究,从热力学定律和化学势角度出发,推导了等温过程的化学-力学耦合本构关系和控制方程,利用变分方法建立了化学-力学耦合系统的能量泛函,得到化学-力学耦合控制方程的等效积分形式和相应的有限元列式. 结合算例,对连续介质的化学-力学耦合行为进行了数值计算,数值结果反映了化学与力学系统的相互耦合作用,即浓度变化能引起介质的变形,同样力学作用也能引起浓度重分布. 从全新的角度建立了描述连续介质的化学-力学耦合行为的基本理论和数值方法,能够较好地反映一类连续介质的化学-力学耦合行为.      

无耦合条件下的多子系统静动态统一本构关系(上)--基本理论

引入了同步性的概念。按照是否同步对所有的内变量进行分组，得到了整个系统的子系统结构。证明了子系统间无耦合时各子系统内部一致性条件成立的判别准则，从而得到了可以统一描述子系统内部的时间相关、时间无关不可逆行为的内变量演化方程和加卸载准则。本构关系既适用于静态加载，也适用于动态加载。经典塑性本构关系、多屈服面本构关系及过应力型粘塑性本构关系只是本文的特例。     

无耦合条件下的多子系统静动态统一本构关系(下)--铝的双子系统静动态统一本构模型

构造了标量形式的无耦合条件下双子系统静动态统一本构模型.推导出第1和第2子系统中加载应变速率临界值ε·c1和ε·c2,当应变速率ε·分别低于和高于某个临界值时,相应子系统中的不可逆行为分别是与时间无关的和与时间相关的.由于当ε·跨越ε·c1和ε·c2时,内变量的求解公式发生变化,所以动态强度随ε·变化的规律发生变化.经与铝的实验结果比较确认,本构模型能够描述材料的多种静动态力学行为.     

Non-linear constitutive relations for magnetostrictive materials

In this paper, non-linear deformation behavior of magnetostrictive materials is studied and three magnetoelastic coupling constitutive models are developed. The standard square (SS) constitutive model is developed by means of truncating the polynomial expansion of the Gibbs free energy. The hyperbolic tangent (HT) constitutive equations, which involve a hyperbolic tangent magnetic-field dependence, are proposed to model the magnetic-field-induced strain saturation of magnetostrictive materials in the region of intense magnetic fields. A new model based on density of domain switching (DDS) is established in terms of the basic truth that magnetic domain switching underlies magnetostrictive deformation. In this model, it is assumed that the relation between density of domain switching, defined by the quantity of magnetic domains switched by per unit magnetic field and magnetic field can be described by a density function with normal distribution. The moduli in these constitutive models can be determined by a material function that is proposed to describe the dependence of the peak piezomagnetic coefficient on the compressive pre-stress for one-dimensional cases based on the experimental results published. The accuracy of the non-linear constitutive relations is evaluated by comparing the theoretical values with experimental results of a Terfenol-D rod operated under both compressive pre-stress and bias magnetic field. Results indicate that the SS constitutive equations can accurately predict the experimental results under a low or moderate magnetic field while the HT model can, to some extent, reflect the trend of saturation of magnetostrictive strain under a high magnetic field. The model based on DDS, which is more effective in simulating the experimental curves, can capture the main characteristics of the mechanism of magnetoelastic coupling deformation of a Terfenol-D rod, such as the notable dependence of magnetoelastic response on external stress and the saturation of magnetostrictive strain under intense magnetic fields. In addition, the SS constitutive relation for a general three-dimensional problem is discussed and an approach to characterize the modulus tensors is proposed.     

Basic Equations of Kelvin's Medium and Analogy with Ferromagnets

A nonlinear elastic Cosserat continuum of special kind whose particles have dynamical spins (Kelvin's medium) is considered. The basic equations for this medium are obtained on the basis of fundamental mechanical laws. Various complete sets of strain tensors are discussed. The analogies between Kelvin's medium and other types of media are discussed. First, we show the analogy between constitutive equations for Kelvin's medium and for elastic shells. The main difference is caused by the dimensionality of the media. Then, for a special case of fast proper rotation we establish the analogy between the dynamic equations of Kelvin's medium and elastic ferromagnetic saturated insulators. The most general way to take into account the coupling between magnetic and elastic subsystems is presented. All results can be interpreted in terms of a mechanical medium as well as in terms of ferromagnets. This revised version was published online in July 2006 with corrections to the Cover Date.     

计及热效应和湿分扩散的耦合黏弹性本构方程

高聚物在电子和航空等领域得到广泛运用,由于高聚物的亲水性,常常发生由于湿热引起的结构失效甚至材料的断裂.近年来,有实验显示,高聚物中的湿热效应及其力学反响是相互影响的,同时,考虑到高聚物的黏弹性,发展新的包含湿热效应的黏弹性本构方程来描述该问题是必要的.本文基于不可逆热力学的基本原理,运用连续介质力学的方法,通过引入内变量表示高聚物的黏性效应,基于 Helmholtz自由能导出一种热、湿分和黏弹性力学性质三场耦合的本构关系和系统控制方程,对实际的分析应用有较强的指导意义.     

铁电陶瓷微观电畴演化的成核率模型

大量的实验已经证实电畴翻转是铁电材料非线性和迟滞性本构曲线的根本原因.根据这些描述电畴运动的物理实验结果,提出了成核率模型来模拟铁电材料的非线性本构行为.进一步由模型建立了基于微观电畴运动的电畴翻转体积分数演化方程.应用该模型得到的理论计算结果与实验数据的比较表明,模型能够描述铁电材料的非线性本构行为.同时模型所揭示的微观电畴成核的演化行为可进一步的指导宏观唯象模型的建立与改进.     

Overall dynamic constitutive relations of layered elastic composites

A method for the homogenization of a layered elastic composite is presented. It allows direct, consistent, and accurate evaluation of the averaged overall frequency-dependent dynamic material constitutive relations without the need for a point-wise solution of the field equations. When the spatial variation of the field variables is restricted by Bloch-form (Floquet-form) periodicity, then these relations together with the overall conservation and kinematical equations accurately yield the displacement or stress mode-shapes and, necessarily, the dispersion relations. The method can also give the point-wise solution of the elastodynamic field equations (to any desired degree of accuracy), which, however, is not required for the calculation of the average overall properties. The resulting overall dynamic constitutive relations are general and need not be restricted by the Bloch-form periodicity.The formulation is based on micromechanical modeling of a representative unit cell of the composite. For waves in periodic layered composites, the overall effective mass-density and compliance (stiffness) are always real-valued whether or not the corresponding unit cell (representative volume element used as a unit cell) is geometrically and/or materially symmetric. The average strain and linear momentum are coupled and the coupling constitutive parameters are always each others' complex conjugates. We separate the overall constitutive relations, which depend only on the composition and structure of the unit cell, from the overall field equations which hold for any elastic composite; i.e., we use only the local field equations and material properties to deduce the overall constitutive relations. Finally, we present solved numerical examples to further clarify the structure of the averaged constitutive relations and to bring out the correspondence of the current method with recently published results.     

Domain switching effect on fracture of piezoelectric solids

Rational design of smart sensors and actuators that consist of piezoelectric solids requires a thorough understanding of the constitutive behavior of this material under mechanical and electrical loading. Domain switching is the cause of significant nonlinearity in the constitutive behavior of piezoelectric solids, which may be enhanced in the presence of cracks. In this paper, the response of piezoelectric solids is formulated by coupling thermal, electrical, and mechanical effects. The corresponding finite element equations are derived and applied in the solution of the piezoelectric center crack problems. The effects of domain switching are evaluated on the near tip stress intensity factors.     

Geometrically nonlinear continuum thermomechanics with surface energies coupled to diffusion

Surfaces can have a significant influence on the overall response of a continuum body but are often neglected or accounted for in an ad hoc manner. This work is concerned with a nonlinear continuum thermomechanics formulation which accounts for surface structures and includes the effects of diffusion and viscoelasticity. The formulation is presented within a thermodynamically consistent framework and elucidates the nature of the coupling between the various fields, and the surface and the bulk. Conservation principles are used to determine the form of the constitutive relations and the evolution equations. Restrictions on the jump in the temperature and the chemical potential between the surface and the bulk are not a priori assumptions, rather they arise from the reduced dissipation inequality on the surface and are shown to be satisfiable without imposing the standard assumptions of thermal and chemical slavery. The nature of the constitutive relations is made clear via an example wherein the form of the Helmholtz energy is explicitly given.     

On the coupling of anisotropic damage and plasticity models for ductile materials

In this contribution various aspects of an anisotropic damage model coupled to plasticity are considered. The model is formulated within the thermodynamic framework and implements a strong coupling between plasticity and damage. The constitutive equations for the damaged material are written according to the principle of strain energy equivalence between the virgin material and the damaged material. The damaged material is modeled using the constitutive laws of the effective undamaged material in which the nominal stresses are replaced by the effective stresses. The model considers different interaction mechanisms between damage and plasticity defects in such a way that two-isotropic and two-kinematic hardening evolution equations are derived, one of each for the plasticity and the other for the damage. An additive decomposition of the total strain into elastic and inelastic parts is adopted in this work. The elastic part is further decomposed into two portions, one is due to the elastic distortion of the material grains and the other is due to the crack closure and void contraction. The inelastic part is also decomposed into two portions, one is due to nucleation and propagation of dislocations and the other is due to the lack of crack closure and void contraction. Uniaxial tension tests with unloadings have been used to investigate the damage growth in high strength steel. A good agreement between the experimental results and the model is obtained.     

Electroelastic waves in dielectrics modeled as polarizable continua

On the basis of a dielectric microcontinuum model, we investigate the problem of bulk wave propagation in a dielectric crystal with hexagonal material symmetry. The present linear micropolar model allows to express electric polarization via mechanical macro and micro-strain measures so that the coupling between acoustic and polarization modes can be described in terms of intrinsic dipole and quadrupole densities. The governing differential systems for different coupled modes are equivalent to some previous results of the classical phenomenological approach to ferroelectrics but also hold for piezoelectric solids with null intrinsic polarization. Resonance couplings between polaritons and acoustic waves arise from the dispersion equations depending on suitable relations among the micropolar constitutive parameters. Exploiting the dynamical representation of polarization for the admitted modes, we obtain piezoelectric coefficients and electromechanical coupling factors as functions of the wavelength (or frequency). As an application, a numerical example is given for the hexagonal phase of zinc sulfide.     

DNA纳米管的扭转刚度模型

刚度是衡量材料弹性变形难易程度的一个定量表征参数,与DNA纳米管静动力学特性及其结构生物功能密切相关．本文致力于研究DNA纳米管的扭转刚度．首先,在六角形均匀封装条件下,考虑到单个DNA杆件弯扭组合问题的静不定特点,我们利用平衡方程、变形协调方程和弹性本构方程,合理预测了DNA纳米管扭转实验中单个DNA杆件的弯扭组合变形,由此给出了DNA纳米管扭转刚度预测的解析模型．最后的结果表明：随着DNA杆数的增加,DNA纳米管的弯曲刚度显著增加,而其扭转刚度却几乎不变,合理解释了扭转实验中发现的现象．有关结论为DNA折叠结构的设计和应用提供了参考．     

On the steady solution of linear advection problems in steady recirculating flows

Numerical modelling of non-Newtonian flows typically involves the coupling between equations of motion characterized by an elliptic behaviour, and the fluid constitutive equation, which is an advection equation linked to the fluid history. In this paper we prove that linear steady advection problems in steady recirculating flows have only one solution when the kinematics differs from a rigid motion. We also give a numerical procedure to determine this steady solution. We will describe this numerical procedure for two linear models the first will be the SFRT flow model and the second will be a simplified linear formulation of the Pom–Pom viscoelastic model.