Nonlinear viscoelastic multi-scale repetitive unit cell model of 3D woven composites with damage evolution

Three-dimensional (3D) textile composites have great potential applications to aircrafts and high speed vehicles because of the high strength/weight ratios and the capabilities of manufacturing complex, net-shape preforms. This paper reports the nonlinear viscoelastic responses and damage mechanisms of one kind of 3D textile composites, named as 3D orthogonal woven composite (3DOWC) under quasi-static tensile loading based on a micro/meso-scale repetitive unit cells (RUCs) model. In the RUCs model, the resin is described with a nonlinear viscoelastic material and the fibers/tows with an elastic material. The damage initiation and propagation in resin are simulated by the post-damage constitutive models with maximum principal theory failure criteria. The fibers/tows impregnated with resin are defined by elastic transverse-isotropic material model with ultimate strengths failure of ‘expanded smeared crack’ both along and perpendicular to fibers/tows axis direction. The engineering parameters and ultimate strengths of homogenized fibers/tows filled with matrix in meso-RUCs model are transferred from the numerical analysis of the micro-RUCs. The results are compared with experimental and theoretical values of RUC deformation and damage initiation and propagation under monotonic axial loading. The methodology of establishing the nonlinear visco-elastic multi-scale model of 3D textile composites without introducing the real fabric architecture in finite element analyses is explained. With the multi-scale RUCs model, the mechanical behaviors of other kinds of 3D textile composites can also be predicted.     

A micromechanical analysis of the nonlinear elastic and viscoelastic constitutive relation of a polymer filled with rigid particles

Micromechanical theory is applied to study the nonlinear elastic and viscoelastic constitutive relations of polymeric matrix filled with high rigidity solid particles. It is shown that Eshelby's method can be extended to the case of nonlinear matrix and Eshelby's tensor still exists provided that Poisson's ratio of the nonlinear matrix assumes constant value in deforming process and the rigidity of elastic filling particles is much higher than that of the matrix. A new method for averaging process is proposed to overcome the difficulty that occured in applying the ordinary equivalent inclusion method or the self-consistant method to nonlinear matrices. A rather simple constitutive equation is obtained finally and the strengthening effect of solid particles to composites is investigated. The work supported by the LNM, Institute of Mechanics, Chinese Academy of Sciences and by the National Natural Science Foundation of China     

A two-dimensional micromechanical model of anisotropic elastic-microplastic damage evolution

Summary A two-dimensional micromechanically based model of anisotropic elastic-microplastic damage evolution is presented. The deterioration of the material is represented by equally oriented Dugdale microcracks. Assuming a physically plausible crack growth law a consequent homogenization gives the macroscopic nonlinear stress-strain behavior during a loading process as well as the evolution of the corresponding anisotropic damage and typical characteristics during a final unloading process.

---

Ein zweidimensionales mikromechanisches Modell anisotroper elastisch-mikroplastischer Schädigungsentwicklung

Übersicht Es wird ein zweidimensionales mikromechanisch begründetes Modell anisotroper elastischmikroplastischer Schädigungsentwicklung dargestellt. Die Materialschädigung wird durch gleich-orientierte Dugdale-Mikrorisse modelliert. Unter der Annahme eines physikalisch plausiblen Rißwachstumsgesetzes liefert eine entsprechende Homogenisierung sowohl das makroskopische nichtlineare Spannungs-Dehnungs-Verhalten während eines Belastungsprozesses als auch die Evolution der entsprechenden anisotropen Schädigung sowie typische Charakteristika beim Entlastungsprozeß.

     

Fatigue damage evolution in silicon films for micromechanical applications

In this paper we examine the conditions for surface topography evolution and crack growth/fracture during the cyclic actuation of polysilicon microelectromechanical systems (MEMS) structures. The surface topography evolution that occurs during cyclic fatigue is shown to be stressassisted and may be predicted by linear perturbation analyses. The conditions for crack growth (due to pre-existing or nucleated cracks) are also examined within the framework of linear elastic fracture mechanics. Within this framework, we consider pre-existing cracks in the topical SiO₂ layer that forms on the Si substrate in the absence of passivation. The thickening of the SiO₂ that is normally observed during cyclic actuation of Si MEMS structures is shown to increase the possibility of stable crack growth by stress corrosion cracking prior to the onset of unstable crack growth in the SiO₂ and Si layers. Finally, the implications of the results are discussed for the prediction of fatigue damage in silicon MEMS structures.     

QUASI-STATIC ANALYSIS FOR VISCOELASTIC TIMOSHENKO BEAMS WITH DAMAGE

Based on convolution-type constitutive equations for linear viscoelastic materials with damage and the hypotheses of Timoshenko beams, the equations governing quasi-static and dynamical behavior of Timoshenko beams with damage were first derived. The quasi-static behavior of the viscoelastic Timoshenko beam under step loading was analyzed and the analytical solution was obtained in the Laplace transformation domain. The deflection and damage curves at different time were obtained by using the numerical inverse transform and the influences of material parameters on the quasi-static behavior of the beam were investigated in detail.     

Nonlinear vibration analysis of viscoelastic cable with small sag

Both the inplane and out-of-plane transverse vibrations of a viscoelastic cable subjected to an initial stress distributing uniform on the cross section are studied. The constitution of the cable material is assumed to be of the hereditary integral type. The partial differential-integral equations of motion are derived first. Then by applying Galerkin's method, the governing equations are reduced to a set of second-order nonlinear differential-integral equations which are solved by finite difference numerical integration procedures. Finally, the effects of the viscosity parameter and the elastic parameter on the transient amplitudes of the first mode are investigated by numerical simulation. Project supported by the National Natural Science Foundation of China (No. 59635140) and the National Postdoctoral Foundation of China.     

复合材料层合板粘弹性固化残余应力分析

基于Ｓｃｈａｐｅｒｙ积分型粘弹性本构关系，推导了考虑横向剪切效应的复合材料层合板线性热粘弹性有限元分析列式，对层合板的粘弹性响应和加工成型过程中的残余应力进行了分析，给出一些有意义的结果     

Nonlinear vibration analysis of isotropic cantilever plate with viscoelastic laminate

The nonlinear vibration of an isotropic cantilever plate with viscoelastic laminate is investigated in this article. Based on the Von Karman’s nonlinear geometry and using the methods of multiple scales and finite difference, the dimensionless nonlinear equations of motion are analyzed and solved. The solvability condition of nonlinear equations is obtained by eliminating secular terms and, finally, nonlinear natural frequencies and mode-shapes are obtained. Knowing that the linear vibration of this type of plate does not have exact solution, Ritz method is employed to obtain semi-analytical nonlinear mode-shapes of transverse vibration of this plate. Airy stress function and Galerkin method are employed to reduce nonlinear PDEs into an ODE of duffing type. Stability of plate and chaotic behavior are investigated by Runge–Kutta method. Poincare section diagrams are in good agreement with results of Lyapunov criteria.     

Nonlinear progressive damage model for composite laminates used for low-velocity impact

In order to effectively describe the progressively intralaminar and interlam- inar damage for composite laminates, a three dimensional progressive damage model for composite laminates to be used for low-velocity impact is presented. Being applied to three-dimensional （3D） solid elements and cohesive elements, the nonlinear damage model can be used to analyze the dynamic performance of composite structure and its failure be- havior. For the intralaminar damage, as a function of the energy release rate, the damage model in an exponential function can describe progressive development of the damage. For the interlaminar damage, the damage evolution is described by the framework of the continuum mechanics through cohesive elements. Coding the user subroutine VUMAT of the finite element software ABAQUS/Explicit, the model is applied to an example, i.e., carbon fiber reinforced epoxy composite laminates under low-velocity impact. It is shown that the prediction of damage and deformation agrees well with the experimental results.     

含界面脱粘压电复合材料层合板的非线性动力稳定性分析

基于Reddy提出的板高阶剪切变形简化理论,研究了含界面脱粘损伤压电复合材料层合板非线性动力稳定性问题.首先,建立了分层模型,推导了考虑几何非线性、阻尼效应、纵向惯性力和力-电耦合效应的Mathieu方程,并且给出了该方程解的解析表达式.其次,通过典型算例讨论了界面脱粘损伤以及反馈控制力对该层合板动力不稳定区域、纵向、横向共振频率和最大"牵引"深度的影响.由典型算例讨论可知:随着层合板界面脱粘损伤的扩大,其动力稳定性能逐渐减弱,其中在损伤较小时,反馈控制力对智能结构几乎没有影响;而在损伤比较大的情况下,反馈控制力将能有效地减少动力不稳定区域重合面积.     

Nonlinear viscoelastic properties and change in entanglement structure of linear polymer

Simultaneous measurements of stress relaxation and differential dynamic modulus were made at 268 K over a time scale of 10 to 10⁴⁵ s for nearly monodisperse polybutadiene (M _w =2.2x10⁵, 1,2-structure 70%, M _e =3600) and also one having coarse cross-linking (M _c =29000). Static shear strain ranged from 0.1 to 2.0. In a long-time region (t> _k ), the relaxation modulus G (; t) could be expressed by the product G (0; t) h (y). The observed h() agreed well with the Doi-Edwards theory without use of IA approximation. Both the cured and uncured samples showed initial drop of the differential storage modulus G (), ; t) followed by gradual recovery, but did not attain the value before shearing G (, ; t) for the uncured sample showed smaller values than that for the cured one in the whole measured time scale at the higher strain, confirming the two origins of nonlinear viscoelasticity of well entangled polymer; induced chain anisotropy and induced decrement in entanglement density. G (, ; t) curves for the cured sample agreed well with the BKZ predictions. But the curves for the uncured sample agreed well with the BKZ prediction only at the time scale of t< _k . BKZ prediction showed significant upward deviations at t> _k . Such the differences are discussed in terms of the two origins.Dedicated to Prof. John D. Ferry on the occasion of his 85th birthday.     

Effects of inter-fibre spacing on damage evolution in unidirectional (UD) fibre-reinforced composites

A three dimensional (3D) micromechanical study has been performed in order to investigate local damage in UD composite materials under transverse and longitudinal tensile loading. In particular, the influence of non-uniform distribution of fibres in RVEs (representative volume element) with a hexagonal packing array and the effects of thermal residual stresses has been investigated. To examine the effect of inter-fibre spacing and residual stress on failure, a study based on the Maximum Principal Stress failure criterion and a stiffness degradation technique has been used for damage analysis of the unit cell subjected to mechanical loading. Results indicate a strong dependence of damage onset and its evolution from the fibres position within the RVE. Predicted mechanical properties, damage initiation and evolution are also clearly influenced by the presence of residual stress.     

Load transfer in fibre-reinforced composites with viscoelastic matrix: an analytical study

In A fibre-reinforced 2D composite material with elastic fibres and viscoelastic, isotropic matrix is studied. Starting from the solution of a reference-problem with elastic matrix material the elastic matrix parameters are substituted by their viscoelastic correspondents in the Laplace domain. For simplification the time-dependent solution is approximated by using limiting value theorems that give information about the time-dependent solution for t → 0 and t → ∞. Then the method of asymptotically equivalent functions is used and illustrated with examples of a steel fibre in a PMMA matrix. The analytical solutions are compared with their numerical counterparts. In summary it can be stated that this paper is a further contribution to the vast literature about the application of the correspondence principle to the solution of special problems of the linear viscoelasticity.     

MICROMECHANICS ANALYSIS ON EVOLUTION OF CRACK IN VISCOELASTIC MATERIALS

ＩｎｔｒｏｄｕｃｔｉｏｎＡｓｓｈｏｗｎｂｙＺＨＡＮＧａｎｄＸＩＯＮＧ ( 1 997) [1],ｕｎｄｅｒｃｏｎｓｔａｎｔｌｏａｄｉｎｇｔｈｅｃｒａｃｋｉｎｖｉｓｃｏｅｌａｓｔｉｃｍａｔｅｒｉａｌｓｈａｓａｎｉｎｃｕｂａｔｉｏｎｔｉｍｅ ,ｗｉｔｈｉｎｉｔｔｈｅｃｒａｃｋｏｐｅｎｉｎｇｄｉｓｐｌａｃｅｍｅｎｔ (ＣＯＤ)ｉｓｃｈａｎｇｉｎｇｌａｒｇｅｒｗｉｔｈｔｉｍｅ,ｗｈｉｌｅｔｈｅｃｒａｃｋｌｅｎｇｔｈｋｅｅｐｓｃｏｎｓｔａｎｔ.Ｉｎｖｅｓｔｉｇａｔｉｏｎｏｎｔｈｅｇｒａｄｕａｌｏｐｅｎｉｎｇｏｆｔｈｅｃｒａｃ…     

Nonlinear dynamic analysis of electrically actuated viscoelastic bistable microbeam system

Nonlinear Dynamics - Viscoelastic phenomena widely exist in MEMS materials, which may have certain effects on transition mechanism of nonlinear jumping phenomena and transient chaotic behaviors....     

Spallation analysis with a closed trans-scale formulation of damage evolution

A closed, trans-scale formulation of damage evolution based on the statistical microdamage mechanics is summarized in this paper. The dynamic function of damage bridges the mesoscopic and macroscopic evolution of damage. The spallation in an aluminium plate is studied with this formulation. It is found that the damage evolution is governed by several dimensionless parameters, i.e., imposed Deborah numbersDe ^* andDe, Mach numberM and damage numberS. In particular, the most critical mode of the macroscopic damage evolution, i.e., the damage localization, is determined by Deborah numberDe ^*. Deborah numberDe ^* reflects the coupling and competition between the macroscopic loading and the microdamage growth. Therefore, our results reveal the multi-scale nature of spallation. In fact, the damage localization results from the nonlinearity of the microdamage growth. In addition, the dependence of the damage rate on imposed Deborah numbersDe ^* andDe, Mach numberM and damage numberS is discussed. The project supported by the National Natural Science Foundation of China (10172084, 10232040, 10232050, 10372012, 10302029) and the Special Funds for Major State Research Project (G200077305)     

Nonlinear analysis of the surface tension driven breakup of viscoelastic filaments

The surface tension driven breakup of viscoelastic filaments into droplets is qualitatively different from that of Newtonian liquid filaments. Disturbances on filaments of dilute polymer solutions often grow to a configuration consisting of nascent droplets connected by a thin ligament; the breakup time for this configuration is much longer than that predicted by extensions of Rayleigh's linear stability theory. We present here a nonlinear analysis of surface tension driven breakup of viscoelastic filaments using two complementary approaches that given equivalent results: a transient finite element solution and a one-dimensional thin filament approximation. We show that significant nonlinear effects lead to the experimentally-observed nascent droplet-ligament configuration, and we predict the entire evolution of the filament profile. Agreement with available experimental data for profile evolution and breakup of jets of Newtonian fluids and dilute polymer solutions is excellent.