Next Article in Journal
OpenAP.top: Open Flight Trajectory Optimization for Air Transport and Sustainability Research
Next Article in Special Issue
A Strong Form Meshless Method for the Solution of FGM Plates
Previous Article in Journal
Design and Analysis of MataMorph-3: A Fully Morphing UAV with Camber-Morphing Wings and Tail Stabilizers
Previous Article in Special Issue
Identification of Time Variations of Moving Loads Applied to Plates Resting on Viscoelastic Foundation Using a Meshfree Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Shape Optimisation of Assembled Plate Structures with the Boundary Element Method

by
Llewellyn Morse
1,*,
Vincenzo Mallardo
1,
Zahra Sharif-Khodaei
2 and
Ferri M.H. Aliabadi
2
1
Department of Architecture, University of Ferrara, Via Quartieri 8, 44121 Ferrara, Italy
2
Department of Aeronautics, Imperial College London, South Kensington Campus, City and Guilds Building, Exhibition Road, London SW7 2AZ, UK
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(7), 381; https://doi.org/10.3390/aerospace9070381
Submission received: 30 May 2022 / Revised: 5 July 2022 / Accepted: 13 July 2022 / Published: 15 July 2022
(This article belongs to the Special Issue Recent Advances in Computational Mechanics)

Abstract

:
A novel methodology is presented for performing sensitivity analyses of assembled plate structures using the Boundary Element Method (BEM). The main novelty of this work is that the exact implicit derivatives of the BEM formulations for assembled plate structures have been derived for the first time and incorporated into a newly developed Implicit Differentiation Method (IDM), enabling sensitivity analyses to be conducted for more complex and realistic structures in a more accurate, robust, and efficient manner than previous approaches. Three numerical examples are investigated to validate the derived exact implicit derivatives and to demonstrate how they could be used for a potential application involving the shape optimisation of a complex X-core structure from the canard of a Eurofighter Typhoon fighter jet. Results show that the newly developed IDM is more accurate, robust, and efficient when compared to alternative methodologies using derivatives obtained from methods such as the Finite Difference Method (FDM) and the Finite Element Method (FEM).

1. Introduction

Sensitivity analysis plays a key role in the design of structures. It is very important for engineers to understand the sensitivity of a structure’s performance to changes in geometrical, material, or loading parameters. Sensitivity analysis also enables engineers to optimise the design of their structures and maximise their safety.
The Boundary Element Method (BEM) is a very useful tool for conducting sensitivity analyses of structures. When modelling a structure, the BEM only requires the outer boundary of the structure to be discretised into elements. This makes the BEM a very efficient tool for sensitivity analyses, because only a relatively small part of the model needs to be re-meshed during the analyses. Furthermore, the BEM can often achieve a similar level accuracy to the Finite Element Method (FEM) while using a courser mesh with fewer elements, indicating that the BEM would be a very efficient method for optimisation problems involving the calculation of many sensitivities over many iterations [1].
Most prior work on structural sensitivity analysis with the BEM amongst the research community has involved 2D structures [2,3,4,5,6,7,8,9,10,11], and 3D structures [12,13,14,15,16,17] to a smaller degree. There have only been a few prior research works on the topic of structural sensitivity analysis with the BEM that have involved plate structures, with some examples being [18,19]. There are many structures in engineering, particularly in aircraft, that are not able to be modelled accurately in 2D. Furthermore, although modelling a structure in 3D can lead to very accurate results, it is often more computationally efficient to model it as a plate structure, without a noticeable loss in modelling accuracy. Previous works on plate structural sensitivity analysis with the BEM [18,19] have only considered simple single-plate structures, which have a limited range of application. Assembled-plate structures, on the other hand, have a much greater range of application and enable more complex plate structures to be investigated. However, there are several challenges associated with the sensitivity analysis of assembled structures—there are significantly more geometric parameters involved compared to single-plate structures and the interaction between these parameters is more complex, increasing the computation time and complexity. The BEM, due to its inherent advantages described earlier, would be very effective at addressing these challenges. Therefore, this current work aims to build upon previous work by the research community and develop a novel methodology for conducting sensitivity analyses for assembled-plate structures—enabling complex plate structures to be investigated both accurately and efficiently.
A common approach to calculating sensitivities is with the Finite Difference Method (FDM). The FDM is a relatively simple method that uses finite differences to calculate sensitivities. Although simple to implement, a drawback of the FDM is that its accuracy is highly dependent on the step-size used. An alternative approach to the FDM would be to derive the exact implicit derivatives of the BEM formulations. This Implicit Differentiation Method (IDM), since it uses the exact derivatives, would be more accurate and robust than the FDM. In addition, by exploiting the inherent advantages of the BEM mentioned earlier, it would also be more computationally efficient, since large parts of the BEM mesh would remain unchanged during the sensitivity analyses. Previous examples of work by the research community involving the derivation of the exact implicit derivatives of BEM formulations include [2,20,21,22] for 2D structures, [23] for 3D structures, and [18,19] for simple single-plate structures. There have not yet been any works involving assembled-plate structures.
In summary, the main novelty of this current work is that the exact implicit derivatives of the BEM assembled-plate structure formulations have been derived for the first time. Previous works have only derived these derivatives for simple single-plate structures. By deriving the derivatives for assembled plate structures, a wider range of more realistic and complex structures can be investigated. The use of these derivatives, along with the inherent advantages of the BEM, can significantly improve the accuracy, robustness, and computational efficiency when performing structural sensitivity analyses of these more realistic and complex structures. Potential applications of this newly developed methodology include the shape optimisation or structural reliability analysis of complex plate structures.
The formulations behind the BEM for assembled plate structures, as well as the exact implicit derivatives of these formulations, are presented in Section 2. Several numerical examples are presented in Section 3 to validate the derived exact implicit derivatives and to demonstrate their accuracy, robustness, and efficiency. A numerical example involving the shape optimisation of a complex assembled plate structure is also presented to demonstrate a potential application of the derived exact implicit derivatives.

2. Methodology

In this section, the formulations of the Boundary Element Method (BEM) and the Implicit Differentiation Method (IDM) for assembled plate structures are presented. In this work, Latin letter indexes (e.g., i, j, k) can take values from 1 to 3, whereas Greek letter indexes (e.g., α , β , ρ , γ ) can take values of either 1 or 2.

2.1. Plate Theory Notation in the Boundary Element Method

Figure 1 shows a plate of thickness h ( h / 2 x 3 + h / 2 ). The membrane of the plate is the x 1 x 2 plane.
The sign convections for the displacements, rotations, tractions, and moments for plates is shown in Figure 2.
The in-plane displacements of the plate are given in Equation (1):
u = u 1 u 2
The rotations and out-of-plane displacement of the plate are given in Equation (2):
w = w 1 w 2 w 3
where w 1 and w 2 are rotations, and w 3 is the out-of-plane displacement of the plate in direction x 3 .
The membrane tractions are given in Equation (3):
t = t 1 t 2 = N 11 n 1 + N 12 n 2 N 21 n 1 + N 22 n 2
where N α β are the membrane stress resultants, and n β are the components of the outward normal vector of the boundary of the plate.
The bending and shear tractions are given in Equation (4):
p = p 1 p 2 p 3 p 1 p 2 = M 11 n 1 + N 12 n 2 M 21 n 1 + N 22 n 2 p 3 = Q 1 n 1 + Q 2 n 2
where M α β are the bending stress resultants, and Q α are the shear stress resultants.

2.2. The Boundary Element Method for Plate Structures

In the approach developed by Dirgantara and Aliabadi [24], the Boundary Element Method (BEM) for plate structures involves using the Reissner theory to model the bending and shear behaviour of plates, whereas 2D plane stress theory is used to model the behaviour of the membrane of plates. The BEM equation for the plate membrane is given in Equation (5):
C m u + n e = 1 N e Γ n e T m ψ e T d Γ n e u n e = n e = 1 N e Γ n e U m ψ e T d Γ n e t n e
The BEM equation for bending and shear is given in Equation (6):
C b w + n e = 1 N e Γ n e T b ψ e T d Γ n e w n e = n e = 1 N e Γ n e U b ψ e T d Γ n e p n e + q 3 n e = 1 N e Γ n e B b d Γ n e
N e denotes the number of elements on the external boundary. T m , T b , U m , U b , and B b are fundamental solutions evaluated along the surface of element n e , Γ n e . q 3 is a uniform load applied over the surface of the plate. The terms C m and C b are free terms and their values can be directly evaluated from a consideration of rigid body motion [1]. Quadratic boundary elements, consisting of three nodes and ten Gaussian quadrature integration points, were used in this work. ψ e T contains the continuous shape functions for the quadratic boundary elements. More information on the numerical integration approach used in this work can be found in [1]. Detailed forms of the above equations can be found in [1,19,22,24].
u n e is a ( 2 × 1 ) vector that contains the in-plane displacements for element n e and is shown in Equation (7):
u n e = u 1 n e u 2 n e
w n e is a ( 3 × 1 ) vector that contains rotations w 1 and w 2 , and vertical displacement w 3 for element n e and is shown in Equation (8):
w n e = w 1 n e w 2 n e w 3 n e
t n e is a ( 2 × 1 ) vector that contains the tractions due to membrane stress resultants for element n e and is shown in Equation (9):
t n e = t 1 n e t 2 n e
p n e is a ( 3 × 1 ) vector that contains bending tractions p 1 and p 2 , and shear traction p 3 for element n e and is shown in Equation (10):
p n e = p 1 n e p 2 n e p 3 n e
By using the above equations, the system of equations shown in Equation (11) can be created:
Hu = Gt
where u contains displacements and rotations, and t contains membrane, bending, and shear tractions. G and H are matrices of coefficients. This system can be rewritten as shown in Equation (12):
AX = F
where X contains all of the unknown displacements, rotations, and tractions. A is a matrix that contains the coefficients from H and G . F contains the known displacements, rotations, and tractions multiplied by the corresponding coefficients from H and G . Details of this procedure are shown in [22].

2.3. The Boundary Element Method for Assembled Plate Structures

The above procedure is straightforward for a structure consisting of a single plate. However, the range of application for a structure consisting of a single plate is very limited. To widen the range of application of the above procedure, it must be extended to structures consisting of multiple plates joined together. This would significantly expand the range of application of the above procedure and enable more complex structures, including many structures found in aircraft, to be investigated in an accurate, robust, and computationally efficient manner.
An example of an assembled plate structure consisting of two plates is shown in Figure 3. Each plate has its own local coordinate system, which is used when calculating the matrix A of each plate. The plates are connected along a junction line.
The A matrix of the assembled plate structure in Figure 3 is shown in Figure 4. The parts of the A matrices corresponding to the edges 1–4 are shown. To account for the connection between the two plates at the junction line, additional columns are added to the A matrix of each plate, and are labelled as ‘J’. The number of elements on the boundary of each plate is assumed to be equal, and is N E . The number of elements along the junction line is N J . The equilibrium equations, which describe the relationships between the displacements and tractions of the connected plates, can be found in the matrices labelled E .
For an assembled plate structure, the plates are joined together by junction lines, as shown in Figure 3. The junction line is along the x 2 -axis. Each plate has its own local set of coordinates ( x 1 i , x 2 i , x 3 i ) with i = 1 , 2 , , N , where N is the number of plates. To account for the fact that neighbouring plates will influence one another in terms of displacements and tractions, equilibrium equations must be derived along junction lines. These equilibrium equations are used in the matrices labelled E in Figure 4.
The displacement equilibrium equations are shown in Equations (13)–(17):
u 1 1 c o s ( θ 1 ) w 3 1 s i n ( θ 1 ) = u 1 2 c o s ( θ 2 ) w 3 2 s i n ( θ 2 ) = = u 1 N c o s ( θ N ) w 3 N s i n ( θ N )
u 1 1 s i n ( θ 1 ) + w 3 1 c o s ( θ 1 ) = u 1 2 s i n ( θ 2 ) + w 3 2 c o s ( θ 2 ) = = u 1 N s i n ( θ N ) + w 3 N c o s ( θ N )
u 2 1 = u 2 2 = = u 2 N
w 1 1 = w 1 2 = = w 1 N
w 2 1 = w 2 2 = = w 2 N = 0
The traction equilibrium equations are shown in Equations (18)–(21):
i = 1 N t 1 i c o s ( θ i ) p 3 i s i n ( θ i ) = 0
i = 1 N t 1 i s i n ( θ i ) + p 3 i c o s ( θ i ) = 0
i = 1 N t 2 i = 0
i = 1 N p α i = 0
where θ i is the angle of the local coordinate axis x 1 i of the i’th plate to the global coordinate axis x 1 , and u 1 i is the displacement in the local coordinate axis x 1 i of the i’th plate, etc. In the equilibrium equations shown above, it is assumed that there is no offset between the membranes of the joined plates. The equations in the case where is there an offset can be found in [1].

2.4. The Implicit Differentiation Method for Plate Structures

In the Implicit Differentiation Method (IDM), the derivatives of Equations (5) and (6) can be calculated with respect to a geometric parameter Z m (such as the length or width of a plate), or with respect to the thickness of the plate h.
The derivative of Equation (5) with respect to Z m is shown in Equation (22):
C m u , m + n e = 1 N e Γ n e T , m m ψ e T d Γ n e + Γ n e T m ψ e T ( d Γ n e ) , m u n e + n e = 1 N e Γ n e T m ψ e T d Γ n e u , m n e = n e = 1 N e Γ n e U , m m ψ e T d Γ n e + Γ n e U m ψ e T ( d Γ n e ) , m t n e + n e = 1 N e Γ n e U m ψ e T d Γ n e t , m n e
The derivative of Equation (5) with respect to h is shown in Equation (23):
C m u , h + n e = 1 N e Γ n e T , h m ψ e T d Γ n e + Γ n e T m ψ e T ( d Γ n e ) , h u n e + n e = 1 N e Γ n e T m ψ e T d Γ n e u , h n e = n e = 1 N e Γ n e U , h m ψ e T d Γ n e + Γ n e U m ψ e T ( d Γ n e ) , h t n e + n e = 1 N e Γ n e U m ψ e T d Γ n e t , h n e
The derivative of Equation (6) with respect to Z m is shown in Equation (24):
C b w , m + n e = 1 N e Γ n e T , m b ψ e T d Γ n e + Γ n e T b ψ e T ( d Γ n e ) , m w n e + n e = 1 N e Γ n e T b ψ e T d Γ n e w , m n e = n e = 1 N e Γ n e U , m b ψ e T d Γ n e + Γ n e U b ψ e T ( d Γ n e ) , m p n e + n e = 1 N e Γ n e U b ψ e T d Γ n e p , m n e + q 3 n e = 1 N e Γ n e B , m b d Γ n e + Γ n e B b ( d Γ n e ) , m
The derivative of Equation (6) with respect to h is shown in Equation (25):
C b w , h + n e = 1 N e Γ n e T , h b ψ e T d Γ n e + Γ n e T b ψ e T ( d Γ n e ) , h w n e + n e = 1 N e Γ n e T b ψ e T d Γ n e w , h n e = n e = 1 N e Γ n e U , h b ψ e T d Γ n e + Γ n e U b ψ e T ( d Γ n e ) , h p n e + n e = 1 N e Γ n e U b ψ e T d Γ n e p , h n e + q 3 n e = 1 N e Γ n e B , h b d Γ n e + Γ n e B b ( d Γ n e ) , h
In the above equations, T , m m , T , m b , U , m m , U , m b , and B , m b are the derivatives of the fundamental solutions with respect to a geometric parameter Z m , and their full expressions can be found in [22]. T , h m , T , h b , U , h m , U , h b , and B , h b are the derivatives with respect to plate thickness h, and their full expressions can be found in [19].
The derivatives of u n e , w n e , t n e , and p n e in Equations (7)–(10) are shown in Equations (26)–(29):
u , m n e = u 1 , m n e u 2 , m n e u , h n e = u 1 , h n e u 2 , h n e
w , m n e = w 1 , m n e w 2 , m n e w 3 , m n e w , h n e = w 1 , h n e w 2 , h n e w 3 , h n e
t , m n e = t 1 , m n e t 2 , m n e t , h n e = t 1 , h n e t 2 , h n e
p , m n e = p 1 , m n e p 2 , m n e p 3 , m n e p , h n e = p 1 , h n e p 2 , h n e p 3 , h n e
Using Equations (22) and (24), a system of equations can be formed as shown in Equation (30):
H , m u + H u , m = G , m t + G t , m
Likewise, using Equations (23) and (25), another system of equations can be formed as shown in Equation (31):
H , h u + H u , h = G , h t + G t , h
where H , G , u , and t are from Section 2.2. H , m , G , m , u , m , t , m , H , h , G , h , u , h , and t , h are their derivatives. Equations (30) and (31) can be rewritten as shown in Equations (32) and (33):
A X , m = F , m A , m X
A X , h = F , h A , h X
In the above system of equations, A and X are outputs from Equation (12). A , m and F , m are found by rearranging Equation (30), whereas A , h and F , h are found by rearranging Equation (31). Therefore, the right-hand sides of Equations (32) and (33) are known, and the only unknowns are X , m and X , h , which can be found via LU decomposition.

2.5. The Implicit Differentiation Method for Assembled Plate Structures

The matrix A , m and the matrix A , h of the assembled plate structure in Figure 3 are shown in Figure 5 and Figure 6, respectively. The matrix A , m and the matrix A , h of each plate were calculated in terms of its local coordinate system. The derivatives of the equilibrium equations, which describe the relationships between the derivatives of the displacements and tractions of the connected plates, can be found in the matrices labelled E , m and E , h in Figure 5 and Figure 6, respectively.
For an assembled plate structure, it is necessary to calculate the derivatives of the displacement equilibrium Equations (13)–(17) and the derivatives of the traction equilibrium Equations (18)–(21). These derivatives are used in the matrices labelled E , m and E , h in Figure 5 and Figure 6, respectively.
The derivatives of the displacement equilibrium Equations (13)–(17) with respect to a geometrical parameter Z m are shown in Equations (34)–(38):
u 1 , m 1 c o s ( θ 1 ) + u 1 1 c o s ( θ 1 ) , m w 3 , m 1 s i n ( θ 1 ) w 3 1 s i n ( θ 1 ) , m = u 1 , m 2 c o s ( θ 2 ) + u 1 2 c o s ( θ 2 ) , m w 3 , m 2 s i n ( θ 2 ) w 3 2 s i n ( θ 2 ) , m = = u 1 , m N c o s ( θ N ) + u 1 N c o s ( θ N ) , m w 3 , m N s i n ( θ N ) w 3 N s i n ( θ N ) , m
u 1 , m 1 s i n ( θ 1 ) + u 1 1 s i n ( θ 1 ) , m + w 3 , m 1 c o s ( θ 1 ) + w 3 1 c o s ( θ 1 ) , m = u 1 , m 2 s i n ( θ 2 ) + u 1 2 s i n ( θ 2 ) , m + w 3 , m 2 c o s ( θ 2 ) + w 3 2 c o s ( θ 2 ) , m = = u 1 , m N s i n ( θ N ) + u 1 N s i n ( θ N ) , m + w 3 , m N c o s ( θ N ) + w 3 N c o s ( θ N ) , m
u 2 , m 1 = u 2 , m 2 = = u 2 , m N
w 1 , m 1 = w 1 , m 2 = = w 1 , m N
w 2 , m 1 = w 2 , m 2 = = w 2 , m N = 0
The derivatives of the traction equilibrium Equations (18)–(21) with respect to a geometrical parameter Z m are shown in Equations (39)–(42):
i = 1 N t 1 , m i c o s ( θ i ) + t 1 i c o s ( θ i ) , m p 3 , m i s i n ( θ i ) p 3 i s i n ( θ i ) , m = 0
i = 1 N t 1 , m i s i n ( θ i ) + t 1 i s i n ( θ i ) , m + p 3 , m i c o s ( θ i ) + p 3 i c o s ( θ i ) , m = 0
i = 1 N t 2 , m i = 0
i = 1 N p α , m i = 0
The derivatives of the displacement equilibrium Equations (13)–(17) with respect to plate thickness h are shown in Equations (43)–(47):
u 1 , h 1 c o s ( θ 1 ) w 3 , h 1 s i n ( θ 1 ) = u 1 , h 2 c o s ( θ 2 ) w 3 , h 2 s i n ( θ 2 ) = = u 1 , h N c o s ( θ N ) w 3 , h N s i n ( θ N )
u 1 , h 1 s i n ( θ 1 ) + w 3 , h 1 c o s ( θ 1 ) = u 1 , h 2 s i n ( θ 2 ) + w 3 , h 2 c o s ( θ 2 ) = = u 1 , h N s i n ( θ N ) + w 3 , h N c o s ( θ N )
u 2 , h 1 = u 2 , h 2 = = u 2 , h N
w 1 , h 1 = w 1 , h 2 = = w 1 , h N
w 2 , h 1 = w 2 , h 2 = = w 2 , h N = 0
The derivatives of the traction equilibrium Equations (18)–(21) with respect to a geometrical parameter Z m are shown in Equations (48)–(51):
i = 1 N t 1 , h i c o s ( θ i ) + t 1 i p 3 , h i s i n ( θ i ) = 0
i = 1 N t 1 , h i s i n ( θ i ) + p 3 , h i c o s ( θ i ) = 0
i = 1 N t 2 , h i = 0
i = 1 N p α , h i = 0

3. Numerical Examples

3.1. Numerical Example 1: L-Shaped Plate

To validate the developed Implicit Differentiation Method (IDM) for assembled plate structures, the L-shaped structure shown in Figure 7 was investigated. This structure is composed of two plates joined together. One plate is clamped at one of its edges, and the other plate is subjected to tension along one of its edges. The two plates are composed of steel with a Young’s modulus E = 206.8 GPa and Poisson’s ratio ν = 0.29 . The geometry of the structure can be defined in terms of the four parameters L, W, θ , and h, as shown in Figure 7 and detailed in Table 1.
The maximum displacement magnitude u m a g m a x occurs at the loaded edge of plate 1. To validate the IDM, the derivatives of u m a g m a x with respect to the parameters in Table 1 were calculated via three different methods:
  • Method 1: The IDM methodology described in Section 2. This method is referred to as ‘IDM-BEM’.
  • Method 2: The BEM with the forward Finite Difference Method (FDM). This method is referred to as ‘FDM-BEM’.
  • Method 3: The Finite Element Method (FEM) with the forward FDM. This method is referred to as ‘FDM-FEM’. The FEM software Abaqus FEA was used in this work.
Before a comparison can be made between these three methods, a mesh convergence study must be carried out for the BEM and the FEM. The results are shown in Figure 8. The BEM mesh is considered to have reached convergence with 80 elements, whereas the FEM mesh is considered to have reached convergence with 800 elements. These two converged meshes are shown as red markers in Figure 8.
Using the converged meshes for the BEM and the FEM, the sensitivities of u m a g m a x with respect to the parameters in Table 1 can be calculated using the three methods introduced earlier. The results are shown in Figure 9.
It is shown in Figure 9 that there is excellent agreement between the three methods ‘IDM-BEM’, ‘FDM-BEM’, and ‘FDM-FEM’. The average difference between IDM-BEM and FDM-FEM was, at most, only 1.9% for W, and only 0.11% for L. The average difference between IDM-BEM and FDM-BEM was even smaller, and was, at most, 0.61% for θ and only 0.073% for W. The excellent agreement between IDM-BEM and FDM-FEM, and between IDM-BEM and FDM-BEM, demonstrates the very high accuracy of the IDM when calculating sensitivities.

3.2. Numerical Example 2: X-Core

To further validate the developed Implicit Differentiation Method (IDM) for assembled plate structures, a significantly more complex structure was investigated with many more geometrical parameters.
The Eurofighter Typhoon has canards to increase its manoeuvrability, as shown in Figure 10. These canards are made of titanium and are manufactured using super-plastic forming and diffusion bonding, which results in a continuously connected structure. A cross-section of one of these canards is shown in Figure 11. It consists of eight X-core sections. An approximation of this cross-section is shown in Figure 12 (front-view) and in Figure 13 (side-view). The cross-section was subjected to the uniform pressure loading q 3 = 1 MPa on the top and bottom surfaces.
The canard cross-section has two planes of symmetry; these are plane x 2 x 3 in Figure 12 (front-view) and plane x 1 x 3 in Figure 13 (side-view). By using these planes of symmetry, a simplified cross-section can be obtained, as shown in Figure 14 (front-view) and Figure 15 (side-view). The geometry of this simplified cross-section can be described in terms of the six variables W 1 , W 2 , L 1 , L 2 , h 1 , h 2 and also the height H and depth D, details of which are shown in Table 2. The structure is made of titanium with a Young’s modulus E = 110 GPa and Poisson’s ratio ν = 0.3 . The structure has symmetry boundary conditions along the planes of symmetry.
The BEM model consists of 14 individual plates joined together. Symmetric boundary conditions were applied along the edges corresponding to the x 2 x 3 and x 1 x 3 planes of symmetry such that the displacement in direction x 1 is zero ( u 1 = 0 ) for edges on the plane x 2 x 3 and the displacement in direction x 2 ( u 2 = 0 ) for edges on the plane x 1 x 3 . For the node corresponding to the centre of the coordinate system (0,0,0), the vertical displacement was set to zero ( w 3 = 0 ).
A displacement contour plot of the X-core structure is shown in Figure 16. The maximum displacement magnitude u m a g m a x occurs at the middle of the top and bottom plates. To validate the IDM, the sensitivities of u m a g m a x with respect to the parameters in Table 1 were calculated via the three different methods introduced in Section 3.1: ‘IDM-BEM’, ‘FDM-BEM’, and ‘FDM-FEM’.
Before a comparison can be made between these three methods, a mesh convergence study must be carried out for the BEM and the FEM. The results of this mesh convergence study are shown in Figure 17. The BEM mesh is considered to have reached convergence with 1120 elements, whereas the FEM mesh is considered to have reached convergence with 19,800 elements. These two converged meshes are shown as red markers in Figure 17. The converged BEM mesh is shown in Figure 18 and it is composed of 14 individual plates joined together.
Using the converged meshes for the BEM and the FEM, the sensitivities of u m a g m a x with respect to the parameters in Table 2 can be calculated using the three methods introduced earlier. The results are shown in Figure 19.
It is shown in Figure 19 that there is excellent agreement between the three methods ‘IDM-BEM’, ‘FDM-BEM’, and ‘FDM-FEM’. The average difference between IDM-BEM and FDM-FEM was, at most, only 1.5% for W 1 , and only 0.29% for h 1 . The average difference between IDM-BEM and FDM-BEM was even smaller, and was, at most, 0.31% for L 1 , and only 0.034% for L 2 . The excellent agreement between IDM-BEM and FDM-FEM and between IDM-BEM and FDM-BEM demonstrates the very high accuracy of the IDM when calculating sensitivities for this very complex example.
To compare the computational efficiency of the three methods, they were each used to calculate 100 sensitivities of the X-core structure. This comparison is shown in Table 3. It is shown that the methods ‘FDM-FEM’ and ‘FDM-BEM’ require a similar amount of CPU time (1120 s vs. 1149 s) to calculate 100 sensitivities. The method ‘IDM-BEM’ is significant more computationally efficient than the other two methods—it is 22% faster than ‘FDM-FEM’ and 24% faster than ’FDM-BEM’.

3.3. Numerical Example 3: Shape Optimisation of the X-Core

The previous two numerical examples validated the newly developed IDM for assembled structures with a simple example and with a significantly more complex example. To demonstrate how the IDM could be used in practice, it is used in this example to optimise the shape of the X-core structure investigated in the previous numerical example.
The shape of the X-core structure was optimised such that its mass was minimised while the maximum displacement magnitude u m a g m a x was less than 0.1 mm. The height H and depth D of the X-core were also fixed at 21.5 mm and 25 mm, respectively. The total width of the X-core W t o t = 3 W 1 + 2 W 2 was also fixed at 129 mm. The optimisation problem is shown in Equation (52):
Minimise : M a s s ( d ) Subject   to : d L d d U , d R n d u m a g m a x 0.1 mm H = 21.5 mm D = 25 mm W t o t = 3 W 1 + 2 W 2 = 129 mm
where the design variables are d = [ W 1 , W 2 , L 1 , L 2 , h 1 , h 2 ] , and n d = 6 . The initial design is d 0 = [ 15 , 42 , 10 , 6 , 2.8 , 0.8 ] . The lower and upper bounds are: d L = 0.75 × d 0 and d U = 1.25 × d 0 . The mass can be calculated using the design variables d and the density of titanium ρ = 4.41 × 10 3 g/mm3.
The IDM was used to solve the above optimisation problem. The optimisation procedure was stopped once the relative change in mass between iterations fell to below 0.1% and u m a g m a x below or equal to 0.1 mm with a tolerance of 0.001 mm. The results of this optimisation procedure can be seen in Figure 20. It is shown in this Figure that the optimisation procedure successfully minimised the mass of the X-core structure—the mass was initially 27 g and was reduced to 23.5 g after 62 iterations. The maximum displacement magnitude u m a g m a x was initially 0.139 mm and was reduced to 0.101 mm after 62 iterations.
The shape and the values for the shape parameters of the X-core structure before and after the shape optimisation procedure are shown in Figure 21 and Table 4. It is shown that the height H and depth D of the X-core structure are unchanged, which was expected. The width of the X-sections have been made wider and the thickness of the inner plates has also been increased, whereas the thickness of the outer plates has been reduced. This had the effect of reducing the mass of the X-core structure and its maximum displacement magnitude u m a g m a x .
In summary, this numerical example demonstrated the effectiveness of the IDM when used in practice for optimising the shape of a complex assembled structure, such as an X-core structure with many geometric design variables.

4. Conclusions

In conclusion, a novel methodology was presented for conducting sensitivity analyses of assembled plate structures using the Boundary Element Method (BEM). The main novelty of this work was that the exact implicit derivatives of the BEM formulations for assembled plate structures were derived for the first time and incorporated into a newly developed Implicit Differentiation Method (IDM) for performing structural sensitivity analyses—enabling sensitivity analyses to be conducted for more complex and realistic structures in a more accurate, robust, and efficient manner. Several numerical examples were presented to validate the derived exact implicit derivatives and to demonstrate how they could be used for a potential application involving the shape optimisation of a complex X-core structure from the canard of a Eurofighter Typhoon fighter jet. The values of six geometrical parameters were optimised to minimise the maximum displacement in the X-core. The newly developed IDM was found to be very similar in terms of accuracy to the Finite Difference Method (FDM), demonstrating percentage differences less than 2%. The IDM was also found to be over 20% more computationally efficient than the FDM. Furthermore, unlike the FDM, the IDM does not require a step-size, meaning that, compared to the FDM, the IDM is significantly more robust. Potential applications include the shape optimisation or structural reliability analysis of complex plate structures.
Future work will aim to develop an IDM for assembled shell structures, enabling the sensitivity analysis of complex aircraft structures such as fuselage panels or wing sections. Sensitivity analyses involving buckling and large-deflections will also be investigated to expand the range of possible applications.

Author Contributions

Conceptualisation, L.M., V.M., Z.S.-K. and F.M.H.A.; formal analysis, L.M.; funding acquisition, V.M.; investigation, L.M.; methodology, L.M., V.M. and F.M.H.A.; project administration, V.M., Z.S.-K. and F.M.H.A.; software, L.M.; supervision, V.M.; validation, L.M. and V.M.; visualisation, L.M.; writing—original draft, L.M. and V.M.; writing—review and editing, V.M., Z.S.-K. and F.M.H.A. All authors have read and agreed to the published version of the manuscript.

Funding

This project has received funding from the Clean Sky 2 Joint Undertaking (JU) under grant agreement No 864154 Project MASCOT. The JU receives support from the European Union’s Horizon 2020 research and innovation programme and the Clean Sky 2 JU members other than the Union.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Aliabadi, M.H. The Boundary Element Method: Applications in Solids and Structures; John Wiley and Sons: Chichester, UK, 2002; Volume 2. [Google Scholar]
  2. Tafreshi, A. Shape optimization of two-dimensional anisotropic structures using the boundary element method. J. Strain Anal. Eng. Des. 2003, 38, 219–232. [Google Scholar] [CrossRef] [Green Version]
  3. Abe, K.; Kazama, S.; Koro, K. A boundary element approach for topology optimization problem using the level set method. Commun. Numer. Methods Eng. 2006, 23, 405–416. [Google Scholar] [CrossRef]
  4. Canelas, A.; Herskovits, J.; Telles, J.C.F. Shape optimization using the boundary element method and a SAND interior point algorithm for constrained optimization. Comput. Struct. 2008, 86, 1517–1526. [Google Scholar] [CrossRef]
  5. Ullah, B.; Trevelyan, J. Correlation between hole insertion criteria in a boundary element and level set based topology optimisation method. Eng. Anal. Bound. Elem. 2013, 37, 1457–1470. [Google Scholar] [CrossRef] [Green Version]
  6. Chang, Y.; Cheng, H.; Chiu, M.; Chien, Y. Shape Optimisation of Multi-Chamber Acoustical Plenums Using BEM, Neural Networks, and GA Method. Arch. Acoust. 2016, 41, 43–53. [Google Scholar] [CrossRef] [Green Version]
  7. Ullah, B.; Trevelyan, J. A boundary element and level set based topology optimisation using sensitivity analysis. Eng. Anal. Bound. Elem. 2016, 70, 80–98. [Google Scholar] [CrossRef] [Green Version]
  8. Liu, C.; Chen, L.; Zhao, W.; Chen, H. Shape optimization of sound barrier using an isogeometric fast multipole boundary element method in two dimensions. Eng. Anal. Bound. Elem. 2017, 85, 142–157. [Google Scholar] [CrossRef]
  9. Ullah, B.; Trevelyan, J.; Sirajul, I. A boundary element and level set based bi-directional evolutionary structural optimisation with a volume constraint. Eng. Anal. Bound. Elem. 2017, 80, 152–161. [Google Scholar] [CrossRef] [Green Version]
  10. Takahashi, T.; Yamamoto, T.; Shimba, Y.; Isakari, H.; Matsumoto, T. A framework of shape optimisation based on the isogeometric boundary element method toward designing thin-silicon photovoltaic devices. Eng. Comput. 2018, 35, 423–449. [Google Scholar] [CrossRef]
  11. Matsushima, K.; Isakari, H.; Takahashi, T.; Matsumoto, T. A topology optimisation of composite elastic metamaterial slabs based on the manipulation of far-field behaviours. Struct. Multidiscip. Optim. 2020, 63, 231–243. [Google Scholar] [CrossRef]
  12. Maduramuthu, P.; Fenner, R.T. Three-dimensional shape design optimization of holes and cavities using the boundary element method. J. Strain Anal. Eng. Des. 2004, 39, 87–98. [Google Scholar] [CrossRef]
  13. Bandara, K.; Cirak, F.; Of, G.; Steinbach, O.; Zapletal, J. Boundary element based multiresolution shape optimisation in electrostatics. J. Comput. Phys. 2015, 297, 584–598. [Google Scholar] [CrossRef] [Green Version]
  14. Ullah, B.; Trevelyan, J.; Ivrissimtzis, I. A three-dimensional implementation of the boundary element and level set based structural optimisation. Eng. Anal. Bound. Elem. 2015, 58, 176–194. [Google Scholar] [CrossRef] [Green Version]
  15. Chen, L.L.; Lian, H.; Liu, Z.; Chen, H.B.; Atroshchenko, E.; Bordas, S.P.A. Structural shape optimization of three dimensional acoustic problems with isogeometric boundary element methods. Comput. Methods Appl. Mech. Eng. 2019, 355, 926–951. [Google Scholar] [CrossRef]
  16. Gaggero, S.; Vernengo, G.; Villa, D.; Bonfiglio, L. A reduced order approach for optimal design of efficient marine propellers. Ships Offshore Struct. 2019, 15, 200–214. [Google Scholar] [CrossRef]
  17. Li, S.; Trevelyan, J.; Wu, Z.; Lian, H.; Wang, D.; Zhang, W. An adaptive SVD–Krylov reduced order model for surrogate based structural shape optimization through isogeometric boundary element method. Comput. Methods Appl. Mech. Eng. 2019, 349, 312–338. [Google Scholar] [CrossRef] [Green Version]
  18. Morse, L.; Sharif Khodaei, Z.; Aliabadi, M.H. A multi-fidelity boundary element method for structural reliability analysis with higher-order sensitivities. Eng. Anal. Bound. Elem. 2019, 104, 183–196. [Google Scholar] [CrossRef]
  19. Morse, L.; Mallardo, V.; Aliabadi, F.M.H. Manufacturing cost and reliability-based shape optimization of plate structures. Int. J. Numer. Methods Eng. 2022, 123, 2189–2213. [Google Scholar] [CrossRef]
  20. Mallardo, V.; Aliabadi, M.H. A BEM sensitivity and shape identification analysis for acoustic scattering in fluid–solid problems. Int. J. Numer. Methods Eng. 1998, 41, 1527–1541. [Google Scholar] [CrossRef]
  21. Mallardo, V.; Alessandri, C. Inverse problems in the presence of inclusions and unilateral constraints a boundary element approach. Comput. Mech. 2000, 26, 571–581. [Google Scholar] [CrossRef]
  22. Morse, L.; Sharif Khodaei, Z.; Aliabadi, M.H. A dual boundary element based implicit differentiation method for determining stress intensity factor sensitivities for plate bending problems. Eng. Anal. Bound. Elem. 2019, 106, 412–426. [Google Scholar] [CrossRef]
  23. Brancati, A.; Aliabadi, M.H.; Mallardo, V. A BEM sensitivity formulation for three-dimensional active noise control. Int. J. Numer. Methods Eng. 2012, 90, 1183–1206. [Google Scholar] [CrossRef]
  24. Dirgantara, T.; Aliabadi, M.H. Crack Growth analysis of plates Loaded by bending and tension using dual boundary element method. Int. J. Fract. 1999, 105, 27–47. [Google Scholar]
Figure 1. Plate geometry. Reprinted with permission from Ref. [1], 2002, Ferri M.H. Aliabadi.
Figure 1. Plate geometry. Reprinted with permission from Ref. [1], 2002, Ferri M.H. Aliabadi.
Aerospace 09 00381 g001
Figure 2. Sign convections for the displacements, rotations, tractions, and moments.
Figure 2. Sign convections for the displacements, rotations, tractions, and moments.
Aerospace 09 00381 g002
Figure 3. An assembled plate structure consisting of two plates joined together along a junction line. Each plate has its own local coordinate system ( x 1 i , x 2 i , x 3 i ). The edges of each plate are labelled 1–4.
Figure 3. An assembled plate structure consisting of two plates joined together along a junction line. Each plate has its own local coordinate system ( x 1 i , x 2 i , x 3 i ). The edges of each plate are labelled 1–4.
Aerospace 09 00381 g003
Figure 4. The A matrix for the assembled plate structure in Figure 3.
Figure 4. The A matrix for the assembled plate structure in Figure 3.
Aerospace 09 00381 g004
Figure 5. The A , m matrix for the assembled plate structure in Figure 3.
Figure 5. The A , m matrix for the assembled plate structure in Figure 3.
Aerospace 09 00381 g005
Figure 6. The A , h matrix for the assembled plate structure in Figure 3.
Figure 6. The A , h matrix for the assembled plate structure in Figure 3.
Aerospace 09 00381 g006
Figure 7. The L-shaped plate investigated in numerical example 1.
Figure 7. The L-shaped plate investigated in numerical example 1.
Aerospace 09 00381 g007
Figure 8. Mesh convergence with the Boundary Element Method (BEM) and the Finite Element Method (FEM). The BEM and FEM meshes considered to have reached convergence are shown as red markers.
Figure 8. Mesh convergence with the Boundary Element Method (BEM) and the Finite Element Method (FEM). The BEM and FEM meshes considered to have reached convergence are shown as red markers.
Aerospace 09 00381 g008
Figure 9. The sensitivities of the maximum displacement magnitude u m a g m a x with respect to (a) length L, (b) width W, (c) plate angle θ , and (d) plate thickness h, for the three methods ‘IDM-BEM’, ‘FDM-BEM’, and ‘FDM-FEM’.
Figure 9. The sensitivities of the maximum displacement magnitude u m a g m a x with respect to (a) length L, (b) width W, (c) plate angle θ , and (d) plate thickness h, for the three methods ‘IDM-BEM’, ‘FDM-BEM’, and ‘FDM-FEM’.
Aerospace 09 00381 g009aAerospace 09 00381 g009b
Figure 10. A canard of the Eurofighter Typhoon (highlighted with a red circle).
Figure 10. A canard of the Eurofighter Typhoon (highlighted with a red circle).
Aerospace 09 00381 g010
Figure 11. A cross-section of one of the canards from the Eurofighter Typhoon.
Figure 11. A cross-section of one of the canards from the Eurofighter Typhoon.
Aerospace 09 00381 g011
Figure 12. An approximation of the canard cross-section (front-view).
Figure 12. An approximation of the canard cross-section (front-view).
Aerospace 09 00381 g012
Figure 13. An approximation of the canard cross-section (side-view).
Figure 13. An approximation of the canard cross-section (side-view).
Aerospace 09 00381 g013
Figure 14. A simplified version of the canard cross-section (front-view).
Figure 14. A simplified version of the canard cross-section (front-view).
Aerospace 09 00381 g014
Figure 15. A simplified version of the canard cross-section (side-view).
Figure 15. A simplified version of the canard cross-section (side-view).
Aerospace 09 00381 g015
Figure 16. Displacement contour plot.
Figure 16. Displacement contour plot.
Aerospace 09 00381 g016
Figure 17. Mesh convergence with the Boundary Element Method (BEM) and the Finite Element Method (FEM). The BEM and FEM meshes considered to have reached convergence are shown as red markers.
Figure 17. Mesh convergence with the Boundary Element Method (BEM) and the Finite Element Method (FEM). The BEM and FEM meshes considered to have reached convergence are shown as red markers.
Aerospace 09 00381 g017
Figure 18. The converged BEM mesh. The 14 plates are labelled P1–P14.
Figure 18. The converged BEM mesh. The 14 plates are labelled P1–P14.
Aerospace 09 00381 g018
Figure 19. The sensitivities of the maximum displacement magnitude u m a g m a x with respect to (a) W 1 , (b) W 2 , (c) L 1 , (d) L 2 , (e) h 1 , and (f) h 2 for the three methods ‘IDM-BEM’, ‘FDM-BEM’, and ‘FDM-FEM’.
Figure 19. The sensitivities of the maximum displacement magnitude u m a g m a x with respect to (a) W 1 , (b) W 2 , (c) L 1 , (d) L 2 , (e) h 1 , and (f) h 2 for the three methods ‘IDM-BEM’, ‘FDM-BEM’, and ‘FDM-FEM’.
Aerospace 09 00381 g019aAerospace 09 00381 g019bAerospace 09 00381 g019c
Figure 20. Results of the shape optimisation of the X-core structure.
Figure 20. Results of the shape optimisation of the X-core structure.
Aerospace 09 00381 g020
Figure 21. The X-core structure before the IDM shape optimisation (top). The X-core structure after the IDM shape optimisation (bottom).
Figure 21. The X-core structure before the IDM shape optimisation (top). The X-core structure after the IDM shape optimisation (bottom).
Aerospace 09 00381 g021
Table 1. The parameters of the L-shaped plate investigated in numerical example 1.
Table 1. The parameters of the L-shaped plate investigated in numerical example 1.
ParameterDescriptionValue
LLength of plate 1 and 2100 mm
WWidth of plate 1 and 250 mm
θ Angle of plate 250
hThickness of plate 1 and 210 mm
Table 2. The parameters of the X-core structure investigated in numerical example 2.
Table 2. The parameters of the X-core structure investigated in numerical example 2.
ParameterDescriptionValue
W 1 Width of outer long plate15 mm
W 2 Width of outer short plate42 mm
L 1 Length of inner inclined plate10 mm
L 2 Length of inner horizontal plate6 mm
h 1 Thickness of outer plate2.8 mm
h 2 Thickness of inner plate0.8 mm
HHeight21.5 mm
DDepth25 mm
Table 3. CPU time to calculate 100 sensitivities using each of the three methods.
Table 3. CPU time to calculate 100 sensitivities using each of the three methods.
MethodCPU Time (s)
FDM-FEM1120
FDM-BEM1149
IDM-BEM870
Table 4. The parameters of the X-core structure before and after shape optimisation.
Table 4. The parameters of the X-core structure before and after shape optimisation.
Shape Parameter (mm)
Design W 1 W 2 L 1 L 2 h 1 h 2 H D W tot Mass (g) u mag max (mm)
Before15.042.010.06.02.80.8021.52512927.00.139
After18.836.410.87.52.10.9821.52512923.50.101
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Morse, L.; Mallardo, V.; Sharif-Khodaei, Z.; Aliabadi, F.M.H. Shape Optimisation of Assembled Plate Structures with the Boundary Element Method. Aerospace 2022, 9, 381. https://doi.org/10.3390/aerospace9070381

AMA Style

Morse L, Mallardo V, Sharif-Khodaei Z, Aliabadi FMH. Shape Optimisation of Assembled Plate Structures with the Boundary Element Method. Aerospace. 2022; 9(7):381. https://doi.org/10.3390/aerospace9070381

Chicago/Turabian Style

Morse, Llewellyn, Vincenzo Mallardo, Zahra Sharif-Khodaei, and Ferri M.H. Aliabadi. 2022. "Shape Optimisation of Assembled Plate Structures with the Boundary Element Method" Aerospace 9, no. 7: 381. https://doi.org/10.3390/aerospace9070381

APA Style

Morse, L., Mallardo, V., Sharif-Khodaei, Z., & Aliabadi, F. M. H. (2022). Shape Optimisation of Assembled Plate Structures with the Boundary Element Method. Aerospace, 9(7), 381. https://doi.org/10.3390/aerospace9070381

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop