Medical, Pharma, Engineering, Science, Technology and Business

**Kostopoulos V ^{*}, Kotzakolios T and Vlachos DE**

Applied Mechanics Laboratory, Department of Mechanical Engineering and Aeronautics, University of Patras, Patras University Campus, GR-26500 Patras, Greece

- *Corresponding Author:
- Kostopoulos V

Applied Mechanics Laboratory, Department of Mechanical Engineering and Aeronautics

University of Patras, Patras University Campus, GR-26500 Patras, Greece

**E-mail:**kostopoulos@mech.upatras.gr

**Received Date:** February 04, 2017;** Accepted Date:** March 06, 2017;** Published Date:** March 10, 2017

**Citation:** Kostopoulos V, Kotzakolios T, Vlachos DE (2017) The Buckling Response of Lattice Fuselage Structures: Validation of Finite Element Models by Using Smeared Unit Cell Analytical Methodology. J Aeronaut Aerospace Eng 6:185. doi: 10.4172/2168-9792.1000185

**Copyright:** © 2017 Kostopoulos V, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Visit for more related articles at** Journal of Aeronautics & Aerospace Engineering

Composite lattice structures are shells that are reinforced by unidirectional helical and hoop ribs. Their main advantage over contemporary composite structures is their superior stiffness to mass ratio. However, their application in industry is still limited. In this paper, the lattice structure concept was applied for the case of a small business aircraft. Emphasis here is given at the initial stages of the design. More specifically, the buckling modes for bending loads were calculated by utilizing a continum unit cell model which was correlated with finite element models for a cylindrical small fuselage structure and two scaled down versions.

Composite lattice strictures are composed of helical and hoop ribs that reinforce a shell. They are characterized by high specific strength and stiffness, making them promising for lightweight applications. The history and development of this kind of structures can be found in literature, mostly published by Vasiliev et al. [1-5].

Lattice structures are using reinforcing ribs in a regular pattern, which allows to be analyzed by smearing the ribs over the skin surface. Thus, the lattice structure can be analyzed as a continuous layer with calculated effective stiffness. Stress and strain equations that are based on theory of orthotropic shells can be used. Such continuum models are published by Vasiliev [6] and Vasiliev and Morozov [7]. Strength and buckling analysis of cylindrical lattice shells based on different continuum models are described in articles from Slinchenko and Verijenko [8], Totaro and Gurdal [9], Buragohain and Velmurugan [10], Paschero and Hyer [11], Totaro [12,13], and Zheng et al. [14].

Additionally, composite lattice structures have been analyzed by using finite element analysis. Results from FE models can be found in articles published by Hou et al., Zhang et al., Frulloni et al., Fan et al., Morozov et al. and Azarov et al. [15-20].

However, in most applications, both methodologies are used, where the basic specifications of the structure are calculated by employing the continuum models. The parameters are then further refined by using the finite element models. Correlation between these methodologies has been done successfully by Azarov [20].

Despite the aforementioned developments in designing composite lattice structures, their application is still limited in some spacecraft applications [3,21].

In this work, the case of a fuselage section of a small business aircraft is examined. The fuselage section has an outer diameter of 1.8 m and a length of 4 m. The fuselage section was assumed to be subjected to buckling loads and its buckling limit was calculated by a smeared unit cell method and finite element analysis. Furthermore, a scaling down approach was followed for estimating the buckling loads for scaled down prototypes, namely 0.5 and 1 m in diameter. The scaling down approach is useful since many times full scale testing is prohibited for this kind of structures since they are integrated with no assembly required, making testing of single subcomponents not feasible. The correlation between the smeared unit cell theory and finite element analysis was very good.

In developing the analytical model, a unit cell of the stiffener structure has to be defined first. The unit cell is chosen such that the whole grid structure can be reproduced by repetition of this unit cell.

The equivalent stiffness parameters of this unit are determined and then applied to the whole cylinder panel. Validation comes from the generation of the panel by repetition of the cell.

1. The transverse modulus of the unidirectional stiffeners is much lower than the longitudinal modulus, and cross sectional dimensions are also very small compared to the length dimension, therefore the stiffeners are assumed to support axial loads only.

2. The strain is uniform across the cross sectional area of the stiffeners. Hence a uniform stress distribution is assumed.

3. Load is transferred through shear forces between the stiffeners and the shell.

Strains used as the matching condition of the stiffener and the shell (inner surface of the shell as interface) is given by:

(1)

(2)

(3)

Strain relationship along the stiffeners’ directions is given by:

(4)

Where c=cos(ϕ), s=sin(ϕ) and ϕ is the stiffener orientation angle from vertical direction.

In this specific case of study, ϕ±26° was selected for the helical ribs and ϕ=90° was selected for the hoop ribs (**Figure 1**).

Resultant Forces of the force diagram is as shown below (**Figure 2**):

(5)

(6)

(7)

Summing up the forces of each opposite sides in both directions, F_{x} and F_{θ} are obtained:

(8)

(9)

(10)

Resultant Forces on the unit cell is given by:

(11)

(12)

(13)

Following the same procedure as the force analysis on the unit cell, resultant moments are computed with respect to the moment diagram (**Figure 3**):

(14)

(15)

(16)

Resultant Moments on the unit cell:

(17)

(18)

(19)

The functions of the mid plane strains of the shell then become:

(20)

The superscript‘s’ denotes the force and moment contributions of the stiffener.

By applying the Rule of Mixtures, the forces and the moments are superimposed according to the volume fractions of the stiffeners (V_{S}) and the shell (V_{Sh}).

(21)

The resultant stiffness parameters obtained from the above equation are the equivalent stiffness parameters of the whole panel.

For the calculation of the critical buckling load, the Ritz method (Przemieniecki, 1968) was used which is based on the minimization of the total energy Π of the cylinder. The total energy of the cylindrical structure is the sum of the strain energy U and the external force work V.

(22)

The strain energy of an orthotropic cylinder with length L is given by the following equation:

(23)

U is dependent from the total stiffness matrix of the unit cell, radius r of the cylinder and axial, circumferential and radial displacements u, v, w respectively.

The dynamic energy term V, due to the external work for length L is given by:

(24)

The displacement field u, v and w are defined by kinematically admissible functions and they are approximated by a double Fourier series with clamped boundary condition:

(25)

(26)

(27)

m, n = 1,2,3…

The total energy expression is a function of the stiffness matrix elements of the equivalent laminate and the unknown displacement field coefficients A_{mn}, B_{mn} and D_{mn}. For the equilibrium to be stable, the total potential energy of the system must be minimum. This can be satisfied by finding the first derivative of the total potential energy with respect to the unknown constants A_{mn}, B_{mn} and D_{mn} and equating to zero. This results in an eigenvalue problem. The resulting Equation is then solved for the unknown in-plane load Nx.

The specifications for the lattice with 1.8 m diameter are summarized in **Table 1**, whereas the material properties in **Table 2** (**Tables 1** and **2**).

Parameter | [90/±26º] ribs orientation |
---|---|

Material | HTS5631 |

Skin Thickness, mm | 1.4 |

No of Reels (each side) | 20 |

No. of Helical ribs | 40 |

No of Circum. Hoops | 6 |

Area of Rib Section, mm2 | 134 |

Mass, Kg | 48.8 |

**Table 1:** Specifications for lattice fuselage structures.

Property | Value |
---|---|

Young's modulus E1 | 128.9GPa |

Young's modulus E2 | 10.4 GPa |

Shear modulus G12 | 4.11 GPa |

v12 | 0.34 |

Thickness t | 0.35 mm |

**Table 2:** HTS-5631 material properties.

By using the above requirements and theory, a MATLAB code was written where the buckling load was calculated and verified by finite element analysis using MSC NASTRAN. The Finite element model consisted of 4400 beam elements for the hoop and helical ribs and 12000 shell elements for the composite skin. The one side of the model was clamped whereas the other one was subjected to a compressive load which was distributed to the nodes via multi-point constraints (MPCs). The first buckling mode from the finite element model of the full scale fuselage is shown in **Figure 3**, whereas the comparison with the unit cell analysis is shown in **Table 3** (**Figures 2** and **3**).

Calculation method | Full Scale prototype 1.8m in diameter, 2m in length |
---|---|

Buckling load calculated by FEM | 467.9 kN |

Buckling load calculated analytically | 463.9 kN |

**Table 3:** Computation of full scale fuselage buckling load by unit cell method and finite element analysis.

In order to scale down the wafer to the desired main dimensions, the analysis was based on the of buckling load/mass ratio (P_{cr}/M) of the scaled down structure versus the cross-sectional area (A_{sc}) of the ribs. Moreover, the number of reels was decreased with respect to manufacturing restrictions, such as keeping the spacing between the reels the same. Finally, the skin thickness remained constant, since further scaling down its thickness below the four layers with total thickness of 1.4 mm was not possible because it could lead to an unsymmetrical fuselage skin layup. Given the above descriptions, a parametric analysis was done with the unit cell theory, in order to extract the required ribs cross section for two scaled down fuselage geometries with 1 m and 0.5 m in diameter. The results for both scaled down structures are summarized in **Figures 4** and **5**, where the ratio of critical buckling load to mass (P_{cr}/M) against ribs cross section is shown. The results were used with final goal to achieve the same P_{cr}/M ratio with the full scale fuselage (**Figures 4 **and **5**).

The final step was the correlation of the unit cell models of the scaled down structures with finite element models. **Figure 6** depicts the first buckling mode for the 1 m diameter prototype whereas **Figure 7** shows the first buckling mode of the 0.5 m diameter scaled down structure. **Table 4** shows the comparison of the critical buckling loads by finite element models and unit cell method (**Figures 6** and **7**).

Scaled down prototype 1m in diameter, 1m in length | Scaled down prototype 0.5m in diameter, 1m in length | |
---|---|---|

Buckling load calculated by FEM | 263.9 kN | 136.5 kN |

Buckling load calculated analytically | 260.6 kN | 132.1 kN |

**Table 4:** The comparison of the critical buckling loads by finite element models and unit cell method.

A smeared unit cell analytical solution has been developed for the case of a small business aircraft fuselage section, which was also correlated with finite element analysis approach. Afterwards, a scaling down approach was followed for two scaled down fuselage structures. By keeping the skin thickness the same and reducing the full scale fuselage diameter, the ribs cross section and distance were calculated in order to achieve the same buckling load to mass ratio. All models have shown an excellent correlation between the unit cell approach and the finite element analysis.

Moreover, a close to linear relationship was found between the buckling load and the diameter of the structures, enabling testing and validation of computational tools in scaled down prototypes.

- Vasiliev VV, Barynin VA, Rasin AF (2001)Anisogrid lattice structures – Survey of development and application. Compos Struct54:361–370.
- Vasiliev VV, Razin AF (2006)Anisogrid composite lattice structures for spacecraft and aircraft applications. Compos Struct 76:182–189.
- Vasiliev VV, Barynin VA, Rasin AF, Petrokovskii SA, Khalimanovich VI (2009)Anisogrid composite lattice structures – development and space applications. Compos Nanostruct3:38–50.
- Vasiliev VV, Barynin VA, Razin AF (2012)Anisogrid composite lattice structures – development and aerospace applications. Compos Struct 94:1117–1127.
- Vasiliev VV, Razin AF, Nikityuk VA (2014) Development of geodesic composite fuselage structure. Int Rev AerospEng7(1):48–54.
- Vasiliev VV (1993) Mechanics of composite structures. Washington: Taylor & Francis.
- Vasiliev VV, Morozov EV (2013) Advanced mechanics of composite materials and structural elements. 3rd ed. Elsevier, Amsterdam.
- Slinchenko D, Verijenko VE (2001) Structural analysis of composite lattice shells of revolution on the basis of smearing stiffness. Compos Struct 54:341–348.
- Totaro G, Gurdal Z (2009)Optimal design of composite lattice shell structures for aerospace applications. Aerosp Sci Technol 13:157–164.
- Buragohain M, Velmurugan R (2009) Buckling analysis of composite hexagonal lattice cylindrical shell using smeared stiffener model. Defence Sci J50(3):230–238.
- Paschero M, Hyer MW (2009) Axial buckling of an orthotropic circular cylinder: Applications to orthogrid concept. Int J Solids Struct46:2151–2171.
- Totaro G (2013) Local buckling modelling of isogrid and anisogrid lattice cylindrical shells with triangular cells. Compos Struct 94:446–452.
- Totaro G (2013)Local buckling modelling of isogrid and anisogrid lattice cylindrical shells with hexagonal cells. Compos Struct 95:403–410.
- Zheng Q, Ju S, Jiang D (2014)Anisotropic mechanical properties of diamond lattice composites structures. Compos Struct 109:23–30.
- Hou A, Gramoll K (1998) Compressive strength of composite latticed structures. JReinforc Plastics Compos17:462–483.
- Zhang Y, Xue Z, Chen L, Fang D (2009) Deformation and failure mechanisms of lattice cylindrical shells under axial loading. Int J Mech Sci 51:213–221.
- Frulloni E, Kenny JM, Conti P, Torre L (2007) Experimental study and finite element analysis of the elastic instability of composite lattice structures for aeronautic applications. Compos Struct 78:519–528.
- Fan H, Jin F, Fang D (2009) Uniaxial local buckling strength of periodic lattice composites. Mater Des 30:4136–4145.
- Morozov EV, Lopatin AV, Nesterov VA (2011) Finite-element modelling and buckling analysis of anisogrid composite lattice cylindrical shells. Compos Struct93:308–323.
- Azarov AV (2012) Continuum and discrete models of lattice composite cylindrical shells. Mech Compos Mater Struct 18:121–130.
- Lopatin AV, Morozov EV, Shatov AV (2015) Deformation of a cantilever composite anisogrid lattice cylindrical shell loaded by transverse inertia forces, Composite Structures 129: 27–35.

Select your language of interest to view the total content in your interested language

- ACS Nano
- AI Algorithm
- Advanced Materials
- Advantages of E-Business
- Aerodynamic Flow Control
- Aerodynamics
- Aeroelasticity
- Aerospace Dynamics
- Aerospace Engineering techniques
- Air Safety
- Aircraft
- Aircraft Flight Mechanics
- Applied Engineering
- Applied Statistics
- Applied Theories on Machines
- Architect
- Architectural Drawing
- Architectural Engineering
- Artificial Intelligence
- Astrodynamics
- Automated Mining
- Automation
- Automation Devices
- Automobile Engineer
- Automotive Engineering
- Automotive Industry
- Aviation
- Avionics
- Behavior-based systems
- Biological Engineering
- Biomechanics
- Biomechanics and Robotics
- Brittle Materials
- Building design
- Business & Management
- Carbon Nanotubes
- Ceramics Engineering
- Cognitive Systems Engineering
- Commercialization of New Techniques
- Composite Materials
- Computational Fluid Dynamics
- Computational Neuroscience
- Computer Simulation
- Concrete
- Construction
- Construction Engineering
- Construction Estimating Software
- Control System
- Cosmoparticle Physics
- Design and Microfabrication
- Detectors and Optical Sensors
- Diesel Engine
- Digital Image Processing
- Dynamic Information
- Dynamical System
- Electronic Material Development
- Engine
- Engine Performance
- Engineering Design
- Engineering Drawing
- Entrepreneurship
- External Force
- Fabric Formwork
- Finance management
- Flight Dynamics
- Fluid Bodies
- Flying Wheel
- Fuel Economy
- Fusion Plasmas
- Fuzzy Logic
- Health care management
- Helmet
- Human-Machine-Interfaces
- Hydraulic Engineering
- Ignititon System
- Industrial Crystallization
- Industrial Engineering
- Industrial Robotics
- Intelligent Robotics
- Interior Design
- Interior Designing
- Landscape Architecture
- Logistics
- Luminescence
- Machine
- Management Cybernetics
- Manufacturing system
- Materials Engineering
- Materials Management
- Mathematical Model
- Mechanical Properties
- Mechanical Systems
- Mechanism
- Mechatronics
- Medical Device
- Medical Robotics
- Metallic Materials (Ferrous & Nonferrous)
- Mobile Device
- Mobile Repots
- Modular Architecture
- Molecular Electronics
- Nano Composites
- Nano Materials
- Nano Particles
- Nano Structures
- Neural Networks
- Neurorobotics
- Nonlinear Dynamics
- Nuclear Physics
- Operations Research
- Parallel Programming
- Polymeric Materials
- Porous Materials
- Process Engineering
- Product Quality
- Production and Operations Management
- Quantum Gravity
- Reliability engineering
- Robotic Rehabilitation
- Robotics Methods
- Rocketry
- Semiconductors
- Simulation
- Social Robots
- Sociology of Architecture
- Soil Physics
- Solid State Plasmas
- Space
- Space Plasmas
- Spark Ignition
- Splitting Method
- Stellar Spectroscopy
- Stochastic control
- Supply Chain Management
- Suprnova Remnants
- Surface Physics
- Technologies Management
- Telerobotics
- Thermodynamics Methods
- Unmanned-Vehicles
- Urban Design
- Urban Planner
- swarm intelligence and robotics

- Journal of Material Sciences & Engineering
- Journal of Electrical Engineering and Electronic Technology
- Journal of Architectural Engineering Technology
- Industrial Engineering & Management
- Advances in Automobile Engineering
- Advances in Robotics & Automation
- Journal of Astrophysics & Aerospace Technology
- Journal of Applied Mechanical Engineering

- Total views:
**429** - [From(publication date):

March-2017 - Aug 18, 2017] - Breakdown by view type
- HTML page views :
**366** - PDF downloads :
**63**

Peer Reviewed Journals

International Conferences 2017-18