复合材料中的渐近均匀化方法

本文将非均质弹性体的渐近均匀化方法应用于复合材料的宏观与细观分析之中。该方法基于平均化的思想，将复合材料视作由周期性的细观结构所构成，其场变量依赖于宏观和细观两个尺度的坐标变量而变化。通过建立位移和应力的渐近表达式，推导出关于周期性基元的细观平衡方程和细观本构关系，并与有限元数值方法相结合，得到材料的宏观等效性能和细观应力分布。对典型算例的分析，反映出该方法的有效性及准确性。     

基于迭代法的非线性弹性均质化研究

微观结构对复合材料的宏观力学性能具有至关重要的影响,通过合理设计复合材料微观结构可以得到期望的宏观性能.均质化方法作为一种有效的设计方法,它从微观结构的角度出发,利用均匀化的概念,实现了对复合材料宏观力学性能的预测和设计.而当考虑非线性因素,均质化的实现就非常困难.本文利用双渐近展开方法,将位移按照宏观位移和微观位移展开,推导了非线性弹性均质化方程.通过直接迭代法,对非线性弹性均质化方程进行了求解,并给出了具体的迭代方法和实现步骤.本文基于迭代步骤和非线性弹性均质化方程编写MATLAB程序,对3种典型本构关系的周期性多孔材料平面问题进行了计算,对比细致模型的应变能、最大位移和等效泊松比,对程序及迭代方法的准确性进行了验证.之后对一种三元橡胶基复合材料进行多尺度均质化,将其分为芯丝尺度和层间尺度.用线弹性的均质化方法得到了芯丝尺度的等效弹性参数,并将其作为层间尺度的材料参数.在层间尺度应用非线性弹性均质化方法对结构进行计算,得到材料的宏观等效性能,并以实验结果为基准进行评价.     

基于迭代法的非线性弹性均质化研究

微观结构对复合材料的宏观力学性能具有至关重要的影响, 通过合理设计复合材料微观结构可以得到期望的宏观性能. 均质化方法作为一种有效的设计方法, 它从微观结构的角度出发, 利用均匀化的概念, 实现了对复合材料宏观力学性能的预测和设计. 而当考虑非线性因素, 均质化的实现就非常困难. 本文利用双渐近展开方法, 将位移按照宏观位移和微观位移展开, 推导了非线性弹性均质化方程. 通过直接迭代法, 对非线性弹性均质化方程进行了求解, 并给出了具体的迭代方法和实现步骤. 本文基于迭代步骤和非线性弹性均质化方程编写MATLAB 程序, 对3种典型本构关系的周期性多孔材料平面问题进行了计算, 对比细致模型的应变能、最大位移和等效泊松比, 对程序及迭代方法的准确性进行了验证. 之后对一种三元橡胶基复合材料进行多尺度均质化, 将其分为芯丝尺度和层间尺度. 用线弹性的均质化方法得到了芯丝尺度的等效弹性参数, 并将其作为层间尺度的材料参数. 在层间尺度应用非线性弹性均质化方法对结构进行计算, 得到材料的宏观等效性能, 并以实验结果为基准进行评价.      

确定复合材料宏观屈服准则的细观力学方法

运用细观力学中的均匀化方法，分析了含周期性微结构复合材料的宏观屈服准则，并对Hill-Tsai准则进行了修正。从基于复合材料细观结构的代表性胞元入手，运用塑性极限理论中的机动分析以及有限元方法，计算了细观结构的极限载荷域。通过宏细观尺度对应关系，得到复合材料的宏观屈服准则。     

复合材料刚度性能表征的协同损伤力学模型

针对复合材料层合板的弥散型损伤,提出一个刚度性能表征的协同损伤力学模型. 该模型兼顾了微观物理损伤响应和宏观材料刚度性能表征. 从微观角度,建立细观RVE 模型求解裂纹表面张开位移和滑开位移,以此定义损伤张量,并在宏观上通过对材料应变和损伤表面位移进行均匀化处理,建立单向板或层合板的损伤刚度矩阵和损伤张量之间的联系. 以基体裂纹为例,详细分析并建立了横向裂纹和纵向裂纹的损伤本构. 计算了[±θ/90₄]_S 铺层层合板中基体横向裂纹对刚度性能的影响,结果表明该方法能够准确地预测复合材料层合板由损伤导致的刚度性能衰减.      

混凝土材料宏观力学特性分析的细观单元等效化模型

提出了一种混凝土材料宏观力学特性分析的新方法—细观单元等效化模型。该方法从描述混凝土材料的细观尺度入手,采用Monte Carlo法生成由骨料颗粒及砂浆基质组成的混凝土试件的随机骨料模型;然后,依据混凝土材料特征单元尺度来剖分有限元网格并投影到建立的随机骨料模型上,各细观单元的有效力学特性则采用复合材料等效化方法来确定。本文方法体现了材料非线性宏观力学特性源于其内在的不均匀性这一认识,而对不均匀性的描述则是以网格剖分是否影响其宏观力学特性为准则。因此,本文方法较其他细观力学方法最大的优点在于极大地减小了体系自由度数目（特别是对于三维问题）,从而提高了计算效率。算例分析初步验证了本文方法的高效性。     

基于细观力学方法的混凝土热膨胀系数预测

建立混凝土材料的有效性质与微结构参数之间的关系,是混凝土材料优化设计的基础。本文用细观力学方法对复合材料宏观有效热膨胀系数进行研究,得到了含有一球形夹杂物的无限大介质在均匀变温作用下的应力场。假定混凝土为由骨料和砂浆基质组成的二相复合材料,根据混凝土宏观体积热膨胀量与组成混凝土的各相介质细观体积热膨胀量相等的原则,采用基于Mori-Tanaka方法的混凝土宏观有效剪切模量,推导出混凝土有效热膨胀系数的解答。对稀疏解法、自洽方法和有限单元数值试验结果的比较说明,本文提出的基于自洽方法的混凝土宏观有效热膨胀系数的理论公式能够较好的描述混凝土的热学特性,该方法可以推广到多相复合材料宏观有效热膨胀系数的预测中。     

多孔材料塑性极限载荷及其破坏模式分析

运用塑性力学中的机动极限分析理论，研究韧性基体多孔材料的塑性极限承载能力和破坏模式。以多孔材料的细观结构为研究对象，将细观力学中的均匀化理论引入到塑性极限分析中，并结合有限元技术，建立细观结构极限载荷的一般计算格式，并提出相应的求解算法。数值算例表明：细观孔洞对材料的宏观强度影响明显；在单向拉伸作用下，孔洞呈现膨胀扩大规律；多孔材料破坏源于基体塑性区的贯通。     

复合材料的宏观性能与参数设计

本文综述了预测复合材料宏观性能──有效刚度的几类方法：自洽模型、单胞模型以及它们的结合──自洽有限元法．阐述了复合材料发生弹塑性变形时的有关力学问题．基于细观力学的定量分析结果，探讨了面向材料宏观刚度的细观结构参数设计的基本原则，以期对建立复合材料细观结构设计的力学和数学模型有所启发．     

基于均匀化理论韧性复合材料塑性极限分析

运用细观力学中的均匀化方法分析了韧性复合材料的塑性极限承载能力.从反映复合材料细观结构的代表性胞元入手,将均匀化理论运用到塑性极限分析中,计算由理想刚塑性、Mises组分材料构成的复合材料的极限承载能力.运用机动极限方法和有限元技术,最终将上述问题归结为求解一组带等式约束的非线性数学规划问题,并采用一种无搜索直接迭代算法求解.为复合材料的强度分析提供了一个有效手段.     

Three-dimensional stochastic analysis using a perturbation-based homogenization method for elastic properties of composite material considering microscopic uncertainty

This paper discusses evaluation of influence of microscopic uncertainty on a homogenized macroscopic elastic property of an inhomogeneous material. In order to analyze the influence, the perturbation-based homogenization method is used. A higher order perturbation-based analysis method for investigating stochastic characteristics of a homogenized elastic tensor and an equivalent elastic property of a composite material is formulated.As a numerical example, macroscopic stochastic characteristics such as an expected value or variance, which is caused by microscopic uncertainty in material properties, of a homogenized elastic tensor and homogenized equivalent elastic property of unidirectional fiber reinforced plastic are investigated. The macroscopic stochastic variation caused by microscopic uncertainty in component materials such as Young’s modulus or Poisson’s ratio variation is evaluated using the perturbation-based homogenization method. The numerical results are compared with the results of the Monte-Carlo simulation, validity, effectiveness and a limitation of the perturbation-based homogenization method is investigated. With comparing the results using the first-order perturbation-based method, effectiveness of a higher order perturbation is also investigated.     

PLASTIC LIMIT ANALYSIS OF DUCTILE COMPOSITE STRUCTURES FROM MICRO-TO MACRO-MECHANICAL ANALYSIS

The load-bearing capacity of ductile composite structures comprised of periodic composites is studied by a combined micro/macromechanicai approach. Firstly, on the microscopic level, a representative volume element （RVE） is selected to reflect the microstructures of the composite materials and the constituents are assumed to be elastic perfectly-plastic. Based on the homogenization theory and the static limit theorem, an optimization formulation to directly calculate the macroscopic strength domain of the RVE is obtained. The finite element modeling of the static limit analysis is formulated as a nonlinear mathematical programming and solved by the sequential quadratic programming method, where the temperature parameter method is used to construct the self-stress field. Secondly, Hill＇s yield criterion is adopted to connect the micromechanicai and macromechanical analyses. And the limit loads of composite structures are worked out on the macroscopic scale. Finally, some examples and comparisons are shown.     

Micromechanical analysis of heterogeneous materials subjected to overall Cosserat strains

In the framework of the computational homogenization procedures, the problem of coupling a Cosserat continuum at the macroscopic level and a Cauchy medium at the microscopic level, where a heterogeneous periodic material is considered, is addressed. In particular, non-homogeneous higher-order boundary conditions are defined on the basis of a kinematic map, properly formulated for taking into account all the Cosserat deformation components and for satisfying all the governing equations at the micro-level in the case of a homogenized elastic material. Furthermore, the distribution of the perturbation fields, arising when the actual heterogeneous nature of the material is taken into account, is investigated. Contrary to the case of the first-order homogenization where periodic fluctuations arise, in the analyzed problem more complex distributions emerge.     

Asymptotic homogenization of three-dimensional thermoelectric composites

Thermoelectric composites are promising for high efficiency energy conversion between thermal flows and electric conduction, though their effective behaviors remain poorly understood due to nonlinear thermoelectric coupling. In this paper, we develop an asymptotic homogenization theory to analyze the effective behavior of three-dimensional (3D) thermoelectric composites, built on the observation that the equations governing microscopic field fluctuations in the composite are actually linear instead of nonlinear after separation of length scales. A set of solutions similar to Green's function method are used to construct the unit cell problem, and appropriate interfacial continuity conditions and boundary conditions are derived. The homogenized governing equations are then developed for thermoelectric composites, and they are further reduced for a special case wherein the heat flow and electric conduction in the composite remains one-dimensional (1D) at macroscopic scale, even though the composite itself is 3D in general. The general homogenization theory is implemented using finite element method, and a key constant in the constructed solutions is determined using the reformulated eigenvalue problem. The algorithm is validated, and is applied for a number of case studies for the effective behavior of thermoelectric composites.     

Micro/macromechanical plastic limit analyses of composite materials and structures

The load-bearing capacities of ductile composite materials and structures are studied by means of a combined micro/macromechanics approach. Firstly, on the microscopic scale, the aim is to get the macroscopic strength domains by means of the homogenization theory of micromechanics. A representative volume element (RVE) is selected to reflect the microstructures of the composite materials. By introducing the homogenization theory into the kinematic limit theorem of plastic limit analysis, an optimization format to directly calculate the limit loads of the RVE is obtained. And the macroscopic yield criterion can be determined according to the relation between macroscopic and microscopic fields. Secondly, on the macroscopic scale, by introducing the Hill's yield criterion into the kinematic limit theorem, the limit loads of orthotropic structures such as unidirectional fiber-reinforced composite structures are worked out. The finite element modeling of the kinematic limit analysis is deduced into a nonlinear mathematical programming with equality-constraint conditions that can be solved by means of a direct iterative algorithm. Finally, some examples are illustrated to show the application of the present approach. Project supported by the National Natural Science Foundation of China (No. 19902007), the National Foundation for Excellent Doctoral Dissertation of China (No. 200025), the Fund of the Ministry of Education of China for Returned Oversea Scholars and the Basic Research Foundation of Tsinghua University.     

Micro-stress prediction in composite laminates with high stress gradients

The objective of this research is to develop a macroscopic theory, which can provide the connection between macro-mechanics and micro-mechanics in characterizing the micro-stress of composite laminates in regions of high macroscopic stress gradients. The micro-polar theory, a class of higher-order elasticity theory, of composite laminate mechanics is implemented in a well-known Pipes–Pagano free edge boundary problem. The micro-polar homogenization method to determine the micro-polar anisotropic effective elastic moduli is presented. A displacement-based finite element method based on micro-polar theory in anisotropic solids is developed in analyzing composite laminates. The effects of fiber volume fraction and cell size on the normal stress along the artificial interface resulting from ply homogenization of the composite laminate are also investigated. The stress response based on micro-polar theory is compared with those deduced from the micro-mechanics and classical elasticity theory. Special attention of the investigation focuses on the stress fields near the free edge where the high macro-stress gradient occurs. The normal stresses along the artificial interface and especially, the micro-stress along the fiber/matrix interface on the critical cell near the free edge where the high macro-stress gradient detected are the focus of this investigation. These micro-stresses are expected to dominate the failure initiation process in composite laminate. A micro-stress recovery scheme based on micro-polar analysis for the prediction of interface micro-stresses in the critical cell near the free edge is found to be in very good agreement with “exact” micro-stress solutions. It is demonstrated that the micro-polar theory is able to capture the micro-stress accurately from the homogenized solutions.     

Stochastic homogenization analysis on elastic properties of fiber reinforced composites using the equivalent inclusion method and perturbation method

This paper describes a methodology for evaluation of influence of microscopic uncertainty in material properties and geometry of a microstructure on a homogenized macroscopic elastic property of an inhomogeneous material. For the analysis of the stochastic characteristics of a homogenized elastic property, the first-order perturbation method is used. In order to analyze the influence of microscopic geometrical uncertainty, the perturbation-based equivalent inclusion method is formulated. In this paper, an analytical form of the perturbation term using the equivalent inclusion method is provided.As a numerical example, macroscopic stochastic characteristics such as an expected value or variance of the homogenized elastic tensor of a unidirectional fiber reinforced plastic, which is caused by microscopic uncertainty in material properties or geometry of a microstructure, are estimated with computing the first order perturbation term of the homogenized elastic tensor. Compared the results of the proposed method with the results of the Monte-Carlo simulation, validity, effectiveness and a limitation of the perturbation-based homogenization method is investigated.