刚性微粒填充的高聚物非线性粘弹性本构关系的微力学分析

在应变具有可加性条件下,本文分析表明Eshelby线性微力学理论,包括其等效包含体方法和相应的平均化方法可推广应用于研究高刚度微粒填充高聚物的非线性粘弹性本构关系。     

智能复合材料的热力学特性

由形状记忆合金丝或颗粒增强的智能复合材料具有特殊的力学性能，本文从理论上预测了智能复合材料的热力学特性。利用Ｅｓｈｅｌｂｙ的等效夹杂原理与Ｍｏｒｉ－Ｔａｎａｋａ的平均场理论导出了本构列式和相变条件，揭示了形状记忆合金在弹性介质约束下的相变机理与过程。     

Material multipole mechanics of elastic dielectric composites

The overall mechanical and electrical behaviors of elastic dielectric composites are investigated with the aid of the concept of material multipoles. In particular, by introducing a statistical continuum material multipole theory, the effects of the electric-elastic interaction and the microstructure (size, shape, orientation,...) of inhomogeneous particles on the overall behaviors of the composites can be obtained. A basic solution for an ellipsoidal elastic inhomogeneity with electric polarization in an infinite elastic dielectric medium is first given, which shows that classical Eshelby ’s elastic solution is modified by the presence of electric-elastic interaction. The overall macroscopic constitutive relations and their overall macroscopic material parameters accounting for electroelastic interaction effect are then derived for the elastic dielectric composites. Some quantitative calculations on the problems with statistical anisotropy, the shape effect and the electric-elastic interaction are finally given for dilute composites.     

A concurrent micromechanical model for predicting nonlinear viscoelastic responses of composites reinforced with solid spherical particles

A concurrent micromechanical model for predicting nonlinear viscoelastic responses of particle reinforced polymers is developed. Particles are in the form of solid spheres having micro-scale diameters. The composite microstructures are idealized by periodically distributed cubic particles in a matrix medium. Each particle is assumed to be fully surrounded by polymeric matrix such that contact between particles can be avoided. A representative volume element (RVE) is then defined by a single particle embedded in the cubic matrix. A spatial periodicity boundary condition is imposed to the RVE. One eighth unit-cell model with four particle and polymer subcells is generated due to the three-plane symmetry of the RVE. The solid spherical particle is modeled as a linear elastic material. The polymeric matrix follows nonlinear viscoelastic behaviors of thermorheologically simple materials. The homogenized micromechanical relation is developed in terms of the average strains and stresses in the subcells and traction continuity and displacement compatibility at the subcells’ interfaces are imposed. A stress–strain correction scheme is also formulated to satisfy the linearized micromechanical and the nonlinear constitutive relations. The micromechanical model provides three-dimensional (3D) effective properties of homogeneous composite responses, while recognizing microstructural geometries and in situ material properties of the heterogeneous medium. The micromechanical formulation is designed to be compatible with general displacement based finite element (FE) analyses. Experimental data and analytical micromechanical models available in the literature are used to verify the capability of the above micromechanical model for predicting the overall composite behaviors. The proposed micromodel is also examined in terms of computational efficiency and accuracy.     

Stress intensity factor of an elliptical crack as influenced by a spherical inhomogeneity

The interaction between an elliptical crack and a spherical inhomogeneity embedded in a three-dimensional solid subject to uniaxial tension is investigated. Both the inhomogeneity and the solid are isotropic but have different elastic moduli. The Eshelby's equivalent inclusion method is applied together with the principle of superposition. An approximate solution for the stress intensity factor is obtained by an approach that expands the distance between the center of the crack and inhomogeneity in series. The local stress field can be increased or decreased depending on the relative modulus of the spherical inhomogeneity and matrix. If the inhomogeneity modulus is larger than that of the matrix, a reduction in the stress intensity factor prevails. Displayed numerically are results to exhibit the influence of inhomogeneity and its distance to the crack.     

Interaction of Anisotropic Ellipsoidal Inclusions in an Isotropic Medium under Mechanical Loads and Uniform Heating

Eshelby's equivalent-inclusion method is extended to a finite number of arbitrarily oriented anisotropic ellipsoidal inclusions in an elastic isotropic matrix under polynomial mechanical loading and heating. The interaction of two identical and two different triaxial ellipsoidal inclusions in an elastic medium is studied as numerical examples. In special cases, the results are compared with those obtained by other authors     

The effective properties of piezocomposites, part I: Single inclusion problem

The problem of a piezoelectric ellipsoidal inclusion in an infinite nonpiezoelectric matrix is very important in engineering. In this paper, it is solved via Green's function technique. The closed-form solutions of the electroelastic Eshelby's tensors for this kind of problem are obtained. The electroelastic Eshelby's tensors can be expressed by the Eshelby's tensors of the perfectly elastic inclusion problem and the perfectly dielectric inclusion problem. Since the closed-form solutions of the Eshelby's tensors of the perfectly elastic inclusion problem and the perfectly dielectric inclusion problem can be given by theory of elasticity and electrodynamics, respectively, the electroelastic Eshelby's tensors can be obtained conveniently. Using these results, the closed-form solutions of the constraint elastic fields and the constraint electric fields inside the piezoelectric ellipsoidal inclusion are also obtained. These expressions can be readily utilized in solutions of numerous problems in the micromechanics of piezoelectric solids, such as the deformation and energy analysis, damage evolution and fracture of the piezoelectric materials. The project supported by the National Natural Science Foundation of China     

一种土的非线性弹性本构模型参数的反演方法

旨在提出一种土的非线性弹性本构模型参数反演的方法。以现今普遍实行的地基载荷试验为基础,依据遗传算法的组合优化理论,采用正演计算和遗传算法优化相结合的方式,建立了土层非线性弹性本构模型参数反演的方法;并依据某黄土场地地基载荷试验数据,实施了黄土土层非线性弹性本构模型参数反演的全过程。计算结果表明,所建立的方法可以实现土层非线性弹性本构模型中相互关联的多个参数的组合优化,并在对初始值要求较低的情况下,可以获得良好的参数反演精度。从而为土的变形特性分析和土与其中及相邻结构的共同作用分析,提供了较好的土体本构模型参数的确定方法。     

Discriminating linear from nonlinear elastic damage using a nonlinear time reversal DORT method

The DORT method is a selective detection and focusing technique originally developed to detect defects and damages which induce linear changes of the elastic moduli. It is based on the time reversal (TR) where a signal collected from an array of transducers is time reversed and then back-propagated into the medium to obtain focusing on selected targets. TR is based on the principle of spatial reciprocity. Attenuation, dispersion, multiple scattering, mode conversion, etc. do not break spatial reciprocity. The presence of defects or damage, may cause materials to show nonlinear elastic wave propagation behavior that will break spacial reciprocity. Therefore the DORT method will not allow focusing on nonlinear elastic scatterers. This paper presents a new method for the detection and identification of multiple linear and nonlinear scatterers by combining nonlinear elastic wave spectroscopy, time reversal and DORT method. In the presence of nonlinear hysteretic elastic scatterers, forcing the solid with a harmonic excitation, the time reversal operator can be obtained not only at the fundamental frequency of excitation, but also at the odd harmonics. At the fundamental harmonic, either inhomogeneities and linear damages can be individually selected but only at odd harmonics nonlinear hysteretic elastic damages exist. A procedure was developed where by decomposing the operator at the odd harmonics, it was possible to focus on nonlinear scatterers and to differentiate them from the linear inhomogeneities. A complete mathematical nonlinear DORT formulation for 1 and 2D structures is presented. To model the presence of nonlinear elastic hysteretic scatterers a Preisach–Mayergoyz (PM) material constitutive model was used. Results relative to 1 and 2 dimensional structures are reported showing the capability of the method to focus and discern selectively linear and nonlinear scatterers. Furthermore, an analysis was conducted to study the influence of the number of sources and their location on the imaging process showing that using a higher numbers of sensors does not automatically bring to a minor uncoupled behaviour between the nonlinear targets.     

Finite element analysis of axisymmetric elastic body problems

Linear form functions are commonly used in a long time for a toroidal volume element swept by a triangle revolved about the symmetrical axis for general axisymmetrical stress problems. It is difficult to obtain the rigidity matrix by exact integration, and instead, the method of approximate integration is used. As the locations of element close to the symmetrical axis, the accuracy of this approximation deteriorates very rapidly. The exact integration have been suggested by various authors for the calculation of rigidity matrix. However, it is shown in this paper that these exact integrations can only be used for those axisymmetric bodies with central hole. For solid axisymmetric body, it can be proved that the calculation fails due to the divergent property of rigidity matrix integration. In this paper a new form function is suggested. In this new form function, the radial displacementu vanishes as radial coordinatesr approach to zero. The calculated rigidity matrix is convergent everywhere, including these triangular toroidal element closed to the symmetrical axis. This kind of element is useful for the calculation of axisymmetric elastic solid body problems.     

The Valanis-Landel Strain-Energy Function

A constitutive equation is derived for the Cauchy stress matrix for arbitrary deformations of an isotropic elastic solid characterized by a Valanis-Landel strain-energy function. A simple example is given of the way in which results for controllable deformations of an incompressible elastic solid, with a Valanis-Landel strain-energy function, can be obtained from the known results for the more general strain-energy function employed by Rivlin.     

The crack tip zone shielding effect for ductile particle reinforced brittle materials

The crack tip zone shielding effect for the ductile particle reinforced brittle materials is analyzed by using a micromechanics constitutive theory. The theory is developed here to determine the elastoplastic constitutive behavior of the composite. The elastoplastic particles, with isotropic or kinematical hardening, are uniformly dispersed in the brittle elastic matrix. The method proposed is based on the Mori-Tanaka's concept of average stress in the composite. The macroscopic yielding condition and the incremental stress strain relation of the composite during plastic deformation are explicity given in terms of the macroscopioc applied stress and the microstructural parameters of the composite such as the volume fraction and yield stress of ductile particles, elastic constants of the two phases, etc. Finally, the contribution of the plastic deformation in the particles near a crack tip to the toughening of the composite is evaluated. The project supported by National Natural Science Foundation of China     

橡胶弹簧非线性刚度的有限元解

本文探讨了用双重非线性有限元手段，求解橡胶弹簧的非线性刚度。橡胶的本构关系采用了广义虎克定律和基于应变能的超弹性材料的本构关系联合使用的方案，用单元内的静水压力ｐ，来修正广义虎克定律的弹性常数Ｅ和Ｇ。还提出了带刚性域的超单元和波前凝聚法两种新的办法，获得橡胶弹簧的刚度，使内存用量和ＣＰＵ时间大幅度下降。因而几何、物理双重非线性橡胶弹簧的计算，能在微机上实现。文内选用英国ＤＵＮＬＯＰ橡胶公司橡胶元件     

The averaged mechanical properties of prismatic lattice formed by ellipsoidal inclusions

Summary The objective of this paper is to evaluate the averaged elastic properties of 3-D grained composites in which identical inclusions form a prismatic network interacting with the matrix material. The inclusions are of ellipsoidal shape with transverse circular sections located at the nodes of a doubly-periodic lattice with an orthogonal elementary cell. When the arrays of inclusions are set at equal spacings in normal directions through the thickness of the matrix, the material formed is an anisotropic composite with tetragonal symmetry at planes transverse to the fiber axis. The longitudinal and transverse elastic and shear moduli as well as the longitudinal Poisson's ratios of such composites are evaluated in this paper. The averaged properties are studied in terms of the aspect ratio and volume fraction of the inclusions as well as the relative rigidity of the constituent phases. Employing the Eshelby's theory for the stress field around a single ellipsoidal inhomogeneity, which is surrounded by the effective anisotropic material, and considering the Mori-Tanaka's concept for the mutual interaction of the neighboring inclusions, we may evaluate the averaged elastic properties of grained composites with aligned ellipsoidal inclusions at finite concentrations. The results provided in a closed-form solution concern the stiffness of 3-D grained composites with parallely dispersed ellipsoidal inclusions forming a prismatic network inside the principal material. It is shown that the stiffness is affected by both the geometry of the inclusions and their concentration. The use of different composite models in the analysis shows that intense variations of stiffness occur mainly in hard composites weakened by soft ellipsoidal inclusions. These findings come in full verification with experimental or theoretical results from the literature. Received 10 February 1998; accepted for publication 27 November 1998     

含有随机夹杂非均匀体的有效弹性模量

在对含有随机夹杂的非均匀体求有效弹性模量时,一般多根据Eshelby的等效夹杂法,但由于该方法没有充分考虑非均匀体内部的微结构,所以其理论具有一定的局限性。本文认为Kunin的微结构理论与Eshelby的等效夹杂法相比更具一般性,因而本文采用了文[9]中一些合理的思想,摒弃了其中不合理的假设,并且建立了一种新的理论模型.最后,本文针对球夹杂的情况给出了非均匀体有效弹性模量依赖于夹杂体积份数的关系,并将该结果与文[10]中的结果进行了比较.     

A micromechanical model of elastoplastic and damage behavior of a cohesive geomaterial

The present study is devoted to the development and validation of a nonlinear homogenization approach of the mechanical behavior of Callovo-Oxfordian argillites. The material is modeled as an heterogeneous composite composed of an elastoplastic clay matrix and of linear elastic or elastic damage inclusions. The macroscopic constitutive law is obtained by adapting the incremental method proposed by Hill [Hill, R., 1965. Continuum micro-mechanics of elastoplastic polycrystals. J. Mech. Phys. Solids 13, 89–101]. The approach consists in formulating the macroscopic tangent operator of the material by considering the nonlinear local behavior of each phase. Due to the matrix/inclusion morphology of the microstructure of the argillite, a Mori–Tanaka scheme is considered for the localization step. The developed model is first compared to Finite Element calculations and then validated and applied for the prediction of the macroscopic stress–strain responses of argillites.     

Meso-mechanical constitutive model for ratchetting of particle-reinforced metal matrix composites

Based on the Eshelby equivalent inclusion theory and Mori-Tanaka averaging method, a meso-mechanical cyclic elasto-plastic constitutive model is proposed to predict the ratchetting of particle-reinforced metal matrix composites. In the proposed model, a Hill-typed incremental formulation is used to simulate the elasto-plastic responses of the composites during cyclic loading with assumptions of elastic particle, elasto-plastic metal matrix and perfect interfacial bond between metal matrix and particles. A new nonlinear kinematic hardening rule extended from the Ohno-Abdel-Karim model [M. Abdel-Karim, N. Ohno, Kinematic hardening model suitable for ratchetting with steady-state, Int. J. Plasticity 16 (2000) 225-240] is employed to describe the ratchetting of metal matrix which dominates the ratchetting of the composites. With further assumption of spherical particles, the proposed meso-mechanical cyclic constitutive model is verified by comparing the predicted uniaxial ratchetting of SiC_P/6061Al composites with corresponding experiments obtained at room temperature [G.Z. Kang, Uniaxial time-dependent ratchetting of SiC_P/6061Al alloy composites at room and high temperature, Comp. Sci. Tech. 66 (2006) 1418-1430]. In the meantime, the effects of different tangent operators employed in the numerical implementation of the proposed model, i.e., continuum (or elasto-plastic) tangent operator C^ep and algorithmic (or consistent) one C^alg, on the predicted ratchetting are also discussed. It is concluded that the proposed model predicts the uniaxial ratchetting of SiC_P/6061Al composites at room temperature reasonably.     

Discretized virtual internal bond model for nonlinear elasticity

The virtual internal bond (VIB) is a micro–macro constitutive model. Although this model is based on a postulated discrete microstructure, it ultimately returns to a continuum constitutive relation through a homogenization process. The homogenization process can reduce the internal degrees of freedom, but it omits the effect of the individual micro bond that may play an important role in the fracture process. The present research develops a discrete system to represent the nonlinear elasticity by discretizing the continuous VIB. This discrete system is composed of unit cells, which can adopt any geometry with any number of bonds. The system is characterized by the force–displacement, not the stress–strain constitutive relationship. The nonlinear properties of this discrete system are governed by the micro-bond potential. The micro bond properties are related to Young’s modulus of the material, the volume and the bond number of the unit cell. For a given material, the unit cell has a certain topological structure and configuration. A discussion of two specific cases (the 2D triangular and 3D tetrahedral unit cells) suggests that the discrete system converges with decreasing unit cell size. In the unstructured unit cell scheme, the discrete system can almost precisely represent the initial Young’s modulus and the Poisson ratio of a nonlinear continuum. A mixed fracture example demonstrates that the present method can efficiently simulate the fracture propagation. The present paper provides a theory for developing a lattice-type mechanical model for nonlinear elasticity and provides new method for the fracture simulation of a nonlinear elastic material.     

基于三维接触有限元分析的管片等效抗弯刚度和错台研究

针对盾构隧道管片接头等效抗弯刚度预测研究中，对梁-弹簧与三维模型整体等效性考虑的不充分，以三维接触有限元计算结果作为测量信息，借助基于挠度等效的反问题求解，提出了一种确定管片接头等效抗弯刚度的新方法；并利用不同轴力-偏心距组合下的反演结果，建立了基于Kriging代理模型的轴力-弯矩-等效抗弯刚度的非线性关系，提出了由此关联的管片结构非线性问题的数值求解方法。与三维有限元结果相比，所提方法可较为准确地预测管片结构的内力和变形，表现出良好的整体等效性。此外，借助三维接触有限元分析，进一步深入探讨了螺栓孔间隙与衬垫本构关系对管片纵向错台量的影响。     

A model for compression-weakening materials and the elastic fields due to contractile cells

We construct a homogeneous, nonlinear elastic constitutive law that models aspects of the mechanical behavior of inhomogeneous fibrin networks. Fibers in such networks buckle when in compression. We model this as a loss of stiffness in compression in the stress–strain relations of the homogeneous constitutive model. Problems that model a contracting biological cell in a finite matrix are solved. It is found that matrix displacements and stresses induced by cell contraction decay slower (with distance from the cell) in a compression weakening material than linear elasticity would predict. This points toward a mechanism for long-range cell mechanosensing. In contrast, an expanding cell would induce displacements that decay faster than in a linear elastic matrix.