(Tag: Visual edit) |
m (Marherna moved page Review 202715644752 to Montes Salmeron 2024a) |
||
(4 intermediate revisions by one other user not shown) | |||
Line 19: | Line 19: | ||
Durante años Airbus ha desarrollado soluciones de pandeo, desde las típicas soluciones cerradas a complejos desarrollos energéticos. Este artículo resume la evolución de dichos métodos con especial énfasis en los desarrollos energéticos y en las prácticas seguidas por Airbus para su desarrollo. Finalmente se incluye un resumen de potenciales mejoras y desarrollos posibles para el futuro. | Durante años Airbus ha desarrollado soluciones de pandeo, desde las típicas soluciones cerradas a complejos desarrollos energéticos. Este artículo resume la evolución de dichos métodos con especial énfasis en los desarrollos energéticos y en las prácticas seguidas por Airbus para su desarrollo. Finalmente se incluye un resumen de potenciales mejoras y desarrollos posibles para el futuro. | ||
+ | |||
+ | '''Palabras clave: '''Métodos analíticos, Rayleigh-Ritz, pandeo, métodos de elementos finitos, MEF | ||
'''ABSTRACT:''' During the last 20 years computational progress has significantly increased the capacity to determine the behavior of structures under complex loading conditions through finite element models analysis. This computational power can be used to build complex detailed finite element models but also to resolve more complex equations which were not possible to be coded before. | '''ABSTRACT:''' During the last 20 years computational progress has significantly increased the capacity to determine the behavior of structures under complex loading conditions through finite element models analysis. This computational power can be used to build complex detailed finite element models but also to resolve more complex equations which were not possible to be coded before. | ||
Line 262: | Line 264: | ||
− | Where Z<sub>n</sub> is a factor depending on the number of waves: | + | Where Z<sub>n</sub> is a factor depending on the number of waves (figure 4): |
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | ||
Line 277: | Line 279: | ||
Therefore, the problem should be adapted to take into account the curvature effect. This increases the complexity of the problem adding two new set of unknowns to the problem. Typically for flat panels it was enough to know the out of plane displacements (''w(x,y)''). Within a curve panel, the axial (''u'') and transverse displacements (''v'') enter now into the problem | Therefore, the problem should be adapted to take into account the curvature effect. This increases the complexity of the problem adding two new set of unknowns to the problem. Typically for flat panels it was enough to know the out of plane displacements (''w(x,y)''). Within a curve panel, the axial (''u'') and transverse displacements (''v'') enter now into the problem | ||
− | definition.[[Image:Draft_Montes Salmeron_393066637-image5.png|center| | + | definition.[[Image:Draft_Montes Salmeron_393066637-image5.png|center|612x612px]]<div class="center" style="width: auto; margin-left: auto; margin-right: auto;"> |
Figure 5. Buckling onset load depending on panel curvature radius</div> | Figure 5. Buckling onset load depending on panel curvature radius</div> | ||
Line 352: | Line 354: | ||
* Comparison with real test data | * Comparison with real test data | ||
− | Obviously non exhaustive validation can be presented in this brief paper but different illustrative examples are enclosed in the following figures, given the adequate confidence that the problem solution is well captured and therefore applicable for sizing and certification: | + | Obviously non exhaustive validation can be presented in this brief paper but different illustrative examples are enclosed in the following figures, from figure 7 to figure 10, given the adequate confidence that the problem solution is well captured and therefore applicable for sizing and certification: |
{| style="width: 100%;border-collapse: collapse;" | {| style="width: 100%;border-collapse: collapse;" | ||
|- | |- | ||
− | | style="border: 1pt solid black;text-align: center;vertical-align: top;width: 50%;" |[[Image:Draft_Montes Salmeron_393066637-image7.png| | + | | style="border: 1pt solid black;text-align: center;vertical-align: top;width: 50%;" |[[Image:Draft_Montes Salmeron_393066637-image7.png|593x593px]] |
− | | style="border: 1pt solid black;text-align: center;vertical-align: top;width: 50%;" |[[Image:Draft_Montes Salmeron_393066637-image8.png| | + | | style="border: 1pt solid black;text-align: center;vertical-align: top;width: 50%;" |[[Image:Draft_Montes Salmeron_393066637-image8.png|728x728px]][[Image:Draft_Montes Salmeron_393066637-image9.png|center|405x405px]] |
|} | |} | ||
Line 365: | Line 367: | ||
{| style="width: 100%;border-collapse: collapse;" | {| style="width: 100%;border-collapse: collapse;" | ||
|- | |- | ||
− | | style="border: 1pt solid black;vertical-align: top;width: 52%;"|[[Image:Draft_Montes Salmeron_393066637-image10.png| | + | | style="border: 1pt solid black;vertical-align: top;width: 52%;"|[[Image:Draft_Montes Salmeron_393066637-image10.png|660x660px]] |
− | | style="border: 1pt solid black;vertical-align: top;width: 47%;"|[[Image:Draft_Montes Salmeron_393066637-image11.png| | + | | style="border: 1pt solid black;vertical-align: top;width: 47%;"|[[Image:Draft_Montes Salmeron_393066637-image11.png|602x602px]] |
|} | |} | ||
Line 387: | Line 389: | ||
===3.6. Implementation=== | ===3.6. Implementation=== | ||
− | The mathematical development of the solutions presented here is really complex and its implementation is not straight forward for an efficient and optimum CPU management | + | The mathematical development of the solutions presented here is really complex and its implementation is not straight forward for an efficient and optimum CPU management (see figure 11 for examples of different performances on computational time vs operations) |
For example, a cubic running time function will be the result of a naive multiplication of two nxn matrices. Matrix chain multiplication with brute force can be time exponential [5]. | For example, a cubic running time function will be the result of a naive multiplication of two nxn matrices. Matrix chain multiplication with brute force can be time exponential [5]. |
RESUMEN: En los últimos 20 años el poder computacional se ha incrementado exponencialmente lo que ha dado lugar a que se puedan construir modelos de elementos finitos muy detallados que representan de manera fiable la estructura que es objeto de dimensionado. Sin embargo, esta potencia de cálculo también se puede usar para resolver ecuaciones complejas que hasta este momento sólo podían resolverse basándose en simplificaciones.
Calcular el pandeo de una estructura es un factor clave si se quiere que esta sea óptima. Es esencial determinar si este pandeo es inestable, conduciendo a un colapso estructural, o si por el contrario existe capacidad de redistribución y post-pandeo. Diseñar una estructura con capacidad de redistribución de carga optimiza el peso incrementando la seguridad y, por tanto es el principal objetivo del dimensionado de algunos componentes estructurales tales como cajones de torsión o revestimiento de fuselajes.
Durante años Airbus ha desarrollado soluciones de pandeo, desde las típicas soluciones cerradas a complejos desarrollos energéticos. Este artículo resume la evolución de dichos métodos con especial énfasis en los desarrollos energéticos y en las prácticas seguidas por Airbus para su desarrollo. Finalmente se incluye un resumen de potenciales mejoras y desarrollos posibles para el futuro.
Palabras clave: Métodos analíticos, Rayleigh-Ritz, pandeo, métodos de elementos finitos, MEF
ABSTRACT: During the last 20 years computational progress has significantly increased the capacity to determine the behavior of structures under complex loading conditions through finite element models analysis. This computational power can be used to build complex detailed finite element models but also to resolve more complex equations which were not possible to be coded before.
Accurate buckling prediction is still a key factor on optimized structural designs. The accurate determination of the onset of buckling is absolutely essential to determine the initial point of instability and the subsequent post-buckling capability. Determining if the buckling will lead to immediate failure or the loading will be able to be redistributed afterwards is a major point to be considered. Designing a structural component with loading redistribution capability optimize the weight and increase the accuracy, so is a key objective of designers on structural components like torque boxes or fuselage skin.
During years Airbus has been focused on developing accurate buckling predictions from close form solutions to complex energy methods development. This paper aims to summarize the evolution of such methods with special emphasis in energy methods and the best practices followed in Airbus for development. On top a summary of potential room for improvement and further evolution will be proposed.
Keywords: Analytical method, Rayleigh-Ritz, buckling, finite element method, FEM
There are typically two ways to approach the problem of the panel stability analysis. One solution will be the use of Finite Element (FE) solution. There are several well-known commercial FE solvers on use in the industry, NASTRAN, ABAQUS, etc. Second solution is the semi-analytical approach development based on energy methods (EM).
Both approaches are complementary and can be used at different maturity status of the design based on following assessments:
In conclusion, Airbus have considered that the development of both approaches leads to the most efficient solution. In this paper we will be limited to the internal developments for energy method approach (analytical solutions).
During the last 40 years the computational power has been developed overriding all the different limits law pronounce over such period of time [1], see figure 1:
This has permitted that methodologies such as described in [2] have become a reality as a computer program, evolving from tabulated graphs, such as can be seen in figure 2, to current computer programs that can be executed as many times as the user requires, calculating each time the specific equations linked to the exact geometry and loads given. The effort, interpretation, interpolations that an engineer should do are reduced to the minimum, increasing the engineering efficiency.
Energy method is based on the general theorem of the equilibrium of mechanical systems: the potential energy of the system reaches an extreme value at its equilibrium position.
Two states of the panel are considered: the equilibrium state that is immediately previous to the buckling one, (at which the middle surface of the panel remains straight under the action of the applied loads); and the equilibrium state that is immediately subsequent to the buckling one (at which the middle surface of the panels is bent). Calling ΠI the potential energy of the first panel state and ΠΙΙ the potential energy of the second one, the relationship between both of them is given by:
|
(1) |
According to the potential energy theorem, the state I shows a stable equilibrium if the next condition is fulfilled:
ΠI < ΠII (δΠ > 0) | (2) |
The equilibrium will be unstable if:
ΠI > ΠII (δΠ < 0) | (3) |
Critical loads are those for which the state I is no longer stable, i.e:
|
(4) |
The total energy of the system can be split in: the strain energy (U) and the potential energy of the applied loads (V).
To illustrate, the strain energy (U) of a simple flat panel can be expressed as follow:
|
(5) |
Where [A], [B] and [D] are the laminate stiffness matrices according to the classical laminate theory.
The potential energy of the applied loads (V) will be:
|
(6) |
Where ʎ indicates the factor applied to the loads which makes the panel instable.
Integrating and equating equations (5) and (6) we can determine the deformed shape function (w(x,y)), and we can obtain the critical buckling load. However this solution cannot be solved mathematically and the energy method propose an approximated solution assuming a deformed shape function as follows:
|
(7) |
|
(8) |
|
(9) |
The Rayleigh-Ritz approach implies defining an approximate solution representing a displacement field that automatically complies with the boundary conditions.
For example, for a simply supported edge a simple Fourier series development can be proposed:
|
(10) |
The greater the number of terms is (nx and ny), the more accurate the solution will be. These functions satisfy the geometrical boundary conditions (w = M = 0 null displacements and moments at the four edges).
Nevertheless, different boundary conditions are sometimes needed, one free edge for example. For those cases typical Fourier series are no more adequate [10]. To overcome this problem, and define various boundary conditions, these have to be imposed as a set of nc linear equations that must be subsequently satisfied by the solution:
|
(11) |
The minimization problem becomes a conditional optimization that will be solved by means of the well-known Lagrange multipliers method.
The Lagrange multiplier method consists in multiplying each of these constraint equations with arbitrary parameters ʎi and adding the result to the minimum potential energy equation, obtaining a new objective function Π2
|
(12) |
The problem is now more complex, new equations are added but the flexibility and accuracy to represent different problems is significantly improved.
With this mathematical approach it is possible to find analytical solutions to buckling problems as shown in figure 3. An important remark is that with this solution only local buckling problems can be solved. For global stability another solution must be used, FE solvers.
In order to represent properly the stiffeners attached to the panels, its contribution need to be added to the strain and potential energy of the system:
|
(13) |
The problem of this idealization is that an additional correction is needed on the flexural inertia respect the panel mid-line following [4]
|
(14) |
Where Zn is a factor depending on the number of waves (figure 4):
Bending inertia is one of the main parameters to guarantee that the buckling waves are contained in the panel in case of buckling and no catastrophic global buckling failure. Therefore, it is important to correctly capture its idealization. A formulation was developed to consider the correct position of the center of gravity without the need of correcting factors.
Many structures are not flat, typically an aircraft fuselage structure is a curved stiffened panel. Curvature effect increases the critical buckling load so an equivalent flat panel analysis will provide conservative results. Nevertheless, this approach leads to a weight increase as it is shown in figure 5. In this figure is shown that a 20% critical buckling onset increase can be expected considering the curvature of the panel for a typical fuselage of 2000 mm radius submitted to pure compression. The effect grows exponentially when the curvature increases.
Therefore, the problem should be adapted to take into account the curvature effect. This increases the complexity of the problem adding two new set of unknowns to the problem. Typically for flat panels it was enough to know the out of plane displacements (w(x,y)). Within a curve panel, the axial (u) and transverse displacements (v) enter now into the problem
definition.The strain at any layer of the laminate (as for flat panels) can be expressed as the addition of the mid-plane strain plus the curvature contribution:
|
(15) |
Where:
|
(16) |
In comparison with flat panels, the existence of panel curvature couples the mid-plane strains {ε}o with the curvatures {κ}, or in other words, the in-plane displacements (u, v) with the out-of-plane displacements w. This coupling is the source of all the differences with the flat panel formulation and forces having to retain the three displacements (u, v, w) throughout the whole calculation process.
Additionally, the displacement functions cannot be only defined for w, but also for u and v. The full development must be done considering now these new equations and the number of new unknowns.
The transverse shear stiffness is a resin dominated property in composite panels. This means that this transverse stiffness is an order of magnitude lower than in-plane properties, leading to the inapplicability of Kirchhoff’s hypothesis in bending for thick composite plates. In figure 6 it is shown how transverse shear stiffness influence grows when the t/b ratio increases.
Kirchhoff’s assumption establishes that a cross section of a plate submitted to bending remains on a plane surface normal to the neutral axis. This could be applicable to isotropic material but it leads to non-conservative results in composite, as it is shown in figure 6. To solve it Mindlin’s plate theory [5] is applied, which considers that plane sections remain plane in bending but non perpendicular to the neutral axis. This hypothesis introduces a rotation due to transverse shear to be added to the usual rotation due to classical plate bending theory. This statement implies that the transverse strains ɣxz and ɣyz are independent of z.
In reality, the transverse shear stresses are not constant in each ply of the laminate, therefore shear correction factors k1 and k2 are used in deriving the transverse shear constitutive relations for the plate. Following the procedures of [5] and [6] ɣxz and ɣyz are replaced by k1·ɣxz and k2·ɣyz
respectively in the strain energy density expression of an orthotropic ply due to transverse shear:
|
(20) |
where [Q]k is stiffness matrix of lamina k in local axes. The value of k1 , k2 is determined such that the coefficient of Qx in the two shear deformation of strain energy expressions is identical, but it will not be developed here, see reference [11].
This is adding again, 2 new unknowns to the problem that can be approximated by the same shape functions described previously. Therefore 5 unknowns are now needed to be solved following procedure indicated in section 3.
Once the development is completed an extensive and intensive validation is done in three different bases:
Obviously non exhaustive validation can be presented in this brief paper but different illustrative examples are enclosed in the following figures, from figure 7 to figure 10, given the adequate confidence that the problem solution is well captured and therefore applicable for sizing and certification:
![]() |
![]() |
![]() |
![]() |
The mathematical development of the solutions presented here is really complex and its implementation is not straight forward for an efficient and optimum CPU management (see figure 11 for examples of different performances on computational time vs operations)
For example, a cubic running time function will be the result of a naive multiplication of two nxn matrices. Matrix chain multiplication with brute force can be time exponential [5].
If A is a 10 × 30 matrix, B is a 30 × 5 matrix, and C is a 5 × 60 matrix, then:
A quicker solution to this problem can be achieved by breaking up the problem into a set of related sub-problems, something that will be extensively used in the energy methods if the adequate performances would like to be achieved.
High complexity problems bring to face quickly the limitation of the technology and therefore adequate technologies need to be understood and managed. Software development skills are important. A poor implementation could ruin a master math development.
During the present paper, the different evolutions of energy methods inside Airbus has been presented. These solutions have the advantage of a high performance when compared with more general FE methods (~10 times faster), which makes them suitable to be used, for example, together with optimization engines. Moreover, they are suitable to be linked to sizing processes that chain different failure modes that can interact between them making prediction reliable, fast and safe as described in the introduction
Nevertheless, not all the effects and geometries are today covered by these solutions and therefore finite element models are used when applicability is superseded. Notwithstanding, the energy methods can be evolved to cover:
The development of these solutions, mainly the post-buckling behavior, will make a step forward on the buckling treatment reducing the need of finite element models and making available to the engineers a wide variety of solutions to support them in the design phase definitions of new composite structures.
[1] Hennessy and Patterson, « Morgan-Grampian Publishers », Computer Architecture, 2019
[2] Handbook of structural stability, Buckling of flat panels, NACA-TN-3781, July 1957
[3] J. N. Reddy, Energy Principles and Variational Methods in Applied Mechanics, Wiley, 2nd edition, 2002.
[4] The effect of longitudinal stiffeners located on one side of a plate on the compressive buckling stress of the plate-stiffener combination, NACA TN 2873, 1953
[5] Mindlin, M. D, Journal of Applied Mechanics, Influence of Rotatory Inertia and Shear on Flexural Motions of Isotropic, Elastic Plates p.31-38. , 1951
[6] Whitney, J. M., Pagano, N. J., Journal of Applied Mechanics, p.1031-1036. Shear Deformation in Heterogeneous Anisotropic Plates, 1970
[7] https://en.wikipedia.org/wiki/Matrix_chain_multiplication
[8] Critical combinations of Shear and Direct Axial Stress for Curved Rectangular Panels, NACA TN 1928, 1949
[9] Flat panels in shear. Buckling of long panels with transverse stiffeners, ESDU 02.03.02, 1983
[10] Martin-Esteban, José Antonio. Energy method for buckling of CFRP interconnected plates with arbitrary boundary conditions, 20th International Conference on Composite Materials, July 2015
[11] Carrasco-Fernández, José. Anisotropic Stiffened Panel Buckling and Bending Analyses using Rayleigh-Ritz Method. International Conference on Advanced Computational Engineering and Experimenting. ACE-X 2011
Accepted on 27/01/24
Submitted on 30/05/23
Licence: Other