On the axial propagation of kink bands in fiber composites : Part ii analysis

A new micromechanical model is presented to simulate the steady-state axial propagation of kink bands investigated experimentally in the accompanying paper (Part I) . The fibers are in a hexagonal array and are assumed to be isotropic and linearly elastic, while the matrix is modeled as an elastic-powerlaw viscoplastic solid. Matrix properties for the model are determined from shear tests on the composite and compression tests on neat PEEK. The model is used to predict the propagation stress (σ_P) of the AS4⧹PEEK composite and to investigate the sensitivity of σ_P to band inclination, matrix properties, and loading rate. A simple model recently reported in the literature is calibrated to the current material system and compared with the present experimental data and model predictions. The micromechanical model is found to predict the propagation stress reasonably well and to capture the rate dependence of the composite. The simple model is found to capture the trends of the behavior.     

Inelastic behavior of an AS4/PEEK composite under combined transverse compression and shear. Part I: experiments

The nonlinear behavior in shear and transverse compression of unidirectional AS4/PEEK and their interaction are investigated experimentally. The composite is rate dependent even at room temperature and its rate exponent is similar to that of neat PEEK. The material is tested under pure shear, pure compression and under biaxial loading histories. The biaxial tests are performed in a custom facility on thin strips of the material. The facility allows freedom to choose the loading path in the biaxial stress and strain spaces of interest. Tests are performed for three biaxial loading paths. In the first, the specimen is sheared then compressed while the shear stress is held constant; in the second, the specimen is compressed then sheared while the compressive stress is held constant; and in the third, the specimen is loaded simultaneously by proportional amounts of compression and shear. It was found that the induced deformation is influenced significantly by the loading history followed. Also, initial loading in shear or compression has only a modest effect on subsequent loading of the other type. An unorthodox yielding behavior for the composite results from this lack of interaction. Finally, the stresses at failure are found to trace an elliptical path in the shear–transverse compression plane, but the failure stress state is not significantly affected by the loading path followed.     

Modeling Transverse Compression of Cylindrical Bodies in Bending

A nonclassical theory of the stress–strain state of composite beams in bending is developed. Both transverse shear and transverse compression are taken into account. The results demonstrate the high accuracy of the theory in investigation of short transversally isotropic beams     

Solutions for effective shear properties in three phase sphere and cylinder models

Solutions are presented for the effective shear modulus of two types of composite material models. The first type is that of a macroscopically isotropic composite medium containing spherical inclusions. The corresponding model employed is that involving three phases: the spherical inclusion, a spherical annulus of matrix material and an outer region of equivalent homogeneous material of unlimited extent. The corresponding two-dimensional, polar model is used to represent a transversely isotropic, fiber reinforced medium. In the latter case only the transverse effective shear modulus is obtained. The relative volumes of the inclusion phase to the matrix annulus phase in the three phase models are taken to be the given volume fractions of the inclusion phases in the composite materials at large. The results are found to differ from those of the well-known Kerner and Hermans formulae for the same models. The latter works are now understood to violate a continuity condition at the matrix to equivalent homogeneous medium interface. The present results are compared extensively with results from other related models. Conditions of linear elasticity are assumed.     

Nonlinear transient response of strain rate dependent composite laminated plates using multiscale simulation

The effects of strain rate dependency and inelasticity on the transient responses of composite laminated plates are investigated. A micromechanics model which accounts for the transverse shear stress effect, the effect of strain rate dependency and the effect of inelasticity is used for analyzing the mechanical responses of the fiber and matrix constituents. The accuracy of the micromechanics model under transverse shear loading is verified by comparing the results with those obtained using a general purpose finite element code. A higher order laminated plate theory is extended to capture the inelastic deformations of the composite plate and is implemented using the finite element technique. A complete micro–macro numerical procedure is developed to model the strain rate dependent behavior of inelastic composite laminates by implementing the micromechanics model into the finite element model. Parametric studies of the transient responses of composite plates are conduced. The effects of geometry, ply stacking sequence, material models, boundary conditions and loadings are investigated. The results show that the strain rate dependency and inelasticity influence the transient responses of composite plates via two significantly different mechanisms.     

Mechanical response of neo-Hookean fiber reinforced incompressible nonlinearly elastic solids

The mechanical behavior of an incompressible neo-Hookean material, directionally reinforced by neo-Hookean fibers, is examined under homogeneous deformations. A composite model for this transversely isotropic material is developed based on a multiplicative decomposition of the deformation gradient which considers interaction between the fiber and the matrix. The so-called standard reinforcing model exhibits non-monotonic behavior in compression. The present composites-based approach leads to a modification of the standard reinforcing model in which monotonic behavior in compression is observed. This stems from the micromechanical basis of the model in which the fiber is treated as a neo-Hookean material. The conditions for loss of monotonicity and positivity in the stress-shear behavior in off-axis simple 2D shear are also obtained.     

Open-cell cellular solids: A constitutive equation for hyperelasticity with deformation induced anisotropy

A constitutive theory is developed for an open-cell flexible cellular solid consisting of a network of struts each connecting two vertex points. A hypothesis is proposed that vertex points move affinely in the large-deformation regime, when the struts buckle, and that the force carried by a strut is a function of the longitudinal and rotational change of its vertex-to-vertex vector. The forces consist of one longitudinal force, parallel with the vertex-to-vertex vector of the strut and one transverse force. The overall stress response is initially dominated by the longitudinal force whilst the addition of the transverse force becomes significant at large deformations. The model contains three parameters: longitudinal stiffness, bending stiffness and critical stretch of a strut. These three parameters are calibrated against a simple compression test. The model is then validated against independent experiments in a simple tension, simple shear and a combined shear-compression test on an isotropic flexible polyether urethane foam. Excellent agreement is obtained between the experiments and the model.     

形状记忆合金纤维复合材料的等效力学行为

在Aboudi提出的胞元模型以及Liu等建立的形状记忆合金的本构模型的基础上，由Legendre多项式，假设每个子胞元的位移场、应变场和应力场，再由子胞元间交界面的应力连续条件和外荷载边界条件推导出基体为弹塑性材料的形状记忆合金纤维复合材料的胞元模型；模拟了呈周期对称的形状记忆合金纤维复合材料受轴向单向拉伸、横向拉伸和横向剪切荷载作用下的等效力学行为，与有限元解进行了比较，结果基本一致。与有限元法比较起来，本文推导出的形状记忆合金纤维复合材料的胞元模型更具高效性。     

Cracks’ closure in 3-D fracture dynamics: the effect of relative location of two coplanar cracks

In this paper, the problem of two equal coplanar cracks with allowance for the crack faces contact interaction was investigated. The problem of the cracks located in homogeneous, isotropic, and linearly elastic solid subjected to normally incident tension–compression wave is solved by the boundary integral equations method. The influence of the distance between two cracks on the stress intensity factors (opening mode and transverse shear mode) is studied for a range of wave numbers. The results are compared with those obtained neglecting cracks’ closure.     

On constitutive modelling of porous neo-Hookean composites

In this paper a hyperelastic constitutive model is developed for neo-Hookean composites with aligned continuous cylindrical pores in the finite elasticity regime. Although the matrix is incompressible, the composite itself is compressible because of the existence of voids. For this compressible transversely isotropic material, the deformation gradient can be decomposed multiplicatively into three parts: an isochoric uniaxial deformation along the preferred direction of the material (which is identical to the direction of the cylindrical pores here); an equi-biaxial deformation on the transverse plane (the plane perpendicular to the preferred direction); and subsequent shear deformation (which includes “along-fibre” shear and transverse shear). Compared to the multiplicative decomposition used in our previous model for incompressible fibre reinforced composites [Guo, Z., Peng, X.Q., Moran, B., 2006, A composites-based hyperelastic constitutive model for soft tissue with application to the human annulus fibrosus. J. Mech. Phys. Solids 54(9), 1952–1971], the equi-biaxial deformation is introduced to achieve the desired volume change. To estimate the strain energy function for this composite, a cylindrical composite element model is developed. Analytically exact strain distributions in the composite element model are derived for the isochoric uniaxial deformation along the preferred direction, the equi-biaxial deformation on the transverse plane, as well as the “along-fibre” shear deformation. The effective shear modulus from conventional composites theory based on the infinitesimal strain linear elasticity is extended to the present finite deformation regime to estimate the strain energy related to the transverse shear deformation, which leads to an explicit formula for the strain energy function of the composite under a general finite deformation state.     

Influence of the loading path on the strength of fiber-reinforced composites subjected to transverse compression and shear

The influence of the loading path on the failure locus of a composite lamina subjected to transverse compression and out-of-plane shear is analyzed through computational micromechanics. This is carried out using the finite element simulation of a representative volume element of the microstructure, which takes into account explicitly fiber and matrix spatial distribution within the lamina. In addition, the actual failure mechanisms (plastic deformation of the matrix and interface decohesion) are included in the simulations through the corresponding constitutive models. Two different interface strength values were chosen to explore the limiting cases of composites with strong or weak interfaces. It was found that failure locus was independent of the loading path for the three cases analyzed (pseudo-radial, compression followed by shear and shear followed by compression) in the composites with strong and weak interfaces. This result was attributed to the fact that the dominant failure mechanism in each material was the same in transverse compression and in shear. Failure is also controlled by the same mechanisms under a combination of both stresses and the failure locus depended mainly on the magnitude of the stresses that trigger fracture rather than in the loading path to reach the critical condition.     

Localization,delocalization and compression fracture in externally pressurized thick cross-ply (very) long cylindrical shells with material non-linearity: A multi-scale and multi-physics analysis

Compression fracture in carbon fiber reinforced plastics (CFRP) involves multiple physical mechanisms operating at multiple scales ranging from angströms to cms and beyond. First, at the macro/meso-scale, combined effects of modal imperfections, transverse shear/normal deformation along with the non-linear hypo-elastic transverse shear (G_TT) material property on the emergence of interlaminar shear crippling type instability modes, related to the localization (onset of deformation softening), delocalization (onset of deformation hardening) and propagation of mode II compression fracture/damage, in thick imperfect cross-ply very long cylindrical shells under applied hydrostatic pressure, are investigated. The primary accomplishment is the (hitherto unavailable) computation of the layer-wise mode II stress intensity factor, energy release rate and kink crack band-width, under hydrostatic compression, from a non-linear finite element analysis (FEA), using Maxwell’s construction and Griffith׳s energy balance approach. Numerical results include the effects of hypoelastic (G_TT only) material property, on localization and delocalization leading to compression fracture.At the micro-scale, a novel three-dimensional eigenfunction expansion technique, based in part on separation of the cylinder length-variable and partly utilizing a modified Frobenius type series expansion in conjunction with an affine transformation to compute the local stress singularity, in the vicinity of a kinked-fiber/matrix trimaterial junction front. Such computed stress singularities represent a measure of the degree of inherent flaw sensitivity of unidirectional CFRP under compression. Finally, dislocation glide in graphite crystallites plays a dominant role in kink band nucleation and propagation at the nano-meter scale.     

Micromechanical modeling of the elasto-viscoplastic behavior of semi-crystalline polymers

A micromechanically based constitutive model for the elasto-viscoplastic deformation and texture evolution of semi-crystalline polymers is developed. The model idealizes the microstructure to consist of an aggregate of two-phase layered composite inclusions. A new framework for the composite inclusion model is formulated to facilitate the use of finite deformation elasto-viscoplastic constitutive models for each constituent phase. The crystalline lamellae are modeled as anisotropic elastic with plastic flow occurring via crystallographic slip. The amorphous phase is modeled as isotropic elastic with plastic flow being a rate-dependent process with strain hardening resulting from molecular orientation. The volume-averaged deformation and stress within the inclusions are related to the macroscopic fields by a hybrid interaction model. The uniaxial compression of initially isotropic high density polyethylene (HDPE) is taken as a case study. The ability of the model to capture the elasto-plastic stress-strain behavior of HDPE during monotonic and cyclic loading, the evolution of anisotropy, and the effect of crystallinity on initial modulus, yield stress, post-yield behavior and unloading-reloading cycles are presented.     

Corotational nonlinear analyses of laminated shell structures using a 4-node quadrilateral flat shell element with drilling stiffness

A new 4-node quadrilateral flat shell element is developed for geometrically nonlinear analyses of thin and moderately thick laminated shell structures. The fiat shell element is constructed by combining a quadrilateral area co- ordinate method （QAC） based membrane element AGQ6- II, and a Timoshenko beam function （TBF） method based shear deformable plate bending element ARS-Q12. In order to model folded plates and connect with beam elements, the drilling stiffness is added to the element stiffness matrix based on the mixed variational principle. The transverse shear rigidity matrix, based on the first-order shear deformation theory （FSDT）, for the laminated composite plate is evaluated using the transverse equilibrium conditions, while the shear correction factors are not needed. The conventional TBF methods are also modified to efficiently calculate the element stiffness for laminate. The new shell element is extended to large deflection and post-buckling analyses of isotropic and laminated composite shells based on the element independent corotational formulation. Numerical re- sults show that the present shell element has an excellent numerical performance for the test examples, and is applicable to stiffened plates.     

On the extension of the Gurson-type porous plasticity models for prediction of ductile fracture under shear-dominated conditions

One of the major drawbacks of the Gurson-type of porous plasticity models is the inability of these models to predict material failure under low stress triaxiality, shear dominated conditions. This study addresses this issue by combining the damage mechanics concept with the porous plasticity model that accounts for void nucleation, growth and coalescence. In particular, the widely adopted Gurson–Tvergaard–Needleman (GTN) model is extended by coupling two damage parameters, representing the volumetric damage (void volume fraction) and the shear damage, respectively, into the yield function and flow potential. The effectiveness of the new model is illustrated through a series of numerical tests comparing its performance with existing models. The current model not only is capable of predicting damage and fracture under low (even negative) triaxiality conditions but also suppresses spurious damage that has been shown to develop in earlier modifications of the GTN model for moderate to high triaxiality regimes. Finally the modified GTN model is applied to predict the ductile fracture behavior of a beta-treated Zircaloy-4 by coupling the proposed damage modeling framework with a recently developed J₂–J₃ plasticity model for the matrix material. Model parameters are calibrated using experimental data, and the calibrated model predicts failure initiation and propagation in various specimens experiencing a wide range of triaxiality and Lode parameter combinations.     

基于t准则的各向异性强度准则及变换应力法

岩土材料通常呈现出成层水平分布特点, 即可将其视为横观各向同性材料, 横观各向同性对于岩土材料的变形以及强度值都会产生显著的影响. 基于已提出的t强度准则, t强度准则是基于各向同性单元体中存在有效滑移面来构建的, 并根据该空间有效滑移面上主剪应力与主法向应力的比值达到一定阈值为破坏条件. 在空间中存在有效滑移面与物理沉积面, 基于上述两个面在空间的位置关系, 用两面夹角作为表征横观各向同性对剪切强度影响程度的参量, 并假定当该夹角值越大, 则各向异性对强度贡献程度越大, 对应更大的应力比强度值, 反之, 则对应更小的应力比强度值. 基于上述思路并类比将其推广为正交三维各向异性准则, 基于三维各向异性材料的三维沉积面, 提出了三维特征沉积面的概念, 并基于空间滑移面与三维特征沉积面之间的夹角作为度量各向异性程度的变量, 提出了基于两面角作为参量考虑原生各向异性的应力比强度公式, 并利用该应力比强度公式来修正已提出的t强度准则, 最终建立了考虑各向异性影响的t准则公式. 在上述准则基础上, 考虑将各向异性应力空间转换为各向同性应力空间的思路, 在各向异性t准则基础上, 推导得到了基于各向异性强度t准则的变换应力公式, 利用变换应力公式可以将传统的以p, q为变量的各向同性本构模型转变为可考虑各向异性的三维本构模型. 通过对岩土材料的强度以及真三轴条件下的应力应变关系试验数据预测, 验证了所提的各向异性t准则及其变换应力公式的有效性及适用性.      

Fundamental solutions of the thermoelastic equations for transversely isotropic plates

The problem of thermoelasticity for transversely isotropic plates acted upon by concentrated heat sources is solved. The {1, 2}-order equations of thermoelasticity that incorporate the transverse shear and normal stresses are used. A bending heat source with symmetric heat transfer is considered. The dependence of thermal stress components on the thermal and thermomechanical parameters of transversely isotropic plates is studied     

ANISOTROPIC STRENGTH CRITERION BASED ON T CRITERION AND THE TRANSFORMATION STRESS METHOD 1)

岩土材料通常呈现出成层水平分布特点, 即可将其视为横观各向同性材料, 横观各向同性对于岩土材料的变形以及强度值都会产生显著的影响. 基于已提出的t强度准则, t强度准则是基于各向同性单元体中存在有效滑移面来构建的, 并根据该空间有效滑移面上主剪应力与主法向应力的比值达到一定阈值为破坏条件. 在空间中存在有效滑移面与物理沉积面, 基于上述两个面在空间的位置关系, 用两面夹角作为表征横观各向同性对剪切强度影响程度的参量, 并假定当该夹角值越大, 则各向异性对强度贡献程度越大, 对应更大的应力比强度值, 反之, 则对应更小的应力比强度值. 基于上述思路并类比将其推广为正交三维各向异性准则, 基于三维各向异性材料的三维沉积面, 提出了三维特征沉积面的概念, 并基于空间滑移面与三维特征沉积面之间的夹角作为度量各向异性程度的变量, 提出了基于两面角作为参量考虑原生各向异性的应力比强度公式, 并利用该应力比强度公式来修正已提出的t强度准则, 最终建立了考虑各向异性影响的t准则公式. 在上述准则基础上, 考虑将各向异性应力空间转换为各向同性应力空间的思路, 在各向异性t准则基础上, 推导得到了基于各向异性强度t准则的变换应力公式, 利用变换应力公式可以将传统的以p, q为变量的各向同性本构模型转变为可考虑各向异性的三维本构模型. 通过对岩土材料的强度以及真三轴条件下的应力应变关系试验数据预测, 验证了所提的各向异性t准则及其变换应力公式的有效性及适用性.     

基于临界状态的砂土本构模型研究

砂土孔隙比及所受压力是其力学特性的重要影响因素. 本文基于砂土临界状态线特性分析,采用以e-(p/p_a)ξ平面内的线性关系描述其等向压缩线. 通过对比分析两种不同压缩线函数 与临界状态线函数之间的关系提出更适合描述砂土在等向压缩下的参考压缩线,并给出了基于参考压缩线的等向硬化规律. 建议了适用于 描述砂土剪切特性的屈服面函数,并给出利用等向压缩和等p路径确定屈服面形状参数μ的方法. 将不同应力比对应的压缩线作为砂土状态参量参考线,以获取潜在强度M_f与特征状态应力比M_c,进而描述砂土压缩与剪切特性;基于等向压缩与等p路径建立了当前应力比与状态参量参考线之间的相关关系,从而实现了砂土状 态参量参考线由参考压缩线向临界状态线平稳过渡. 建立的砂土本构模型共11个参数,均能够通过常规土工试验或经验获取. 基于模型预测与Toyoura砂的等向压缩、三轴不排水剪切试验及排水剪切试验的对比结果,本文建立的砂土本构模型很好地描述了Toyoura 砂在不同孔隙比和不同压力下的压缩与剪切特性.