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

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

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.     

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

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

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

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

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.     

Assessment of existing and introduction of a new and robust efficient definition of the representative volume element

Accurate numerical homogenization necessitates the thorough determination of the Representative Volume Element (RVE). There exists several seminal works on the notion of the RVE in homogenization, its definitions and methods of determination for efficient computation of composite effective properties. The objective of the current work is to assess the ability of numerical RVE determination methods to deliver accurate effective properties of composite materials. This paper demonstrates that common and well-established RVE determination methods, based on studying the convergence rate of the effective properties with respect to the volume element size, are invalid for the case of composites reinforced by randomly oriented fibers and yield erroneous estimates of their effective properties. Following the failure of traditional RVE determination methods, we proposed a new RVE determination criterion that is not based on the average property stability, but its statistical variations. Our new proposed criterion has been shown to be more accurate than other criteria in computing the effective properties of composites for aspect ratios up to 60. Moreover, the proposed criterion does not necessitate a convergence study over the volume element size, hence reducing considerably the RVE determination cost. Finally, our work questions the validity of many published works dealing with composites including heterogeneities of high aspect ratios.     

A computational homogenization approach for the yield design of periodic thin plates. Part II : Upper bound yield design calculation of the homogenized structure

In the first part of this work (Bleyer and de Buhan, 2014), the determination of the macroscopic strength criterion of periodic thin plates has been addressed by means of the yield design homogenization theory and its associated numerical procedures. The present paper aims at using such numerically computed homogenized strength criteria in order to evaluate limit load estimates of global plate structures. The yield line method being a common kinematic approach for the yield design of plates, which enables to obtain upper bound estimates quite efficiently, it is first shown that its extension to the case of complex strength criteria as those calculated from the homogenization method, necessitates the computation of a function depending on one single parameter. A simple analytical example on a reinforced rectangular plate illustrates the simplicity of the method. The case of numerical yield line method being also rapidly mentioned, a more refined finite element-based upper bound approach is also proposed, taking dissipation through curvature as well as angular jumps into account. In this case, an approximation procedure is proposed to treat the curvature term, based upon an algorithm approximating the original macroscopic strength criterion by a convex hull of ellipsoids. Numerical examples are presented to assess the efficiency of the different methods.     

A mathematical programming algorithm for limit analysis

This paper deals with the limit analyses of perfect rigid-plastic continua. Based on the kinematic theorem of the limit analysis theory, a mathematical programming finite element formula for determining the upper bound load multiplier has been established, and an iteration algorithm proposed accordingly. In this algorithm the plastic and rigid zones are distinguished for every iteration step, and the goal function is modified gradually. The difficulties caused by the nonsmoothness of the goal function are overcome. Some examples solved by this algorithm are presented. The project supported by National Natural Science Foundation of China.     

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.     

A nonlinear programming approach to kinematic shakedown analysis of frictional materials

This paper develops a novel nonlinear numerical method to perform shakedown analysis of structures subjected to variable loads by means of nonlinear programming techniques and the displacement-based finite element method. The analysis is based on a general yield function which can take the form of most soil yield criteria (e.g. the Mohr–Coulomb or Drucker–Prager criterion). Using an associated flow rule, a general yield criterion can be directly introduced into the kinematic theorem of shakedown analysis without linearization. The plastic dissipation power can then be expressed in terms of the kinematically admissible velocity and a nonlinear formulation is obtained. By means of nonlinear mathematical programming techniques and the finite element method, a numerical model for kinematic shakedown analysis is 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 shakedown load can be calculated. An effective, direct iterative algorithm is then proposed to solve the resulting nonlinear programming problem. The calculation is based on the kinematically admissible velocity with one-step calculation of the elastic stress field. Only a small number of equality constraints are introduced and the computational effort is very modest. The effectiveness and efficiency of the proposed numerical method have been validated by several numerical examples.     

HOMOGENIZATION—BASED TOPOLOGY DESIGN FOR PURE TORSION OF COMPOSITE SHAFTS

In conjunction with the homogenization theory and the finite element method, the mathematical models for designing the corss-section of composite shafts by maximizing the torsion rigidity are developed in this paper. To obtain the extremal torsion rigidity, both the cross-section of the macro scale shaft and the representative microstructure of the composite material are optimized using the new models. The micro scale computational model addresses the problem of finding the periodic microstructures with extreme shear moduli. The optimal microstructure obtained with the new model and the homogenization method can be used to improve and optimize natural or artificial materials. In order to be more practical for engineering applications, cellular materials rather than ranked materials are used in the optimal process in the existence of optimal bounds for the elastic properties. Moreover, the macro scale model is proposed to optimize the cross-section of the torsional shaft based on the tailared composites. The validating optimal results show that the models are very effective in obtaining composites with extreme elastic properties, and the cross-section of the composite shaft with the extremal torsion rigidity. The project supported by the National Natural Science Foundation of China (10172078 and 10102018)     

A kinematic shakedown theorem considering external loading and temperature variation

A kinematic shakedown theorem including the external loading and temperature variation is derived based on the concept of an assumed yield surface. An example is given to show the efficiency of the theorem obtained.The Project Supported by National Natural Science Foundation of China.     

Constitutive response of idealized granular media under the principal stress axes rotation

This paper is dedicated to the understanding of the phenomena, which give rise to anisotropy and non-coaxiality in granular materials. In achieving three-dimensional numerical simulation under static condition of granular media, granular element method (GEM) is adopted in this study. The method has been incorporated into the so-called mathematical homogenization theory for quasi-static equilibrium problems, which enables us to obtain the macroscopic/phenomenological inelastic deformation response of a representative volume element (RVE). To examine the anisotropic macroscopic deformation properties of the assumed RVE, which is solved by granular element method (GEM), a series of numerical experiments involving the pure rotation of the principal stress axes are carried out, and its results are discussed in relation to induced anisotropy and non-coaxiality.     

A Micromechanics-based Finite Element Model for the Constitutive Behavior of Polycrystalline Ferromagnets

A micromechanics-based finite element model for the constitutive behavior of polycrystalline ferromagnets is developed. In the model, the polycrystalline solid is assumed to comprise numerous single crystals with randomly distributed crystallographic orientations, and the single crystals, in turn, consist of ferromagnetic domains, each of which is represented by a cubic element. The dipole directions of the domains are randomly assigned to simulate the crystallographic nature of ferromagnetic polycrystals. A switching criterion for the domains is specified at the microscopic level. The macroscopic constitutive behavior is obtained by averaging the microscopic/local behavior of each domain. The developed model has been applied to the simulation of a ferromagnetic material. With appropriate material parameters adopted, hysteresis loops of the predicted magnetic induction versus magnetic field and those of the strain versus magnetic field are shown to agree well with experimental observations.The project supported by the National Natural Science Foundation of China (90205030, 10472088, 10425210), the National Basic Research Program of China (2006CB601202) and the State Administration of the Foreign Experts Affairs Through the “111” Project (B06024) The English text was polished by Yunming Chen.     

A plasticity model for cellular materials with open-celled structure

Based on a rigid-plastic material model that obeys the von Mises yield criterion, the plastic behavior of foams with an open-celled structure is studied in this paper using a single unit cell. An approximate continuum plasticity model is developed within the framework of the upper bound theorem of plasticity to describe the yield behavior of foams. The microscopic velocity fields are derived for the unit cell, which satisfy the incompressibility and the kinematic boundary conditions, and expressed in macroscopic rate of deformation. From the microscopic velocity fields, a macroscopic yield function is developed for foams under multi-axial stresses and includes the effects of the hydrostatic stress due to the void presence and growth. The dependency of the derived yield surfaces of foams on their relative densities is studied. The plastic behavior of foams is also studied numerically using the finite element method. The newly developed plasticity model is compared with the finite element analysis results and other available foam models and then correlated with the finite element results.     

A computational homogenization approach for the yield design of periodic thin plates. Part I: Construction of the macroscopic strength criterion

The purpose of this paper is to propose numerical methods to determine the macroscopic bending strength criterion of periodically heterogeneous thin plates in the framework of yield design (or limit analysis) theory. The macroscopic strength criterion of the heterogeneous plate is obtained by solving an auxiliary yield design problem formulated on the unit cell, that is the elementary domain reproducing the plate strength properties by periodicity. In the present work, it is assumed that the plate thickness is small compared to the unit cell characteristic length, so that the unit cell can still be considered as a thin plate itself. Yield design static and kinematic approaches for solving the auxiliary problem are, therefore, formulated with a Love–Kirchhoff plate model. Finite elements consistent with this model are proposed to solve both approaches and it is shown that the corresponding optimization problems belong to the class of second-order cone programming (SOCP), for which very efficient solvers are available. Macroscopic strength criteria are computed for different type of heterogeneous plates (reinforced, perforated plates,…) by comparing the results of the static and the kinematic approaches. Information on the unit cell failure modes can also be obtained by representing the optimal failure mechanisms. In a companion paper, the so-obtained homogenized strength criteria will be used to compute ultimate loads of global plate structures.     

Thermomechanical Multiscale Constitutive Modeling: Accounting for Microstructural Thermal Effects

In this work we present a thermomechanical multiscale constitutive model for materials with microstructure. In these materials thermal effects at microscale have an impact on the effective macroscopic stress. As a result, it turns out that the homogenized stress depends upon the macroscopic temperature and its gradient. In order to allow this interplay to be thermodynamically valid, we resort to a macroscopic extended thermodynamics whose elements are derived from the microscopic behavior using homogenization concepts. Hence, the thermodynamics implications of this new class of multiscale models are discussed. A variational approach based on the Hill–Mandel Principle of Macro-homogeneity, and which makes use of the volume averaging concept over a local representative volume element (RVE), is employed to derive the thermal and mechanical equilibrium problems at the RVE level and the corresponding homogenization expressions for the effective heat flux and stress. The material behavior at the RVE level is described through standard phenomenological constitutive models. To sum up, the novel contribution of the model presented here is that it allows to include the microscopic temperature fluctuation field, obtained from the multiscale thermal analysis, in the micro-mechanical problem at the RVE level while keeping thermodynamic consistency.     

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

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