Using the displacement potential approach of orthotropic composite materials for the plane stress conditions, an orthotropic panel subjected to a combined loading at its right lateral edge is solved. Effects of fiber orientations and material isotropy on the elastic field of an orthotropic panel subjected to a combined loading are discussed. The analytical elastic solutions at different sections of the panel with fiber orientation θ=90° are compared with those of finite element predictions to ensure the reliability of our present solutions.
1. Introduction
Recently in lightweight structures, the use of orthotropic composite material is increasing due to its higher strength weight ratio compared to its corresponding isotropic material. Aircraft wings which are generally made of aluminium are stiffened. Nowadays, boron/epoxy, graphite/epoxy orthotropic composite materials are used for manufacturing of aircraft panels. In structures, the stiffened panel experiences tension and bending at the same time. Effect of fiber orientation and material isotropy on the exact elastic field is necessary for designing purpose because strength, toughness, and all other mechanical properties of any orthotropic composite depend on its fiber orientation. Under loading, the stress and displacement components of any orthotropic composite vary with the variation of fiber orientation. This type of analysis is carried out by the finite element method. Exact analytical solution of this type of problem is hardly possible to carry out till now due to lack of effective elasticity formulations of orthotropic composite materials and boundary management techniques.
The limited progress in solution of the two-dimensional problems of elasticity theory for anisotropic bodies was achieved by means of Lekhnitskii formalism [1, 2] and Stroh's formalism [3, 4]. Some solutions of the two-dimensional problems for anisotropic bodies were obtained by these approaches. An exposition of Stroh’s formalism and its application to a limited number of two-dimensional problems was given by Ting [5]. The proof of the equivalence between the Lekhnitskii and Stroh formalisms was obtained by Barnett and Kirchner [6]. The stress analysis problems are still suffering from a lot of shortcomings in spite of becoming the fundamental subject in the field of elasticity and thus are being constantly investigated [7–15].
Elasticity problems are usually formulated either in terms of deformation parameters or stress parameters. Among the existing mathematical models of plane boundary-value stress problems, the stress function approach [15] and the displacement formulation [16] are noticeable. The application of the stress function formulation in conjunction with finite-difference technique has been reported for the solution of plane elastic problems where all of the boundary conditions are prescribed in terms of stresses only [7–15, 17, 18]. Further, Conway and Ithaca [18] modified the stress function formulation in the form of Fourier integrals to the case where the material is orthotropic and obtained analytical solutions for a number of ideal problems. The shortcoming of the stress function approach is that it accepts boundary conditions only in terms of loadings of the structure. Boundary restraints specified in terms of the displacement components of the structure cannot be satisfactorily imposed on the stress function approach.
Stress analysis of composite structures is carried out using different analytical, semianalytical approaches where loading is applied along the perpendicular plane of the lamina of laminated composite [19–29] with limited mixed boundary conditions but along the directions of the plane of the composite laminas of the composite should be analyzed for the perfect design of composite structures with mixed boundary conditions. Kulikov and Plotnikova [19] focused on the use of efficient approach to exact elastic solutions of cross-ply and angle-ply composite plates. Interlaminar stresses resulting from bending of thick rectangular laminated plates with arbitrary laminations and boundary conditions are analyzed analytically based on a three-dimensional multiterm extended Kantorovich method (3DMTEKM) [20]. Using the principle of minimum total potential energy, three systems of coupled ordinary differential equations with nonhomogeneous boundary conditions are obtained [20]. Zhang and Hoa [21] proposed a limit-based approach to provide solution to the composite tube having cylindrical orthotropy. The Taylor series expansion method was used to drive new equations to replace those that are identically satisfied, while the other equations are kept to be same [21]. Singh and Shukla [22] presented the nonlinear response of laminated composite plates. The mathematical formulation of the actual physical problem of the laminated composite plate subjected to mechanical loading is presented utilizing higher order shear deformation theory and von-Karman nonlinear kinematics [10]. These nonlinear governing differential equations of equilibrium are linearized using quadratic extrapolation technique [22]. A meshfree technique based on multiquadric RBF is used for analysis of the problems, and isotropic, orthotropic, and laminated composite plates with immovable simply supported and clamped edges are analyzed [22]. A semianalytical method for bending of corrugated-core, honeycomb-core, and X-core sandwich panel is presented [23]. In the global displacement field, the governing equations of three sandwich panels are derived using energy variation principle and solved by employing Fourier series and Galerkin approach [23]. Shahbazi et al. [24] introduced a meshfree approach for analysis of isotropic/orthotropic cross-ply laminated plates with symmetric/non-symmetric layers. A series of exponential basis functions, satisfying the partial differential equations were used to approximate the solution on the whole domain [24]. The boundary conditions are enforced through a collocation approach on a set of boundary points [24]. Aifantis strain gradient elasticity theories and Zhang’s two-variable method are used to study elastic bending problems of bilayered microcantilever beams, containing a gradient layer, subjected to a transverse concentrated load [25]. Based on three-dimensional theory of elasticity, axisymmetric static analysis of functionally graded circular and annular plates imbedded in piezoelectric layers is investigated using differential quadrature method [26]. Liew et al. [27] discussed mainly on the developments of element free or meshless methods and their applications in the analysis of composite structures. Rodrigues et al. [28] proposed to use the Murakami’s zig-zag theory for the static and vibration analysis of laminated plates, by local collocation with radial basis functions in a finite differences framework. Stürzenbecher and Hofstetter [29] presented an axiomatic equivalent single layer plate theory for cross-ply laminated composites, which is based on the work of Lekhnitskii and Ren and delivers accurate deformation and stress prognoses at the cost six solution variables.
Deb Nath et al. [30, 31] analyzed the elastic field of different stiffened composite structures analytically, with the verifications of the solution of FE (finite element) and FD (finite difference) methods. Deb Nath et al. [30, 31] studied the effect of fiber orientation on the elastic field of a stiffened orthotropic panel under tension. For the case of combined effects of bending and tension at the ends of the stiffened panel, effects of fiber orientation and material isotropy are not studied analytically till now which are necessary for the optimum design of orthotropic members in any loaded structure because during loading members of any structure experience tension and bending together. The present study presents the effects of fiber orientation and material isotropy on the elastic field of an orthotropic composite panel subjected to a combined loading at its tip. This analysis is carried out analytically using Fourier integral and the displacement potential formulation of orthotropic composite materials. Two types of fiber orientation are considered in the present analysis such as θ=0° and θ=90°. Appendix A illustrates the boundary and governing equations for isotropic panel and Appendix B illustrates the equations to derive the isotropic mechanical properties from orthotropic composite materials. The soundness and reliability of the present solution are shown comparing the present analytical solution at some particular sections with those of finite element predictions.
2. Analytical Model of the Present Problem
The geometry of a stiffened boron/epoxy panel subjected to a combined loading whose fiber orientation θ=90° is shown in Figure 1. The supporting edge of the panel is rigidly fixed, two opposite edges are stiffened, and at the tip of the panel, uniform tension and shear are applied. Mechanical properties of the boron/epoxy composite used are given by Table 1. The applied tensile and shear stress on the tip of the panel is 41.4 MPa.
Properties of composites used to obtain numerical results.
Material
Property
Boron/Epoxy
Fiber
Ef (103 MPa)
414
νf
0.20
Resin
Er (103 MPa)
3.45
νr
0.35
Composite
E1 (103 MPa)
282.9
E2 (103 MPa)
24.2
G12 (103 MPa)
10.4
ν12
0.27
ν21
0.023
Model of an orthotropic composite panel subjected to a combined loading (uniform tension and shear).
3. Displacement Potential Formulation of Orthotropic Materials for Plane Stress Conditions
With reference to a rectangular coordinate system (x,y), the differential equations of equilibrium for the plane stress problems of orthotropic composite materials having fiber orientation θ=90° [32] are as follows:
(1)E1E2E1-ν122E2∂2ux∂x2+(ν12E1E2E1-ν122E2+G12)∂2uy∂x∂y+G12∂2ux∂y2=0,E12E1-ν122E2∂2uy∂y2+(ν12E1E2E1-ν122E2+G12)∂2ux∂x∂y+G12∂2uy∂x2=0.
The governing equation of the orthotropic material with fiber orientation 90° in terms of the function ψ is as follows [31]:
(2)E2G12∂4ψ∂x4+E2(E1-2ν12G12)∂4ψ∂x2∂y2+E1G12∂4ψ∂y4=0.
The corresponding governing equation for isotropic materials can be obtained when the respective conditions [E1=E2=E,ν12=ν21=ν,G12=G=E/2(1+ν)] are substituted in (2) which is as follows:
(3)∂4ψ∂x4+2∂4ψ∂x2∂y2+∂4ψ∂y4=0.
The components of different displacement components with fiber orientation 90° [31] are as follows:
(4)ux(x,y)=∂2ψ∂x∂y,uy(x,y)=-1Z11[E1E2∂2ψ∂x2+G12(E1-ν122E2)∂2ψ∂y2].
The stress-strain relations of orthotropic composite lamina with fiber orientation θ=90° in terms of the potential function ψ are obtained as follows [31]:
(5)σxx(x,y)=E1E2G12Z11[∂3ψ∂x2∂y-ν12∂3ψ∂y3],σyy(x,y)=E1Z11[E2(ν12G12-E1)∂3ψ∂x2∂y-E1G12∂3ψ∂y3],σxy(x,y)=E1E2G12Z11[ν12∂3ψ∂x∂y2-∂3ψ∂x3].
4. Solution of the Present Problem4.1. Solution of the Panel (θ=90°)
A stiffened rectangular boron/epoxy orthotropic composite panel is considered. Its fiber orientation is θ=90°. Its supporting edge is rigidly fixed and the two opposing edges are stiffened. The panel is considered to be of unit thickness, and its configuration with respect to coordinate axes is illustrated in Figure 1. In this case, if the potential function ψ is assumed to be
(6)ψ=∑m=1∞Ymsinαx+Bx3+Cxy2,
where Ym is a function of y only, and α=mπ/a, then Ym has to satisfy the ordinary differential equation
(7)Ym′′′′-(E2G12-2ν12E2E1)α2Ym′′+E2E1α4Ym=0.
The general solution of the differential equation (7) can be given by
(8)Ym=Amem1y+Bmem2y+Cmem3y+Dmem4y,
where
(9)m1,m2=α2[{K1±K12-4E2E1}]1/2,m3,m4=-α2[{K1±K12-4E2E1}]1/2,K1=E2G12-2ν12E2E1,
where Am, Bm, Cm, and Dm are arbitrary constants.
Now, combining (4) to (5) and (6), the expressions of displacement and stress components are obtained as follows:
(10)ux(x,y)=∑m=1∞αYm′cosαx+2Cy,uy(x,y)=1Z11∑m=1∞[-Z22Ym′′+E1E2α2Ym]sinαx+(6BE1E2+2CZ33)x,σxx(x,y)=-E1E2G12Z11∑m=1∞[ν12Ym′′′+α2E1Ym′]sinαx,σyy(x,y)=E1E2Z11∑m=1∞[Ym′α2(E1-ν12G12)E1E2Z11∑m=1∞-G12Ym′′′]sinαx,σxy(x,y)=E1E2G12Z11∑m=1∞[αν12Ym′′+α3Ym]cosαx-6BE1E2G12Z11+2Cν12E1E2G12Z11,
where Z22=G12(E1-ν122E2), Z33=(ν12E2Z11-E12E2)/(E1-ν122E2).
Substituting the different derivatives of Ym in the expressions of the displacement and stress components equations (10), we get
(11)ux(x,y)=∑m=1∞[(m1Amem1y+m2Bmem2y+m3Cmem3y+m4Dmem4y)]αcosαx+2Cy,uy(x,y)=-1Z11∑m=1∞[(m12Z22-α2E1E2)Amem1y+(m22Z22-α2E1E2)Bmem2y+(m32Z22-α2E1E2)Cmem3y+(m42Z22-α2E1E2)Dmem4y]×sinαx+(6BE1E2+2CZ33)x,σxx(x,y)=-E1E2G12Z11×∑m=1∞[(m1α2+m13ν12)Amem1y+(m2α2+m23ν12)Bmem2y+(m3α2+m33ν12)Cmem3y+(m4α2+m43ν12)Dmem4y]sinαx,σyy(x,y)=-E1E2Z11×∑m=1∞[{m1α2(ν12G12-E1)+m13G12}Amem1y+{m2α2(ν12G12-E1)+m23G12}Bmem2y+{m3α2(ν12G12-E1)+m33G12}Cmem3y+{m4α2(ν12G12-E1)+m43G12}Dmem4y]sinαx,σxy(x,y)=E1E2G12Z1∑m=1∞[(α3+αm12ν12)Amem1y+(α3+αm22ν12)Bmem2y+(α3+αm32ν12)Cmem3y+(α3+αm42ν12)Dmem4y]cosαx-6BE1E2G12Z11+2Cν12E1E2G12Z11.
4.2. Boundary Conditions and Its Application
For the present problem, it is observed that the boundary conditions on the two stiffened edges as observed in Figure 1 are
(12)uy(0,y)=0,σxx(0,y)=0,uy(a,y)=0,σxx(a,y)=0,
satisfied automatically.
The associated boundary conditions of the supporting edge are
(13)ux(x,0)=0,uy(x,0)=0.
Now, the axial loading on the right lateral edge of the panel y=b can be expressed mathematically as follows:
(14)σyy(x,b)=P=∑m=1∞Emsinαx,σxy(x,b)=P=E0+∑m=1∞Em¯cosαx,
where Em=(2/a)∫0aPsinαxdx=-4P/mπ, E0=(1/a)∫0aPdx=P; Em¯=(2/a)∫0aPcosαxdx=0.
Substituting the boundary conditions of (13) to (14) into the general expressions of (11), we get the following six equations in terms of the six unknown coefficients:
(15)m1Am+m2Bm+m3Cm+m4Dm=0,(16)(m12Z22-α2E1E2)Am+(m22Z22-α2E1E2)Bm+(m32Z22-α2E1E2)Cm+(m42Z22-α2E1E2)Dm=0,(17)(m1α2Z33+m13G12E1)Amem1b+(m2α2Z33+m23G12E1)Bmem2b+(m3α2Z33+m33G12E1)Cmem3b+(m4α2Z33+m43G12E1)Dmem4b=-EmZ11E1,(18)(α3+αm12ν12)Amem1b+(α3+αm22ν12)Bmem2b+(α3+αm32ν12)Cmem3b+(α3+αm42ν12)Dmem4b=0,(19)-3B+Cν12=PZ112E1E2G12,(20)3BE1E2+CZ33=0.
Solution of the previous algebraic equations (15) to (18) yields the unknown constants Am, Bm, Cm, and Dm.
From (19) and (20), we get the values of B and C as
(21)B=-PZ11Z336G12E1E2(ν12E1E2+Z33),C=PZ112G12(ν12E1E2+Z33).
Once the values of the unknowns are known, they are directly substituted into (11) to obtain the explicit expressions for the different parameters of interest, namely, the two displacement and the three stress components. Following the previous procedure, the panel having the orientation of fibers θ=0° and isotropic material are solved.
5. Results and Discussions5.1. Elastic Analysis of of the Stiffened Panel (θ=90°)
Here, the deformed shape of the orthotropic beam panel subjected to a combined loading having fiber orientation θ=90° under uniform tension and shear at its tip is shown in Figure 2. The effect of fiber orientation and material isotropy on the elastic field of the stiffened orthotropic panel is analyzed. Finally, the present elastic solutions of the supporting edge and the loaded tip are compared with those of finite element predictions.
Deformed shape of the panel (deformation is magnified by multiplying factor of 1000).
5.2. Effect of Fiber Orientation and Material Isotropy
Effects of fiber orientation and material isotropy on the displacement and stress components are analyzed in this paper. Two different fiber orientations, such as θ=0° and 90°, are taken into account for the purpose of the present analysis. Here, a serious attention is taken on the effect of fiber orientation and material isotropy of a boron/epoxy-stiffened panel subjected to a combined loading (tension and shear) at the tip. The aspect ratio of the panel b/a, as observed in Figure 1, is 3. The effective materials properties are given in Table 1. Effects of fiber orientation θ=90° and θ=0° are considered as case 1 and case 2, respectively, and material isotropy is considered in case 3.
The distribution of the lateral displacement component at the section y/b=1.0 for three different cases is shown in Figure 3(a). In this figure, there is observed slight effect on fiber orientations. Over the range 0<x/a<1, for case 1, the lateral displacement with x/a at the section y/b=1 remains same with x/a. For case 2, the lateral displacement with x/a at the section y/b=1 increases. At the neutral axis, x/a=0.5, the lateral displacement for case 1 and case 2 is equal. The lateral displacement at the section y/b=1 of the isotropic panel increases with x/a. The magnitude of lateral displacement component ux of the isotropic panel at the section y/b=1 is the lowest. Figure 3(b) illustrates the effect of fiber orientation and material isotropy on the normalized axial displacement component uy at the section y/b=1.0. At two stiffened edges, the axial displacement component uy is zero, which satisfies the physical boundary conditions of the problem. The axial displacement component uy of the orthotropic panel is the highest for case 2 because in this case, the stiffness of the panel in y-direction is the lowest, and for case 1, the opposite phenomenon occurs due to having the highest stiffness in the same directions of the panels among three cases, and the axial displacement component uy of the isotropic panel remains in between them due to having intermediate stiffness. From Figure 4, it is observed that at the section y/b=1, for case 2, case 1, and case 3, the panel shows the highest, lowest, and intermediate lateral stress σxx, and the variation of the lateral stress occurs due to having different stiffness in x-direction of the panel; for example, for case 2, case 1, and case 3, in x-direction, the panel shows the highest, lowest, and intermediate stiffness. From Figure 5(a), it is observed that the lateral displacement at the section y/b=0.5 with x/a for case 2 remains constant with x/a but it decreases with x/a for case 1. At the neutral axis, x/a=0.5, the lateral displacement at the section y/b=0.5 for case 1 and 2 is same. For case 3, the lateral displacement at the section y/b=0.5 is the lowest in compared to other two cases, 1 and 2. Figure 5(b) illustrates the comparison of the axial displacement components uy at the section y/b=0.5 at different fiber orientation and material isotropy of the panel. The axial displacement component uy at the section y/b=0.5 is the highest for the 90° fiber orientation (case 1) and the lowest for the isotropic material (case 3), and the displacement component of 90° fiber orientation (case 2) remains in between them. From Figure 6(a), it is shown that at the section y/b=0.5, the isotropic panel shows the highest lateral stress, σxx but effect of fiber orientation on the stress distribution is little. The panel whose fiber orientation is 90° shows highest axial stress σyy at the sectiony/b=0.5 as shown in Figure 6(b). From Figure 6(c), it is observed that at the section y/b=0.5, the shear stress is the highest for case 1 and the lowest for case 2. The shear stress for case 3 at the section y/b=0.5 remains in between that of other two cases. Figure 7 shows the effect of fiber orientation and material isotropy on the shear stress distribution at the section x/a=0.0. There is also observed effect of fiber orientation and material isotropy on the shear stress distribution.
Effect of fiber orientation and material isotropy on the normalized displacement components at the section y/b=1.0 of the panel: (a) lateral displacement; (b) axial displacement.
Effect of fiber orientation and material isotropy on the normalized lateral stress component at the section y/b=1.0 of the panel.
Effect of fiber orientation and material isotropy on the normalized displacement components at the section y/b=0.5 of the panel: (a) lateral displacement; (b) axial displacement.
Effect of fiber orientation and material isotropy on the normalized stress components at the section y/b=0.5 of the panel: (a) lateral stress component; (b) axial stress component; (c) shear stress component.
Effect of fiber orientation and material isotropy on the normalized shear stress component at the section x/a=0 of the panel.
6. Comparison with the Available Results
To compare the present solutions with those of FEM (finite element modeling), the panel made of boron/epoxy having fiber orientation θ=90° is solved using the commercial software ANSYS. Rectangular orthotropic plane element is used to mesh the present model. To get higher accuracy of solutions, the numbers of elements are decreased or increased.
Figure 8 shows the comparison of solutions of lateral and axial displacement components at the section y/b=1 by the analytical and FE predictions. The lateral displacement component of the section y/b=1 obtained by FE and analytical method coincide each other as shown in Figure 8(a). Away from the stiffened edge to the neutral axis of the panel, the deviation between the solutions of axial displacement obtained by FEM and analytical method increases, and axial displacement obtained by analytical method is higher than those of FEM as shown in Figure 8(b). From Figure 9, it is observed that the lateral stress of the section y/b=1 obtained by analytical and FEM is almost same in most of the regions of this section. Figure 10 illustrates the comparison of different stress components of the section y/b=0 by analytical and FE methods. From Figures 10(a) and 10(b), it is observed that the deviation between the solution of lateral stress obtained by the analytical and FE methods is significant; for the lateral stress, finite element method gives higher values than that of analytical method, but for the axial stress, the opposite phenomenon occurs. From Figure 10(c), it is clearly observed that the shear of the section y/b=0 obtained by analytical and FE methods is almost same. From the comparative analysis of the solutions between those of analytical and FEM methods, it is observed that in most of the cases, the solutions obtained by the analytical and FE methods quantitatively agree well, but in some cases, the solutions do not agree quantitatively but agree well qualitatively. So, the present analytical solutions are highly reliable for designing such type of structures.
Comparison of solutions of different displacement components at the section y/b=1 of the beam by analytical and FE methods: (a) lateral displacement; (b) axial displacement.
Comparison of solutions of the lateral stress component at the section, y/b=1 by analytical and FE methods.
Comparison of solutions of different stress components at the section y/b=0 by analytical and FE methods: (a) lateral stress; (b) axia stress and (c) shear stress.
7. Conclusions
The effects of fiber orientation and material isotropy on the stress and displacement components of the panels under uniform tension and shear at the tip of the panel are discussed. Effect of stiffeners is strongly observed in the solution. Due to stiffeners, maximum shear stress is observed at the stiffened edges. Significant effect of fiber orientation and material isotropy are observed on the displacement and stress components at different important sections of the panel. To check the reliability and soundness of the present solution, the present analytical solutions of the panel are compared with those of Finite element predictions.
AppendicesA. Solution of the Isotropic Panel
The mathematical background for the solution of isotropic stiffened panel using the present ψ-formulation is summarized in this appendix for ready reference of interested readers. Combining (3) and (6), the governing ordinary differential equation for the isotropic panel problem is obtained, which is
(A.1)Ym′′′′-2α2Ym′′+α4Ym=0.
The general solution of (A.1) is as follows:(A.2)Ym=Amcoshαy+Bmαysinhαy+Cmsinhαy+Dmαycoshαy.
The corresponding expressions for the displacement and stress components in terms of the four arbitrary constants are as follows:(A.3)ux(x,y)=∑m=1∞α2[(Am+Bm)sinhαy+(Cm+Dm)coshαy+Bmαycoshαy+Dmαysinhαy(Am+Bmh)]cosαx+2Cy,(A.4)uy(x,y)=∑m=1∞α2[(Am-2tBm)coshαy+(Cm-2tDm)sinhαy+Bmαysinhαy+Dmαycoshαy]sinαx+{6B+C(1-ν)}E2x,(A.5)σxx(x,y)=-∑m=1∞Eα3(1+ν)[(Am+t′Bm)sinhαy+(Cm+t′Dm)coshαy+Bmαycoshαy+Dmαysinhαy(Am+t′Bm)]sinαx,(A.6)σyy(x,y)=-∑m=1∞Eα3(1+ν)[(Am-tBm)sinhαy∑m=1∞Eα3(1+ν)+(Cm-tDm)coshαy∑m=1∞Eα3(1+ν)+Bmαycoshαy∑m=1∞Eα3(1+ν)+Dmαysinhαy]sinαx,(A.7)σxy(x,y)=∑m=1∞Eα3(1+ν)2[{Am(1+ν)+2νBm}coshαy∑m=1∞Eα3(1+ν)2+Bm(1+ν)αysinhαy∑m=1∞Eα3(1+ν)2+{Cm(1+ν)+2νDm}sinhαy∑m=1∞Eα3(1+ν)2+Dm(1+ν)αycoshαy]×cosαx+E(1+ν)2(-6B+2Cν),
where t=(1-ν)/(1+ν) and t′=(1+3ν)/(1+ν).
Now, as before, solving (A.3), (A.4), (A.6), and (A.7) with the appropriate boundary conditions given by (12)–(14), the values of the four constants are evaluated.
B. Effect of Material Isotropy on the Solution
It can be noted here that the corresponding isotropic material properties of the boron/epoxy orthotropic composite material are obtained by solving the following set of equations [33]:
(B.1)VmKm-K3Km+4G+VfKf-K3Kf+4G=0,VmGm-GGm+β+VfGf-GGf+β=0,E=9KG3K+G,G=E2(1+ν),
where β=G(9K+8G)/(6K+12G).
The solutions of the previous equations are carried out for a composite composition of 60% boron fiber and 40% epoxy resin. The supporting mathematical treatments required to solve the isotropic stiffened panel problem are given in Appendix A.
NomenclatureE1
:Elastic modulus of the material in x-direction
E2:
Elastic modulus of the material in y-direction
ν12:
Major Poisson’s ratio
ν21:
Minor Poisson’s ratio
θ:
Fiber orientation
ux:
Displacement components in the x direction
uy:
Displacement components in the y direction
σxx:
Normal stress components in the x direction
σyy:
Normal stress components in the y direction
σxy:
Shearing stress component in the xy plane
ψ:
Displacement potential function
a:
Width of the panel
b:
Length of the panel
σyy0:
Uniformly distributed axial loading on the panel
σxy0:
Uniform shear stress
E:
Elastic modulus of isotropic material
ν:
Poison’s ratio of isotropic material.
LekhnitskiiS. G.1963San Francisco, Calif, USAHolden-DayLekhnitskiiS. G.1968New York, NY, USAGordon and BreachStrohA. N.Dislocations and cracks in anisotropic elasticity1958330625646StrohA. N.Steady-state problems in anisotropic elasticity19624177103TingT. C. T.1996New York, NY, USAOxford University PressBarnettD. M.KirchnerH. O. K.A proof of the equivalence of the Stroh and Lekhnitskii sextic equations for plane anisotropic elastostatics19977612312392-s2.0-0031184684ConwayH. D.ChowL.MorganG. W.Analysis of deep beams195118163172ConwayL. Chow H. D.WinterG.Stresses in deep beams19522557HorganC. O.KnowelsJ. K.Recent developments concerning Saint Venant’s principle198323179269ParkerD. F.The role of Saint-Venant’s solutions in rod and beam theories19794648618662-s2.0-0018702114HardyS. J.PipelzadehM. K.Static analysis of short beams1991261529Krishna MurtyA. V.Towards a consistent beam theory19842268118162-s2.0-0021440858DurelliA. J.RanganayakammaB.On the use of photoelasticity and some numerical methods814Photomechanics and Speckle Metrology198718Proceedings of the SPIEDurelliA. J.RanganayakammaB.Parametric solution of stresses in beams198911524014142-s2.0-0024607275TimoshenkoS.GoodierV. N.19793rd New York, NY, USAMcGraw- HillUddinM. W.1966Ottawa, CanadaCarleton UniversityChapelR.SmithH. W.Finite difference solutions for plane stresses1968611561157ConwayH. D.IthacaN. Y.Some problems of orthotropic plane stress195352-A-47276KulikovG. M.PlotnikovaS. V.Exact 3D stress analysis of laminated composite plates by sampling surface method20129436543663TahaniM.AndakhshidehA.Interlaminar stresses in thick rectangular laminated plates with arbitrary laminations and boundary conditions under transverse loads20129417931804ZhangC.HoaS. V.A limit-based approach to the stress analysis of cylindrically orthotropic composite cylinders (0/90) subjected to pure bending20129426102618SinghJ.ShuklaK. K.Nonlinear flexural analysis of laminated composite plates using RBF based meshless method20129417141720HeL.-S ChengY.LiuJ.Precise bending stress analysis of corrugated-core, honey comb-core and X-core sandwich panels20129416561668ShahbaziM.BoroomandB.SoghratiS.A mesh-free method using exponential basis functions for laminates modeled by CLPT, FSDT and TSDT - Part I: formulation201193311231192-s2.0-7996035405510.1016/j.compstruct.2011.06.023ZhangN. H.MengW. L.AifantisE. C.Elastic bending analysis of bilayered beams containing a gradient layer by an alternative two-variable method20119331303139AlibeiglooA.SimintanV.Elasticity solution of functionally graded circular and annular plates integrated with Sensor and actuator layers using differential quadrature20119324732486LiewK. M.ZhaoX.FerreiraA. J. M.A review of meshless methods and functionally graded plates and shells20119320312041RodriguesJ. D.RoqueC. M. C.FerreiraA. J. M.CarreraE.CinefraM.Radial basis functions-finite differences collocation and a unified formulation for bending, vibration and buckling analysis of laminated plates, according to Murakami's zig-zag theory2011937161316202-s2.0-7995548545010.1016/j.compstruct.2011.01.009StürzenbecherR.HofstetterK.Bending of cross-ply laminated composites: an accurate and efficient plate theory based upon models of Lekhnitskii and Ren2011933107810882-s2.0-7995489540710.1016/j.compstruct.2010.09.020Deb NathS. K.AfsarA. M.Analysis of the effect of fiber orientation on the elastic field in a stiffened orthotropic panel under uniform tension using displacement potential approach20091643003072-s2.0-7044958908810.1080/15376490802664281Deb NathS. K.AhmedS. R.Investigation of elastic field of a short orthotopic composite column by using finite-difference technique20082228116111692-s2.0-5734910292910.1243/09544100JAERO374JonesR. M.1975McGraw-HillNanC. W.YuanR. Z.ZhangL. M.HoltJ. B.The physics of metal/ceramic functionally gradient materials199334Westerville, Ohio, USAAmerican Ceramic society7582