A modified couple stress model for bending analysis of composite laminated beams with first order shear deformation
Wanji Chen,Lì Lì,Ma Xu
DOI: https://doi.org/10.1016/j.compstruct.2011.05.032
IF: 6.3
2011-01-01
Composite Structures
Abstract:Based on a modified couple stress theory, a model for composite laminated beam with first order shear deformation is developed. The characteristics of the theory are the use of rotation–displacement as dependent variable and the use of only one constant to describe the material’s micro-structural characteristics. The present model of beam can be viewed as a simplified couple stress theory in engineering mechanics. An example as a cross-ply simply supported beam subjected to cylindrical bending loads of f w = q 0 sin ( πx / L ) is adopted and explicit expression of analysis solution is obtained. Numerical results show that the present beam model can capture the scale effects of microstructure, and the deflections and stresses of the present model of couple stress beam are smaller than that by the classical beam mode. Additionally, the present model can be reduced to the classical composite laminated Timoshenko beam model, Isotropic Timoshenko beam model of couple stress theory, classical isotropic Timoshenko beam, composite laminated Bernoulli–Euler beam model of couple stress theory and isotropic Bernoulli–Euler beam of couple stress theory. Keywords Composite laminated Timoshenko bending beam Modified couple stress theory Material length parameter Scale effect 1 Introduction Since the 1960s, experiments have shown that micro-structure has scale effects due to impurities, crystal lattice mismatch and micro cracks at micro scales. With the material size scaling down to the order of micro scales, the stiffness and the strength of metal materials can increase with the size decreasing, which is called size effects. The size effects have been proved by many experiments in the recent two decades. For example, Fleck et al. [1] observed that the scaled shear strength increases by a factor of three as the wire diameter decreases from 170 μm to 12 μm in the twisting of thin copper wires; Stolken and Evan [2] reported a significant increase in the normalized bending hardening with the beam thickness decreasing in bending of ultra thin beams. Sun et al. [3] put forward a alternative view of the size effects in the nano-scale structures. As conventional continuum theory cannot explain or solve the problems of the scale effects, theories for microstructures need to be developed. Theories for microstructures include couple stress theory and strain gradient theory. A series of research in the couple stress/strain gradient theories have been made. For example, in the 1960s, Toupin [4] , Koiter [5] and Mindlin proposed strain gradient theory by introducing strain gradient into constitutive equation [6] . Between the 1980s and 1990s, Aifantis [7] , Fleck and Hutchinson [8,9] developed the strain gradient theory in plasticity. Gao et al. [10] further improved the strain gradient theory in plasticity. A modified couple stress theory has recently been proposed by Yang et al. in which the couple stress tensor is symmetric and only one internal material length scale parameter is considered [11] . The couple stress theory can be viewed as a special format of strain gradient theory which uses rotation as a variable to describe curvature, while the strain gradient theory uses strain as variable to describe curvature. Though both theories can describe the scale defects at micro-scale, the couple stress theory contains fewer rotation variables than the strain gradient theory does for the strain variables. In the couple stress theory, the variables related to micro-scale impurities or defects are formulated into rotation equilibrium equations. In the strain gradient theory, these variables are formulated into higher order strain terms in geometric equations. In both cases, new parameters which describe the material scale characteristics are introduced as higher order term (4th order) into the partial differential governing equation. Yet, in conventional continuum mechanics, this partial differential governing equation is a 2nd order equation. Generally speaking , the couple stress/strain gradient theory for microstructures can be classified into two respective theories, C 1 theory and C 0 theory. For C 1 theory the displacements and rotations/strains are dependent variables. For C 0 theory, the displacements and rotations/strains are independent variables. In the application in the engineering, the micro-structures such as sensors and actuators in micro-electromechanical systems (MEMS) and nano-electromechanical systems (NEMS) are often consist in the components of beam, plate and membrane et al. According to the application in engineering, the beam, plate and shell theories based on couple stress/strain gradient theory should be developed. The researchers have focused on the beam theory on micro-scale in recent years. A number of papers have been published for attempting to develop microstructure-dependent non-local Timoshenko beam models and apply them to analyze nanotubes and other small beam-like members/devices. All of these models are based on a C 0 theory in which the rotation–displacement as dependent variables. For example, the model for pure bending proposed by Anthoine [12] is based on the classical C 0 couple stress elasticity theory, which includes two additional internal material length scale parameters. The higher-order Bernoulli–Euler beam model developed by Papargyri-Beskou et al. [13] is based on the C 0 gradient elasticity theory, which involves two internal material length scale parameters. The non-local Bernoulli–Euler beam model by Peddieson et al. [14] , in the formulation the constitutive equation suggested by Eringen [15] contains two additional material constants. More background related to the couple stress beam based on the C 0 couple stress theory, especially Cosserat-type theories which contain more than two additional material constants, can be found in the review by Altenbach al. [16] . Recently, due to the difficulty of determining more than one microstructure-dependent length scale parameters and the approximate nature of beam theories, C 1 non-classical beam models involving only one material length scale parameter are getting many attentions. One model, as a simpler Bernoulli–Euler beam model based on modified couple stress theory with only one material length parameter, has recently been developed by Park and Gao [17] . Ma et al. proposed a microstructure-dependent Timoshenko beam model based on a modified couple stress theory with only one material length parameter [18] . Tsiatas proposed a new Kirchhoff plate model based on a modified couple stress theory [19] . Metin developed a general nonlocal beam theory based on C 0 theory [20] , where the nonlocal constitutive equations proposed by Eringen [15] are adopted. The nonclassical R–L beam model based on the higher order shear deformation theory and C 1 couple stress theory was developed by Ma et al. [21] . The non-classical R–L model can be reduced to the existing classical elasticity-based R–L model by using the material length scale parameter and Poisson’s ratio are both taken to be zero. The classical R–L beam model [22] is a third-order beam model satisfied the condition of shear stress equal zero on the upper and lower surfaces of the beam. For moderate thickness beam, the accuracy is higher than first-order shear beam model. Furthermore the R–L beam model can be reduce the non-classical Bernoulli–Euler beam model when the normality assumption is introduced. Composite laminate beam and plate are widely used in engineering. Due to the microscale such as fiber, impurities and micro cracks at micro matrix are involved in a laminated composite structure, it results in classical laminate theory invalid in some problems related to the miro-scale of laminate composites. The objective of this paper is to develop a microstructure-dependent model for the laminated Timoshenko beam based on a modified couple stress theory with only one material length scale parameter. 2 Formulations for modified couple stress theories Unlike the conventional continuum mechanics, the rotation vector ω i is introduced to kinematic relation of the classical couple stress theory, as well as the curvatures tensor χ ij and couple stress tensor m ij . Unlike the classical couple stress theory, Yang et al. [11] developed a modified couple stress theory in which the part of rotation gradient in the strain tensor is symmetric. 2.1 Modified coupled stress theory According to the symmetric couple stress theory proposed by Yang et al., the strain tensor and curvature tensor can be defined as ε ij = 1 2 ( u i , j + u j , i ) , χ ij = 1 2 ( ω i , j + ω j , i ) respectively, where ω = 1 2 curl u , u = ( u i ) is the displacement vector and ω ( ω ι ) is the rotation vector. The main differences of modified couple stress theory with standard couple stress theory are that for modified couple stress theory the couple stress tensor is symmetric and only one internal material length scale parameter is considered [11] , however, for standard couple stress theory, the couple stress tensor is asymmetric and number of internal material length scale parameters is one not always. The beam theory is a special plane problem of the plane elasticity, so the related 2-D couple stress theories can be given as follows. Considering conventional representation in the engineering, the component representation for the couple stress theory is adopted. 2.2 Formulations of plane modified couple stress theory (C 1 theory) The displacements are represented by u and v , which are displacements along x and y directions. Consider the strain tensor and curvature tensor can be defined respectively as ε ij = 1 2 ( u i , j + u j , i ) , χ ij = 1 2 ( ω i , j + ω j , i ) , we introduce γ xy = γ 12 + γ 21 , χ x = χ 13 + χ 31 , χ y = χ 23 + χ 32 . The geometric equations can be written as: (2-1) ε x = ∂ u ∂ x ε y = ∂ v ∂ y γ xy = ∂ v ∂ x + ∂ u ∂ y χ x = 1 2 ∂ 2 v ∂ x 2 - ∂ 2 u ∂ x ∂ y χ y = 1 2 ∂ 2 v ∂ x ∂ y - ∂ 2 u ∂ y 2 Strain : { ε x , ε y , γ xy , χ x , χ y } . Stress : { σ x , σ y , τ xy , m x , m y } . where ε x , ε y , andγ xy are normal and shear strains in continuum mechanics. σ x , σ y , τ xy are normal and shear stresses in continuum mechanics. χ x , χ y are curvatures and torsional shear strain for microstructures. m x , m y are bending momentums and torsional shear stress for microstructures. The constitutive equations can be written as: (2-2) σ = [ σ x σ y τ xy ] T = D ε where D = D 1 μ D 1 μ D 1 D 1 G , and (2-3) m = m x m y = 2 ℓ 2 G 2 ℓ 2 G χ x χ y ℓ are an internal material length scale parameter, and E , μ are constants of elasticity, D 1 = E ( 1 - μ 2 ) , G = E 2 ( 1 + μ ) . The strain energy is expressed as (2-4-1) U = ∫ V ( ( σ x ε x + σ y ε y + τ xy γ xy ) + ( m x χ x + m y χ y ) ) dx dy Substituting (2-1) into (2-4) , we have, (2-4-2) U = ∫ V σ x ∂ u ∂ x + σ y ∂ v ∂ y + τ xy ∂ u ∂ y + ∂ v ∂ x + 1 2 m x ∂ 2 v ∂ x 2 - ∂ 2 u ∂ x ∂ y + 1 2 m y ∂ 2 v ∂ x ∂ y - ∂ 2 u ∂ y 2 dxdy Integrating by parts of Eq. (2-4-2) , we obtain (2-5) U = ∫ S σ x n x + τ xy n y u + ( σ y n y + τ xy n x ) v + 1 2 ( m x n x + m y n y ) ∂ v ∂ x - ∂ u ∂ y ds - ∫ V ∂ σ x ∂ x + ∂ τ xy ∂ y - 1 2 ∂ 2 m x ∂ x ∂ y + ∂ 2 m y ∂ y 2 u + ∂ σ y ∂ y + ∂ τ xy ∂ x + 1 2 ∂ 2 m x ∂ x 2 + ∂ 2 m y ∂ x ∂ y v dx dy = ∫ S σ x n x + τ xy n y - 1 2 ∂ m x ∂ x + ∂ m y ∂ y n y u + σ y n y + τ xy n x + 1 2 ∂ m x ∂ x + ∂ m y ∂ y n x v ds + ∫ S 1 2 ( m x n x + m y n y ) ∂ v ∂ x - ∂ u ∂ y ds From (2-5) , following governing equations can be obtained. The equilibrium equations (no body forces) are (2-6) ∂ σ x ∂ x + ∂ τ xy ∂ y - 1 2 ∂ 2 m x ∂ x ∂ y + ∂ 2 m y ∂ y 2 = 0 ∂ τ xy ∂ x + ∂ σ y ∂ y + 1 2 ∂ 2 m x ∂ x 2 + ∂ 2 m y ∂ x ∂ y = 0 The boundary forces are (2-7) T = T mx T my T ω = σ x n x + τ s n y - 1 2 ∂ m x ∂ x + ∂ m y ∂ y n y σ y n y + τ s n x + 1 2 ∂ m x ∂ x + ∂ m y ∂ y n x m x n x + m y n y The boundary displacements are (2-8) u ˜ = u v 1 2 ∂ v ∂ x - ∂ u ∂ y Substituting (2-1), (2-2), (2-3) into (2-6) , The couple stress equilibrium equation in terms of displacements: (2-9) E 1 - μ 2 ∂ 2 u ∂ x 2 + ( 1 - μ ) 2 ∂ 2 u ∂ y 2 + ( 1 + μ ) 2 ∂ 2 v ∂ x ∂ y - 1 2 ℓ 2 G ∂ 4 u ∂ x 2 ∂ y 2 - ∂ 4 v ∂ x 3 ∂ y + ∂ 4 u ∂ y 4 - ∂ 4 v ∂ x ∂ y 3 = 0 E 1 - μ 2 ( 1 + μ ) 2 ∂ 2 u ∂ x ∂ y + ( 1 - μ ) 2 ∂ 2 v ∂ x 2 + ∂ 2 v ∂ y 2 + 1 2 ℓ 2 G ∂ 4 u ∂ x 3 ∂ y - ∂ 4 v ∂ x 4 + ∂ 4 u ∂ x ∂ y 3 - ∂ 4 v ∂ x 2 ∂ y 2 = 0 3 Basic equations of composite laminated beam of modified couple stress theory In the point of view of theory of elasticity, the beam theory can be described by introducing the hypothesis of the cross-section into the plane elasticity. It is also true for the composite laminated beam for the couple stress theory. Considering conventional representation of beam theory in the engineering, the x – y coordinate of the plane is replaced by x – z coordinate shown in Fig. 1 . Based on the couple stress theory, only ω y is included among the rotations are ω x = 0 and ω z = 0. 3.1 hypothesis of composite laminated beam of modified couple stress theory The hypothesis of the cross-section of the classical beam can be adopted in the couple stress theory of the Timoshiko beam [18] . In order to avoid distortion and warping of beam section under pure bending, the fiber orientation of the composite laminated beam should be orthogonal. Assumed displacements in a section of composite laminated beam can be described by (3-1) u ( x , z ) = u 0 ( x ) - z θ ( x ) v = 0 w = w ( x ) where θ is the angle of rotation around the y -axis of the cross-section (see Fig. 1 ). Substituting (3-1) into the expression of the rotation as ω = 1 2 curl u , we have, (3-2) ω x = 1 2 ( w , y - v , z ) = 0 ω y = 1 2 ( u , z - w , x ) = - 1 2 ( θ + w , x ) ω z = 1 2 ( v , x - u , y ) = 0 3.2 Strain of composite laminated beam of modified couple stress theory Consider the strain tensor and curvature tensor can be defined respectively as ε ij = 1 2 ( u i , j + u j , i ) , χ ij = 1 2 ( ω i , j + ω j , i ) ∘ According to the engineering conventional representation, the strain tensor and curvature tensor for the beam can be expressed in the vector form as follows: ε x = u , x , γ xz = γ zx = 2 γ 13 , χ xy = χ 12 , χ yx = χ 21 , and ε y = γ xy = γ yz = χ x = χ y = χ xz = χ yz = 0 . The strain can be written as (3-3) ε = ε x γ xz γ zx χ xy χ yx Substituting (3-1), (3-2) into (3-3) , we have (3-4) ε x = ∂ u ∂ x = du 0 dx - z d θ dx γ xz = γ zx = 1 2 ∂ u ∂ z + dw dx = 1 2 dw dx - θ χ yx = χ xy = 1 2 d ω y dx + d ω x dy = 1 2 d ω y dx = - 1 4 d θ dx + d 2 w dx 2 3.3 Constitutive relations of composite laminated beam of modified couple stress theory The constitutive relations of composite laminated beam are defined in layer-by-layer. The stress–strain relations of k th layer in the local coordinate ( x ′, z ′) can be expressed as follows (3-5) σ k = C k ε where (3 - 6) σ k = σ x ′ k τ x ′ z ′ k τ z ′ x ′ k m x ′ y ′ k m y ′ x ′ k T (3-7) ε = ε x ′ γ x ′ z ′ γ z ′ x ′ χ x ′ y ′ χ y ′ x ′ T (3-8) C k = c 11 k c 44 k c 44 k 2 ℓ 2 c 44 k 2 ℓ 2 c 44 k where x ′ aligns with the direction of the fiber in k th layer, C 11 k = E 1 k 1 - v 12 k 2 , C 44 k = G 12 k , v 21 = E 2 k v 12 k E 1 k , E 1 k is elastic constant of k th layer, G 12 k is shear elastic constant of k th layer, v 12 k is Poisson ratio of k th layer, in which subscripts 1 and 2 represent the direction of fiber and matrix, respectively. After coordinate transformation, the stress–strain relations of k th layer in the global coordinate ( x , z ) can be written as follows (3-9) σ k = Q k ε where (3 - 10) σ k = σ x k τ xz k τ zx k m xy k m yx k T (3-11) ε = ε x γ xz γ zx χ xy χ yx T (3-12) Q k = Q 11 k Q 44 k Q 44 k 2 ℓ 2 Q ∼ 44 k 2 ℓ 2 Q ∼ 44 k . The components of Q k are expressed as (3-13) Q 11 k = m 4 C 11 k + n 4 C 22 k Q 44 k = C 44 k m 2 + C 44 k n 2 = C 44 k (3-14) Q ∼ 44 k = C 44 k H ( ϕ k ) where m = conϕ k , n = sin ϕ k and ϕ k is angle of ply, H ( ϕ k ) = 0 when ϕ k = 0 ° 1 when ϕ k = 90 ° . In order to avoid distortion and warping of beam section under pure bending, the effect of couple stress can be ignored when angle of ply is ϕ k = 0. 3.4 principle of virtual work for composite laminated beam of modified couple stress theory It is well known that the principle of virtual work can be used to derive the equilibrium equation and the boundary condition. The principle of virtual work for composite laminated beam of couple stress theory ca be given by (3-16) δ U - δ W = 0 where (3-17) δ U = ∑ k - 1 n δ U k = ∑ k - 1 n ∫ Ω k ε T Q k δ ε dv (3-18) δ W = ∫ Ω f ¯ T δ u dv + ∫ δ Ω T ¯ T δ u ds where f ¯ T and T ¯ T are body force and boundary force respectively. Substituting equation (3-4) and (3-5) and (3-15) into the equation (3-17) , by the integration of the y and z coordinates in the section of beam, the equation of beam becomes (3-19) δ U = ∑ k = 1 n δ U k = ∑ k = 1 n ∫ V k σ kT δ ε dx dy dz = ∫ Ω ∑ k = 1 n ∫ - h k 2 h k 2 σ kT δ ε dz dx dy = ∑ k = 1 n ∫ Ω k σ x k δ ε x + τ xz k δ γ xz + τ yz k δ γ yz + m xy k δ χ xy + m yx k δ χ yx dv = ∑ k = 1 n ∫ Ω k σ x k δ ε x + 2 τ xz k δ γ xz + 2 m xy k δ χ xy dv = ∫ 0 L - dN dx δ u 0 - dQ dx + 1 2 d 2 Y dx 2 δ w + dM dx - Q + 1 2 dY dx δ θ dx + N δ u 0 + Q + 1 2 dY dx δ w - Y 2 δ dw dx - M + Y 2 δ θ x = 0 x = L where N , M , Q are the classical tractions of the beam, Y is the traction of couple stress moment of the beam. They are (3-20) N = ∑ k = 1 n ∫ A k σ x k dA , M = ∑ k = 1 n ∫ A k ( σ x k z ) dA Q = ∑ k = 1 n ∫ A k τ xz k dA , Y = ∑ k = 1 n ∫ A k m xy k dA The expression of the work by the external forces on the beam in the modified couple stress theory can be expressed as (3-21) δ W = ∫ l ( f u δ u 0 + f w δ w + f c δ ω ) dx + [ N ¯ δ u 0 + V ¯ δ w + M ¯ δ θ ] | x = 0 x = L where f u and f w are, respectively, the x - and z -components of the body force per unit length along the x -axis, f c is the y -component of body force per unit length along the x -axis, and N ¯ , V ¯ and M ¯ are the applied axial force, transverse force, and bending moment at the two ends of the beam respectively, and (3-22) ∫ ℓ f c δ ω dx = - 1 2 ∫ ℓ f c δ ( θ + w , x ) dx = - 1 2 ∫ ℓ f c δ θ dx + df c dx δ w x = 0 x = L - ∫ ℓ f c δ θ dx Substituting Eqs. (3-19) and (3-21) into the equation (3-16) we have (3-23) ∫ 0 L - dN dx + f u δ u 0 - dQ dx + 1 2 d 2 Y dx 2 + 1 2 df c dx + f w δ w + dM dx - Q + 1 2 dY dx + 1 2 f c δ θ dx + ( N - N ¯ ) δ u 0 + Q + 1 2 dY dx + f c 2 - V ¯ δ w + - Y 2 - Y ¯ δ dw dx + - M - Y 2 - M ¯ δ θ x = 0 x = L = 0 From (4-8) , the equilibrium equations are obtained as (3-24) dN dx + f u = 0 dQ dx + 1 2 d 2 Y dx 2 + 1 2 df c dx + f w = 0 dM dx - Q + 1 2 dY dx + 1 2 f c = 0 and traction boundary conditions at x = 0 and x = L are (3-25) N = N ¯ Q + 1 2 dY dx + 1 2 f c = V ¯ Y = 0 M x + Y 2 = M ¯ and displacement boundary conditions are (3-26) u 0 = u ¯ 0 w = w ¯ dw dx = d w ¯ dx θ = θ ¯ 4 Composite laminated beam of modified couple stress theory 4.1 Equilibrium equations in terms of displacements for the composite laminated beam of modified couple stress theory with first order shear deformation Substituting geometric (3-4) and stress–strain relation (3-9) into (3-20) , we have (4-1) N = Q ¯ 11 du 0 dx - J ¯ 11 d θ dx M = J ¯ 11 du 0 dx - I ¯ 11 d θ dx Y = - ℓ 2 Q ¯ ¯ 44 2 d 2 w dx 2 + d θ dx Q = k s Q ¯ 44 dw dx - θ the equilibrium equations in terms of displacements of Timoshenko’s beam of couple stress theory can be obtained as follows (4-2) Q ¯ 11 d 2 u 0 dx 2 - J ¯ 11 d 2 θ dx 2 + f u = 0 k s Q ¯ 44 d 2 w dx 2 - d θ dx - ℓ 2 Q ¯ ¯ 44 2 d 4 w dx 4 + d 3 θ dx 3 + 1 2 df c dx + f w = 0 J ¯ 11 d 2 u 0 dx 2 - I ¯ 11 d 2 θ dx 2 - k s Q ¯ 44 dw dx - θ - ℓ 2 Q ¯ ¯ 44 2 d 3 w dx 3 + d 2 θ dx 2 + 1 2 f c = 0 where k s is the Timoshenko shear coefficient, which depends on the geometry of beam cross-section, and (4-3) Q ¯ ¯ 44 = ∑ k = 1 n Q ∼ 44 k b ( z k + 1 - z k ) Q ¯ jj = ∑ k = 1 n Q jj k b ( z k + 1 - z k ) ( j = 4 , 5 ) J ¯ ii = ∑ k = 1 n Q ii k b z k + 1 2 - z k 2 2 ( i = 1 ) I ¯ ii = ∑ k = 1 n Q ii k b z k + 1 3 - z k 3 3 ( i = 1 ) 4.2 Degradation of the composite laminated beam of modified couple stress theory 4.2.1 Classical composite laminated Timoshenko beam Substituting ℓ = 0 into the Eq. (4-1) , the equilibrium equations in terms of displacement of classical composite laminated Timoshenko beam can be given as (4-4) Q ¯ 11 d 2 u 0 dx 2 - J ¯ 11 d 2 θ dx 2 + f u = 0 k s Q ¯ 44 d 2 w dx 2 - d θ dx + 1 2 df c dx + f w = 0 J ¯ 11 d 2 u 0 dx 2 - I ¯ 11 d 2 θ dx 2 - k s Q ¯ 44 dw dx - θ + 1 2 f c = 0 4.2.2 Isotropic Timoshenko beam of modified couple stress theory For the isotropic Timoshenko beam ( H ( ϕ k ) ≡ 1), the elastic constants become Q ¯ 11 = EA , Q ¯ 44 = GA , Q ¯ 55 = GA / 2 , Q ¯ 66 = 0 , J ¯ 11 = 0 , I ¯ 11 = EI and I = bh 3 12 . Substituting these constants into the (4-4) the equilibrium equations in terms of displacements classical isotropic Timoshenko beam of couple stress theory can be obtained as (4-5) EA d 2 u 0 dx 2 + f u = 0 k s GA d 2 w dx 2 - d θ dx - ℓ 2 GA 4 d 4 w dx 4 + d 3 θ dx 3 + 1 2 df c dx + f w = 0 EI d 2 θ dx 2 + k s GA dw dx - θ + ℓ 2 GA 4 d 3 w dx 3 + d 2 θ dx 2 - 1 2 f c = 0 These are identical to results in the reference [18] . 4.2.3 Classical isotropic Timoshenko beam Substituting ℓ = 0 and f c = 0 into the (4-5) , the equilibrium equations of the Classical isotropic Timoshenko beam in terms of displacements can be obtained as follows (4-6) EA d 2 u 0 dx 2 + f u = 0 k s GA d 2 w dx 2 - d θ dx + f w = 0 EI d 2 θ dx 2 + k s GA dw dx - θ = 0 These are identical to results in the reference [18] . 4.2.4 Composite laminated Bernoulli–Euler beam of modified couple stress theory Substituting θ = dw dx into the (4-5) , the equilibrium equations in terms of displacements of composite laminated Bernoulli–Euler beam of couple stress theory can be obtained as follows (4-7) Q ¯ 11 d 2 u 0 dx 2 - J ¯ 11 d 3 w dx 3 + f u = 0 - ℓ 2 Q ¯ ¯ 44 2 d 4 w dx 4 + df c dx + f w = 0 J ¯ 11 d 3 u 0 dx 3 - I ¯ 11 + ℓ 2 2 Q ¯ ¯ 44 d 4 w dx 4 + 1 2 f c = 0 4.2.5 Isotropic Bernoulli–Euler beam of couple stress theory Substituting equations θ = dw dx , u 0 = 0 , f c = 0 and H ( ϕ k ) = 1 into the (4-7) , considering only bending deformation of beam, the equilibrium equation in terms of displacements of classical isotropic Bernoulli–Euler beam of couple stress theory can be obtained as (4-8) ( I ¯ 11 + ℓ 2 Q ¯ 44 ) d 4 w dx 4 = f w This is identical to result in the reference [17] . 5 Numerical example for scale effect: simply supported beam subjected to cylindrical bending A cross-ply simply supported beam shown in Fig. 2 is analyzed here. The beam is only subjected to cylindrical bending loads of f w = q 0 sin π x L , i.e., f u = f c = 0. The thickness and the material properties of each layer are uniform. Boundary conditions: (5 - 1) u | x = 0 = u | x = L = 0 , w | x = 0 = w | x = L = 0 , (5-2) Y | x = 0 = Y | x = L = 0 , M | x = 0 = M | x = L = 0 . Substituting (4-1) into the (4-2) , in terms of θ and w , the (5-2) become (5-3) d 2 w dx 2 x = 0 = d 2 w dx 2 x = L = 0 , d θ dx x = 0 = d θ dx x = L = 0 5.1 Solution of the composite laminated beam of modified couple stress theory The trial function is assumed as (5-4) u 0 ( x ) = 0 , w ( x ) = w 0 sin π x L , θ ( x ) = θ 0 cos π x L Substituting (5-4) into (5-2) , we have (5 - 5) π 2 L 2 k s Q ¯ 44 + π 2 ℓ 2 Q ¯ ¯ 44 2 L 2 w 0 + π L π 2 ℓ 2 Q ¯ ¯ 44 2 L 2 - k s Q ¯ 44 θ 0 - q 0 = 0 (5-6) π L π 2 ℓ 2 Q ¯ ¯ 44 2 L 2 - k s Q ¯ 44 w 0 + k s Q ¯ 44 + π 2 Q ¯ ¯ 44 ℓ 2 2 L 2 + π 2 I ¯ 11 L 2 θ 0 = 0 The solution of the displacements in the center of beam is obtained as follows (5 - 7) w 0 = q 0 L 4 2 k s Q ¯ 44 L 2 + 2 π 2 I ¯ 11 + π 2 Q ¯ ¯ 44 ℓ 2 π 4 2 k s I ¯ 11 Q ¯ 44 L 2 + ℓ 2 Q ¯ ¯ 44 4 k s Q ¯ 44 L 2 + π 2 I ¯ 11 (5-8) θ 0 = q 0 L 3 2 k s Q ¯ 44 L 2 - π 2 ℓ 2 Q ¯ ¯ 44 π 3 2 k s I ¯ 11 Q ¯ 44 L 2 + ℓ 2 Q ¯ ¯ 44 4 k s Q ¯ 44 L 2 + π 2 I ¯ 11 where k s = 5 + 5 ν 12 6 + 5 ν 12 for the rectangular section. The stress in the beam is (5-9) σ x k = π Q 11 k θ 0 z L sin π x L 5.2 Solution of Bernoulli–Euler beam of modified couple stress theory By using the trial functions of w ( x ) = w 0 sin π x L , the solution of the displacement in the center of beam can be obtained (5-10) w 0 = q 0 L 4 π 4 I ¯ 11 + ℓ 2 Q ¯ 44 σ x k = zQ 11 k π L 2 w 0 sin π x L 5.2.1 Numerical examples for the scale effects of microstructure In order to test characteristics of the scale effects of microstructure, models of simply supported laminated cross-ply beam are adopted. The sizes of the beam model are width b = 25 μm, thickness h = 25 μm, length L = 200 μm. Cylindrical bending load is q 0 = 1 N mm. The material constants [23] : E 2 = 6.98 GPa, E 1 = 25 E 2 , G 12 = 0.5 E 2 , G 22 = 0.25 E 2 , ν 12 = ν 22 = 0.25, in which subscripts 1 and 2 represent the direction of fiber and matrix, respectively. We choose the next two types of cross-ply laminated beam with three-layer as follows. The parameters of the first one [0°/90°/0°] are identical to Q ¯ 44 E 2 bh = 0.4 , Q ¯ 55 E 2 bh = 0.2 , 12 I ¯ 11 E 2 bh 3 = 2.014 , Q ¯ ¯ 44 [ 0 ° / 90 ° / 0 ° ] = (0/ 0.4/0). The parameters of second one [90°/0°/90°] are identical to Q ¯ 44 E 2 bh = 0.3 , Q ¯ 55 E 2 bh = 0.15 , 12 I ¯ 11 E 2 bh 3 = 0.158 , Q ¯ ¯ 44 [ 90 ° / 0 ° / 90 ° ] = ( 0.3 / 0 / 0.3 ) . Next, keep the thickness of the beam constant and change the material constant l to examine the scale effect. Numerical results of the deflection of the beam are given in Fig. 3 , which show that the deflection of the beam in couple stress theory is smaller than that in the classical elasticity as the material constant l increases. Numerical results of the deflection of the beam are given in Fig. 3 which show that the deflection of the beam in couple stress theory is smaller than that in the classical elasticity as the material constant l increases. Numerical results of the angle of rotation of the beam are given in Fig. 4 , which show that the angle of rotation of the beam in couple stress theory is smaller than that in the classical elasticity as the material constant l increases. Numerical results of the stress in section of the beam are given in Fig. 5 , which show that the stress in the section of the beam in couple stress theory is smaller than that in the classical elasticity as the material constant l increases. 5.2.2 Numerical examples to compare Euler–Bernoulli and Timoshenko beam couple stress theories for microstructures In order to compare Euler–Bernoulli and Timoshenko beam couple stress theories for microstructures, aforementioned models of simply supported laminated cross-ply beam are adopted. However, various sizes of the beam are chosen firstly as length L = 200 μm and thickness h = 25 μm, and secondly for a slender beam with a large aspect ratio, length L = 2000 μm and thickness h = 25 μm. We choose the cross-ply laminated beam with three-layer of [0°/90°/0°] and [90°/0°/90°], respectively, change the material constant as l = (0, h /2, h ), respectively, to examine the scale effect. Numerical results of the deflection of the beam are given in Figs. 6 and 7 which show that the deference of the Timoshenko beam in couple stress theory is less than Euler–Bernoulli beam of couple stress theory for Medium thickness beam with h / L = 0.125 and for thin beam with h / L = 0.0125 under the material constant l as the same. Numerical results of the stress in section of the beam are given in Figs. 8 and 9 , which show that the stress σ x of the Timoshenko beam in couple stress theory is smaller than Euler–Bernoulli beam of couple stress theory for Medium thickness beam with h / L = 0.125 and for thin beam with h / L = 0.0125 under the material constant l as the same. 6 Conclusion A new model for composite laminated beam with first order shear deformation on the couple stress theory is developed. The characteristics of the couple stress theory are the use of rotation–displacement as dependent variables and the use of only one constant to describe the material’s micro-structural characteristics. By introducing the hypothesis of the cross-section of beam, the governing equations of the composite laminated beam of couple stress theory are established by the principle of virtual work. In order to avoid distortion and warping of beam section under pure bending, the fiber orientation of the composite laminated beam should be orthogonal. The present model of beam can be viewed as a simplified couple stress theory in engineering mechanics. A cross-ply simply supported beam subjected to cylindrical bending loads of f w = q 0 sin ( πx / L ) is solved by directly applying the newly developed beam model. Numerical results show that the present beam model can capture the scale effects of microstructure. The deflections and stresses of the present model of beam of couple stress theory are always smaller than that by the classical beam model. Additionally, the present model can be reduced directly to the classical composite laminated Timoshenko beam, Isotropic Timoshenko beam of couple stress theory, classical isotropic Timoshenko beam, composite laminated Bernoulli–Euler beam of couple stress theory and isotropic Bernoulli–Euler beam of couple stress theory. Conflict of interest None declare. Acknowledgement The work in this paper was supported by the National Natural Sciences Foundation of China (No. 11072156 ). This support is gratefully acknowledged. References [1] N.A. Fleck G.M. Muller M.F. Ashby J.W. Ashby Strain gradient plasticity: theory and experiment Acta Metall Mater 42 2 1994 475 487 [2] J.S. Stolken A.G. Evans A microbend test method for measuring the plasticity length-scale Acta Mater 46 14 1998 5109 5115 [3] Z.H. Sun X.X. Wang A.K. Soh H.A. Wu Y. Wang Bending of nanoscale structures: Inconsistency between atomistic simulation and strain gradient elasticity solution Comput Mater Sci 40 1 2007 108 113 [4] R.A. Toupin Elastic materials with couple-stresses Arch Rational Mech Anal 11 1962 385 414 [5] W.T. Koiter Couple stresses in the theory of elasticity I & II Proc K Ned Akad Wet (B) 67 1964 17 44 [6] R.D. Mindlin Microstructure in linear elasticity Arch Rational Mech Anal 16 1964 51 78 [7] E.C. Aifantis On the microstructural origin of certain inelastic models Trans ASME J Eng Mater Tech 106 1984 326 330 [8] N.A. Fleck J.W. Hutchinson A phenomenological theory for strain gradient effects in plasticity J Mech Phys Solids 41 1993 1825 1857 [9] N.A. Fleck J.W. Hutchinson Strain gradient plasticity Adv Appl Mech 33 1997 295 361 [10] H. Gao Y. Huang W.D. Nix J.W. Hutchinson Mechanism-based strain gradient plasticity – I. Theory J Mech Phys Solids 47 1999 1239 1263 [11] F. Yang A.M. Chong D.C.C. Lam P. Tong Couple stress based strain gradient theory of elasticity Int J Solids Struct 39 2002 2731 2743 [12] A. Anthoine Effect of couple-stresses on the elastic bending of beams Int J Solids Struct 37 2000 1003 1018 [13] S. Papargyri-Beskou K.G. Tsepoura D. Polyzos D.E. Beskos Bending and stability analysis of gradient elastic beams Int J Solids Struct 40 2003 385 400 [14] J. Peddieson G.R. Buchanan R.P. McNitt Application of nonlocal continuum models to nanotechnology Int J Eng Sci 41 2003 305 312 [15] A.C. Eringen On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves J Appl Phys 54 1983 4703 4710 [16] J. Altenbach H. Altenbach V.A. Eremeyev On generalized Cosserat-type theories of plates and shells: a short review and bibliography Arch Appl Mech 80 2010 73 92 [17] S.K. Park X.L. Gao Bernoulli–Euler beam model based on a modified couple stress theory J Micromech Microeng 16 2006 2355 2359 [18] H.M. Ma X.L. Gao J.N. Reddy A microstructure-dependent Timoshenko beam model based on a modified couple stress theory J Mech Phys Solids 56 2008 3379 3391 [19] G.C. Tsiatas A new Kirchhoff plate model based on a modified couple stress theory Int J Solids Struct 46 2009 2757 2764 [20] Aydogdu Metin A general nonlocal beam theory: its application to nanobeam bending, buckling and vibration Phys E: Low-Dimens Syst Nanostruct 41 9 2009 1651 1655 [21] H.M. Ma X.L. Gao J.N. Reddy A nonclassical Reddy–Levinson beam model based on a modified couple stress theory Int J Multi-Scale Computat Eng 8 2 2010 167 180 [22] J.N. Reddy a simply higher-order theory for laminated composite plates J Appl Mech 51 12 1984 745 752 [23] N.J. Pagano Exact solutions for rectangular bi-directional composites J Compos Mater 4 1970 20 34
What problem does this paper attempt to address?
-
Elastoplastic Analysis with Fine Beam Model of Vector Form Intrinsic Finite Element
X. F. Yuan,C. Chen,Y. F. Duan,R. J. Qian
DOI: https://doi.org/10.1177/1369433217718984
2017-01-01
Advances in Structural Engineering
Abstract:A new fine beam model of vector form intrinsic finite element is proposed to conduct elastoplastic analysis of beams with different depth-span ratios. Vector form intrinsic finite element is a new structural analysis method to compute large deformation, collision, and fracture problems. The presented fine beam model is different from Euler beam or Timoshenko beam. The coupled effects of multifactor are considered in this new model, such as shear deformation, additional axial displacement induced by additional rotation of cross section, and additional transverse displacements considering bending stiffness reduction due to shear deformation. Applying to vector form intrinsic finite element analysis, the formula of strain and internal force of this new beam model is derived. The plane cross-sectional assumption is also modified when the material is at plastic stage. Using the proposed model, computer programs are developed by Fortran. Comparison between the computational results by present method and numerical or experimental results in the literature shows that the present fine beam model of vector form intrinsic finite element is more efficient and accurate when the analyzed beam has a large depth-span ratio and material is at plastic stage.
-
A nonlinear couple stress-based sandwich beam theory
Bruno Reinaldo Goncalves,Anssi T. Karttunen,Jani Romanoff
DOI: https://doi.org/10.48550/arXiv.1903.07805
2019-03-19
Abstract:A geometrically nonlinear sandwich beam model founded on the modified couple stress Timoshenko beam theory with Kármán kinematics is derived and employed in the analysis of periodic sandwich structures. The constitutive model is based on the mechanical behavior of sandwich beams, with the bending response split into membrane-induced and local bending modes. A micromechanical approach based on the structural analysis of a unit cell is derived and utilized to obtain the stiffness properties of selected prismatic cores. The model is shown to be equivalent to the classical thick-face sandwich theory for the same basic assumptions. A two-node finite element interpolated with linear and cubic shape functions is proposed and its stiffness and geometric stiffness matrices are derived. Three examples illustrate the model capabilities in predicting deflections, stresses and critical buckling loads of elastic sandwich beams including elastic size effects. Good agreement is obtained throughout in comparisons with more involved finite element models.
Classical Physics
-
Nonlinear dynamics of cantilevered microbeams based on modified couple stress theory
h l dai,yan kun wang,liang wang
DOI: https://doi.org/10.1016/j.ijengsci.2015.05.007
IF: 7.155
2015-01-01
International Journal of Engineering Science
Abstract:The aim of this paper is to develop a new nonlinear theoretical model for cantilevered microbeams and to explore the nonlinear dynamics based on the modified couple stress theory, taking into account of one single material length scale parameter. The full nonlinear equation of motion, which is valid when the motion is large, is derived using the Hamilton’s principle. The governing partial differential equation is further discretized with the aid of Galerkin’s method. The numerical results, in which the existence of primary resonances of the first mode of the microbeam due to base excitations is demonstrated, are presented in the form of frequency–response curves, phase portraits and time histories. For a cantilevered microbeam subjected to harmonic base excitations, it is found that the frequency–response curve exhibits a clear softening-type behavior. For the same system but with an intermediate linear spring support, it is shown that the linear spring is capable of increasing the resonance frequency and decreasing the resonance amplitudes of the microbeam. Interestingly, it is found that the softening behavior could be changed to a hardening one if an intermediate nonlinear spring is added somewhere along the microbeam’s length.
-
A Transversely Isotropic Magneto-Electro-elastic Timoshenko Beam Model Incorporating Microstructure and Foundation Effects
G. Y. Zhang,Y. L. Qu,X. -L. Gao,F. Jin
DOI: https://doi.org/10.1016/j.mechmat.2020.103412
IF: 4.137
2020-01-01
Mechanics of Materials
Abstract:A new model is developed for transversely isotropic magneto-electro-elastic Timoshenko beams by using a variational formulation based on Hamilton's principle. An extended modified couple stress theory and a two-parameter Winkler-Pasternak elastic foundation model are applied to describe the microstructure and foundation effects, respectively. The equations of motion and complete boundary conditions for the beam are simultaneously obtained. The new model includes the models for transversely isotropic piezoelectric and piezomagnetic Timoshenko beams and the model for transversely isotropic magneto-electro-elastic Bernoulli–Euler beams, all incorporating the microstructure and foundation effects, as special cases. To illustrate the newly developed beam model, the static bending and free vibration problems of a simply supported transversely isotropic magneto-electro-elastic beam subjected to a uniformly distributed load are analytically solved by directly applying the general formulas derived. Numerical results are also provided to qualitatively show the microstructure, foundation and magneto-electro-elastic coupling effects on the responses of the simply supported beam.
-
Nonlinear Non-Classical Microscale Beams: Static Bending, Postbuckling and Free Vibration
W. Xia,L. Wang,L. Yin
DOI: https://doi.org/10.1016/j.ijengsci.2010.04.010
IF: 7.155
2010-01-01
International Journal of Engineering Science
Abstract:This paper initiates the theoretical analysis of nonlinear microbeams and investigates the static bending, postbuckling and free vibration. The nonlinear model is conducted within the context of non-classical continuum mechanics, by introducing a material length scale parameter. The nonlinear equation of motion, in which the nonlinear term is associated with the mean axial extension of the beam, is derived by using a combination of the modified couple stress theory and Hamilton’s principle. Based on this newly developed model, calculations have been performed for microbeams simply supported between two immobile supports. The static deflections of a bending beam subjected to transverse force, the critical buckling loads and buckled configurations of an axially loaded beam, and the nonlinear frequencies of a beam with initial lateral displacement are discussed. It is shown that the size effect is significant when the ratio of characteristic thickness to internal material length scale parameter is approximately equal to one, but is diminishing with the increase of the ratio. Our results also indicate that the nonlinearity has a great effect on the static and dynamic behaviors of microscale beams. To attain accurate and reliable characterization of the static and dynamic properties of microscale beams, therefore, both the microstructure-dependent parameters and the nonlinearities have to be incorporated in the design of microscale beam devices and systems.
-
Modified Timoshenko Beam Model for Bending Behaviors of Layered Materials and Structures
Huasong Qin,Yingbo Yan,Huichao Liu,Jingran Liu,Yong-Wei Zhang,Yilun Liu
DOI: https://doi.org/10.1016/j.eml.2020.100799
IF: 4.728
2020-01-01
Extreme Mechanics Letters
Abstract:Layered materials and structures (LMS), such as van der Waals two-dimensional (2D) layered materials and nacre-like layered structures, often exhibit highly anisotropic mechanical properties, i.e., strong in in-plane directions but weak in out-of-plane direction. Despite the strong anisotropy in their mechanical properties, Timoshenko beam model (TBM) is usually used to describe the bending deformation of LMS. We note, however, that there are two fundamental issues in using TBM to describe LMS: First, the stiffness of LMS approaches zero when the interlayer shear modulus G approaches zero; and second, the first derivative of deflection becomes discontinuous at the point of concentrated force. Clearly, both are not true for LMS. In this work, by introducing the bending energy of monolayer into the potential energy of TBM, we develop a modified Timoshenko beam model (MTBM), which is able to not only address these two issues, but also correctly predict the bending stiffness of LMS without any fitting parameters. Our analysis shows that the bending behaviors of LMS are determined by a dimensionless parameter λL, where L is the length of the beam and λ=kGA∕D0+kGA∕(nDbend), where, kGA and D0 are, respectively, the shear and bending rigidity of the beam cross-section, Dbend is the bending rigidity of monolayer, and n is the number of layer. When λL→0, the MTBM degenerates to the multi-beam model with bending stiffness of nDbend; while it degenerates to the TBM whenλL→∞. Furthermore, if kGA is much larger than D0, both MTBM and TBM degenerate to the classical Euler–Bernoulli beam model. We further perform molecular dynamics simulations, finite element simulations and experiments to validate the MTBM. Based on the MTBM, a couple of interesting applications of LMS are also demonstrated. Hence, the MTBM presented here captures the necessary intrinsic deformation modes of LMS and provides an accurate tool for the prediction and optimization of the mechanical properties of LMS.
-
Non-classical Timoshenko beam element based on the strain gradient elasticity theory
Bo Zhang,Yuming He,Dabiao Liu,Zhipeng Gan,Lei Shen
DOI: https://doi.org/10.1016/j.finel.2013.10.004
IF: 2.618
2014-01-01
Finite Elements in Analysis and Design
Abstract:This paper presents a novel Timoshenko beam element based on the framework of strain gradient elasticity theory for the analysis of the static bending, free vibration and buckling behaviors of Timoshenko microbeams. The element proposed is a two-node element which has 6-DOF (degrees of freedom) at each node considering both bending and stretching deformations, and 4-DOF considering only bending deformation. Unlike the classical Timoshenko beam element, the current element satisfies the C"0 continuity and C"1 weak continuity and contains three material length scale parameters to capture the size effect. Finite element formulations are derived by utilizing the corresponding weak form equations. Convergence, shear locking and comparison studies are carried out to examine the reliability and accuracy of the numerical solutions. The shear locking study shows that the present beam element is free of shear locking. Besides, it is established that there is a good agreement between the present results with the results in existing literature. To further illustrate the applicability and accuracy of the new Timoshenko beam element, the static bending, free vibration and buckling problems of microbeams with various boundary conditions are covered by the analysis. The results show that such small size effects are significant when the beam thickness is small, but become negligible with increasing beam thickness. Some results are believed to be the first known in the open literature and can be used as a benchmark for further studies.
-
On Nonlocal Vertical and Horizontal Bending of a Micro-Beam
Chengxiu Zhu,Yingting Chen,Jingbo Zhao,Cheng Li,Zuxiang Lei
DOI: https://doi.org/10.1155/2022/5121377
IF: 1.43
2022-06-04
Mathematical Problems in Engineering
Abstract:The vertical and horizontal bending of micro-beams subjected to axial compressive and transverse concentrated loadings is a common lateral deformation in micro-/nano-engineering that plays a significant role in the design and optimization of micro-/nano-devices. The present study aims to investigate it using the nonlocal theory. For this purpose, the simplified mathematical model is developed, and the nonlocal differential constitutive equation is applied. Since the mechanical properties of micro-beams are different from those of macro-beams, the non-classical nonlocal bending moment is introduced to improve the classical bending formulation in order to adapt to the vertical and horizontal bending of micro-beams. The effects of the external load, external size, structural stiffness, and internal characteristic scale on the vertical and horizontal bending deformation including the midpoint deflection and critical compression are presented. The present analytical model and results are validated by the finite element method. It is shown that the critical compression decreases with increasing the internal characteristic scale. Moreover, the midpoint deflection varies remarkably with respect to the axial and transverse loadings, structural stiffness, internal characteristic scale, and external size. An obvious nonlocal scale effect is found, in which the internal characteristic scale cannot be neglected compared with the external size. Besides, a threshold value of the structural stiffness is determined, and the connotation of the nonlocal interaction requires that the structural stiffness shall not be lower than that threshold. A mutual restriction between structural stiffness and external loadings is observed in the vertical and horizontal bending. In particular, it is further proved that the classical continuum mechanics can not be used in micro-/nano-scaled mechanics through a strange phenomenon that is contrary to mechanical common sense in the calculation example. The study is expected to be beneficial to the design and application of micro-beams subjected to the vertical and horizontal bending.
engineering, multidisciplinary,mathematics, interdisciplinary applications
-
A general dynamic theoretical model of elastic micro-structures with consideration of couple stress effects and its application in mechanical analysis of size-dependent properties
Yilin Qu,Peng Li,Feng Jin
DOI: https://doi.org/10.1007/s00707-019-02534-4
IF: 2.645
2019-10-29
Acta Mechanica
Abstract:A general and systematic theoretical framework of elastic micro-structures is established with the aid of modified couple stress theory for investigating the size-dependent property in small scale, in which the size-dependence is considered by introducing a material length scale parameter. Mathematically, dynamic governing equations and corresponding boundary conditions are derived and simplified by using single power series expansion for a micro-plate and double power series expansion for a micro-beam. It is demonstrated that this method exhibits extraordinary superiority, i.e., different vibration modes can be extracted easily from artificial truncations. This theoretical model can be reduced to some classical cases, including the Bernoulli–Euler beam, Timoshenko beam, Kirchhoff plate and Mindlin plate, if some specific assumptions are made. After validation, a systematic numerical investigation is carried out, which focuses on the couple stress effect on shear resonance of a cantilever micro-plate. Finally, a methodology for proposing the critical size that distinguishes micro-scale from macro-scale is illustrated in detail.
mechanics
-
A novel beam model suitable at the microscale
Young-Rok Kwon
DOI: https://doi.org/10.1080/15376494.2024.2303374
2024-02-14
Mechanics of Advanced Materials and Structures
Abstract:A novel beam model is proposed to address the limitations of existing microscale beam models by introducing a higher-order displacement field assumption for the modified couple stress theory. The governing equations and boundary conditions are derived using the variational principle. Validation through finite element analysis for a simply supported beam problem demonstrates the beam model successfully captures the width-dependent behavior of microbeams, which existing beam models cannot explain. Additionally, the beam model accurately represents the plane strain behavior as the beam width increases and the classical beam behavior at a small length scale parameter.
materials science, multidisciplinary,mechanics, composites, characterization & testing
-
Hierarchy of beam models for lattice core sandwich structures
Anssi T. Karttunen,J.N. Reddy
DOI: https://doi.org/10.1016/j.ijsolstr.2020.08.020
2020-10-18
Abstract:A discrete-to-continuum transformation to model 2-D discrete lattices as energetically equivalent 1-D continuum beams is developed. The study is initiated in a classical setting but results in a non-classical two-scale micropolar beam model via a novel link within a unit cell between the second-order macrorotation-gradient and the micropolar antisymmetric shear deformation. The shear deformable micropolar beam is reduced to a couple-stress and two classical lattice beam models by successive approximations. The stiffness parameters for all models are given by the micropolar constitutive matrix. The four models are compared by studying stretching- and bending-dominated lattice core sandwich beams under various loads and boundary conditions. A classical 4th-order Timoshenko beam is an apt first choice for stretching-dominated beams, whereas the 6th-order micropolar model works for bending-dominated beams as well. The 6th-order couple-stress beam is often too stiff near point loads and boundaries. It is shown that the 1-D micropolar model leads to the exact 2-D lattice response in the absence of boundary effects even when the length of the 1-D beam (macrostructure) equals that of the 2-D unit cell (microstructure), that is, when L=l.
Applied Physics
-
Size-dependent couple stress Timoshenko beam theory
Ali R. Hadjesfandiari,Arezoo Hajesfandiari,Haoyu Zhang,Gary F. Dargush
DOI: https://doi.org/10.48550/arXiv.1712.08527
2017-12-21
Abstract:In this paper, a new size-dependent Timoshenko beam model is developed based on the consistent couple stress theory. In the present formulation, the governing equations and corresponding boundary conditions are obtained. Afterwards, this formulation is used to investigate size-dependency for several elementary beam problems. Analytical solutions are obtained for pure bending of a beam and for a cantilever beam with partially and fully clamped boundary conditions. These analytical results are then compared to the numerical results from a two-dimensional finite element formulation for the corresponding couple stress continuum problem.
General Physics
-
A refined quasi-3D zigzag beam theory for free vibration and stability analysis of multilayered composite beams subjected to thermomechanical loading
Bin Han,Wei-Wei Hui,Qiancheng Zhang,Zhenyu Zhao,Feng Jin,Qi Zhang,Tian Jian Lu,Bingheng Lu
DOI: https://doi.org/10.1016/j.compstruct.2018.08.005
IF: 6.3
2018-01-01
Composite Structures
Abstract:A refined four-unknown quasi-3D zigzag beam theory is developed to model the free vibration and buckling behaviors of multilayered composite beams subjected to axial mechanical loading (e.g., distributed load and terminal force) and uniform temperature variation. Types of the composite beams considered include laminated composite beams, sandwich beams with composite face sheets, and fiber metal laminates. The proposed theory accounts for not only thickness stretching but also interlaminar continuity of transverse shear stresses and displacements. Associated eigenvalue problems for various boundary conditions are derived using the Ritz method. Accuracy and effectiveness of the theoretical predictions are verified by comparison with existing results and present finite element simulations. The theory is employed to quantify the effects of axial distributed load/terminal force and temperature variation on free vibration and buckling for different boundary conditions, geometric parameters and material properties. The present theory could produce sufficiently accurate predictions of natural frequencies and buckling capacities of multilayered beams at a very low computational cost.
-
A microstructure-dependent Kirchhoff plate model based on a reformulated strain gradient elasticity theory
Gongye Zhang,Chenyi Zheng,Changwen Mi,Xin-Lin Gao
DOI: https://doi.org/10.1080/15376494.2020.1870054
2021-01-18
Mechanics of Advanced Materials and Structures
Abstract:<span>A new microstructure-dependent non-classical model for Kirchhoff plates is developed by using a reformulated strain gradient elasticity theory that incorporates both the strain gradient and couple stress effects. The equation of motion and the boundary conditions are simultaneously obtained through a variational formulation based on Hamilton's principle. The new plate model contains one material constant to account for the strain gradient effect and one material length scale parameter to capture the couple stress effect. The newly developed non-classical plate model includes the plate model incorporating the couple stress effect alone and the plate model based on the classical elasticity as two special cases. To illustrate the new model, the buckling, static bending and free vibration problems of a simply supported rectangular plate are analytically solved by directly applying the general formulas derived. The numerical results reveal that the presence of the strain gradient and couple stress effects leads to reduced plate deflections, enlarged critical buckling loads and increased natural frequencies. These microstructure effects are significant when the plate is very thin, but they are diminishing as the plate thickness increases. These predicted trends of the size effects at the micron scale agree with those observed experimentally.</span>
materials science, multidisciplinary,mechanics, composites, characterization & testing
-
Microbeam Model and Related Differential Quadrature Finite Elements
LIU Songzheng,ZHANG Bo,SHEN Huoming,ZHANG Xu
DOI: https://doi.org/10.21656/1000-0887.410260
2021-01-01
Abstract:A size-dependent quasi-3D functionally graded (FG) microbeam model was presented within the combined framework of the modified couple stress theory and a 4-unknown higher-order shear and normal deformation theory. Then the model was applied to analyze the static bending and free vibration of FG microbeams. With the 2nd Lagrange equation, the corresponding motion equations and the appropriate boundary conditions were obtained. A 2-node 16DOF differential quadrature finite element combining the Gauss-Lobatto quadrature rule with the differential quadrature rule was constructed to handle the general static/dynamic boundary value problems of FG microbeams. A comparison study was performed to show the efficacy of the proposed theoretical model and solution method. Finally, the effects of the gradient index, the intrinsic length scale parameter, the geometrical parameters and the boundary conditions on the static and dynamic characteristics of FG microbeams were examined. Numerical results reveal that the developed beam model and element are applicable to the analysis of mechanical behaviors of FG microbeams with various slenderness ratios. Besides, introduction of the couple stress effect can significantly change the static and dynamic characteristics of FG microbeams.
-
A unified high-order model for size-dependent vibration of nanobeam based on nonlocal strain/stress gradient elasticity with surface effect
Weidong Yang,Shuo Wang,Wenbing Kang,Tao Yu,Yan Li
DOI: https://doi.org/10.1016/j.ijengsci.2022.103785
IF: 7.155
2023-01-01
International Journal of Engineering Science
Abstract:In this work, a unified high-order nanobeam model considering various high-order shear deformation beam theories is established to investigate the vibration response of nanobeam on the basis of two-phase local/nonlocal strain and stress gradient theory, as well as surface elasticity theory. The unified model also includes the Euler-Bernoulli beam model and Timoshenko beam model. Herein, the elastic dynamics governing equations and boundary conditions are derived using Hamilton's principle, and the analytical solutions, such as exact formulas for natural frequencies, are obtained by employing the Navier method for simply supported boundary conditions. The effects of local volume fraction, nonlocal parameter, material length scale parameter, shear deformation and surface energy in stress and strain-driven models are analyzed in detail, respectively. The parametric studies reveal that the two scale parameters (nonlocal parameters and material length characteristic parameters) have opposite effects on the stiffness of the nanobeams in the two driving models, while the surface parameters have the same effect on the stiffness of the two driving models. The influence of the slenderness ratio on the surface effect and scale effect is opposite, meaning that the increase of the slenderness ratio deepens the influence of the surface effect but weakens the influence of the scale effect. There are also differences in the effects of higher-order modes on the two effects. Higher modes lead to more significant scale effects, but the effect of higher modes on surface effects depends on the surface elastic properties of the material. We also find that the introduction of surface elasticity increases the gap between the TBT and other higher-order beams, which indicates the prediction results of the higher-order beam model are more accurate when both surface and nonlocal effects are considered. In addition, it is represented that the surface elasticity makes aluminum nanobeams exhibit a stiffness softening effect, while the effect of surface elasticity on the stiffness of silicon nanobeams is significantly dependent on the slenderness ratio and the number of modes.
engineering, multidisciplinary
-
Nonlinear postbuckling analysis of magneto-electro-thermo-elastic laminated microbeams based on modified couple stress theory
Yu-fang Zheng,Li-Chuan Liu,De-yong Qu,Chang-ping Chen
DOI: https://doi.org/10.1016/j.apm.2023.01.021
IF: 5.336
2023-01-23
Applied Mathematical Modelling
Abstract:In this article, the nonlinear postbuckling behavior of the magneto-electro-thermo-elastic (METE) laminated microbeams is presented. According to the modified couple stress theory (MCST) and Reddy's three-order shear deformation theory (RTSDT), in conjunctions with the von Karman geometric nonlinearity, the nonlinear static model of METE laminated microbeam is established. The nonlinear governing equations and the corresponding boundary conditions are derived by using the principle of virtual work principle. Afterwards, using an analytical method, the nonlinear postbuckling behavior of METE laminated microbeam with simple supported boundary conditions is described. Moreover, numerical examples are exhibited to reveal the effects of the material length scale parameter, slenderness ratio, temperature rise, magneto-electric potential on the critical buckling load and the nonlinear postbuckling response. Also the distribution law of magneto-electric potential through the thickness direction of the microbeam with various lay-up modes is discussed.
mechanics,engineering, multidisciplinary,mathematics, interdisciplinary applications
-
The Flapwise Bending Free Vibration Analysis of Micro-rotating Timoshenko Beams Using the Differential Transform Method
Hadi Arvin
DOI: https://doi.org/10.1177/1077546317736706
IF: 2.633
2017-10-17
Journal of Vibration and Control
Abstract:The flapwise bending free vibration analysis of isotropic rotating Timoshenko microbeams, including the size effects, is presented in this paper. A nonclassical theory, i.e. the Modified Couple Stress theory, has been employed to include the size effect in the presented formulation. By consideration of the Timoshenko beam assumptions the shear deformation and the rotary inertia effects are taken into account. The Hamilton’s principle is applied to the obtained strain and kinetic energy relations to derive the nonlinear equations of motion and the associated boundary conditions. After nondimensionalization of the equations of motion and the corresponding boundary conditions, the linearized form of the equations of motion and the accompanied boundary conditions are developed. A semi-analytical approach, i.e. the differential transform method, is implemented to achieve the flapping and axial frequencies. The achieved results are validated via comparison with the available results in the literature. The material length scale, shear deformation consideration, rotating speed and the slenderness ratio influences on the natural frequencies are examined. The results demonstrate that the slenderness ratio and the thickness to the material length scale parameter quotient are the dominant indicators in determining the usage of the nonclassical theories against classical theories. On the other hand, the precision in determination of the higher modes frequencies motivates us to implement the Timoshenko beam model instead of the Euler–Bernoulli beam model.
engineering, mechanical,mechanics,acoustics
-
Size-dependent and microinertia effects on statics and dynamics based on the modified couple stress theory with arbitrary boundary conditions
Zuokun Yang,Fuqi Zhou,Liangliang Chu,Guansuo Dui
DOI: https://doi.org/10.1007/s00707-023-03836-4
IF: 2.645
2024-01-25
Acta Mechanica
Abstract:In this work, we discuss size-dependent and microinertia effects on the static and dynamic performances of a microscale model based on the microinertia-based modified couple stress theory, a non-classical continuum theory capable of capturing the behavior of size dependence and frequency dispersion characteristics. In the framework of the variational statement, a microscale structure model is developed and the governing equations of equilibrium as well as all boundary conditions for statics and dynamics are reformulated. The developed theory is imposed to tackle microstructure-dependent Timoshenko beam model in two distinct scale parameters: the material length scale parameter is utilized to determine the size dependence and the microinertia length scale parameter is employed to describe the higher-order microrotation relation. The generally valid closed-form analytic expressions are obtained and suitable for various formats of boundaries and mechanical loads. As case studies, the predicted trends agree with those observed within the framework of the modified couple stress theory. Results indicate that the material microlength scale parameter strengthens the static deformations, while the microinertia length scale parameter weakens the dynamic frequencies. In addition, boundary conditions are also an important aspect in statics and dynamics as well as the mechanical response predicted by non-classical continuum theories.
mechanics
-
Transversely heterogeneous nonlocal Timoshenko beam theory: A reduced-order modeling via distributed-order fractional operators
Wei Ding,Sansit Patnaik,Fabio Semperlotti
DOI: https://doi.org/10.1016/j.tws.2024.111608
IF: 5.881
2024-01-20
Thin-Walled Structures
Abstract:Nonlocal effects are widely identified in micro/nano systems and have been studied extensively in recent decades. Nevertheless, existing studies mostly focus on the modeling of homogeneous nonlocal systems and as a result, cannot model micro/nanostructures with heterogeneously distributed nonlocal effects. This study introduces a reduced-order nonlocal Timoshenko beam model that facilitates precise modeling of micro-elastic beams with layer-wise heterogeneous nonlocal effects. First, a fully-resolved two-dimensional (2D) plane-strain fractional-order nonlocal elastic framework is developed to provide reference solutions for 2D nonlocal beam problems. Building upon this general nonlocal elasticity theory , an equivalent one-dimensional (1D) nonlocal Timoshenko model is then developed. The layer-wise heterogeneous nonlocal information is captured by using distributed-order (DO) operators, and the impacts arising from the nonlocal heterogeneity are further characterized by introducing two auxiliary parameters. Both 1D and 2D approaches are applied to simulate the mechanical responses of nonlocal beams. Direct comparisons of numerical simulations produced by either the 1D or the fully-resolved 2D model confirm that the DO Timoshenko beam formulation (together with the two auxiliary parameters) can capture not only the overall beam deflection but also the additional shear effect induced by the heterogeneous nonlocality. Computational cost assessment also indicates the proposed approach's superior performance. The proposed DO Timoshenko model enables model-order reduction without compromising the heterogeneous nonlocal description of the material hence leading to an efficient and accurate reduced-order nonlocal modeling approach that addresses the limitations of prior nonlocal microstructural theories and extends applicability to a broader range of heterogeneous nano/microstructures.
engineering, mechanical,mechanics, civil