An analysis of axisymmetric buckling and growth of circular-shaped delamination

The buckling and post-buckling of clamped circular plate subjected to distributed radial compressed load is presented by using the high-order perturbation analysis and shooting method. The sixth-order solution shows good agreement with the FEM results in [11]. The results in this paper are applied to investigate the buckling and growth of pressed thin film delamination in the film/substrate system. Under a certain residual pressure in the thin film, two characteristic blister radiiR _c andR _g, the critical radius and growing radius respectively, are obtained. The numerical result shows that the growth criterion of delamination in [9,10] is not perfect. In variant residual stress or interface toughness, the conditions of no growth, stable growth and unstable growth of the delamination are obtained by comparing the driving force at the interface crack tip with the interface toughness.This project is supported by National Natural Science Foundation of China.     

圆形脱层的轴对称屈曲及扩展分析

利用高阶摄动结合打靶法，分析了固支圆板在均匀径向力作用下的轴对称屈曲和过屈曲，所得结果与文（６）的ＦＥＭ结构吻合得很好，应用于薄膜－基底结构，研究了受压薄膜脱层的屈曲、扩展问题，得到了在一定的残余压应力作用下，脱层屈曲的临界尺寸Ｒｃ和扩展尺寸Ｒｇ。     

The effects of bridging in a 3D composite on buckling, postbuckling and growth of delamination

Summary Buckling and postbuckling solutions to circular delamination constrained by transversal restoring forces, which occur extensively in stitched or woven composites with three-dimensional (3D) reinforcement, are obtained by using von Karman's geometrically nonlinear thin plate theory by means of Taylor's series expansion. The through-thickness tows are assumed to provide continuous and linear restoring tractions, opposing the deflection of the annular delaminated region adjacent to a penny-shaped crack. When the end of the delaminated layer is clamped, and the deflection is permitted in the positive direction of the z-axis only, there exists a characteristic delamination radius a _* for initial buckling. In the case that the initial delamination radius a ₀ exceeds a _*, it will consist of waves whose sizes decrease gradually, as they are apart from the delamination center with larger distances, and will usually not span the whole crack region. Therefore, buckling profiles can be divided into two types: (1) lacking contact phenomena between the delaminated layer and the base plate; (2) having contact surfaces inside the delamination region. In this paper, growth laws of buckling, postbuckling and growth of delamination at lacking contact surface are discussed. The corresponding stability of the delamination growth under fixed boundary load is studied, and the dependence of stable scope upon the fracture toughness of the composite and the elastic constant of bridging fiber is summarized. It follows from the analysis that bridging can increase the load-bearing capacity of composite structure, improve its mechanical performances and restrain the growth of delamination. Received 23 November 1998; accepted for publication on 13 January 1999     

Creep buckling and postbuckling of circular cylindrical shells under axial compression

An infinitely long, axially compressed, circular cylindrical shell with an imperfection in the shape of the axisymmetric classical buckling mode, undergoing steady or non-steady creep, is analyzed. The axisymmetric problem is solved incrementally using nonlinear shell equations The ratio of the applied stress to the classical buckling stress determines if the shell will collapse axisymmetrically or if it will bifurcate into a nonaxisymmetric mode, and whether or not bifurcation will result in instantaneous collapse. The bifurcation problem is formulated exactly and the initial postbuckling behavior is investigated via an asymptotic elastic analysis, based on Koiter's general theory Numerical results are compared with available experimental data.     

Dynamics of delamination buckling

The dynamics of delamination growth in layered plates as a result of delamination buckling is studied in the context of the model problem of a two layer circular plate with an initial penny-shaped delaminalion subject to a time dependent radial compressive load. The interaction of delamination buckling and normal impact is modeled through the application of a uniform transverse tension pulse, simulating the effect of a wave front passing through the layer interface, in conjunction with the radial loading. In addition three types of radial loading functions are considered independently and characcteristic features of the response are observed. Comparison is made with growth predictions made using a quasi-static analysis.     

Optimizing strength and toughness of nanofiber-reinforced CMCs

Nanofibers used in current ceramic matrix composites (CMCs) are usually wavy and of finite length. Here, a shear-lag model for predicting the tensile strength and work of fracture of a composite containing a “single matrix crack”, as a function of all the material parameters including constant plus Coulomb interfacial friction, is presented for a CMC containing wavy, finite-length nanofibers having a statistical distribution of strengths. Literature results are recovered for straight infinite fibers, with strength and toughness depending on a characteristic strength σ_c and a characteristic length δ_c. For nanofibers of finite length L, radius of curvature R, and with Coulomb friction coefficient μ, the strength and toughness are found to depend only on σ_c, δ_c, and two new dimensionless parameters μδ_c/R and L/δ_c. Parametric results for a typical carbon nanotube CMC show an optimal region of morphology (L and R) that maximizes composite strength and toughness, exceeding the properties of composites with straight (R=∞) and/or long (L=∞) reinforcements. Therefore, while factors such as the nanofiber strength distribution and the nanofiber–matrix interface sliding resistance may not be easily controlled, the tuning, via processing, of purely geometrical properties of the nanofibers (L and R) alone can aid in maximizing composite properties. These results have important application in hybrid CMCs such as “fuzzy fiber” CMCs, where micron-scale fibers are covered with a forest of nanofibers such that the nanofiber array can bridge longitudinal cracks and thus improve delamination resistance.     

Post-buckled propagation model for compressive fatigue of impact damaged laminates

An analytical model for prediction of compressive fatigue threshold strains in composite plates with barely visible impact damage (BVID) is presented. The model represents the complex damage morphology as a single circular delamination at a critical level and calculates the strain at which thin-film buckling of the circular delaminated region occurs. The threshold strain is defined as the strain at which the strain energy release rate for the fracture of post-buckled delaminated plies along the delamination is equal to the critical Mode I value (G_1C) for the resin. The model predicts the critical through-thickness level for delamination, the stability of delamination growth and also the sensitivity to experimental error in geometric measurements of the damage area. Results obtained using the model show an agreement of fatigue strain to within 4% of experimental values for four sets of data reported in the literature.     

A review of the volume-based strain energy density approach applied to V-notches and welded structures

A large bulk of experimental data from static tests of sharp and blunt V-notches and from fatigue tests of welded joints are presented in an unified way by using the mean value of the Strain Energy Density (SED) over a given finite-size volume surrounding the highly stressed regions. When the notch is blunt, the control area assumes a crescent shape and R₀ is its width as measured along the notch bisector line. In plane problems, when cracks or pointed V-notches are considered, the volume becomes a circle or a circular sector, respectively. The radius R₀ depends on material fracture toughness, ultimate tensile strength and Poisson’s ratio in the case of static loads; it depends on the fatigue strength Δσ_A of the butt ground welded joints and the Notch Stress Intensity Factor (NSIF) range ΔK₁ in the case of welded joints under high cycle fatigue loading (with Δσ_A and ΔK₁ valid for 5 × 10⁶ cycles).Dealing with welded joints characterised by a plate thickness greater than 6 mm, the final synthesis based on SED summarises nine hundred data taken from the literature while a new synthesis from spot-welded joints under tension and shear loading, characterised by a limited thickness of the main plate, is presented here for the first time (more than two hundred data).Dealing with static tests, about one thousand experimental data as taken from the recent literature are involved in the synthesis. The strong variability of the non-dimensional radius R/R₀, ranging from about zero to about 1000, makes the check of the approach based on the mean value of the SED severe.     

Stability of an annular plate of a shape memory alloy

Analytical solutions using various hypotheses are obtained for the problem of buckling of an annular circular plate of a shape memory alloy in the forward or reverse thermoelastic phase transformations under uniform radial compression. It is established that, despite the axial symmetry of the body geometry, loads, and boundary conditions, the minimum critical loads for annular plates correspond to nonaxisymmetric buckling modes, in contrast to continuous plates.     

Three-dimensional numerical simulation of detonations in coaxial tubes

Three-dimensional numerical simulation of detonations in both a circular tube and a coaxial tube are simulated to reveal characteristics of single spinning and two-headed detonations. The numerical results show a feature of a single spinning detonation which was discovered in 1926. Transverse detonations are observed in both tubes, however, the single spinning mode maintains the complex Mach reflection whereas the two-headed mode develops periodically from the single Mach reflection to the complex one. The calculated cell aspect ratio for the two-headed mode changes from 1.09 to 1.34 as the radius of axial insert increases from r ₁/R = 0.1 to 0.9. The calculated cell aspect ratio for r ₁/R = 0.1 is close to the experimental results without an axial insert. The formation of an unreacted gas pocket behind the detonation front was not observed in the single spinning mode; however, the two-headed mode has unreacted gas pocket behind the front near the axial insert.      

Two-dimensional analyses of delamination buckling of symmetrically cross-ply rectangular laminates

The conventional approach to analysis the buckling of rectangular laminates containing an embedded delamination subjected to the in-plane loading is to simplify the laminate as a beam-plate from which the predicted buckling load decreases as the length of the laminate increases. Two-dimensional analyses are employed in this paper by extending the one-dimensional model to take into consideration of the influence of the delamination width on the buckling performance of the laminates. The laminate is simply supported containing a through width delamination. A new parameter β defined as the ratio of delamination length to delamination width is introduced with an emphasis on the influence of the delamination size. It is found that (i) when the ratio β is greater than one snap-through buckling prevails, the buckling load is determined by the delamination size and depth only; (ii) as the ratio β continues to increase, the buckling load will approach to a constant value. Solutions are verified with the well established results and are found in good agreement with the latter.     

Nonlinear bending and post-buckling of a functionally graded circular plate under mechanical and thermal loadings

Based on the classical nonlinear von Karman plate theory, axisymmetric large deflection bending of a functionally graded circular plate is investigated under mechanical, thermal and combined thermal–mechanical loadings, respectively, and axisymmetric thermal post-buckling behavior of a functionally graded circular plate is also investigated. The mechanical and thermal properties of functionally graded material (FGM) are assumed to vary continuously through the thickness of the plate, and obey a simple power law of the volume fraction of the constituents. Governing equations for the problem are derived, and then a shooting method is employed to numerically solve the equations. Effects of material constant n and boundary conditions on the temperature distribution, nonlinear bending, critical buckling temperature and thermal post-buckling behavior of the FGM plate are discussed in details.     

ANALYSIS ON BUCKLING AND POSTBUCKLING OF DELAMINATION IN 3D COMPOSITES

In this paper, the problem of axisymmetric buckling and postbuckling of a circular thin-film delamination bridged by through-thickness fiber tows in 3D composites is presented. An iterativeprocedure based on Taylor's series expansion is used to generate a family of nondimensionalized post-buckling solutions of the above problem by yon Karman's nonlinear plate theory. Attention is fo-cused, herein, on the effects of the bridge force of through-thickness fibers on the buckling and post-buckling behavior of the delamination. It is found that fiber bridge not only increases the ability of re-sisting delamination buckling and postbuckling, but also brings on the jump of the delamination deflec-tion mode during the postbuckling phase. Consequently the behavior of the composite structure with de-lamination is greatly improved, such as increasing the residual strength and prolonging the service life.     

Mode III fracture behavior of laminated composite with edge crack in torsion

A three-dimensional (3-D) finite element analysis was performed on a [90,(+45/−45)_n,(−45/+45)_n,90]_s class of laminated composites under the edge crack torsion (ECT) test configuration. Finite element delamination models were established and formulas for calculating the Mode III fracture toughness from 3-D finite element models were developed. The relations between the interlaminar fracture behavior and various configuration parameters were investigated and the effects of point loads, ends, geometry, Mode II interference, and friction were evaluated. Results showed that with proper selection of ECT specimen configuration and layup, the delamination could grow in pure Mode III in the middle region of the specimen. Specimen end effect played an important role in the ECT test. A Mode II component occurred in the end regions but it did not interfere significantly with the Mode III delamination state. Specimen dimension ratio, layup, and crack length exhibited significant effect on the interlaminar fracture behavior and the calculated strain energy release rates. However, friction between crackfaces was found to have negligible effect on the interlaminar properties.     

The characterization of telephone cord buckling of compressed thin films on substrates

The topology of the telephone cord buckling of compressed diamond-like carbon films (DLC) on glass substrates has been characterized with atomic force microscopy (AFM) and with the focused ion beam (FIB) imaging system. The profiles of the several buckles have been measured by AFM to establish the symmetry of each repeat unit, revealing similarity with a circular buckle pinned at its center. By making parallel cuts through the buckle in small, defined locations, straight-sided buckles have been created on the identical films, enabling the residual stress in the film to be determined from the profile.It has been shown that the telephone cord topology can be effectively modeled as a series of pinned circular buckles along its length, with an unpinned circular buckle at its front. The unit segment comprises a section of a full circular buckle, pinned to the substrate at its center. The model is validated by comparing radial profiles measured for the telephone cord with those calculated for the pinned buckle, upon using the residual stress in the film, determined as above. Once validated, the model has been used to determine the energy release rate and mode mixity, G(ψ).The results for G(ψ) indicate that the telephone cord configuration is preferred when the residual stress in the DLC is large, consistent with observations that straight-sided buckles are rarely observed, and, when they occur, are generally narrower than telephone cords. Telephone cords are observed in many systems, and can be regarded as the generic morphology. Nevertheless, they exist subject to a limited set of conditions, residing within the margin between complete adherence and complete delamination, provided that the interface has a mode II toughness low enough to ensure that the buckle crack does not kink into the substrate.     

受弯脱层层板后屈曲有限元分析

脱层是复合材料层板结构中主要的缺陷形式之一。当脱层层板受到压力载荷的作用会造成脱层的局部屈曲和扩展,从而使结构的强度和刚度大为降低。含脱层层板的弯曲问题包含了脱层的压缩问题,却比压缩问题更加复杂。本文对含穿透脱层层板在纯弯载荷作用下的后屈曲问题进行了基于一阶剪切层板理论的几何非线性有限元分析,运用虚裂纹闭合技术求解了纯弯载荷作用下的脱层尖端的能量释放率各型分量,并用脱层扩展判据求解了脱层起始扩展载荷。     

Exact solutions for radial deformations of a functionally graded isotropic and incompressible second-order elastic cylinder

We find closed-form solutions for axisymmetric plane strain deformations of a functionally graded circular cylinder comprised of an isotropic and incompressible second-order elastic material with moduli varying only in the radial direction. Cylinder's inner and outer surfaces are loaded by hydrostatic pressures. These solutions are specialized to cases where only one of the two surfaces is loaded. It is found that for a linear through-the-thickness variation of the elastic moduli, the hoop stress for the first-order solution (or in a cylinder comprised of a linear elastic material) is a constant but that for the second-order solution varies through the thickness. The radial displacement, the radial stress and the hoop stress do not depend upon the second-order elastic constant but the hydrostatic pressure and hence the axial stress depends upon it. When the two elastic moduli vary as the radius raised to the power two or four, the radial and the hoop stresses in an infinite space with a pressurized cylindrical cavity equal the pressure in the cavity. For an affine variation of the elastic moduli, the hoop stress in an internally loaded cylinder made of a linear elastic isotropic and incompressible material at the point is the same as that in a homogeneous cylinder. Here R_in and R_ou equal, respectively, the inner and the outer radius of the undeformed cylinder and R the radial coordinate of a point in the unstressed reference configuration.     

Bifurcation buckling of shells of revolution including large deflections,plasticity and creep

A summary is first presented of the conceptual difficulties and paradoxes surrounding plastic bifurcation buckling analysis. Briefly discussed are nonconservativeness, loading rate during buckling, and the discrepancy of buckling predictions with use of J₂ flow theory vs J₂ deformation theory. The axisymmetric prebuckling analysis, including large deflections, elastic-plastic material behavior and creep is summarized. Details are given on the analysis of nonsymmetric bifurcation from the deformed axisymmetric state. Both J₂ flow theory and J₂ deformation theory are described. The treatment, based on the finite-difference energy method, applies to layered segmented and branched shells of arbitrary meridional shape composed of a number of different elastic-plastic materials. Numerical results generated with a computer program based on the analysis are presented for an externally pressurized cylinder with conical heads.     

Analysis on small scale bridging in fibrous monolithic ceramics

Based on shear lag model of interface between fiber and matrix, a new formula that relates the crack opening displacement and bridging force in fibrous monolithic ceramics was constructed under the framework of small scale bridging. This formula was applied to predict the fracture resistance orR-curve response of a three-point bending specimen made of fibrous monolithic ceramics. A parametric investigation on the influences of fiber volume fraction, fiber radius, characteristics of constituents, BN's fracture toughness and specimen's geometry on the bridging forces and fracture resistance in Si₃N₄/BN composite was carried out. The upper and lower limits of theR-curve of Si₃N₄/BN in small scale bridging were derived. This research revealed the role played by the above parameters in the fracture toughness of fibrous monolithic ceramics. Supported by National Natural Science Foundation of China (59632090).     

Deformations and stability of spherical caps under centrally distributed pressures

In this paper the deformations and stability in large axisymmetric deflection of spherical caps under centrally distributed pressures are investigated. The Newton- spline method for solving the nonlinear equations governing large axisymmetric deflection of spherical caps is presented. The buckling behavior is studied for a cap with fixed geometry when the size of the loaded radius is allowed to vary, and for a fixed loaded radius when the shell geometry is allowed to vary. The influence of the buckling modes on the critical loads is analysed. Numerical results are given for =0.3.Project Supported by National Science Foundation of China.