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.     

Optimum design of laminated composite foam-filled sandwich plates subjected to strength constraint

Optimum design of laminated composite sandwich plates with both continuous (core thickness) and discrete (layer group fiber angles and thicknesses) design variables subjected to strength constraint is studied via a two-level optimization technique. The strength of a sandwich plate is determined in a failure analysis using the Tsai–Wu failure criterion and the finite element method which is formulated on the basis of the layerwise linear displacement theory. In the first level optimization of the design process, the discrete design variables are temporarily treated as continuous variables and the corresponding minimum weight of the sandwich plate is evaluated subject to the strength constraint using a constrained multi-start global optimization method. In the second level optimization, the optimal solution obtained in the first level optimization is used in the branch and bound method for solving a discrete optimization problem to determine the optimal design parameters and the final weight of the plate. Failure test of laminated composite foam-filled sandwich plates with different lamination arrangements are performed to validate the proposed optimal design method. A number of examples of the design of laminated composite foam-filled sandwich plates are given to demonstrate the feasibility and applications of the proposed method.     

Étude du critère de plasticité des matériaux poreux

The ductile failure of porous metallic materials is studied here using both Limit Analysis (LA) methods, a problem treated by Gurson with his famous kinematic approach in 1977. The present work is devoted to determining the strength of porous materials with long circular cylindrical cavities in the case of plane stress. The numerical methods developed here use the Hill–Mandel method based on the homogenization theory of heterogeneous media within the LA framework. The use of kinematic and static approaches gave an excellent estimation of the yield criterion for all the cases studied. The numerical results based on LA methods have been compared with analytical and semi-analytical yield domain expressions proposed by different authors. The results show that the Richmond model was the most accurate in terms of our predictions.     

Limit analysis of orthotropic plates

The limit analysis problem for plates in bending is considered. The failure criterion for the material is assumed as orthotropic, with possible non-symmetric strength properties. According to Kirchhoff’s hypothesis, the plate is conceived as a superposition of layers, individually in plane stress situation, and continuity is enforced by means of a kinematic assumption. By exploiting previous results, recently established by the authors, the expression of the dissipation power per unit plate area is defined on this basis and the kinematic (upper bound) theorem of limit analysis is cast in a form suitable for numerical solutions. To this purpose, efficient algorithms successfully employed in the isotropic case can be used with minor modifications. The effectiveness of the procedure is demonstrated by solving some homogeneous plate examples. Results permit the assessment of the influence of different aspects, such as the ratio between strengths along the orthotropy directions, the tensile to compressive strength differential and the inclination of the orthotropy axes with respect to the sides. The effects of in-plane edge constraints are also discussed and it appears that they are emphasized considerably by anisotropy. Even if referred to specific cases, some conclusions can be regarded as fairly general.     

A Reissner–Mindlin limit analysis model for out-of-plane loaded running bond masonry walls

Earthquake surveys have demonstrated that the lack of out-of-plane strength is a primary cause of failure in many traditional forms of masonry. Moreover, bearing walls are relatively thick and, as a matter of fact, many codes of practice impose a minimal slenderness for them, as for instance the recent Italian O.P.C.M. 3431 [2005. Ulteriori modifiche ed integrazioni all’OPCM 3274/03 (in Italian) and O.P.C.M. 3274, 20/03/2003, Primi elementi in materia di criteri generali per la classificazione sismica del territorio nazionale e di normative tecniche per le costruzioni in zona sismica (in Italian)], in which the upper bound slenderness is fixed respectively equal to 12 for artificial bricks and 10 for natural blocks masonry. In this context, a formulation at failure for regular assemblages of bricks based both on homogenization and Reissner–Mindlin theory seems particularly attractive. In this paper a kinematic limit analysis approach under the hypotheses of the thick plate theory is developed for the derivation of the macroscopic failure surfaces of masonry out-of-plane loaded. The behavior of a 3D system of blocks connected by interfaces is identified with a 2D Reissner–Mindlin plate. Infinitely resistant blocks connected by interfaces (joints) with a Mohr–Coulomb failure criterion with tension cut-off and compressive cap are considered. Finally, an associated flow rule for joints is adopted. In this way, the macroscopic masonry failure surface is obtained as a function of the macroscopic bending moments, torsional moments and shear forces by means of a linear programming problem in which the internal power dissipated is minimized, once that a subclass of possible deformation modes is a priori chosen. Several examples of technical relevance are presented and comparisons with previously developed Kirchhoff–Love static [Milani, G., Lourenço, P.B., Tralli, A., 2006b. A homogenization approach for the limit analysis of out-of-plane loaded masonry walls. J. Struct. Eng. ASCE (in press)] and kinematic [Sab, K., 2003.Yield design of thin periodic plates by a homogenisation technique and an application to masonry walls. C.R. Mech. 331, 641–646] failure surfaces are provided. Finally, two meaningful structural examples are reported, the first concerning a masonry wall under cylindrical flexion, the second consisting of a rectangular plate with a central opening out-of-plane loaded. For both cases, the influence of the shear strength on the collapse load is estimated.     

超高强度平头圆柱形弹体对低碳合金钢板的高速撞击实验

为分析不同组分低碳合金钢板抗超高强度低碳合金钢弹体的高速撞击性能及破坏模式,以两种典型防弹特种钢SS、AS以及常见的Q235A钢为研究对象,通过静态拉伸、静态压缩及动态压缩测试,获得静态拉伸和压缩性能参数以及1 000~6 000 s-1应变率范围内的力学行为,分析了材料组分与力学性能的相关性。采用弹道枪加载撞击方法,获得了两种超高强度合金钢平头圆柱形弹体对3种钢板(14.5~15.9 mm厚)的弹道极限速度,通过分析获得了不同工况下的极限比吸收能,讨论了合金钢板在弹体高速撞击下破坏模式的差异,分析了材料力学性能与破坏模式的相关性。研究表明:3种合金钢板抗弹体撞击性能与材料屈服强度正相关,但其性能间的差异远小于屈服强度间的差异;在超高强度合金钢平头圆柱形弹体的高速撞击下,3种钢板的失效机制与其力学性能密切相关,Si和Mn含量高的AS钢呈硬脆性特征,其断裂失效主要取决于材料的剪切强度,而Si和Mn含量较低的SS钢和Q235A钢具有良好的塑性,其断裂失效主要取决于材料的压缩强度和剪切强度。     

Minimal plate design for singular regions

A method for solving the non-linear problem of minimal volume design of sandwich plates obeying the Mises criterion is presented for the class of plates for which the orientation of the principal directions depends solely upon the position of the stress point on the yield surface. The principal bending coordinates, the yield moment, and the deflection rate at plastic collapse for the optimum design are determined in their most general form. An example is presented for a plate supported over two edges but free on its remaining edge.     

Failure criterion for reinforced concrete beams and plates subjected to membrane force,bending and shear

In this paper, a failure criterion for reinforced concrete plates is derived through the kinematic method in the framework of the limit analysis theory. This criterion is expressed in terms of the stress resultant variables: membrane force, shear force and bending moment at once. The aim of the authors is to be able to predict the failure of reinforced concrete plate structures in statics or in slow dynamics using directly the internal forces (membrane and shear forces and moment) resulting from a finite-element computation.In a first step, a beam criterion is derived. The closed form expression of the criterion shows that it is made up of two parts, one independent of the moment (i.e. depending only on the normal force and the shear force) and one depending on the normal force, the shear force and the bending moment. This structure of the criterion allows to determine two failure modes: shear failure and bending failure.Then in a second step, the beam criterion is extended to the case of reinforced concrete plates. The obtained criterion is partly numerical and partly a close form expression. It gives an upper bound of the load, and when this limit load is reached, the criterion is able to supply, on one hand, the failure mode (as seen in the beam case) and, on the other hand, the angles of the failure plane in the reinforced concrete plate section.Thirdly, the criterion is implemented in the finite element software Europlexus and validated with respect to punching experimental tests. We show that the criterion must be used with an effectiveness factor applied on the concrete compressive strength.     

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.     

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.     

Ultimate strength of postbuckling for simply supported rectangular composite thin plates under compression

1.IntroductionAsiswellthnown,thepostbucklingmetalplatescanbesuccessivelysubjectedtotheload.Forsomemetalplates,theultimatestrengthcanreachthirtytimesofthebucklingstressll].Inordertostudythebearingcapacityofthecompositeplatesthroughtestsandtheoreticalanalysis,f'ttiluretestsofthe283specimensofGFRPthinrectangularplatesundercompres.sionarecarriedout.ItisprovedthattheGFRPplatesinpostbucklingstillhavehigherbearingcapacity,forthetestedspecimensshowedultimatestrengthwhichisfourtotwenty-livetiniesth…     

非均质材料动力分析的广义多尺度有限元法

自然界和工程中的大部分材料都具有多尺度特征,当考察尺度小到一定程度后,都将表现出非均质性.针对非均质材料的动力问题,提出了一种广义多尺度有限元方法,其基本思想是利用静态凝聚法以及罚函数法构造能够反映单元内部材料非均质特性的多尺度位移基函数.与传统扩展多尺度有限元法中的基函数构造方式不同,广义多尺度有限元法的基函数无需通过在子网格域上多次求解椭圆问题得到,而可直接通过矩阵运算获得.其主要步骤如下:利用数值基函数将一个非均质单胞等效为一个宏观单元,进而形成整个结构的等效刚度矩阵,并得到宏观网格的节点位移,最后再次利用数值基函数得到微观尺度上的位移结果.该广义多尺度有限元法是扩展多尺度有限元法的一种新的拓展,可模拟具有更加复杂几何的非均质单胞的力学行为.通过数值算例,模拟了非均质材料的静力问题、广义特征值问题以及瞬态响应问题,计算结果表明:在边界条件一样的情况下,广义多尺度有限元法的计算结果与传统有限元的计算结果保持高度一致.与传统有限元相比,该方法在保证计算精度的同时极大地提高了计算效率.研究结果表明,广义多尺度有限元法能够很好地模拟非均质单胞的力学行为,具有良好的工程应用潜力.      

Elasto-plastic buckling analysis of laminated plates including interfacial damage

Elasto-plastic buckling of orthotropic laminated plates, which include interfacial damage, is analyzed in detail. Firstly, a novel mixed hardening yield criterion, as an improvement of Hill’s counterpart, is proposed for the orthotropic materials on the basis of the plastic theory. And differing from Hill’s theory, the present yield criterion is related to the spherical tensor of stress. Then, the incremental elasto-plastic constitutive relations of the mixed hardening orthotropic materials are presented. Secondly, the incremental static equilibrium equations for laminated plates including interfacial damage are established based on Von-Karman type theory and the principle of minimum potential energy. Finally, the elasto-plastic buckling of laminated plates are solved by adopting the Galerkin method and iteration scheme. The numerical results show that buckling of the plate occurs easier due to the existence of interfacial damage, and the critical load trends to constant when the interfacial damage approaches a certain degree. Also, the effect of anisotropy on buckling is obvious and the analysis of elasto-plastic buckling is necessary.     

Yield design of thin periodic plates by a homogenization technique and an application to masonry wallsCalcul à la rupture des plaques minces périodiques par une technique d'homogénéisation et une application aux murs en maçonnerie

A homogenization method for determining overall yield strength properties of thin periodic plates from their local strength properties is proposed within the framework of the yield design theory. The proposed method is applied to the determination of the in-plane and out of plane strength criterion for masonry described as a regular assemblage of infinitely resistant bricks separated by Coulomb interfaces. To cite this article: K. Sab, C. R. Mecanique 331 (2003).     

Asymptotic expansions and the possibilities to drop the hypotheses in the prandtl problem

The plane problem on the quasistatic compression of a thin perfectly plastic layer between undeformable rough plates (the Prandtl problem) has a well-known analytic solution at all points sufficiently far from the midsection and endpoints of the layer. Both the static and the kinematic component of this solution were obtained on the basis of the Prandtl hypothesis [1] stating that the tangential stress is linear along the layer thickness and is maximal in absolute value on the plate surfaces. (If the plates are perfectly rough, then this maximum value coincides with the shear yield stress.) The Prandtl hypothesis was widely confirmed in experiments carried out after the paper [1] had been published.At the same time, it is natural to ask whether one can construct a classical solution of this problem without imposing any static or kinematic hypotheses on the unknown variables and whether there exist any other mathematical solutions in which these hypotheses do not hold and which themselves are not observed in experiments.In the present paper, we use asymptotic analysis with a natural small geometric parameter and uniquely determine an exact solution (in the sense of finiteness of the number of terms in the asymptotic expansion), which coincides with the Prandtl solution generalized to the case of an arbitrary roughness coefficient of the plates. We rigorously show that such asymptotics cannot hold near the layer midsection, where we construct another, internal asymptotic expansion. In the abovementioned sense, the solution corresponding to the internal expansion is also exact and models the compression of a thin vertical strip in the middle of the layer. We realize two possible versions of matching of the two expansions in the cross-section whose distance from the midsection is equal to the layer thickness.     

Strength homogenization of matrix-inclusion composites using the linear comparison composite approach

A homogenization procedure to estimate the macroscopic strength of nonlinear matrix-inclusion composites with different strength characteristics of the matrix and inclusions, respectively, is presented in this paper. The strength up-scaling is formulated within the framework of the yield design theory and the linear comparison composite (LCC) approach, introduced by Ponte Castaneda (2002) and extended to frictional models by Ortega et al. (2011), which estimates the macroscopic strength of composite materials in terms of an optimally chosen linear thermo–elastic comparison composite with a similar underlying microstructure. In the paper various combinations for the underlying material behavior for the individual phases of the composite are considered: The matrix phase can be a quasi frictional material characterized either by a Drucker–Prager-type (hyperbolic) or an elliptical strength criterion, which predicts a strength limit also in hydrostatic compression, while the inclusion phase either may represent empty pores, pore voids filled with a pore fluid, rigid inclusions, or solid inclusions, whose strength characteristics also maybe described by a Drucker–Prager-type or an elliptical strength criterion. For generating the homogenized strength criterion efficiently in such general cases of matrix-inclusion composites, a novel algorithm is proposed in the paper. The validation of the proposed strength homogenization procedure for selected combinations of strength characteristics of the matrix material and the inclusions is conducted by comparisons with experimental results and alternative existing strength homogenization models.     

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.     

Non-quadratic Kelvin modes based plasticity criteria for anisotropic materials

Novel (non-quadratic) plasticity criteria based on Kelvin modes are formulated here for anisotropic materials. As an example, such a macroscopic criterion is applied with success to the case of FCC nickel-base single crystals. Indeed, relying on the cubic symmetry of the material, the Kelvin decomposition of elasticity tensor easily allows for the definition of an objective and loading independent criterion. The criterion identification is performed from different loading cases for CMSX2 single crystal superalloy. Tension-torsion yield surfaces at room temperature and yield stress dependence on crystal orientation are modeled. The Kelvin modes based criterion is compared to experimental data, to Hill and Barlat and coworkers macroscopic criteria and to Schmid law predictions. The results show that a simple three-parameter yield function built thanks to von Mises equivalent Kelvin stresses accounts for a satisfying plasticity criterion for such alloys.Non-quadratic norm ∥·∥_a plasticity framework is addressed. Intrinsic generalizations of Hershey-Hosford criterion are proposed for cubic material symmetry.     

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

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

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.