Inversion for Constitutive Model Parameters of Bird in Case of Bird Striking

Yu, Peng; Yu, Siliu; He, Sheng; Huang, Xinheng; Yun, Weijing

doi:https://doi.org/10.1155/2022/2456777

Shock and Vibration

On this page

Abstract Introduction Results and Discussion Conclusion Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Responses and Protection of Structures against Multiple Dynamic Hazards

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 2456777 | https://doi.org/10.1155/2022/2456777

Inversion for Constitutive Model Parameters of Bird in Case of Bird Striking

Peng Yu,¹Siliu Yu,¹Sheng He,^1,2Xinheng Huang,¹and Weijing Yun¹

Academic Editor: Jingxuan Wang

Received19 Nov 2021

Revised31 Dec 2021

Accepted07 Jan 2022

Published09 Feb 2022

Abstract

The strike between flying birds and airplanes results in the unacceptable losses of aircraft structures. The normal approach to evaluate the bird impact resistance is the combined full-scale experimental-numerical study. However, the simulation results from the current available bird constitutive models are usually not in good agreement with the experimental data. Establishing a reasonable bird constitutive model is difficult and significant to the simulation of the bird striking process. In this paper, based on the displacement measurements of an aluminum plate subjected to soft body impact, an inversion of the bird constitutive model is conducted by using backpropagation (BP) neural network. A comparative evaluation of this inversion model and other constitutive models is carried out, indicating that the proposed inversion model is more reasonable.

1. Introduction

The strike between flying birds and airplanes frequently results in catastrophic damage to the aircraft. According to the International Bird Strike Committee, 55 fatal accidents, in which 108 aircraft were damaged and 277 people died, were reported between 1912 and 2009 [1]. The serious events of birds striking can lead to an unacceptable loss of aircraft and the death of crewmembers. In early studies, the common method to determine the capability against bird strikes is to carry out a full-scale bird-strike experiment [2]. However, various uncertainties may lead to diversities of experiments, and the experiments cost much. The development of advanced numerical techniques and computing technologies makes it possible to perform the coupling analysis and simulate the whole process of bird strike problems. However, the numerical results of current models are not in good agreement with the experimental data. Establishing a reasonable bird constitutive model is significant to the finite element analysis of bird strikes on aircraft structures.

The anatomic structure of birds includes several internal cavities like pneumatic bones, lungs, and peculiar air sacs, which is a complex nonlinear system. Moreover, the identification of a bird constitutive model and the corresponding material parameters is difficult. Several authors used an elastic-plastic bird constitutive model with a defined failure strain [3–5], while others highlighted the limitations of this simplified material law [6]. It was observed that no fluid-like flow response can be achieved with such an elastic-plastic model, only if the shear modulus G is set very low [7]. In Ref. [8] a rubber-like hyper-elastic Mooney–Rivlin model was adopted to simulate the bird's behavior. A more common method is an equation of state (EoS) for the bird constitutive model, by defining the pressure-volume relationship with parameters of water at room temperature. A polynomial form of the EoS and a simpler Murnaghan EoS were used in most studies [2, 9–14]. A further approach of the Grüneisen EoS was adopted in Refs. [15, 16]. However, in Ref. [17] it was reported that the Grüneisen EoS is only valid for solid materials that remained in the solid-state during the impact process and should therefore be carefully used for bird strike simulations and the usage of Grüneisen EoS dependent on the finite element code, since there are only a few of the abovementioned algorithms available in most commercial software codes. The material constants of all EoS have to be defined, which cannot be measured directly.

A regular technique is the parameter calibrations, for example, the elastic, elastic-plastic, or nonrevolving flow model is assumed to be the bird constitute model, and then the plate response during the bird-plate impact experiment is used to confirm the parameters of the constitutive model. This methodology was performed in studies such as [16, 18] using optimization software to determine the parameters in conjunction with bird strike experiments on instrumented plates. Nishikawa et al. [19] impacted gelatin birds on rigid targets and confirmed that the deformation of the bird model did behave like fluid and numerical simulation results correlated well with experiment data.

In recent years, the artificial neural network (ANN) has provided a fundamentally practical and powerful approach to material modeling. The basic advantage of ANN is that it does not explicitly need physical knowledge of the deformation mechanism and mathematical model. Darryl et al. [20] used a contact algorithm based on the Lagrange multiplier method to predict appropriate impact force, this method solved the severe contact-induced deformation due to a large contact force and highly deformable projectile. Previous studies usually adopted the Lagrangian formulation to model bird strikes. However, due to element instability, refining mesh elements and narrowing the time-step method had to be used to alleviate the problem, to solve this problem, nowadays more researchers use arbitrary Lagrange–Euler (ALE) [19, 21–23] and smoothed particle hydrodynamics (SPH) [19, 24–27] method in bird strike problems [25]. However, there are few reports about the application of ANN to model the bird constitutive model by using SPH Formulation during the impact process. In this paper, based on the displacement measurement results of bird impact on the aluminum plate, an ANN model with an error backpropagation learning algorithm has been applied to an inversion on bird constitutive model parameters and the parameters are obtained by inversion. The validity and reliability of the inversion model are discussed. Furthermore, a comparative evaluation of three other constitutive equations and the trained network model is carried out.

2. Introduction of BP Neural Network

The artificial neural network is an intelligent information-treatment system with the characteristics of adaptive learning and treating complex relationships [28]. An artificial neural network (ANN) consists of a large number of highly connected artificial units (neurons or nodes). Each neural network comprises an input layer, an output layer, and one or more hidden layers, which are connected by the processing units, as can be seen in Figure 1. The input layer is used to receive data, while the output layer sends the output values to users. The role of the hidden layer is to provide the necessary complexity (which denotes to abstract the characteristics of the input data to another dimension space to show its more abstract features) for nonlinear problems. In the present model, sigmoid functions are employed due to the processing units for computational convenience.

Among various types of ANN models, an iterative gradient algorithm called backpropagation (BP) algorithm is popular in the field in engineering applications and it is used for the present work. The training of a network by backpropagation involves three stages: the feedforward of the input training pattern, the training and backpropagation with the associated error, and the adjustment of the weights. When the artificial neural network is established for the soft body constitutive model, a feedforward backpropagation algorithm is selected to train the network.

3. An Example

3.1. Flow of Soft Body Parameter Optimization

First of all, determine the soft body parameters to be optimized, and then according to the experience to determine the range of parameters, the values determined according to the experimental data of the soft body impact aluminum sample points. And then the sample parameters in LS-DYNA in the calculation, can get aluminum midpoint displacement curve of time, choose corresponding time point in the experimental data of the displacement, training sample is completed. After the improved BP neural network training, the required soft body model parameters are output and the stable network unit connection weights are obtained. Then network training is carried out to obtain the soft body model parameters. Lastly, the obtained parameters are put into LS-DYNA for calculation, and the calculation results are compared with the experimental results. The flow chart of soft body parameter optimization as shown in Figure 2.

3.2. Construction of Bird Body Constitutive Model

Newtonian fluid is selected to describe the constitutive model of the bird body, since the real bird body is a structure with flesh and blood [29]. Newtonian fluid is a viscous fluid with a linear proportional relationship between shear stress and deformation gradient, and its stress-strain relationship is determined by the following equation:where σ_ij is the stress tensor, V_ij is the deformation rate tensor, D_ijkl is the viscosity coefficient tensor of fluid, and P is the static pressure. The term pδ_ij represents the state of stress in the fluid at rest (when V_ij = 0). Assume that static pressure P depends on the density and temperature of the fluid according to some equation of state. For Newtonian fluids, we assume that the elements of the tensor D_ijkl can be temperature-dependent and the stress or deformation rate is independent. Tensor D_ijkl is of order 4 with 3⁴ = 81 elements.

If the fluid is isotropic, that is, the tensor D_ijkl has the same component matrix in any Cartesian rectangular system, then D_ijkl can be expressed by two independent constants λ and μ as the following equation:

Then,

Combining equations

If the mean normal stress 1/3σ_kk is assumed to be independent of the volume expansion rate V_kk, then there must be

Thus, the constitutive equation becomes

This formula was put forward by Stokes, so the fluid subject to the above equation is called Stokes fluid, for this kind of fluid only a material constant μ (i.e., viscosity coefficient) is enough to determine its properties.

If the fluid is incompressible, then V_kk = 0, and the constitutive equation of incompressible viscous fluid can be obtained:

If μ = 0, the constitutive equation of inviscid fluid is obtained:

The appearance of hydrostatic pressure P marks the fundamental difference between hydrodynamics and elasticity. To accommodate the needs of this new variable, it is usually assumed that there exists some equation of state, which gives the relationship between density ρ, pressure P, and absolute temperature T, that is

3.3. Finite Element Model

The bird strike setup is conducted based on the experiment of the bird-impact plate in Figure 3 by using LS-DYNA element code.

(a)

(b)

The flat plate fixed along the four edges is impacted by a soft body at the center of the flat plate at 91.1 m/s with an impact angle of 30°. The dimensions of the flat plate are 525 mm × 430 mm × 10 mm. The geometric model of a bird is simulated by a cylinder with a hemisphere at the two ends, the diameter of the hemisphere is 53 mm and the length of the cylinder is 156 mm. The mass of the bird is 1.8 kg.

In the finite element model, the 8-node solid and 164-element type in the LS-DYNA element library is employed for the bird and flat plate. In addition, hourglass control and single-point integration are also used. The numbers of the elements for the bird and flat plate model are 4504 and 2000.

3.4. Materials Properties

3.4.1. Materials Properties of Flat Plate

The plate is made from LY12-CZ aluminum alloy. The material parameters of the flat plate are mentioned in Table 1. The material properties of the flat plate are modeled by the MAT_PLASTIC_KINEMATIC in LS-DYNA code where the values are set in Table 1.

3.4.2. Material Properties of Bird

Generally speaking, a trend is visible towards the utilization of SPH models with an equation of state for the material of bird using the properties of water–air-mixture at room temperature. A very good correlation to experimental results could be obtained with these models. For the fluid model, in LS-DYNA, the material model of NULL (empty material) is used to describe its material properties. The parameters to be input are shown in Table 2.

According to the property of the bird in the experiment, using fluid constitutive and Gruneisen state equation to describe soft body is closer to the reality. The bird model is defined by MAT_NULL and the equation of state (EoS) is specified. The inputs of parameters are shown in Table 3.

During the impact process, the deformation of the bird is described by the equation of state. Elastic modulus and passion ratio are 0, and the density is 938 kg/m³. The four parameters of the constitutive model of the bird to be determined are shown in Table 4.

3.5. Identification of the Bird Model and the Inversion of the Parameters

The displacements at 12 different times are chosen as the inputs, and the outputs of the model are four constitutive parameters of the soft body. The 12 different times are based on the experimental displacement-time data of the center of the flat plate. The experimental displacement-time data of the center of the flat plate is given in Figure 4.

Before the network training, both input and output variables are normalized within the range 0 to 1 in order to obtain a reliable form for the neural network ^[31]. If the number of hidden layers is too small, the trained network might not have sufficient ability to learn the process correctly. As a result, in order to choose the appropriate number of hidden layers, the trial-and-error procedure is started with seven neurons in the hidden layer and further carried out with up to fourteen neurons. It is found that when the number of neurons is 9, the mean square error reaches the minimum.

In the bird constitutive model, the values of the four parameters have their ranges based on the prior knowledge, specific ranges (see Table 5).

The parameters to be optimized are uniformly divided into a number of values with the space. According to the orthogonal design and uniform design principle, these values are divided into 25 groups, and each sample group contains four parameters, which is shown in Table 6.

Substituting these parameters of 25 groups into the bird constitutive model, 25 displacement-time curves of the center of the flat plate can be obtained in LS-DYNA simulation. Selecting the displacements of the center of the flat plate, the training of the network can be finished, which is shown in Table 7.

When developing the ANN model, the samples of 1, 5, 10, 15, 20, and 25 are removed to test the generalization capability of the network. The optimized ANN model consists of 12 input neurons, 4 output neurons, and a single hidden layer (with 9 neurons) and the transfer functions of the optimized ANN model are ‘tan sigmoid’ and ‘pure linear’. A feedforward backpropagation algorithm is selected to train the network. The setting of other training parameters for the neural network is listed in Table 8.

4. Results and Discussion

The performance of ANN in the training stage is shown in Figure 5. It can be seen that the trend of error reduction remains unchanged after 2000 training cycles. After being trained, the ANN model can map the nonlinear relationship between the parameters of the bird constitutive model and the displacements of the center of the flat plate.

The optimized parameters can be obtained, as can be seen in Table 9. The comparison of the numerical and experimental displacement of the center of the flat plate using the optimized constitutive model of bird is shown in Figure 6, indicating the simulation result using the parameters obtained from the developed network is in good agreement with experimental data. In order to verify the identification of the constitutive model after optimization, a comparative evaluation of three other constitutive equations and the present bird constitutive model is carried out. The parameters and stress distributing graphs of elastic model, plastic kinematic model, and elastic-plastic hydrodynamic Model are shown in Tables 10–12 and Figure 7, respectively, and the parameters of EOS_LINEAR_POLYNOMIA are shown in Table 13. As can be seen in Figure 8, the model of the bird remains in good condition while the flat plate with obvious deformation in the elastic model. The result of the plastic kinematic model and the elastic-plastic hydrodynamic model show a similar phenomenon that the model of bird has large deformation while the flat plate with little, and the fringe levels drop significantly.

(a)

(b)

(c)

During the impact process, different bird constitutive models result in different displacement and force-time histories of the center of the flat plate. The maximum displacement and the force of the center of the flat plate of the elastic constitutive model are the largest among the four constitutive models, in which the deviations from experimental results are also the largest. Figure 9 shows that the present bird constitutive model developed from the present network is in better agreement with experimental data than the other three constitutive models.

5. Conclusion

In this paper, the artificial neural network model with a backpropagation learning algorithm was used to develop the constitutive relationship for the soft body during the impact process. The feasibility and practicability of this identification method were verified through a specific example. From the results of this present study, the main conclusions are drawn as follows:(1)Different flexible body constitutive models are used, and the impact of the soft body presents different modes, and the displacement curve and impact force curve of the aluminum plate center is also different. Among them, the displacement and impact force of the linear elastic constitutive model is the largest, but the deviation from experimental results is also the largest.(2)The calculated results of the soft body constitutive parameters optimized by the improved BP neural network algorithm are in good agreement with the experimental results. At the same time, it can be seen that the calculation curve optimized by BP neural network algorithm almost passes all the test points, which proves the accuracy and reliability of the optimized flexible body parameters, and the calculation convergence accuracy also meets the requirements.(3)As a calculation method, BP neural network can be used in the optimization calculation of constitutive parameters of soft body impact. However, it should be pointed out that each constitutive model needs to build its own neural network when it is used for calculation, and the optimization results and errors will be different due to the difference of network structure.

Data Availability

The data used to support the findings of this study are included within the paper.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This study was funded by National Natural Science Foundation (12102095), the Scientific Research Project for Young Innovative Talents of Guangxi Province (20200312), Research grant for Talent of Guangxi Plan (100), the starting research grant for High-level Talents from Guangxi University, the Science and Technology Major Project of Guangxi Province (AA18118055), the Guangxi Natural Science Foundation (20181096), and application of key technology in building construction of prefabricated steel structure (20190528).

References

J. Thorpe, “Update on fatalities and destroyed civil aircraft due to bird strikes with appendix for 2008 & 2009,” in Proceedings of the 29th International Bird Strike Committee Meeting, Cairns, Australia, 2010.
View at: Google Scholar
D. Zhang and Q. Fei, “Effect of bird geometry and impact orientation in bird striking on a rotary jet-engine fan analysis using SPH method,” Aerospace Science and Technology, vol. 54, no. 54, pp. 320–329, 2016.
View at: Publisher Site | Google Scholar
S. Zhu, M. Tong, and Y. Wang, “Experiment and numerical simulation of a full-scale aircraft windshield subjected to bird impact,” in Proceedings of the 50th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, Palm Springs, CA, USA, May 2009.
View at: Google Scholar
R. Doubrava and V. Strnad, “Bird strike analyses on the parts of aircraft structure,” in Proceedings of the 27th Congress of the International Council of the Aeronautical Sciences, Nice, France, September 2010.
View at: Google Scholar
Y. K. Zhang and Y. L. Li, “Analysis of the anti-bird impact performance of typical beam-edge structure based on ANSYS/LS-DYNA,” Advanced Materials Research, vol. 33-37, pp. 395–400, 2008.
View at: Publisher Site | Google Scholar
J. L. Yang, X. J. Cai, and C. H. Wu, “Experimental and FEM study of windshield subjected to high speed bird impact,” Acta Mechanica Sinica, vol. 19, no. 6, pp. 543–550, 2003.
View at: Publisher Site | Google Scholar
F. Stoll and R. A. Brockman, “Finite elelment simulation of fflgh-speed soft-body impacts,” 1997.
View at: Google Scholar
Y. N. Gong and S. Q. Xu, “Bird impact analysis of aircraft windshield transparency,” Chinese Journal of Aeronautics, vol. 5, no. 2, 1992.
View at: Google Scholar
Y. Zhou, Y. Sun, and W. Cai, “Bird-striking damage of rotating laminates using SPH-CDM method,” Aerospace Science and Technology, vol. 84, no. 84, pp. 265–272, 2019.
View at: Publisher Site | Google Scholar
S. A. Meguid, R. H. Mao, and T. Y. Ng, “FE analysis of geometry effects of an artificial bird striking an aeroengine fan blade,” International Journal of Impact Engineering, vol. 35, no. 6, pp. 487–498, 2008.
View at: Publisher Site | Google Scholar
A. Khan, R. Kapania, and E. Johnson, “A review of soft body impact on composite structure,” in Proceedings of the 51st AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference 18th AIAA/ASME/AHS Adaptive Structures Conference, Orlando, FL, USA, April 2010.
View at: Google Scholar
R. H. Mao, S. A. Meguid, and T. Y. Ng, “Effects of incidence angle in bird strike on integrity of aero-engine fan blade,” International Journal of Crashworthiness, vol. 14, no. 4, pp. 295–308, 2009.
View at: Publisher Site | Google Scholar
S. Q. Jia, F. S. Wang, J. J. Zhou, Z. Jiang, and B. Xu, “Study on the mechanical performances of carbon fiber/epoxy composite material subjected to dynamical compression and high temperature loads,” Composite Structures, vol. 258, Article ID 113421, 2021.
View at: Publisher Site | Google Scholar
M. Anghileri, L. L. Castelletti, and D. Molinelli, “A strategy to design bird-proof spinners,” in Proceedings of the 7th European LS-DYNA Users Conference, Salzburg, Austria, May 2009.
View at: Google Scholar
M. Anghileri, L. Castelletti, and V. Mazza, “Birdstrike: approaches to the analysis of impacts with penetration,” Impact loading of lightweight structures, vol. 49, pp. 63–74, 2005.
View at: Google Scholar
L. S. Nizampatnam, “Models and methods for bird strike load predictions,” 2007.
View at: Google Scholar
M. A. Mccarthy, J. R. Xiao, C. T. Mccarthy et al., “Modelling of bird strike on an aircraft wing leading edge made from fibre metal laminates-part 2: modelling of impact with SPH bird model,” Applied Composite Materials, vol. 11, no. 5, pp. 317–340, 2004.
View at: Publisher Site | Google Scholar
M. A. Lavoie, A. Gakwaya, M. N. Ensan, D. G. Zimcik, and D. Nandlall, “Bird's substitute tests results and evaluation of available numerical methods,” International Journal of Impact Engineering, vol. 36, pp. 1276–1287, 2009.
View at: Publisher Site | Google Scholar
M. Nishikawa, K. Hemmi, and N. Takeda, “Finite-element simulation for modeling composite plates subjected to soft-body, high-velocity impact for application to bird-strike problem of composite fan blades,” Composite Structures, vol. 93, no. 5, pp. 1416–1423, 2011.
View at: Publisher Site | Google Scholar
L. Y. Darryl, L. C. John, M. E. Richard, and L. W. Kate, “Hazing and movements of Canada geese near Elmendorf air force base in Anchorage, Alaska,” International Biodeterioration & Biodegradation, vol. 45, pp. 103–110, 2000.
View at: Publisher Site | Google Scholar
A. F. Johnson and M. Holzap, “Modelling soft body impact on composite structures,” Composite Structures, vol. 61, no. 1-2, pp. 103–113, 2003.
View at: Publisher Site | Google Scholar
M. A. Lavoie, A. Gakwaya, and M. N. Ensan, “Validation of available approaches for numerical bird strike modeling tools,” International Review of Mechanical Engineering, vol. 1, 2013.
View at: Google Scholar
F. Allaeys, G. Luyckx, and W. V. Paepegem, “Characterization of real and substitute birds through experimental and numerical analysis of momentum, average impact force and residual energy in bird strike on three rigid targets: a flat plate, a wedge and a splitter,” International Journal of Impact Engineering, vol. 99, pp. 1–13, 2017.
View at: Publisher Site | Google Scholar
R. Hedayati, M. Sadighi, and M. Mohammadi-Aghdam, “On the difference of pressure readings from the numerical, experimental and theoretical results in different bird strike studies,” Aerospace Science and Technology, vol. 32, pp. 260–266, 2014.
View at: Publisher Site | Google Scholar
R. Hedayati and S. Ziaei-Rad, “A new bird model and the effect of bird geometry in impacts from various orientations,” Aerospace Science and Technology, vol. 28, pp. 9–20, 2013.
View at: Publisher Site | Google Scholar
J. Liu and Y. L. Li, “Numerical simulation of a rotary engine primary compressor impacted by bird,” Chinese Journal of Aeronautics, vol. 4, no. 26, pp. 926–934, 2013.
View at: Publisher Site | Google Scholar
J. M. Zurada, Introduction to Artificial Neural Systems, Jaico Publishing House, Chennai, India, 1992.
Y. Zhang and Y. Li, “Identification of the bird constitutive model,” in Proceedings of the 3rd International Conference on Mechanical Engineering and Mechanics, Beijing China, June 2009.
View at: Google Scholar
K. Swingler, Applying Neural Networks: A Practical Guide, Morgan Kaufmann, SanFrancisco, CA, USA, 1996.

Copyright

Copyright © 2022 Peng Yu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

357

Downloads

586

Citations