Constraint effects in adhesive joint fracture

Constraint effects in adhesive joint fracture are investigated by modelling the adherents as well as a finite thickness adhesive layer in which a single row of cohesive zone elements representing the fracture process is embedded. Both the adhesive and the adherents are elastic-plastic with strain hardening. The bond toughness Γ (work per unit area) is equal to Γ₀+Γ_p, where Γ₀ is the intrinsic work of fracture associated with the embedded cohesive zone response and Γ_p is the extra contribution to the bond toughness arising from plastic dissipation and stored elastic energy within the adhesive layer. The parameters of the model are identified from experiments on two different adhesives exhibiting very different fracture properties. Most of the tests were performed using the wedge-peel test method for a variety of adhesives, adherents and wedge thicknesses. The model captures the constraint effects resulting from the change in Γ_p: (i) the plastic dissipation increases with increasing bond line thickness in the fully plastic regime and then decreases to reach a constant value for very thick adhesive layers; (ii) the plastic dissipation in the fully plastic regime increases drastically as the thickness of the adherent decreases. Finally, this model is used to assess a simpler approach which consists of simulating the full adhesive layer as a single row of cohesive elements.     

Asymptotic analysis of layered elastic beams

We consider an elastic beam formed by three layers, fixed at one end and loaded at the free end. We call adherents the upper and lower layers ${\Omega }_{+}^{?}$ and ${\Omega }_{?}^{?}$ and an adhesive layer ${\Omega }_m^{?}$. We denote by $?h_{\pm ,m}$ the thickness of each layer and we suppose that the stiffness of the adhesive layer is ${?}^2$, with respect to that of the adherents. By an asymptotic analysis we obtain the zeroth order limit problem and the form of the second order displacements. To cite this article: M. Serpilli, C. R. Mecanique 333 (2005).     

Asymptotic modelling of the linear dynamics of laminated beams

We study the linear dynamics of a layered elastic beam by means of the asymptotic expansion method. The beam consists of three linearly elastic isotropic layers: the middle layer is considered to be thinner and softer than the upper and lower ones. We characterize the limit models by distinguishing three cases of natural frequencies: the low frequencies associated with flexural vibrations, the mean frequencies associated with axial vibrations and the high frequencies, associated with transversal shear and pinching vibrations.     

Elastic field due to an edge dislocation in a multilayered composite

A numerical procedure is presented for the analysis of the elastic field due to an edge dislocation in a multilayered composite. The multilayered composite consists of n perfectly bonded layers having different material properties and thickness, and two half-planes adhere to the top and bottom layers. The stiffness matrices for each layer and the half-planes are first derived in the Fourier transform domain, then a set of global stiffness equations is assembled to solve for the transformed components of the elastic field. Since the singular part of the elastic field corresponding to the dislocation in the full-plane has been extracted from the transformed components, regular numerical integration is needed only to evaluate the inverse Fourier transform. Numerical results for the elastic field due to an edge dislocation in a bimaterial medium are shown in fairly good agreement with analytical solutions. The elastic field and the Peach–Kohler image force are also presented for an edge dislocation in a single layered half-plane, a two-layered half-plane and a multilayered composite made of alternating layers of two different materials.     

Junction of thin plates

The elasticity problem of a thin plate with an edge is considered using asymptotic methods. The small parameter ε describes the relative thickness of the plate. In the case when the elasticity coefficients are everywhere of the same order of magnitude, the asymptotic behaviour of the plate is such that the angle of the edge remains constant under the deformation (the junction is called `rigid'). We also consider a junction mode of a narrow filet (the order of its width is O(ε)) of a `soft' elastic material, the elasticity coefficients being O(ε) with respect to those of the plates. In this case (called `elastic junction'), the asymptotic modelling contains an energy bilinear form associated with the filet which involves the variation of the angle and the sliding along the junction.     

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.     

Conjugate natural convection heat transfer between two fluids separated by a horizontal wall: steady-state analysis

The conjugate heat transfer across a thin horizontal wall separating two fluids at different temperatures is investigated both numerically and asymptotically. The solution for large Rayleigh numbers is shown to depend on two nondimensional parameters;α/ε ², withα being the ratio of the thermal resistance of the boundary layer in the hot medium to the thermal resistance of the wall andε the aspect ratio of the plate, andβ, the ratio of the thermal resistances of the boundary layers in the two media. The overall Nusselt number is an increasing function ofα/ε ² taking a finite maximum value forα/ε ² → ∞ and tending to zero forα/ε ² → 0.     

Boundary Perturbations Due to the Presence of Small Linear Cracks in an Elastic Body

In this paper, Neumann cracks in elastic bodies are considered. We establish a rigorous asymptotic expansion for the boundary perturbations of the displacement (and traction) vectors that are due to the presence of a small elastic linear crack. The formula reveals that the leading order term is ε ² where ε is the length of the crack, and the ε ³-term vanishes. We obtain an asymptotic expansion of the elastic potential energy as an immediate consequence of the boundary perturbation formula. The derivation is based on layer potential techniques. It is expected that the formula would lead to very effective direct approaches for locating a collection of small elastic cracks and estimating their sizes and orientations.     

Étude théorique et numérique du comportement d'un assemblage de plaquesTheoretical and numerical study of the behaviour of bonded plates

We consider the flexural static behavior of two coplanar Kirchhoff–Love plates bonded by a thin third one, softer and narrower. In the case of perfect adhesion or unilateral conditions between the adherents and the adhesive, we study the asymptotic behavior of the assembly when the stiffness and the narrowness of the intermediate plate tend to zero. To cite this article: F. Zaittouni et al., C. R. Mecanique 330 (2002) 359–364.     

The axisymmetric double contact problem for a frictionless elastic layer

The general axisymmetric double contact problem for an elastic layer pressed against a half space by an elastic stamp is considered. The problem is solved under the assumptions that the three materials have different elastic properties, the contact along the interfaces is frictionless and only compressive normal tractions can be transmitted across the interfaces, and, in the case of the elastic stamp, the local radius of curvature of the stamp is large compared to the stamp-layer contact radius. The problem is reduced to a system of singular integral equations in which the contact pressures are the unknown functions. The solution is obtained and extensive numerical results are given for three stamp geometries, namely, rigid and elastic spherical stamps, and a flat-ended rigid cylindrical stamp. The results show that in the case of a flat-ended rigid cylindrical stamp the radius b of the contact area between the layer and the subspace is independent of the magnitude P of the total transmitted load and in all other cases b will depend on P.     

Bounds on the distribution of extreme values for the stress in composite materials

Suitable macroscopic quantities beyond effective elastic properties are used to assess the distribution of stress within a composite. The composite is composed of N anisotropic linearly elastic materials and the length scale of the microstructure relative to the loading is denoted by ε. The stress distribution function inside the composite λ^ε(t) gives the volume of the set where the norm of the stress exceeds the value t. The analysis focuses on the case when 0<ε?1. A rigorous upper bound on lim_ε→0λ^ε(t) is found. The bound is given in terms of a macroscopic quantity called the macro stress modulation function. It is used to provide a rigorous assessment of the volume of over stressed regions near stress concentrators generated by reentrant corners or by an abrupt change of boundary loading.     

The validity range of CPT and Mindlin plate theory in comparison with 3-D vibrational analysis of circular plates on the elastic foundation

The validity and the range of applicability of the classical plate theory (CPT) and the first-order shear deformation plate theory, also called Mindlin plate theory (MPT), in comparison with three-dimensional (3-D) p-Ritz solution are presented for freely vibrating circular plates on the elastic foundation with different boundary conditions. In order to achieve this purpose, a study of the 3-D elasticity solution is carried out to determine the free vibration frequencies of clamped, simply supported and free circular plates resting on an elastic foundation. The Pasternak model with adding a shear layer to the Winkler model is used for describing the elastic foundation. In addition to being employed the p-Ritz algorithm, the analysis is based on the linear, small strain and 3-D elasticity theory. In this analysis method, a set of orthogonal polynomial series in a cylindrical polar coordinate system is used to arrive eigenvalue equation yielding the natural frequencies for the circular plates. The accuracy of these results is verified by appropriate convergence studies and checked with the available literature and the MPT. Furthermore, the effect of the foundation stiffness parameters, thickness-radius ratio, and different boundary conditions on the ill-conditioning of the mass matrix as well as on the vibration behavior of the circular plates is investigated. Afterwards, the validity and the range of applicability of the results obtained on the basis of the CPT and MPT for a thin and moderately thick circular plate with different values of the foundation stiffness parameters are graphically presented through comparing them with those obtained by the present 3-D p-Ritz solution. Finally, the phenomenon of mode shape switching is investigated in graphical forms for a wide range of the Winkler foundation stiffness parameters.     

Adaptive and anisotropic mesh strategy for thin shell problems. Case of inhibited parabolic shells

The aim of this paper is to show the reliability of an adaptive and anisotropic mesh procedure for thin shell problems. We consider singular perturbation problems only for parabolic shells whose behavior is described by the Koiter model. The corresponding system of equations, which depends on the relative thickness ε of the shell, is elliptic except at the limit for ε = 0 where it is parabolic. In a first part of this paper, we study theoretically the phenomena of internal layers appearing during the singular perturbation process, when the loading is somewhat singular. These layers have very different structures either they are along or across the asymptotic lines of the middle surface of the shell. In a second part, numerical computations are performed using a finite element software coupled with an adaptive anisotropic mesh generator. This technique enables to approach accurately the singularities and the layers predicted by the theory especially for very small values of the thickness. The efficiency of such a procedure in comparison with uniform meshes is put in a prominent position.     

Skewed Poiseuille-Couette flows of sPTT fluids in concentric annuli and channels

Analytical solutions have been derived for the helical flow of PTT fluids in concentric annuli, due to inner cylinder rotation, as well as for Poiseuille flow in a channel skewed by the movement of one plate in the spanwise direction, which constitutes a simpler solution for helical flow in the limit of very thin annuli. Since the constitutive equation is a non-linear differential equation, the axial and tangential/spanwise flows are coupled in a complex way. Expressions are derived for the radial variation of the axial and tangential velocities, as well as for the three shear stresses and the two normal stresses. For engineering purposes expressions are given relating the friction factor and the torque coefficient to the Reynolds number, the Taylor number, a nondimensional number quantifying elastic effects (εDe²) and the radius ratio. For axial dominated flows fRe and C_M are found to depend only on εDe² and the radius ratio, but as the strength of rotation increases both coefficients become dependent on the velocity ratio (ξ) which efficiently compacts the effects of Reynolds and Taylor numbers. Similar expressions are derived for the simpler planar case flow using adequate non-dimensional numbers.     

Homogenization of a contact problem for a system of densely situated punches

A linear contact problem of an elastic half space with rigid punches ε-periodically situated on a bounded part of the boundary of the elastic solid is investigated. Using the method of homogenization theory and the method of matched asymptotic expansions, the leading terms of the asymptotic solution are constructed as ε→0. The general capacity of the contact spot is introduced and some its properties are described.     

Vertical correlations in superlattices using the Grinfeld method

In this paper, we analyse the energetics of a multilayered structures like, for instance, B/A/B/A_substrate. It is well-known that a coherent pre-strained B layer on an A substrate will generally results in a corrugation of the free-surface of the B layer. This behavior is the result of stress relaxation in the B-layer and the phenomenon is known as the Asaro-Tiller-Grinfeld instability. We extend the methods used for a two-layer structure to a multilayered structure and the main application is the vertical correlation in superlattices. We analyse the energetics of a corrugated B layer which is grown on a A/B/A_substrate, where the A layers are flat but the intermediate B layer is already corrugated. We show that the self-organization of the second B layer, due to elastic interactions in the bulk, depends on the corrugation of the first B layer and the generic best situation is that of a top-on-top (also called correlated layers) vertical alignment. We also prove that the interaction energy between two successive B layers attains a maximum at a critical thickness of the intermediate A layer. This interaction energy has the same order of magnitude as the elastic energy release due to free-surface corrugation at each upper surface of a B layer.     

Axially compressed buckling of pressured multiwall carbon nanotubes

This paper studies axially compressed buckling of an individual multiwall carbon nanotube subjected to an internal or external radial pressure. The emphasis is placed on new physical phenomena due to combined axial stress and radial pressure. According to the radius-to-thickness ratio, multiwall carbon nanotubes discussed here are classified into three types: thin, thick, and (almost) solid. The critical axial stress and the buckling mode are calculated for various radial pressures, with detailed comparison to the classic results of singlelayer elastic shells under combined loadings. It is shown that the buckling mode associated with the minimum axial stress is determined uniquely for multiwall carbon nanotubes under combined axial stress and radial pressure, while it is not unique under pure axial stress. In particular, a thin N-wall nanotube (defined by the radius-to-thickness ratio larger than 5) is shown to be approximately equivalent to a single layer elastic shell whose effective bending stiffness and thickness are N times the effective bending stiffness and thickness of singlewall carbon nanotubes. Based on this result, an approximate method is suggested to substitute a multiwall nanotube of many layers by a multilayer elastic shell of fewer layers with acceptable relative errors. Especially, the present results show that the predicted increase of the critical axial stress due to an internal radial pressure appears to be in qualitative agreement with some known results for filled singlewall carbon nanotubes obtained by molecular dynamics simulations.     

Crack in functionally graded piezoelectric strip bonded to elastic surface layers under electromechanical loading

Solved is the problem of a crack in a functionally graded piezoelectric material (FGPM) bonded to two elastic surface layers. It is assumed that the elastic stiffness, piezoelectric constant, and dielectric permittivity of the FGPM vary continuously along the thickness of the strip. The outside layers are under antiplane mechanical loading and in-plane electric loading. The solution involves solving singular integral equations by application of the Gauss–Jacobi integration formula. Numerical calculations are carried out to obtain the energy density factors. Their variations with the geometric, loading and material parameters are shown graphically.     

A theoretical study of the coupling effects in piezoelectric ceramics

The objective of this study is to delineate electro-mechanical coupling in piezoceramic materials. The model system investigated is a two-dimensional linear piezoceramic strip polarized in the thickness direction, and it is subjected to local symmetric pressures on the upper and lower edges, traction-free boundary conditions on both end surfaces, and voltages on portions of the upper and lower edges. Under a simplifying assumption of the gradient of the electric potential, closed form solutions of the elastic field have been obtained. It is noticed that instead of the nine constants (including the elastic compliance constants, s_{i j}, the piezoelectric constants, d_{i j}, and the dielectric permittivity constants, ε_{i j}) , the elastic and piezoelectric characteristics of the material can be represented by three parameters, β₁, β₂ and β₃. β₁ consists of elastic compliance constants only. β₂ and β₃ signify the piezoelectric effect. Furthermore, higher values of β₂ imply a more pronounced piezoelectric effect on the elastic field. The identification of these parameters greatly facilitates the study of coupling effects in piezoelectric ceramics.