Shape Sensing of Plate and Shell Structures Undergoing Large Displacements Using the Inverse Finite Element Method

+e inverse Finite Element Method (iFEM) is applied to reconstruct the displacement field of a shell structure which undergoes large deformations using discreet strain measurements as the prescribed data. +e iFEM computations are carried out using an incremental procedure where at each load step, the incremental strains are used to evaluate the incremental displacements which in turn update the geometry of the deformed structure. +e efficacy of the proposed approach to predict large displacements is examined using two case studies involving a cantilevered wing-shaped plate and a clamped plate. +e incremental iFEM procedure is demonstrated to be sufficiently accurate in terms of reproducing the correct nonlinear character of the loaddisplacement curve even when a reduced number of strain sensors is used. +erefore, this approach may have important implications for real-time monitoring of aerospace structures that undergo large displacements.


Introduction
Efficient structural health monitoring (SHM) systems are enabling technologies for the maintenance and control of future aerospace vehicles.When implemented, such systems provide real-time monitoring of the vehicle's structural integrity.
is technology would thus enable improved structural safety, reduced maintenance operations, and significant cost savings.A SHM system uses a set of strain sensors embedded on the surface of the structure.Using real-time sensor measurements, the structural deformations need to be reconstructed, thus requiring the solution of an inverse problem known as "shape sensing".
Existing shape-sensing methods can be classified into four categories based on different approaches [1]: (1) numerical integration of experimental strains, (2) use of global or piecewise continuous basis functions to approximate the displacement field, (3) application of neural networks, and (4) use of variational approaches.Most of the methods based on the integration of measured strains deal with beam problems and make use of the classical beam equations [2,3].For the full-field estimation of the deformed shape of more complex two-dimensional structures, basis-function approaches utilize an a priori set of spatial functions and unknown coefficients to fit discretely measured strains [4][5][6][7][8].A major drawback of the methods based on neural networks is that their accuracy depends on the choice of the training load cases [9].
For many large-scale practical applications, the variational methods appear to be most attractive.In these methods, a suitable error functional is minimized over the entire problem domain, by comparing the estimated and measured strains, commonly in the least-squares sense.One such method is the inverse Finite Element Method (iFEM), advanced by Tessler and Spangler [10][11][12] and initially developed for shear-deformable plate and shell structures.e method employs a weighted least-square variational formulation by discretizing the problem domain using C 0continuous finite elements.
e approach is suitable for modelling complex geometries and allows both static and dynamic responses to be reconstructed without any prior knowledge of material properties or the loading conditions of the structure.e inverse Finite Element Method was adapted by Gherlone [13] for shape sensing of truss, beam, and frame structures.Since the initial introduction of iFEM, the approach has been successfully applied on a variety of simple and complex structures, where the focus has been primarily on the small deformation, linear load-displacement response.e existing iFEM formulations use the simplifying assumption of "small displacements"; hence, applying these methods to large-displacement problems will generally result in significant errors in the prediction of the deformed shape.
e first real attempt at addressing the problem of geometric nonlinearity was made by Tessler et al. [14].e method involves a modification of the standard linear iFEM, where the linear iFEM is used over a number of load steps and the incremental variation in deformation over each step is summed to determine the final deformed shape of the structure.is method bears resemblance to the nonlinear finite element method, which also uses a series of linear displacement increments to recreate the nonlinear nature of the load-deflection curve.Initial results obtained by employing this technique to predict deformations in a complex wing profile were seen as promising.
Shape-sensing techniques capable of predicting large deformations have a wide range of applications for future aircraft and spacecraft missions.Many aircraft components, such as the fuselage skin and certain parts of the wing, are subjected to compressive loading during aircraft operation.An aircraft structural health monitoring system capable of reconstructing the nonlinear structural response should facilitate improved predictions of the onset of buckling and related failure.Such technologies can also be applied to shape morphing structures in aircraft and membrane space structures such as solar sails, spacecraft deorbiters, and sun shields (as in the James Webb Space Telescope), where large deformations are normally observed.e predicted shape can be used as feedback for the actuation system to facilitate control over the deformed structural shape.
is paper aims to describe the basic methodology and procedure of the linear incremental iFEM formulation and to demonstrate the method's suitability to reconstruct large deformations in plate and shell structures.In Section 2, the formulation of the standard linear iFEM is briefly described, followed by the linear incremental iFEM procedure.In Section 3, the iFEM approach for predicting large deformations is demonstrated on two example problems, one involving a cantilevered wing-shaped plate and the other, a clamped square plate under transverse loading.e effect of the number of strain sensors on the accuracy of the reconstructed displacement field is also investigated.Finally, concluding remarks and future directions of this research are presented in Section 4.

The Inverse Finite Element Method
In this section, a brief review of the iFEM approach for applications to beam, plate, and shell structures is presented.
e displacement field experienced by the structure is completely defined by the vector (u x (x), u y (x), u z (x)) in the Cartesian coordinate system x ≡ (x, y, z).Depending on the particular structural theory, the displacement vector can be expressed in terms of a set of kinematic variables u(x).ese kinematic variables can be interpolated within each finite element of the discretization by a set of element-shape functions: where N is the shape-function matrix and q e the nodal degrees-of-freedom.
Corresponding to the adopted structural theory, the strain field is completely described by a set of K independent quantities, ε ≡ ε k  (k � 1, . . ., K), usually referred to as strain measures.Using Equation (1) and the linear straindisplacement relations, the strain measures can be expressed in terms of the nodal degrees-of-freedom: where the matrix L contains the derivatives of the shape functions.e iFEM approach is based on the minimization of a functional defined as the least-squares error between the analytic strain measures (Equation ( 2)) and the corresponding experimental strain measures, e latter are evaluated at n discrete locations using strain gauges, strain rosettes, or other types of strain sensors.For a single element, the error functional, Φ e , is defined as the sum of the products of the least-squares component, Φ e k , referred to the k − th strain measure, and the corresponding weighting coefficient, w e k , Φ e u e  ≡ with Φ e k defined as e weighting coefficients w e k are the products of dimensional parameters and dimensionless coefficients, λ e k .e dimensional parameters guarantee that the terms in Equation ( 3) have the same mathematical units [11][12][13], whereas the λ e k coefficients are assigned either high (1, unity) or relatively low values (10 −4 ) to enforce a stronger or weaker correlation between the analytic strain measures and their experimentally measured counterparts.
e inverse-element error functional, Φ e , is minimized with respect to the unknown nodal degrees-of-freedom, q e , leading to the element matrix equation a e q e � b e .Applying the appropriate coordinate transformations, the element contributions are assembled into a global system of equations, and upon enforcing problem-dependent displacement boundary conditions, the global equations in matrix form are Aq � b. (5)

Shock and Vibration
In Equation ( 5), the matrix A depends on the shape functions and strain-sensor locations, whereas the vector b is a function of the experimental strain measures.Equation ( 5) can be readily solved for the unknown nodal degrees-offreedom, q by inverting the matrix A.
us, for small displacements, the strain-sensor locations are assumed to remain unchanged, and A is inverted only once, whereas the vector b is updated at each strain-data acquisition increment.
Since only strain-displacement relations are used in the de nition of Φ e , iFEM does not require the knowledge of any material properties or the applied loading.us, it is applicable for both static and dynamic conditions, without requiring inertial or damping material properties.
2.1.Plate Formulation.Consider a plate of thickness 2t for which the (x, y) plane of the Cartesian coordinate system represents the midplane (Figure 1).e components of the displacement vector are given according to Mindlin plate theory, as where the ve kinematic variables are u ≡ u, v, w, θ x , θ y T ; u and v are average uniform displacements in the x and y directions, respectively; w is the average transverse deection; θ x and θ y are rotations of the normal about the negative x and positive y axes, respectively.e strain eld is identi ed by eight strain measures, three membrane (k 1−3), three bending (k 4−6), and two transverse shear (k 7, 8), ε k (k 8), given by If Reference [12], a three-node inverse plate/shell element, iMIN3, has been developed and will be applied in the present analysis.
e element is formulated using C 0continuous anisoparametric shape functions to avoid shear locking.
e de ection variable is interpolated with a quadratic polynomial, whereas the other four kinematic variables vary linearly over the element.e strain measures de ned in Equation ( 7) can then be expressed in terms of the nodal degrees-of-freedom, as in Equation ( 2), and the error functional can be written as in Equations ( 3 e membrane and bending curvature strain measures, e ε (i) and k ε (i) , can be evaluated at the location (x (i) , y (i) ) from surface strain measurements on the top (+) and bottom (−) surfaces (Figure 2) as e ε (i) e transverse shear strain measures, g, cannot be evaluated experimentally.us, for the corresponding terms in the functional (Equation (3)), the following expression is adopted (k 7, 8): where A e is the element area.Moreover, λ e k (k 7, 8) are set to small values, e.g., 10 −4 , relative to 1 (unity) which is used for the λ e k (k 1, . . ., 6) coe cients.When only few strain sensors are available, some elements may not have any strain data.For the elements without the strain data, the norms in Equation ( 9) are used for all eight strain measures, where the λ e k (k 1, . . ., 8) coe cients are assigned small values, e.g., 10 −4 , relative to 1 (unity) which is used for the elements that possess strain data.

iFEM Incremental Algorithm for Nonlinear Deformations.
Herein, it is assumed that the structure under consideration is instrumented with a high-speed strain-sensor measurement system which is capable of measuring strain levels in

Shock and Vibration
real time as the structure undergoes deformations under quasistatic or dynamic loading.e standard linear iFEM formulation can be applied at each incremental load level to reconstruct the current deformed shape, using the current strain increments, and then updating the geometry of the iFEM mesh and repeating the shape-sensing analysis based on the next strain increment.us, geometrically nonlinear deformations ("large displacements") can be reconstructed in real time as the strain history becomes available to iFEM from a strain-sensor measurement system [13].
(i) At each load increment, i, evaluate (experimentally or with a nonlinear direct FEM analysis) the incremental section strains that represent measured strain increments.(ii) Perform iFEM analysis using the strain increments to obtain the nodal-degrees-of-freedom, q i .(iii) Update the structural geometry (i.e., the iFEM nodal coordinate positions) using the iFEM determined displacements, i.e., x i+1 � x i + N(x i )q i (iv) Update the orientation of the measured strains corresponding to the current geometry, x i+1 (v) Repeat the above-stated incremental procedure until the strain history is complete.
e iFEM incremental procedure just outlined is thus in complete agreement with the standard incremental procedures used for geometrically nonlinear analysis of the direct FEM.Unlike the direct FEM, the iFEM incremental procedure does not invoke equilibrium iterations at different load increments.
is is because iFEM has no access to equilibrium equations.Instead, the iFEM equations smooth the measured strain data and integrate the straindisplacement relations.Since the in-situ measured strains are constantly updated at each load level, the straindisplacement relations are constantly updated and corrected to the current state of deformation.us, within this incremental formulation, the strain-level updates provide the necessary corrections at each load increment.Consequently, as will be demonstrated in Section 3, highly accurate geometrically nonlinear deformations can be reconstructed by this simple and computationally efficient incremental iFEM procedure.

Numerical Examples
e application of the incremental iFEM method is presented using two example problems for elastostatic deformation of plates subjected to transverse loading.Nonlinear direct FE models are used both to provide the strain measurements as input data and displacements as reference results for the iFEM analysis.e FEM analyses are carried out using the finite-strain three-node S3R shell element in ABAQUS/Standard 6.13 [15].e iFEM reconstruction is performed using the same triangular mesh, where the strain sensors are positioned on both the top and bottom surfaces of the iMIN3 elements.
e accuracy of the iFEM prediction is then assessed by calculating the Root Mean Square (RMS) of the nodal transverse deflection error with respect to the direct FEM results, defined as where m indicates the total number of nodes, which is the same for both the direct and the inverse analyses.
e plate is subjected to a uniformly distributed transverse pressure providing a resultant load that corresponds to 30 times the plate's own weight.Figure 3(b) shows the direct and inverse triangular mesh.Single-strain component measurements, simulating linear strain gauges, are provided to a subset of inverse elements of the iFEM discretization.In particular, 44 elements have linear strain gauges positioned on the top and bottom surfaces (Figure 3(b)).For all the inverse elements, λ e k (k � 7, 8) � 10− 4 has been used for the contributions to the error functional due to shear strain measures, whereas λ e k � 10 −4 has been used for all of the contributions related to "non-measured" strain components or to elements without strain data.
Figure 4 shows the comparison between the loaddeflection curves as obtained using the nonlinear direct FEM and iFEM analyses.e maximum deflection experienced at the end of the loading process is comparable to the wing span, and the direct solution exhibits significant nonlinearity in the load-deflection curve.By using the standard linear iFEM procedure based on the final value of the "measured" strains, the reconstructed load-deflection response cannot reproduce the correct structural response.When adopting the incremental iFEM procedure, the loaddeflection curve is more accurate (Figure 5 shows that the RMS percent error never exceeds 1%).
Figure 6 compares the deformed shapes of the plate obtained using the nonlinear direct FEM analysis and the linear incremental iFEM reconstruction.
In the second example, an isotropic square plate which is clamped along the four edges is subjected to a static pressure load of magnitude 5 MPa on the top surface.e plate's side span is 5 m and the thickness is 10 mm. e elastic modulus and Poisson's ratio are, respectively, E � 200000 MPa and v � 0.2.e iFEM analysis used the same triangular mesh as the direct FEM analysis, as depicted in Figure 7, and was supplied with a complete set of strain values at the element centroids.
Figure 8 shows the load-deflection curve for the transverse deflection at the center of the square plate.e linear iFEM was employed incrementally across a total of 21 load steps.e resulting highly nonlinear load-deflection curve is primarily due to stress stiffening.Even though the straindisplacement relations are kept linear at any load step, the gradual updating of the deformed geometry over a number of steps allows the incremental iFEM procedure to reproduce the highly accurate nonlinear deflections.In 4 Shock and Vibration comparison, the standard single-step linear iFEM application is represented by a straight line from the initial to the nal de ection point, completely missing any intermediate nonlinear response.Note that the load-de ection curve of the nonlinear direct FEM analysis matches closely with the incremental iFEM results.
e results are also examined further by computing the RMS error at each load step.
Figure 9 depicts the RMS error of the maximum deection for each load step.roughout the loading cycle, the linear incremental iFEM procedure results in the RMS error of less than 2%.In contrast, the linear iFEM procedure results in substantially higher RMS error, reaching values of around 40% during the loading process.
Figure 10 shows the nal deformed shape obtained by ABAQUS and the one predicted by the iterative iFEM procedure.
To ascertain the e ect of the number of strain sensors and their locations on the predictive capability of the linear incremental iFEM procedure, a parametric study is performed for the clamped square plate problem.Six di erent strain-sensor con gurations are selected (refer to Figure 11) to perform the parametric study.Each strain sensor is indicated by a red circle and shows the elements in which the strains are "measured" at the element centroids, on both top and bottom surfaces.e elements without the strain sensors are formulated according to Equation (9), using λ e k 10 −4 (k 1, . . ., 8). e total number of strain sensors used in each conguration is di erent, but care is taken to ensure that each  Shock and Vibration con guration has a uniform distribution of strain sensors throughout the surface area of the plate to avoid any errors due to highly localised placement of sensors.Each con guration is characterised using the parameter β, which is de ned as the average number of strain sensors in a 1 m 2 area of the plate.e parameter β is proposed with the intention of nding the minimum sensor density required for accurate shape sensing of nonlinear de ection of the considered plate structure.
e results from the linear incremental iFEM analysis are shown in Figure 12.As expected, higher values of β result in lower RMS errors.As the value of β increases, the RMS errors reduces and is seen to saturate starting at β 4. For β values as low as 2, the error magnitude is seen to be around 2%, which is not a signi cant variation from the minimum error value of 1.61%.ese results indicate that even when using only 50 strain sensors, i.e., one quarter of the number of elements, the linear incremental iFEM procedure is capable of reproducing the nonlinear displacement eld with good accuracy.is serves as further proof of the robustness of the linear incremental iFEM procedure.

Conclusions
e inverse problem of reconstructing the displacement eld from the discrete strain-sensor data for plate and shell structures undergoing large deformations is addressed using the inverse Finite Element Method.e formulation is based on an incremental procedure, where at each load step, iFEM is applied using the incremental strains to obtain the   Shock and Vibration incremental displacements that are subsequently used to update the current deformed geometry.
To demonstrate the accuracy of the incremental iFEM approach, two numerical case studies have been carried out: (1) a cantilevered wing-shaped plate and (2) a clamped square plate subjected to transverse loading.To provide the incremental strain measurements as input data and displacements as reference data for the iFEM analysis, geometrically nonlinear direct FEM models were used.e iFEM accuracy was assessed by calculating the RMS of the nodal transverse de ection error.
e incremental iFEM approach was demonstrated to be highly accurate and leads to precise load-displacement curves.Furthermore, the incremental iFEM procedure was demonstrated to be robust and accurate even when fairly sparse strain-sensor data is available, with many inverse elements without any strain data, and with only a single "measured" strain component used in the elements that had strain data.e results of this paper clearly indicate that the incremental iFEM approach can be readily implemented as a viable real-time structural health monitoring tool for aerospace structures undergoing nonlinear deformations.Our future research will focus on the iFEM applications to geometrically nonlinear dynamic problems.Shock and Vibration

Figure 3 :
Figure 3: Wing-shaped plate model.(a) Geometry and boundary conditions (thickness 3 mm) and (b) triangular mesh for both FEM and iFEM analysis (132 elements), 44 inverse elements with linear strain gauges oriented along two sensing lines (for a total of 88 strain measurements).

Figure 4 :Figure 5 :
Figure 4: Wing-shaped plate.Load-de ection curves corresponding to the nonlinear direct FEM and linear iFEM solutions.emaximum de ection (w max ) is located at the intersection between the tip chord and the trailing edge.

Figure 12 :
Figure12: RMS percent error using linear incremental iFEM analysis for di erent sensor arrangements.