材料非线性对复合材料层合板热自由边界效应的影响

本文用准三维有限元法研究了材料非线性对复合材料层合板热自山边界效应的影响,给出了修正型Hahn-Tsai非线性应力-应变关系的三维形式。由本文非线性分析方法得到的层间应力与以往由线性分析方法得到的层间应力做了比较,结果表明:材料非线性能显著降低层间剪应力的集中程度,但对层间正应力影响不太明显。     

Interlaminar hygrothermal stresses in laminated plates

Within the elasticity formulation the most general displacement field for hygrothermal problems of long laminated composite plates is presented. The equivalent single-layer theories are then employed to determine the global deformation parameters appearing in the displacement fields of general cross-ply, symmetric, and antisymmetric angle-ply laminates under thermal and hygroscopic loadings. Reddy’s layerwise theory is subsequently used to determine the local deformation parameters of various displacement fields. An elasticity solution is also developed in order to validate the efficiency and accuracy of the layerwise theory in predicting the interlaminar normal and shear stress distributions. Finally, various numerical results are presented for edge-effect problems of several cross-ply, symmetric, and antisymmetric angle-ply laminates subjected to uniform hygrothermal loads. All results indicate high stress gradients of interlaminar normal and shear stresses near the edges of laminates.     

Three-dimensional photoelastic analysis of a fiber-reinforced composite model

A three-dimensional photoelastic analysis using the “stress-freezing” technique was conducted to determine the stress distributions in the matrix of a unidirectionally fiber-reinforced composite model subjected to matrix shrinkage and normal transverse loading. The model, consisting of a square array of polycarbonate rods in an epoxy matrix, simulated a boron-filament-reinforced plastic composite with a fiber-volume fraction of 0.50 at the critical temperature of the matrix epoxy. The effects of matrix shrinkage were separated from those of external loading by analyzing two identical models, one loaded and the other unloaded. The Lamé-Maxwell equations of equilibrium were used to separate stresses along axes of symmetry on interior transverse slices. Axial stress components were obtained by subslicing. Results are presented in dimensionless form by dividing the stresses by the average stress through the section. A comparison with theoretical results for a boron-epoxy composite shows excellent agreement, although Poisson's ratio of the model matrix is appreciably different from that of the prototype (0.5 compared to 0.35). One significant result was that the maximum stress occurs in the middle of the matrix section between fibers which is at variance with the theoretical prediction of maximum stress at the interface. Stress-concentration factors vary from 1.80 at the interface to 2.0 at the midpoint of the matrix section between fibers.     

Fracture process zone concept and delamination of composite laminates

Delaminations occur at the free edges of composite laminates because of high interlaminar stresses, which develop due to a mismatch in elastic properties between adjacent plies and the limited toughness of the matrix. The process of delamination is modelled using the concept of a fracture process zone. A finite element analysis that can be used with a constant or linear stress-displacement relationship within the fracture process zone is presented. It is found that the simple constant stress model gives an accurate prediction of the delamination on the free edge of a tensile specimen.     

3-D numerical study on the bending of symmetric composite laminates

３－ＤＮＵＭＥＲＩＣＡＬＳＴＵＤＹＯＮＴＨＥＢＥＮＤＩＮＧＯＦＳＹＭＭＥＴＲＩＣＣＯＭＰＯＳＩＴＥＬＡＭＩＮＡＴＥＳＣｈｉｅｎＷｅｉ－ｚａｎｇ（钱伟长）;ＨｕａｎｇＱｉａｎ（黄默）;ＦｅｎｇＷｅｉ（冯伟）（ＳｈａｎｇｈａｉＵｎｉｖｅｒｓｔｙｏｆＴｅｃ...     

An axisymmetric stress analysis in a single fibre composite of finite length under a thermal expansion mismatch

In this paper, an exact solution is derived for the characterization of thermal stresses in a single-fibre composite of finite length. All the required boundary and interfacial conditions of the thermo-elastic problem are thus satisfied exactly. The proposed method involves a particular solution that is added to a three-dimensional (3D) complementary displacement field which satisfies automatically the Navier's equations. Based on experimental data provided by fibre Bragg grating sensors, an axisymmetric analysis is used then to determine the residual stress field inside the composite due to matrix shrinkage. The numerical results clearly indicate that all stress components vary significantly near the ends. An abrupt change of the shear stresses is thus predicted close to the edges. The results of the model are also found to be in good agreement with those obtained from finite element simulations. A comparison of the proposed approach with three other published theoretical models is also presented.     

Thermal stress analysis of a silicon carbide/aluminum composite

Thermal deformations and stresses were studied in a silicon-carbide/aluminum filamentary composite at temperatures up to 370°C (700°F). Longitudinal and transverse thermal strains were measured with strain gages and a dilatometer. An elastoplastic micromechanical analysis based on a one-dimensional rule-of-mixtures model and an axisymmetric two-material composite cylinder model was performed. It was established that beyond a critical temperature thermal strains become nonlinear with decreasing longitudinal and increasing transverse thermal-expansion coefficients. This behavior was attributed to the plastic stresses in the aluminum matrix above the critical temperature. An elastoplastic analysis of both micromechanical models was performed to determine the stress distributions and thermal deformation in the fiber and matrix of the composite. While only axial stresses can be determined by the rule-of-mixtures model, the complete triaxial state of stress is established by the composite cylinder model. Theoretical predictions for the two thermal-expansion coefficients were in satisfactory agreement with experimental results.     

Interlaminar fracture toughness of graphite/epoxy composite under mixed-mode deformations

An antisymmetric test fixture is employed to investigate interlaminar fracture behavior in graphite/epoxy composite material under mixed-mode deformations. Finite correction factors for the graphite/epoxy fracture specimen with various crack lengths are used to determine the interlaminar fracture toughness by finite-element stress analysis. Interlaminar fracture characteristics of graphite/epoxy composite material under mode-I, mode-II and mixed-mode deformations are evaluated experimentally. A mixed-mode fracture criterion is also investigated to obtain information on mixedmode interlaminar fracture behavior of graphite/epoxy composite material.Paper was presented at the 1988 SEM Spring Conference on Experimental Mechanics held in Portland, OR on June 5–10.     

层合圆柱厚壳热应力分析的状态方程及其半解析法

通过对Ｈｅｌｌｉｎｇｅｒ－Ｒｅｉｓｓｎｅｒ变分原理的修正，导出了变温作用下层合圆柱厚壳的状态方程及其半解析法，该方法在ｚ－θ曲面内采用通常的有限元离散，而沿壳厚（ｒ）方向采用状态空间法给出解析解答，且通过采用状态转移矩阵，建立了层合圆柱壳内外表面应力和位移之间的关系式，然后利用打靶法进行求解，从而大大降低了计算中的未知量数目。     

压电智能复合材料层合梁的力电耦合模型及层间应力分析

由于非凡的物理性能,石墨烯纳米片(GPL)被认为是最有吸引力的复合材料增强材料之一.GPL增强材料可以明显提高聚偏氟乙烯(PVDF)压电性能和力学性能.在力电载荷作用下,对含均匀石墨烯薄片增强(GSR)智能压电复合材料层合梁层间应力预测至关重要.若对受到力电耦合作用且层与层之间材料性能突变的压电层合梁层间剪切变形预测有误,则其层间应力过大可能导致层间失效.因此,论文提出一种适于分析此类问题且满足层与层之间相容性条件的有效力电耦合模型,用于含GSR致动器的复合材料层合梁层间应力分析.应用Reissner混合变分原理(RMVT),可以提高考虑力电耦合效应的横向剪应力预测精度.三维(3D)弹性理论和所选模型计算结果将用于评估所提梁模型性能.此外,还从力电载荷、压电层厚度、石墨烯体积分数和长厚比等方面对含GSR致动器复合材料层合梁力学响应特性进行了系统的研究.     

Inelastic behavior of an AS4/PEEK composite under combined transverse compressionand shear.Part II: modeling

A micromechanical model is presented for simulating the nonlinearities exhibited by AS4/PEEK composites in shear and transverse compression, their interaction, and their rate dependence at room temperature. The fibers are assumed to be transversely isotropic and to be distributed in a hexagonal pattern in the matrix. The PEEK matrix is modeled as an elastic–powerlaw viscoplastic, isotropic solid through two related models. Model I is the simple J₂–type viscoplasticity; Model II is a rate dependent version of the non-associative Drucker-Prager model. Both models are calibrated so that they reproduce the shear response of the composite. Model II is also calibrated to its transverse compression response. Both models capture the rate dependence of the composite well. Model I is significantly less stiff in transverse compression than the experimental data. However, it does a reasonable job of predicting other aspects of the biaxial experiments and captures the important trends of behavior observed. Model II does better in transverse compression, but shearing in the presence of transverse compression is found to be stiffer than the measured responses. The unit cell model allows us to examine the stresses in the composite, providing an explanation for the lack of interaction between the constant stress and the increasing stress observed experimentally for certain loading paths.     

Steady-state aeroelasticity of fluid flow over a laminated composite plate

A stress-function-based variational approach is used to determine the interlaminar stresses in a multilayered strip of laminate subject to arbitrary combinations of axial extension, bending, and the steady-state aerodynamic loading of fluid flow over the upper surface of laminated composite plate. Symmetric four-layer, cross-ply and angle-ply laminates are considered in details. Some numerical solutions by using a personal computer are obtained. The present results for four-layer laminates show that the aerodynamic loading has significant effects on the interlaminar stresses near the free-edge regions.     

Theory of elastic and photoelastic isodynes. Samples of application in composite structures

The theoretical foundations of a system of new methods which the authors call isodyne photoelasticity are presented. It is shown that the plane-elastic isodynes related to two-dimensional stress fields can be obtained experimentally using particular scattered-light techniques. The experimentally obtained photoelastic isodynes allow non-destructive analysis of stresses in structures with local three-dimensional stress states.Isodyne photoelasticity is used to determine stresses in composite structures, in particular the stresses at internal discontinuities and the delaminating stresses. The results are used to assess the reliability of some analytical solutions. The theories of isodynes and of the developed experimental techniques of isodyne photoelasticity are presented in a manner compatible with contemporary concepts of the models of reality—physical models, mathematical models, and experimental models constructed of matter and energy.Paper was presented at V International Congress on Experimental Mechanics held in Montreal, Quebec, Canada on June 10–15, 1985.     

粘弹层合板的稳态振动和层间应力

应用混合分层理论和Ressiner混合变分原理，在板厚方向取二次位移插值函数和三次、四次横向应力插值函数推导出粘弹层合板的动力学方程，得出简支粘弹层合板稳态振动的解。不仅得出与三层弹性板精确的自振频率吻合良好的解，而且对于粘弹层合板，所计算的自振频率和结构损耗因子也与三维结果吻合较好。计算了自由阻尼层合板对应的低阶法向位移响应幅值和层问横向应力的幅值。结果表明，较高的层间横向正应力是低频稳态振动中引起粘弹层合板分层破坏的主要因素，采用适当模量和厚度的粘弹性材料将有效地降低粘弹层合板的层间横向正应力的幅值。     

基于树脂层强度预测分层萌生载荷

 轴向拉伸下的层合板在自由边缘附近存在奇异应力场,容易产生分层萌生进而导致结构破坏．而大多数基于材料力学方法的分层萌生研究对强度参数的确定有较强的主观性,缺乏合理解释．本文通过在层间插入一定厚度的树脂层,将分层萌生视为层间树脂层在三维应力状态下的强度破坏,根据有明确试验标准的树脂强度和Mohr 判据判断是否发生分层萌生．采用Pipes-Pagano 的均匀轴向拉伸模型计算应力场,对层间应力与树脂层面外应力进行了对比．取临界长度为4 个单层板厚度对T800/914 层合板的分层萌生进行了预测,结果表明预测值与实验值吻合良好．     

Three-dimensional scattered-light stress analysis of discontinuous fiber-reinforced composites

This paper reports an investigation of three-dimensional stresses in models of composite materials with discontinuous fibers using the scattered-light photoelastic method. A special scattered-light polariscope with loading system was designed and constructed for this research. Two models were used in this investigation. The first was a cylindrical tension model with a single discontinuous axially aligned aluminum fiber surrounded by a polyester matrix. The second was a five-fiber cylindrical compression model with a central discontinuous fiber surrounded by four continuous ones and polyester matrix. The experimentally determined stress distribution from each model along chosen lines was presented. The stress distribution on the same lines in the single-fiber case was calculated using the finite-element method. The calculated result showed fair agreement with the experimentally determined results.     

An analytic solution for interlaminar stresses in a fiber-reinforced double-layer cylindrical shell

The present paper investigates the basic characteristics of the interlaminar stresses in a double-layer cylindrical shell with both ends simply-supported under uniform external or internal pressure. The double layer shell is composed of a 0° fiber-reinforced composite layer and an isotropic layer. In this paper, this axisymmetric problem is solved exactly with the three-dimensional theory of elasticity. Both the displacement field and the stress field of each layer of the shell are expressed in Fourier series and Fourier-Bessel series. Then we illustrate the effects of the various parameters, such as geometry, material constants, loading conditions and stacking sequence, on the interlaminar stresses.     

Numerical simulation of elasto-plastic deformation of composites: evolution of stress microfields and implications for homogenization models

The deformation of a composite made up of a random and homogeneous dispersion of elastic spheres in an elasto-plastic matrix was simulated by the finite element analysis of three-dimensional multiparticle cubic cells with periodic boundary conditions. “Exact” results (to a few percent) in tension and shear were determined by averaging 12 stress-strain curves obtained from cells containing 30 spheres, and they were compared with the predictions of secant homogenization models. In addition, the numerical simulations supplied detailed information of the stress microfields, which was used to ascertain the accuracy and the limitations of the homogenization models to include the nonlinear deformation of the matrix. It was found that secant approximations based on the volume-averaged second-order moment of the matrix stress tensor, combined with a highly accurate linear homogenization model, provided excellent predictions of the composite response when the matrix strain hardening rate was high. This was not the case, however, in composites which exhibited marked plastic strain localization in the matrix. The analysis of the evolution of the matrix stresses revealed that better predictions of the composite behavior can be obtained with new homogenization models which capture the essential differences in the stress carried by the elastic and plastic regions in the matrix at the onset of plastic deformation.     

考虑界面损伤层合梁板的层间应力分析

基于精确应力分析的广义六自由度板理论，应用变分原理和损伤力学中 的应变等效原理，考虑复合材料铺设层内和层间界面处的损伤效应，建立了具两种损伤模式 的复合材料层合板的三维非线性平衡微分方程，且运用有限差分法对考虑损伤简支层合梁板 的层间应力进行了求解.     

非理想界面对三相复合材料应力场的影响

对非理想界面的三相复合材料,提出了计算弹性应力场的微观力学模型,在适当的简化假设下,对带界相的颗粒增强和纤维增强复合材料,得到了应力场的计算公式。以剪切载荷为例给出了数值例子。给出的数值结果表明非理想界面对三相复合材料应力场的影响。