Bonded lap joints of composite laminates with tapered edges

This study presents a semi-analytical solution method to analyze the geometrically nonlinear response of bonded composite lap joints with tapered and/or non tapered adherend edges under uniaxial tension. The solution method provides the transverse shear and normal stresses in the adhesives and in-plane stress resultants and bending moments in the adherends. The method utilizes the principle of virtual work in conjunction with von Karman’s nonlinear plate theory to model the adherends and the shear lag model to represent the kinematics of the thin adhesive layers between the adherends. Furthermore, the method accounts for the bilinear elastic material behavior of the adhesive while maintaining a linear stress–strain relationship in the adherends. In order to account for the stiffness changes due to thickness variation of the adherends along the tapered edges, the in-plane and bending stiffness matrices of the adherents are varied as a function of thickness along the tapered region. The combination of these complexities results in a system of nonlinear governing equilibrium equations. This approach represents a computationally efficient alternative to finite element method. The numerical results present the effects of taper angle, adherend overlap length, and the bilinear adhesive material on the stress fields in the adherends, as well as the adhesives of a single- and double-lap joint.     

Fully-coupled nonlinear analysis of single lap adhesive joints

This paper presents novel closed-form and accurate solutions for the edge moment factor and adhesive stresses for single lap adhesive bonded joints. In the present analysis of single lap joints, both large deflections of adherends and adhesive shear and peel strains are taken into account in the formulation of two sets of nonlinear governing equations for both longitudinal and transverse deflections of adherends. Closed-form solutions for the edge moment factor and the adhesive stresses are obtained by solving the two sets of fully-coupled nonlinear governing equations. Simplified and accurate formula for the edge moment factor is also derived via an approximation process. A comprehensive numerical validation was conducted by comparing the present solutions and those developed by Goland and Reissner, Hart-Smith and Oplinger with the results of nonlinear finite element analyses. Numerical results demonstrate that the present solutions for the edge moment factor (including the simplified formula) and the adhesive stresses appear to be the best as they agree extremely well with the finite element analysis results for all ranges of material and geometrical parameters.     

On the estimation of Kt at a bottom of a grind-out cavity after one-sided patching

An unsupported plate containing an embedded grind-out cavity repaired with a reinforcement bonded on one side may experience a considerable out-of-plane bending near a cavity due to the load-path eccentricity even when the geometrically nonlinear effect is taken into account. This out-of-plane bending causes stresses in the plate near the bottom of the cavity vary significantly through a remaining plate’s thickness with inner wall fiber stresses in some cases at 20–30% higher than the corresponding plane-stress type results, i.e., those without considering the out-of-plane deflections. A plane-stress type analysis of a repair over a corrosion cavity had been presented in a previous paper by Duong et al. [Theoretical and Applied Fracture Mechanics 36 (2001a) 187] for a constant depth cavity repaired with an elliptical patch and most recently by Duong and Yu [International Journal of Engineering Science 40 (2002a) 347] for a spherical depth cavity repaired with a polygonal patch. Extension of these methods to include the effect of out-of-plane bending therefore will be delineated in the present paper.     

Stress distribution in a lap joint under tension-shear

This paper deals with the elastostatic load transfer of a tensile load in a model of an adhesive lap joint (tension-shear problem). The adhesive layer is regarded as infinitesimally thin and the displacement and traction vectors in the adherends are assumed to be continuous across the bond. The problem is reduced to a pair of Fredholm integral equations of the second kind which involve the mean angle between the deformed bond line and the tensile load. This angle, in turn, is determined by means of a scheme due to Goland and Reissner. Numerical results for the bond line stresses and the stress intensity factors at the ends of the bonded region are presented.     

EXPERIMENTAL AND ANALYTICAL STUDY OF GEOMETRIC NONLINEAR BENDING OF A CANTILEVER BEAM UNDER A TRANSVERSE LOAD

Journal of Applied Mechanics and Technical Physics - Exact and approximate analytical solutions are compared with experimental data on the geometrically nonlinear bending of a...     

Analytical solutions for adhesive composite joints considering large deflection and transverse shear deformation in adherends

This paper presents a novel formulation and analytical solutions for adhesively bonded composite single lap joints by taking into account the transverse shear deformation and large deflection in adherends. On the basis of geometrically nonlinear analysis for infinitesimal elements of adherends and adhesive, the equilibrium equations of adherends are formulated. By using the Timoshenko beam theory, the governing differential equations are expressed in terms of the adherend displacements and then analytically solved for the force boundary conditions prescribed at both overlap ends. The obtained solutions are applied to single lap joints, whose adherends can be isotropic adherends or composite laminates with symmetrical lay-ups. A new formula for adhesive peel stress is obtained, and it can accurately predict peel stress in the bondline. The closed-form analytical solutions are then simplified for the purpose of practical applications, and a new simple expression for the edge moment factor is developed. The numerical results predicted by the present full and simplified solutions are compared with those calculated by geometrically nonlinear finite element analysis using MSC/NASTRAN. The agreement noted validates the present novel formulation and solutions for adhesively bonded composite joints. The simplified shear and peel stresses at the overlap ends are used to derive energy release rates. The present predictions for the failure load of single lap joints are compared with those available in the literature.     

Propagation of Bending Waves in a Beam the Material of Which Accumulates Damage During Its Operation

Journal of Applied Mechanics and Technical Physics - A self-consistent mathematical model is stated in linear and nonlinear formulations, which includes the equation of bending vibrations of a beam...     

Solution for bending of a polygonal plate with cracks by the linear conjugate method

Archive of Applied Mechanics - Stress–strain state of a polygonal multiply connected plate under bending is considered. The plate under the action of different force factors (bending moment M...     

A crack bridging model for bonded plates subjected to tension and bending

A crack bridging model is presented for analysing the tensile stretching and bending of a cracked plate with a patch bonded on one side, accounting for the effect of out-of-plane bending induced by load-path eccentricity inherent to one-sided repairs. The model is formulated using both Kirchhoff–Poisson plate bending theory and Reissners shear deformation theory, within the frameworks of geometrically linear and nonlinear elasticity. The bonded patch is represented as distributed springs bridging the crack faces. The springs have both tension and bending resistances ; their stiffness constants are determined from a one-dimensional analysis for a single strap joint, representative of the load transfer from the cracked plate to the bonded patch. The resulting coupled integral equations are solved using a Galerkin method, and the results are compared with three-dimensional finite element solutions. It is found that the formulation based on Reissners plate theory provides better agreement with finite element results than the classical plate theory.     

Simulation of Rheological Properties of Polyethylene Melts under Uniaxial Tension

Journal of Applied Mechanics and Technical Physics - The modified Vinogradov–Pokrovskii rheological model is used to describe the occurrence of stresses in a polymer melt under uniaxial...     

A Bending-Gradient model for thick plates,Part II: Closed-form solutions for cylindrical bending of laminates

In the first part (Lebée and Sab, 2010a) of this two-part paper we have presented a new plate theory for out-of-plane loaded thick plates where the static unknowns are those of the Kirchhoff–Love theory (3 in-plane stresses and 3 bending moments), to which six components are added representing the gradient of the bending moment. The new theory, called Bending-Gradient plate theory is an extension to arbitrarily layered plates of the Reissner–Mindlin plate theory which appears as a special case when the plate is homogeneous. Moreover, we demonstrated that, in the general case, the Bending-Gradient model cannot be reduced to a Reissner–Mindlin model. In this paper, the Bending-Gradient theory is applied to laminated plates and its predictions are compared to those of Reissner–Mindlin theory and to full 3D (Pagano, 1969) exact solutions. The main conclusion is that the Bending-Gradient gives good predictions of deflection, shear stress distributions and in-plane displacement distributions in any material configuration. Moreover, under some symmetry conditions, the Bending-Gradient model coincides with the second-order approximation of the exact solution as the slenderness ratio L/h goes to infinity.     

WENTZEL–KRAMERS–BRILLOUIN SOLUTIONS TO AN EQUATION OF INTERNAL GRAVITY WAVES IN A STRATIFIED MEDIUM WITH SLOWLY VARYING SHEAR FLOWS

Journal of Applied Mechanics and Technical Physics - Model buoyancy frequency distribution and the Wentzel–Kramers–Brillouin method are used to obtain an asymptotic solution to a...     

A Bending-Gradient model for thick plates. Part I: Theory

This is the first part of a two-part paper dedicated to a new plate theory for out-of-plane loaded thick plates where the static unknowns are those of the Kirchhoff–Love theory (3 in-plane stresses and 3 bending moments), to which six components are added representing the gradient of the bending moment. The new theory, called the Bending-Gradient plate theory is described in the present paper. It is an extension to arbitrarily layered plates of the Reissner–Mindlin plate theory which appears as a special case of the Bending-Gradient plate theory when the plate is homogeneous. However, we demonstrate also that, in the general case, the Bending-Gradient model cannot be reduced to a Reissner–Mindlin model. In part two (Lebée and Sab, 2011), the Bending-Gradient theory is applied to multilayered plates and its predictions are compared to those of the Reissner–Mindlin theory and to full 3D Pagano’s exact solutions. The main conclusion of the second part is that the Bending-Gradient gives good predictions of both deflection and shear stress distributions in any material configuration. Moreover, under some symmetry conditions, the Bending-Gradient model coincides with the second-order approximation of the exact solution as the slenderness ratio L/h goes to infinity.     

Linear and higher order displacement theories for adhesively bonded lap joints

This work presents an adhesive model for stress analysis of bonded lap joints, which can be applied to model thin and thick adhesive layers. In this theory, linear variations of displacement components along the adhesive thickness are firstly assumed, and the longitudinal strain and the Poisson's effect of the adhesive are modeled. A differential form of the equilibrium equations for the adherends is analytically solved by means of compatible relations of the adhesive deformation. The derived shear and peel stresses are compared with the classical adhesive model of continuous springs with constant shear and peel stresses, and validated with two-dimensional finite element results of the geometrically nonlinear analysis using a commercial package. The numerical results show that the present linear displacement theory can be applied to both thin and moderately thick adhesive layers. The present formulation of the linear displacement theory is then extended to the higher order displacement theory for stress analysis of a thick adhesive, whose numerical results are also compared with those of the finite element computation.     

Stress-function variational method for interfacial stress analysis of adhesively bonded joints

High interfacial stresses at the free edges of adherends are responsible for the debonding failure of adhesively bonded joints (ABJs). In this paper, a general stress-function variational method is formulated to determinate the interfacial shear and normal (peeling) stresses in ABJs in high accuracy. By extending authors’ prior work in stress analysis of bonded joints (Wu and Jenson, 2011), all the planar stress components in the adherends and adhesive layer of an ABJ are expressed in terms of four unknown interfacial stress functions, which are introduced at the upper and lower surfaces of the adhesive layer. A set of governing ordinary differential equations (ODEs) of the four interfacial stress functions is obtained via minimizing the complimentary strain energy of the ABJ, which is further solved by using eigenfunctions. The obtained semi-analytic stress field can satisfy all the traction boundary conditions (BCs) of the ABJ, especially the stress continuity across the bonding lines and the shear-free condition at the ends of adherends and adhesive layer. As an example, the stress field in an adhesively single-sided strap joint is determined by the present method, whose numerical accuracy and reliability are validated by finite element method (FEM) and compared to existing models in the literature. Parameter studies are performed to examine the dependencies of the interfacial stresses of the exemplified ABJ upon the geometries, moduli and temperature change of the adherends and adhesive layer, respectively. The present method is applicable for scaling analysis of joint strength, optimal design of ABJs, etc.     

Nonlinear analysis of a composite panel with a cutout repaired by a bonded tapered composite patch

This study presents a solution method to analyze the geometrically nonlinear response of a patch-repaired flat panel (skin) with a cutout under general loading conditions. The effect of induced stiffening due to tensile loading on the in-plane and, particularly, the out-of-plane behaviors of the patch-repaired skin are investigated. The damage to the skin is represented in the form of a cutout under the patch. The patch with tapered edges is free of external tractions. The skin is subjected to general boundary and loading conditions along its external edge. The solution method provides the transverse shear and normal stresses in the adhesive between the skin and the patch, and in-plane and bending stresses in the patch and skin. Both the patch and skin are made of linearly elastic composite laminates, and the adhesive between them is homogeneous and isotropic, exhibiting a bi-linear elastic behavior. Modified Green’s strain–displacement relations in conjunction with von Karman assumptions are employed in determining the in-plane strains in the skin and patch; however, the transverse shear strains in the adhesive are determined based on the shear-lag theory. The present solution method utilizes the principle of virtual work in conjunction with complex potential functions.     

BENDING OF THIN ELECTROMAGNETOELASTIC PLATES

Journal of Applied Mechanics and Technical Physics - Main correlations in the theory of bending of thin electromagnetoelastic plates are obtained, in which complex potentials are used. Exact...     

Exact static solutions to piezoelectric smart beams including peel stresses. II. Numerical results,comparison and discussion

This part presents the numerical results, comparisons and discussion for the exact static solutions of smart beams with piezoelectric (PZT) actuators and sensors including peel stresses presented in Part I. (International Journal of Solids and Structures, 39, 4677–4695) The actuated stress distributions in the adhesive and the adhesive edge stresses varying with the thickness ratios are firstly obtained and presented. The actuated internal stress resultants and displacements in the host beam are then calculated and compared with those predicted by using the shear lag model. The stresses in the adhesive caused by an applied axial force, bending moment and shear force are calculated, and then used to compute the sensing electric charges for comparison with those predicted using the shear lag model. The numerical results are given for the smart beam with (a) one bonded PZT and (b) two symmetrically bonded PZTs, with a comparison to those predicted using the shear lag model. Novel, simple and more accurate formulas for the equivalent force and bending moment induced by applied electric field are also derived for the host beam with one PZT or two symmetrically bonded PZTs. The symmetric shear stress and the anti-symmetric peel stress components caused by a shear force are discussed. In addition, in the case of PZT edge debonding, the stress redistribution in the adhesive and the self-arresting mechanism are also investigated.     

On internal dissipation inequalities and finite strain inelastic constitutive laws: Theoretical and numerical comparisons

Internal dissipation always occurs in irreversible inelastic deformation processes of materials. The internal dissipation inequalities (specific mathematical forms of the second law of thermodynamics) determine the evolution direction of inelastic processes. Based on different internal dissipation inequalities several finite strain inelastic constitutive laws have been formulated for instance by Simo [Simo, J.C., 1992. Algorithms for static and dynamic multiplicative plasticity that preserve the classical return mapping schemes of the infinitesimal theory. Computer Methods in Applied Mechanics and Engineering 99, 61–112]; Simo and Miehe [Simo, J.C., Miehe, C., 1992. Associative coupled thermoplasticity at finite strains: formulation, numerical analysis and implementation. Computer Methods in Applied Mechanics and Engineering 98, 41–104]; Lion [Lion, A., 1997. A physically based method to represent the thermo-mechanical behavior of elastomers. Acta Mechanica 123, 1–25]; Reese and Govindjee [Reese, S., Govindjee, S., 1998. A theory of finite viscoelasticity and numerical aspects. International Journal of Solids and Structures 35, 3455–3482]; Lin and Schomburg [Lin, R.C., Schomburg, U., 2003. A finite elastic–viscoelastic–elastoplastic material law with damage: theoretical and numerical aspects. Computer Methods in Applied Mechanics and Engineering 192, 1591–1627]; Lin and Brocks [Lin, R.C., Brocks, W., 2004. On a finite strain viscoplastic theory based on a new internal dissipation inequality. International Journal of Plasticity 20, 1281–1311]; and Lin and Brocks [Lin, R.C., Brocks, W., 2005. An extended Chaboche’s viscoplastic law at finite strains: theoretical and numerical aspects. Journal of Materials Science and Technology 21, 145–147]. These constitutive laws are consistent with the second law of thermodynamics. As the internal dissipation inequalities are described in different configurations or coordinate systems, the related constitutive laws are also formulated in the corresponding configurations or coordinate systems. Mathematically, these constitutive laws have very different formulations. Now, a question is whether the constitutive laws provide identical constitutive responses for the same inelastic constitutive problems. In the present work, four types of finite strain viscoelastic and endochronically plastic laws as well as three types of J₂-plasticity laws are formulated based on four types of dissipation inequalities. Then, they are numerically compared for several problems of homogeneous or complex finite deformations. It is demonstrated that for the same inelastic constitutive problem the stress responses are identical for deformation processes without rotations. In the simple shear deformation process with large rotation, the presented viscoelastic and endochronically plastic laws also show almost identical stress responses up to a shear strain of about 100%. The three laws of J₂-plasticity also produce the same shear stresses up to a shear strain of 100%, while different normal stresses are generated even at infinitesimal shear strains. The three J₂-plasticity laws are also compared at three complex finite deformation processes: billet upsetting, cylinder necking and channel forming. For the first two deformation processes similar constitutive responses are obtained, whereas for the third deformation process (with large global rotations) significant differences of constitutive responses can be observed.     

Viscoplastic approach for rate-dependent failure analysis of concrete joints and interfaces

In this work, a new rate-dependent interface model for computational analysis of quasi-brittle materials like concrete is presented. The model is formulated on the basis of the inviscid elastoplastic model by [Carol, I., Prat, P.C., López, C.M., 1997. “A normal/shear cracking model. Interface implementation for discrete analysis”. Journal of Engineering Mechanics, ASCE, 123 (8), pp. 765–773.]. The rate-dependent extension follows the continuous form of the classical viscoplastic theory by [Perzyna, P., 1966. “Fundamental problems in viscoplasticity”. Advances in Applied Mechanics, 9, pp. 244–368.]. According to [Ponthot, J.P., 1995. “Radial return extensions for viscoplasticity and lubricated friction”. In: Proceedings of International Conference on Structural Mechanics and Reactor Technology SMIRT-13, Porto Alegre, Brazil, (2), pp. 711–722.] and [Etse, G., Carosio, A., 2002. “Diffuse and localized failure predictions of Perzyna viscoplastic models for cohesive-frictional materials”. Latin American Applied Research (32), pp. 21–31.] it includes a consistency parameter and a generalized yield condition for the viscoplastic range that allows an straightforward extension of the full backward Euler method for viscoplastic materials. This approach improves the accuracy and stability of the numerical solution. The model predictions are tested against experimental results on mortar and concrete specimens that cover different stress paths at different strain rates. The results in this work demonstrate, on one hand, the capabilities of the proposed elasto–viscoplastic interface constitutive formulation to predict the rate-dependency of mortar and concrete failure behavior, and, on the other hand, the efficiency of the numerical algorithms developed for the computational implementation of the model that include the consistent tangent operator to improve the convergence rate at the finite element level.