均匀化方法在粘弹性多层复合材料中的应用

主要研究了由各向同性线弹性加强体和各向同性线粘弹性基体组成的多层复合材料的问题,在已有的线弹性多层材料的均匀化方法的基础上,应用弹性一粘弹性对应原理,在Carson域中求解粘弹性多层材料的问题。通过Burgers模型表示线粘弹性基体材料,反演得到了多层材料的有效松弛模量和有效泊松比在时间域中的表达式,并且与实验结果和其他结果进行了比较。     

M积分与夹杂/缺陷弹性模量的显式关系简

M积分在材料构型力学中表征着缺陷自相似扩展的能量释放率,而有效弹性模量下降量在传统损伤力学中是一个具有内变量属性的损伤参数. 探讨了两者之间的特定关系,以此为材料构型力学与损伤力学搭建桥梁.借助穆斯海里什维利（Muskhelishvili）复势函数方法获取无限大弹性平面含圆形夹杂的弹性场解,根据M 积分的复势函数解析表达式得到M 积分与夹杂弹性模量的显式表达式. 随后通过有限元分析,对含复杂缺陷群的弹塑性材料进行数值模拟,结果表明内部缺陷区域的有效弹性模量下降与M 积分存在着特定关系. 基于此,提出利用材料构型力学中的外变量参数（M 积分）来替代损伤力学中的内变量参数（弹性模量下降量）描述材料的缺陷演化.     

M积分与夹杂/缺陷弹性模量的显式关系

M积分在材料构型力学中表征着缺陷自相似扩展的能量释放率,而有效弹性模量下降量在传统损伤力学中是一个具有内变量属性的损伤参数. 探讨了两者之间的特定关系,以此为材料构型力学与损伤力学搭建桥梁.借助穆斯海里什维利（Muskhelishvili）复势函数方法获取无限大弹性平面含圆形夹杂的弹性场解,根据M 积分的复势函数解析表达式得到M 积分与夹杂弹性模量的显式表达式. 随后通过有限元分析,对含复杂缺陷群的弹塑性材料进行数值模拟,结果表明内部缺陷区域的有效弹性模量下降与M 积分存在着特定关系. 基于此,提出利用材料构型力学中的外变量参数（M 积分）来替代损伤力学中的内变量参数（弹性模量下降量）描述材料的缺陷演化.      

含微裂纹弹性体的应力应变关系

本义建立了考虑裂纹闭合和裂纹表面摩擦影响的含微裂纹弹性体的应力应变关系，给出了柔度张量增量的显式表达式。对于二维平面应力和平面应变状态，给出了等效工程弹性系数。数值计算结果表明，裂纹闭合和裂纹面摩擦对裂纹体的应力应变关系和等效工程弹性系数有重要影响。     

含微裂纹弹性体的应力应变关系

本义建立了考虑裂纹闭合和裂纹表面摩擦影响的含微裂纹弹性体的应力应变关系，给出了柔度张量增量的显式表达式。对于二维平面应力和平面应变状态，给出了等效工程弹性系数。数值计算结果表明，裂纹闭合和裂纹面摩擦对裂纹体的应力应变关系和等效工程弹性系数有重要影响。     

含正交排列夹杂和缺陷材料的等效弹性模量和损伤

研究含正交排列夹杂和缺陷材料的等效弹性模量和损伤,推导了以Eshelby－Mori－Tanaka方法求解多相各向异性复合材料等效弹性模量的简便计算公式,针对含三相正交椭球状夹杂的正交各向异性材料,得到了由细观参量（夹杂的形状、方位和体积分数）表示的等效弹性模量的解析表达式．在此基础上,提出了一个宏细观结合的正交各向异性损伤模型,从而建立了以细观量为参量的含损伤材料的应力应变关系．最后,对影响材料损伤的细观结构参数进行了分析．     

Size-dependent effective elastic moduli of particulate composites with interfacial displacement and traction discontinuities

This work aims at estimating the size-dependent effective elastic moduli of particulate composites in which both the interfacial displacement and traction discontinuities occur. To this end, the interfacial discontinuity relations derived from the replacement of a thin uniform interphase layer between two dissimilar materials by an imperfect interface are reformulated so as to considerably simplify the characteristic expressions of a general elastic imperfect model which is adopted in the present work and include the widely used Gurtin–Murdoch and spring-layer interface models as particular cases. The elastic fields in an infinite body made of a matrix containing an imperfectly bonded spherical particle and subjected to arbitrary remote uniform strain boundary conditions are then provided in an exact, coordinate-free and compact way. With the aid of these results, the elastic properties of a perfectly bonded spherical particle energetically equivalent to an imperfectly bonded one in an infinite matrix are determined. The estimates for the effective bulk and shear moduli of isotropic particulate composites are finally obtained by using the generalized self-consistent scheme and discussed through numerical examples.     

Stresses in rotating heterogeneous viscoelastic composite cylinders with variable thickness

An analytical solution is presented for the rotation problem of a two-layer composite elastic cylinder under a plane strain assumption. The external cylinder has variable-thickness formulation, and is made of a heterogeneous orthotropic material. It contains a fiber-reinforced viscoelastic homogeneous isotropic solid core of uniform thickness. The thickness and elastic properties of the external cylinder are taken as power functions of the radial direction. By the boundary and continuity conditions, the radial displacement and stresses for the rotating composite cylinder are determined. The effective moduli and Illyushin’s approximation methods are used to obtain the viscoelastic solution to the problem. The effects of heterogeneity, thickness variation, constitutive, time parameters on the radial displacement, and stresses are investigated.     

瞬态波在层合板中传播的粘弹比拟法

本文发展了粘弹比拟理论,并将之用于求解半无限空间三层复合材料在垂直层合方向传播的瞬态波问题。对于层合板中应力波的传播问题,寻找到了一等效粘弹体,并用一种较好的Laplace变换的数值反演法求得了等效松弛函数和其它一些必要的辅助函数。用特征线法求得了等效粘弹体的应力和速度,进而得到了三层复合材料中心的应力、速度,进一步就得到了层中任意点的应力和速度。对于一个可由精确理论(射线理论)给出计算结果的问题,将粘弹比拟理论的结果和射线理论的结果进行了比较,结果表明,粘弹比拟理论对三层复合材料的瞬态波传播问题是相当成功的。     

Thermomagnetic viscoelastic responses in a functionally graded hollow structure

This paper presents an analytical solution for the interaction of electric potentials,electric displacements,elastic deformations,and thermoelasticity,and describes electromagnetoelastic responses and perturbation of the magnetic field vector in hollow structures(cylinder or sphere),subjected to mechanical load and electric potential.The material properties,thermal expansion coefficient and magnetic permeability of the structure are assumed to be graded in the radial direction by a power law distribution.In the present model we consider the solution for the case of a hollow structure made of viscoelastic isotropic material,reinforced by elastic isotropic fibers,this material is considered as structurally anisotropic material.The exact solutions for stresses and perturbations of the magnetic field vector in FGM hollow structures are determined using the infinitesimal theory of magnetothermoelasticity,and then the hollow structure model with viscoelastic material is solved using the correspondence principle and Illyushin’s approximation method.Finally,numerical results are carried out and discussed.     

Comparison between the relaxation spectra obtained from homogenization models and finite elements simulation for the same composite

The objective of this work is to determine the relaxation spectrum of spherical particles reinforced viscoelastic and isotropic composites from 3D Finite Elements (FE) simulations of the microstructure. The matrix and the reinforcements are assumed to be incompressible and Maxwellian. The spectra obtained from the FE simulations are compared with those obtained from analytical homogenization models. This paper presents the procedure used for generating the FE models as well as the procedure used for obtaining relaxation spectra meeting the thermodynamics requirements imposed on linear viscoelastic materials. It seems that the relaxation spectrum for the microstructure studied in this paper is composed of a negligible continuous part and a discrete part of higher intensity. In any case, the resulting material does not have a Maxwellian behavior.     

变温粘弹性杆的轴向位移的解析解

关于粘弹性材料的解析解不多,而关于变温粘弹性材料的解析解则更少．本文在三维变温粘弹性本构方程［1］的基础上,将构件材料的终态温度等效松弛曲线表示成负幂指数的形式,推导出了粘弹性材料的杆在轴力作用下处于变温状态下的位移的解析解．     

The elastic–viscoelastic correspondence principle for non-homogeneous materials with time translation non-invariant properties

This paper revisits the elastic–viscoelastic correspondence principle for non-homogeneous materials. Several recent publications discussed this principle for functionally graded materials (FGMs) with time translation invariant viscoelastic properties. It was demonstrated that the correspondence principle is valid only for the FGMs with separable relaxation moduli (moduli in separable form in space and time). This paper reconsiders this issue. It shows that the correspondence principle is valid even for non-homogeneous materials with separable relaxation moduli even if the time-dependences of the relaxation moduli in shear and dilatation are not necessarily time translation invariant. The property of similarity of Volterra operators is used to obtain the corresponding elastic solution. The correspondence is established between the elastic solution and the operator-transformed viscoelastic solution. The transformation operators are combinations of the Laplace transform operator and additional integral operators.     

A continuum theory for isotropic two-phase elastic composites

The paper presents the effective stiffness theory for isotropie two-phase elastic composites. The theory predicts dispersion of longitudinal and transverse plane time-harmonic travelling waves. The limiting phase velocities at vanishing wave numbers serve in the determination of the elastic moduli of the equivalent homogeneous isotropic medium. These elastic moduli are compared with the effective moduli defined statically.     

Load transfer in fibre-reinforced composites with viscoelastic matrix: an analytical study

In A fibre-reinforced 2D composite material with elastic fibres and viscoelastic, isotropic matrix is studied. Starting from the solution of a reference-problem with elastic matrix material the elastic matrix parameters are substituted by their viscoelastic correspondents in the Laplace domain. For simplification the time-dependent solution is approximated by using limiting value theorems that give information about the time-dependent solution for t → 0 and t → ∞. Then the method of asymptotically equivalent functions is used and illustrated with examples of a steel fibre in a PMMA matrix. The analytical solutions are compared with their numerical counterparts. In summary it can be stated that this paper is a further contribution to the vast literature about the application of the correspondence principle to the solution of special problems of the linear viscoelasticity.     

Determination and analysis of the effective relaxation properties of a composite with viscoelastic components

The relaxation properties of a two-component material are determined depending on time, volume fraction, and type of reinforcement, and the relationship among them. The type of reinforcement is determined by the aspect ratio of the ellipsoid of revolution that models the inclusion. The effective moduli of the composite are determined from the relaxation properties of the components. It is assumed that the composite components are made of isotropic viscoelastic materials with volume expansion and shear characteristics described by two Rabotnov’s fractional-exponential functions with different orders of fractionality. To obtain the solution in the time domain, its fractional rational representation in the frequency domain is used. Optimizing the parameters of this representation and transforming the parameters of the solution to the time domain make it possible to obtain solutions in compact form in terms of relaxation kernels     

A Note on Perturbation Formulae for the Surface-Wave Speed Due To Perturbations in Material Properties

In a recent paper by Tanuma and Man, a two-term asymptotic formula was derived for the speed of surface waves propagating in an anisotropic elastic half-space whose elastic moduli differ only slightly from those for a (base) isotropic elastic material. This formula disagrees with that derived by Delsanto and Clark in an earlier paper using a different method. In this short note, we use a simple procedure to derive another two-term asymptotic formula for the surface-wave speed. Our formula takes the same compact form even if the base material is generally anisotropic. We show that when an error in the work of Delsanto and Clark is corrected, the three different methods do give equivalent results.      

Determination of material response functions for prestrained rubbers

The mechanical response of two natural rubber compounds is examined in order to determine relevant material parameters by non-linear finite element analysis. The materials are subjected to (a) combined static torsion and extension, and (b) small, steady-state torsional oscillations superposed on a large static simple extension. The materials are assumed to be incompressible and isotropic in their undeformed state and a time-strain separable relaxation modulus tensor is employed in order to characterize the steady-state harmonic viscoelastic response. The combined static torsion and extension experiments are used to determine the basic delayed elastic response functions. A Rivlin-type strain energy expression of third-order accuracy is used for the purpose. The two-constant, Mooney-Rivlin form is found to be adequate for both materials in the relatively limited range of deformation magnitudes considered.The torsional storage and loss moduli are determined under quasistatic conditions as functions of frequency and axial static pre-strain. The time-strain separability is found to be a resonable approximation in a relatively limited range of static prestrain magnitudes and frequencies considered for the natural gum rubbers investigated. The experimental methodology is discussed in some detail.     

A continuum micromechanical theory of overall plasticity for particulate composites including particle size effect

Classical continuum micromechanics cannot predict the particle size dependence of the overall plasticity for composite materials, a simple analytical micromechanical method is proposed in this paper to investigate this size dependence. The matrix material is idealized as a micropolar continuum, an average equivalent inclusion method is advanced and the Mori–Tanaka's method is extended to a micropolar medium to evaluate the effective elastic modulus tensor. The overall plasticity of composites is predicted by a new secant moduli method based on the second order moment of strain and torsion of the matrix in a framework of micropolar theory. The computed results show that the size dependence is more pronounced when the particle's size approaches to the matrix characteristic length, and for large particle sizes, the prediction coincides with that predicted by classical micromechanical models. The method is analytical in nature, and it can capture the particle size dependence on the overall plastic behavior for particulate composites, and the prediction agrees well with the experimental results presented in literature. The proposed model can be considered as a natural extension of the widely used secant moduli method from a heterogeneous Cauchy medium to a micropolar composite.     

Bounds and perturbation series for incompressible elastic composites with transverse isotropic symmetry

The effective elastic behavior of a transversely isotropic composite made from two incompressible elastic materials is examined. The set of all effective elasticity tensors for transversely isotropic finite rank laminar microstructures is described. The extremal property of this class of microstructures is used to derive a new more precise characterization of the set of effective shear moduli.The perturbation series for the effective elasticity tensor is considered. An explicit formula for the second order perturbation tensor is derived. We describe precisely the set of tensors that correspond to all second order perturbations consistent with transverse isotropy. We apply analytic methods [cf. 27] to show that all second order perturbation tensors are realized by finite rank laminar microstructures.Supported by NSF through Grant DMS-3907658.