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

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

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

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

基于遗传算法的复合材料细观结构拓扑优化设计

利用高精度通用单胞模型将复合材料的细观拓扑结构与宏观力学性能结合起来,采用遗传算法对复合材料的细观结构进行优化,发展了基于遗传算法的复合材料细观结构拓扑优化设计方法.以材料的宏观力学性能为优化目标,从随机的初始细观结构出发,对复合材料纤维体积百分比进行约束,经过迭代获得满足设计要求的代表性体积单元.在优化过程中,对遗传算法的交叉过程作了较大的改进,实现了复合材料细观拓扑结构的任意变化,提高了对可行域的搜索效率.分别以极限剪切模量和泊松比为优化目标,验证了所提出优化方法的正确性和有效性.     

反映材料微观结构效应的宏观动力特性分析

探讨了实现复合材料细观结构与构件宏观动力特性直接关联的细观元法。细观元法先进行宏观单元剖分，然后分析细观结构，将微结构力学量转化为宏观结点量计算。该方法不增加自由度，但使得微结构层次上任意铺设方式和任意微小变化均能在宏观动力特性上得到反映。     

近片层γ-TiAl基合金宏观屈服的微观表征

将近片层-γTiAl基合金视为以等轴γ颗粒为基体,PST颗粒为夹杂的两相复合材料,基于细观力学自洽理论,对合金的有效弹性模量及基体和夹杂中的应力和应变场进行了解析分析计算,并结合细观力学的宏细观关联方法,确定了近片层-γTiAl基合金的宏观屈服的微观表征.结果表明:夹杂颗粒中的应力和应变场与外载及夹杂的体积分数f和椭球长细比ρ有关,软取向PST夹杂颗粒的微变形屈服导致近片层-γTiAl基合金材料的整体宏观屈服.     

界面强度对玻璃微珠填充聚丙烯力学性能的影响

本文针对玻璃微珠填充聚丙烯这一刚性粒子填充聚合物复合体系进行了实验研究。通过偶联剂改性对比,研究了该聚合物复合材料在不同界面粘结状态下的宏观拉伸、冲击力学性能。此外,根据冲击破坏断面的电镜观测结果,发现复合体系的断裂和增韧机制随界面粘结强度不同而发生改变,界面改性使得材料抗冲击破坏能力得到增强。本文还采用在位拉伸过程中的细观观测方法,观测到材料在一维应力作用下,刚性粒子和基体界面的脱粘、开裂过程,分析了该复合体系细观结构和宏观力学性能之间的关系,发现界面改性对于材料细观结构的界面脱粘和宏观屈服现象的重要影响,为发展新型复合材料提供了实验依据。     

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

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

三维五向编织复合材料渐进损伤分析及强度预测

基于材料连续体细观结构单胞,提出了材料的三维渐进损伤分析模型,采用非线性有限元方法并结合均匀化平均思想,首次建立了三维五向编织复合材料的强度预测模型。经研究典型编织角材料在拉伸载荷作用下细观损伤的发生及演化过程,分析了材料的细观失效机理,获得了材料的宏观拉伸应力应变曲线和极限破坏强度,并详细探讨了主要工艺参数编织角对材料宏观力学性能的影响规律。     

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

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

复合材料应力分析的均匀化方法

建立了基于均匀化理论的确定复合材料结构应力场的方法．其实质是用均质的宏观结构和非均质的具有周期性分布的细观结构描述原结构；将力学量表示成关于宏观坐标和细观坐标的函数，并用细观和宏观两种尺度之比为小参数展开，用摄动技术将原问题化为一细观均匀化问题和一宏观均匀化问题．这两个问题的解确定了包含等效位移和一阶近似位移的位移场，由此获得应力场．利用该方法给出了圆柱形孔隙材料和单向纤维复合材料在单向拉伸时的应力场以及空隙材料简支梁的局部应力场，说明了该方法的有效性     

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.     

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.     

Limit analysis of composite materials based on an ellipsoid yield criterion

Employing repeating unit cell (RUC) to represent the microstructure of periodic composite materials, this paper develops a numerical technique to calculate the plastic limit loads and failure modes of composites by means of homogenization technique and limit analysis in conjunction with the displacement-based finite element method. With the aid of homogenization theory, the classical kinematic limit theorem is generalized to incorporate the microstructure of composites. Using an associated flow rule, the plastic dissipation power for an ellipsoid yield criterion is expressed in terms of the kinematically admissible velocity. Based on nonlinear mathematical programming techniques, the finite element modelling of kinematic limit analysis is then developed as a nonlinear mathematical programming problem subject to only a small number of equality constraints. The objective function corresponds to the plastic dissipation power which is to be minimized and an upper bound to the limit load of a composite is then obtained. The nonlinear formulation has a very small number of constraints and requires much less computational effort than a linear formulation. An effective, direct iterative algorithm is proposed to solve the resulting nonlinear programming problem. The effectiveness and efficiency of the proposed method have been validated by several numerical examples. The proposed method can provide theoretical foundation and serve as a powerful numerical tool for the engineering design of composite materials.     

Limit analysis and homogenization of porous materials with Mohr–Coulomb matrix. Part I: Theoretical formulation

The present two-part study aims at investigating the specific effects of Mohr–Coulomb matrix on the strength of ductile porous materials by using a kinematic limit analysis approach. While in the Part II, static and kinematic bounds are numerically derived and used for validation purpose, the present Part I focuses on the theoretical formulation of a macroscopic strength criterion for porous Mohr–Coulomb materials. To this end, we consider a hollow sphere model with a rigid perfectly plastic Mohr–Coulomb matrix, subjected to axisymmetric uniform strain rate boundary conditions. Taking advantage of an appropriate family of three-parameter trial velocity fields accounting for the specific plastic deformation mechanisms of the Mohr–Coulomb matrix, we then provide a solution of the constrained minimization problem required for the determination of the macroscopic dissipation function. The macroscopic strength criterion is then obtained by means of the Lagrangian method combined with Karush–Kuhn–Tucker conditions. After a careful analysis and discussion of the plastic admissibility condition associated to the Mohr–Coulomb criterion, the above procedure leads to a parametric closed-form expression of the macroscopic strength criterion. The latter explicitly shows a dependence on the three stress invariants. In the special case of a friction angle equal to zero, the established criterion reduced to recently available results for porous Tresca materials. Finally, both effects of matrix friction angle and porosity are briefly illustrated and, for completeness, the macroscopic plastic flow rule and the voids evolution law are fully furnished.     

A LOWER BOUND LIMIT ANALYSIS OF DUCTILE COMPOSITE MATERIALS

The plastic load-bearing capacity of ductile composites such as metal matrix composites is studied with an insight into the microstructures. The macroscopic strength of a composite is obtained by combining the homogenization theory with static limit analysis, where the temperature parameter method is used to construct the self-equilibrium stress field. An interface failure model is proposed to account for the effects of the interface on the failure of composites. The static limit analysis with the finite-element method is then formulated as a constrained nonlinear programming problem, which is solved by the Sequential Quadratic Programming （SQP） method. Finally, the macroscopic transverse strength of perforated materials, the macroscopic transverse and off-axis strength of fiber-reinforced composites are obtained through numerical calculation. The computational results are in good agreement with the experimental data.     

Micromechanical approach to the strength properties of frictional geomaterials

The present paper describes a micromechanics-based approach to the strength properties of composite materials with a Drucker–Prager matrix in the situation of non-associated plasticity. The concept of limit stress states for such materials is first extended to the context of homogenization. It is shown that the macroscopic limit stress states can theoretically be obtained from the solution to a sequence of viscoplastic problems stated on the representative elementary volume. The strategy of resolution implements a non-linear homogenization technique based on the modified secant method. This procedure is applied to the determination of the macroscopic strength properties and plastic flow rule of materials reinforced by rigid inclusions, as well as for porous media. The role of the matrix dilatancy coefficient is in particular discussed in both cases. Finally, finite element solutions are derived for a porous medium and compared to the micromechanical predictions.     

A bipotential-based limit analysis and homogenization of ductile porous materials with non-associated Drucker–Prager matrix

In Gurson's footsteps, different authors have proposed macroscopic plastic models for porous solid with pressure-sensitive dilatant matrix obeying the normality law (associated materials). The main objective of the present paper is to extend this class of models to porous materials in the context of non-associated plasticity. This is the case of Drucker–Prager matrix for which the dilatancy angle is different from the friction one, and classical limit analysis theory cannot be applied. For such materials, the second last author has proposed a relevant modeling approach based on the concept of bipotential, a function of both dual variables, the plastic strain rate and stress tensors. On this ground, after recalling the basic elements of the Drucker–Prager model, we present the corresponding variational principles and the extended limit analysis theorems. Then, we formulate a new variational approach for the homogenization of porous materials with a non-associated matrix. This is implemented by considering the hollow sphere model with a non-associated Drucker–Prager matrix. The proposed procedure delivers a closed-form expression of the macroscopic bifunctional from which the criterion and a non-associated flow rule are readily obtained for the porous material. It is shown that these general results recover several available models as particular cases. Finally, the established results are assessed and validated by comparing their predictions to those obtained from finite element computations carried out on a cell representing the considered class of materials.     

Kinematic limit analysis of frictional materials using nonlinear programming

In this paper, a nonlinear numerical technique is developed to calculate the plastic limit loads and failure modes of frictional materials by means of mathematical programming, limit analysis and the conventional displacement-based finite element method. The analysis is based on a general yield function which can take the form of the Mohr–Coulomb or Drucker–Prager criterion. By using an associated flow rule, a general nonlinear yield criterion can be directly introduced into the kinematic theorem of limit analysis without linearization. The plastic dissipation power can then be expressed in terms of kinematically admissible velocity fields and a nonlinear optimization formulation is obtained. The nonlinear formulation only has one constraint and requires considerably less computational effort than a linear programming formulation. The calculation is based entirely on kinematically admissible velocities without calculation of the stress field. The finite element formulation of kinematic limit analysis is developed and solved as a nonlinear mathematical programming problem subject to a single equality constraint. The objective function corresponds to the plastic dissipation power which is then minimized to give an upper bound to the true limit load. An effective, direct iterative algorithm for kinematic limit analysis is proposed in this paper to solve the resulting nonlinear mathematical programming problem. The effectiveness and efficiency of the proposed method have been illustrated through a number of numerical examples.