Finite Element Analysis based on A Parametric Model by Approximating Point Clouds

Abstract

Simplified models are widely applied in finite element computations regarding mechanical and structural problems. However, the simplified model sometimes causes many deviations in the finite element analysis (FEA) of structures, especially in the non-designed structures which have undergone unknowable deformation features. Hence, a novel FEA methodology based on the parametric model by approximating three-dimensional (3D) feature data is proposed to solve this problem in the present manuscript. Many significant and effective technologies have been developed to detect 3D feature information accurately, e.g., terrestrial laser scanning (TLS), digital photogrammetry, and radar technology. In this manuscript, the parametric FEA model combines 3D point clouds from TLS and the parametric surface approximation method to generate 3D surfaces and models accurately. TLS is a popular measurement method for reliable 3D point clouds acquisition and monitoring deformations of structures with high accuracy and precision. The B-spline method is applied to approximate the measured point clouds data automatically and generate a parametric description of the structure accurately. The final target is to reduce the effects of the model description and deviations of the FEA. Both static and dynamic computations regarding a composite structure are carried out by comparing the parametric and general simplified models. The comparison of the deformation and equivalent stress of future behaviors are reflected by different models. Results indicate that the parametric model based on the TLS data is superior in the finite element computation. Therefore, it is of great significance to apply the parametric model in the FEA to compute and predict the future behavior of the structures with unknowable deformations in engineering accurately.

1. Introduction

The present manuscript proposes a novel finite element analysis (FEA) methodology which is based on the accurate parametric surface model by fitting the actual three-dimensional (3D) feature data. 3D feature data can be detected by many efficient measurement methods, e.g., terrestrial laser scanning (TLS), digital photogrammetry, and radar technology. This parametric model is realized by approximating the point clouds data from TLS which is an accurate and real-time measurement method. The final goal is to improve the accuracy of finite element (FE) computation results due to geometric model errors, especially in the analysis of structures which have undergone unknowable deformation features.

1.1. Background

It is widely recognized today that FE methods permit the analysis of engineering systems to be increasingly accurate and intelligent. The influence of the computational uncertainties regarding the system behavior leads the scientific community to recognize the importance of accuracy in engineering application problems [1]. An accurate FE computation can improve the design efficiency and the reliability of the analysis and prediction in the actual engineering applications. Engineering experiences show that the accuracy problems do not only include the assessment of the material [2,3,4,5] but also the geometry model of the engineering structures [6]. The use of the geometric constitution has been verified in the field of the virtual modelling application [7,8].

Geometric simplifications generally have to be adapted in the FE computation process. Losses of geometric information and connection inconsistencies between numerical models and true objects are obvious during the simplification. Consequently, issues related to the geometric model generation and its accuracy in the computational models are continuously gaining significance [9]. The geometric accuracy and structural optimization are still highlights and challenging topics in the field of the FEA. Therefore, accurate and fast measurement technologies play important roles in the object acquisition and geometric description. This is a reasonable and feasible solution to improve the accuracy of the geometry in this FEA process.

1.2. Geometric Model Generation Technologies

Computational model generation techniques play significant roles in the accuracy, the reliability, and usability of the FEA. Computer-aided design (CAD) describes the object in a visualization model, while computer-aided engineering (CAE) focuses on the analysis of the model designed. Simplification in the CAD model generation is commonly applied in describing actual problems, which leads to inaccuracies in the final FEA. Hence, the bottlenecks between model description in the design of CAD and the analysis of CAE still exist [10]. Therefore, structural model generation methods in CAD are powerfully developed and have raised a lot of interest of researchers in recent years, especially regarding the accurate surface geometry which is also necessary in the CAE computation. The robustness, higher accuracy, and lower simulation cost are required in the geometric model generation [9].

Unknowable deformed surfaces are commonly challenging objects when describing complex structures in 3D environments [11]. The B-splines and non-uniform rational B-splines (NURBS) methods in [12] are advanced choices and solutions in terms of processing free-form geometries, especially for complex deformed shapes.

The B-splines approximation enables the geometry surface control to be efficient by adjusting the control points and basis functions [12]. The B-spline-based approximation normally contains three main steps, which includes the parameterization of the data measured, the knot adjustment, and the control point determination. Advantages regarding the spatial approximation robustness of the B-spline method are affirmed in [13,14]. The basic theory about B-splines is discussed in Section 2.

The application of the B-spline fitting in the field of rail track model generation is introduced in [13]. Their B-spline algorithm is robust against outliers when determining the knot vector of B-spline with a mixture of Monte-Carlo methods and an evolutionary algorithm. Researchers in [15] present a B-spline-based procedure to construct a topological surface model from a point set. This reconstruction method shows obvious effectiveness in complex models with corner structures. An efficient generation method is applied in 3D human bone reconstruction based on the B-spline interpolation scheme [16]. The basic surface is the implicit surface which is interpolated to the Computed Tomography image data reliably. The B-spline surfaces are subsequently adopted by extracting the smooth curve net and nodes from the implicit surfaces. The final bone integration of the B-spline model is obtained by means of the projection of all nodes onto the continuous implicit surfaces [16].

When noisy point clouds are processed, a B-spline fitting model is calculated by least-squares through the singular value decomposition and LU methods [17]. This procedure offers a flexible way and extreme accurate results when the object contains complicated shapes, concavities, and singularities. Another reconstruction method based on the B-spline approximation is proposed in [18]. The reconstruction basis is the deformed mesh which is computed from the numerical simulation. It applies the weighted displacement estimation to solve the unorganized 3D points from the calculation results. The triangular B-spline-based modelling technique is introduced in [19]. The innovation is that knots and control points are determined by means of minimizing a linear combination of interpolation and fairness functional. High accuracies of the modelling technique are ensured when applied in several examples.

A rotation-free model is developed in [20] using the NURBS theory and Kirchhoff–Love shell theory. The advantage of both NURBS and B-splines is that curvatures of the model surface can be evaluated in the exact geometric description. This is a popular and preferred method in shape optimization. When aiming at the geometric shape optimization, the NURBS method is an efficient choice [21]. Both the weights and positions of NURBS are helpful and beneficial to adjust the surface shape. This allows an accurate calculation of the surface. An optimized reconstructed NURBS model is successfully obtained. Researchers show the powerful potential of the NURBS model construction ability in a medical application [22]. Hexahedral solid mesh models based on NURBS are generated to analyze the patient-specific blood flow geometry. NURBS application also causes the link between the CAD and the polyhedral model [23]. However, the choice regarding the weight of NURBS is of great challenge when the approximated shape is complicated, especially in the complex and deformed structure. The weight part in the description of B-splines is omitted, which simplifies the approximation process.

The isogeometric analysis (IGA) concept is introduced to perform FE computation about exact geometries of solid, fluid, and structure models to analyze the problems [24] when integrating the accurate geometry design and analysis of the object. IGA constructs CAD geometries to approximate the solution fields by means of NURBS in MATLAB [24]. It bridges the gap between the design in the CAD geometry and the analysis in the FEA [25]. One of the most significant contributions described by Hughes et al. is to apply the NURBS method not only as a geometric but also an analytical discretization tool. On the basis of the IGA concept, a new algorithm is proposed with the application of bivariate B-spline spaces [26,27]. B-spline-based spaces form a De Rham diagram and they are subsequently applied to solve the Maxwell equation and eigenvalue problems in the electromagnetics aspect. Other methods concerning the IGA concept which have been developed have been made very successful by many researchers in recent decades [10,26,27,28,29,30,31]. However, patching multiple IGA parametrizations to form complex topologies is far from trivial when maintaining certain continuity requirements in the complex structures is required [25].

1.3. Geometric Information Acquisition from Sensors

Advanced sensor technologies are becoming increasingly useful tools for on-site measuring and modelling. The sensor system allows the geometric generation of the surveyed structure in its current condition to be realized in a short time. TLS is a modern measuring instrument which allows the real-time geometry to be measured rapidly with high accuracy and spatial resolution.

The detection of geometric deformation is realized by TLS which rapidly capture high precision 3D point clouds of hillside areas [32]. It is proposed to use a novel registration algorithm in the deformation detection which registers TLS stations to the unmanned aerial vehicle dense image points accurately. A method to construct accurate structure models is presented by means of TLS in [33]. This method takes the point clouds from TLS as the input, while final results reveal an efficient and accurate extraction of various structural models regarding forest trees. It is pointed out in [34] that newly available platforms based on mobile laser scanning and multispectral laser scanner make great efforts to the highly qualified acquisition of the forest geometry. 3D virtual individual geometries are created with the use of highly detailed 3D TLS data [35]. The modeling framework in [35] is highly generic, transferable, and adjustable to data acquired from other TLS sensors. TLS is now widely applied in historic building information modelling [36]. It aims at constructing a novel prototype library of parametric historic building models. A system of cross platform programs is obtained to collect historic building information. Furthermore, a real mediaeval bridge geometry is described with benefits from the advanced TLS technique to obtain more reliable results [37]. Based on the resulting geometry surveyed, the aim is to carry out the limit analysis to obtain the collapse load positions and experimental values.

The integration of different multi-sensors is reported in [38] which proposes the integration of TLS and the radar technology. A 3D radar interferometric model is obtained in [38] to facilitate the spatial interpretation of displacements affecting archaeological monuments. Finally, early health warning procedures are suggested. An accurate and efficient 3D model is described according to specific requirements and visibility analysis combining digital photogrammetry and laser scanning [39]. A practical framework for the integration of photogrammetry and TLS is reported in [40]. The proposed method obtains a decimeter-level accuracy for the generated digital surface models and maps [40]. Modern integration innovations and technologies in works published previously can also be found in other references [41,42,43]. Correspondingly, measuring technologies are developing at an increasingly advanced speed, therefore, the reliability and accuracy of 3D geometry acquisition are also increasing.

1.4. Finite Element Analysis Model with High Accuracy

The authors of [44] point out that the structural analysis in the FE computation requires highly accurate spatial models. Traditional model acquisition methods rely on manual measurements. They are prone to obtain lots of errors in the results. Therefore, geometries with unknown accuracies and many geometric simplifications are applied into numerical computation, which results in distortions of the computation results [45]. Combining TLS and FEA was proposed previously, aiming at reducing the errors in the health monitoring of structures [46]. Complex FEA models with more detailed accuracies and less errors have become available by employing remote sensing technologies, such as TLS and photogrammetry. This is also a main trend and method to obtain more accurate FE computation results.

The CLOUD2FEM procedure is presented and validated to transform 3D point clouds of complex and irregular objects to 3D FEA models [47,48]. The flowchart of this method is shown in Figure 1. The procedure includes the point cloud as a stacking of point sections from TLS. Voxel elements are constructed by transforming the points into the eight-node hexahedral finite elements required. In this way, a discretized geometry is constructed with voxel elements and ready to be used in numerical computation. If both inner and outer surfaces of the object are detected by TLS, the whole 3D numerical structure can be obtained efficiently. Detailed parts with greater curvature in the model are processed with smaller discretization to obtain smoother features. This procedure helps to improve the final accuracy of the FE computation.

A pixelated fusion method regarding solid model generation via actual points-based voxelization to carry out FEA is proposed by [49]. The flowchart can be found in Figure 2a. The input of this research is the point cloud which is segmented from light detection and ranging (LiDAR) data. The asterisk in Figure 2a means the preparation of LiDAR and detection of the building object involve multiple steps outside the scope of this reference’s scientific contribution. Most previous voxelization methods focus on the surface representations of structures. Hence, the key step of these methods must identify through which voxels the surfaces pass [49]. By contrast, as it is shown in Figure 2b, the innovation is in the voxel generation step. The voxelization based on point clouds in this research operates on the points data directly. Taking the measured points as internal points, the hexahedral elements are created around them. In this process, the surface derived is not required. The pixelated model is obtained. The resulting pixelated model maps directly to the FE discretization model to realize the whole FE computation process.

A new fusion model strategy for the FE computation is presented in [48] which integrates TLS, deviation analysis, and FE computation. The reference objects with respect to the plane, sphere, cylinder, and cone are optimized through the deviation analysis. The computation of the distance between the point clouds and reference shapes are obtained to provide local deviations. Hence, well-described geometric deformations can be obtained by this comparison in the deviation maps. A 3D FE model is constructed in this way with high accuracy.

In [50] the authors introduce an automatic method which is similar to the method in [47,48] but performs better. The mesh model transformed from point clouds is constructed by four-node tetrahedral elements. Nodes from the bottom part are projected to the horizontal plane for the easy operation of the boundary conditions. The re-topology step is carried out to obtain optimized solid meshes. The mesh model is characterized by a large reduced number of solid finite elements.

An updated numerical mesh model is proposed to obtain the knowledge of the actual structural behavior in the current damaged objects [51]. Damaged structure features are described in the mesh model. This is helpful to gather further in-depth insights into current conditions. Crack behaviors are identified in the final results [51]. This plays a critical role to guarantee the good accuracy of the FE computation.

There are two different disciplines in the development, CAD and CAE, while different model representations are requirements [52]. Consequently, the integration of both CAD and CAE is a key step. Most methods introduced in this section focus on meshing the CAD geometry and combining the meshes with a CAE discretization model. Other branches, B-spline and NURBS approximation introduced in Section 1.2, representing the CAD model, are developing rapidly. Many researchers have made great contributions in this topic.

The IGA is a significant product used for solving problems regarding structural mechanics, solid mechanics, and the fluid field [25]. The NURBS-based IGA method outlined in [24] integrates the accurate geometric design with analysis. This is discussed mainly in Section 1.2. However, major shortcomings are obvious in this method. The greatest difficulty is that IGA makes mesh generation complicated in the computation [25]. This is still a highlighted topic to be developed in the FE computation.

In recent years, a methodology framework has been suggested to integrate the highly accurate two-dimensional (2D) model into the commercial simulation software ANSYS, which shows possibilities to provide a suitable and accurate CAD model for FE computations [6]. This framework also offers a brief to construct a 2D FE geometry with TLS. The integrated geometry is based on B-splines. After the extraction of the point clouds measured by TLS, B-spline is applied to approximate the external structure of the arch concerned, as is shown in Figure 3. Obtained B-spline points are converged to describe the 2D geometry. The B-spline-based geometry captures the deformed details and gaps accurately. The error rate regarding the average displacement is obviously improved when the computational model applies the B-spline-based 2D model.

1.5. Motivation and Framework

The common mechanical and structural models contain unknowable deformations after the unexpected mechanical process. It is necessary to apply advanced measurement methods to detect the accurate deformed surface information. This research focuses on the FEA based on the parametric model by approximating point clouds from TLS. It is expected to apply the parametric approximation method to improve the accuracy in describing the measured data of 3D structural surfaces. Therefore, the entire FE model of composite structures can be generated in a novel parametric method. The parametric model provides a continuous and parametric expression that can be applied in both CAD and CAE. This method makes up for the shortcomings of the FE geometric description regarding the discontinuity of the mesh-based model in Section 1.4. The final goal is to reduce computation errors in FEA due to the geometric simplification of structural models.

A reliable FE computation can decrease the number of actual experiments and reasonably reduce the expensive cost of the experiments. Therefore, improving the computation accuracy has become one of the key goals in the FE computation process of structural problems [53,54,55,56].

The general method to carry out the structural analysis is to apply the simplified model mainly. This method is separated into three steps, as is shown in the blue parts of Figure 4. The first step is to idealize the detailed structures into simplified CAD geometries or designed geometries. The main features following the geometric details in the simplification process are normally the height, the width, the length, and the thickness. The second step is to transfer the CAD model to the CAE model. The third is to compute the CAE model. This is a scientific and significant method if the task is to design and improve the object structure in the computations of structural problems [57,58]. However, there is an inaccurate simplification when the structure contains detailed gaps and irregular parts. It is essential to describe the detailed structures in a more scientific and accurate way.

Accurate model approximation of the structural problem leads to more physically accurate results in the final computation process. The mesh discretization of the geometry based on point clouds is a significant method to obtain accurate model reconstruction, which was confirmed in the previous Section 1.4. Their main advantage lies in their continuity characteristic describing surfaces and curves. Attempts to apply B-spline curves in a 2D arch model indicates and validates the significance of the parametric model for the FEA [6]. However, it should be noticed that the main limits lie in the discontinuity of the discretized model description. This can lead to incompatibilities in processing grids using different commercial FE computation software. Therefore, if an effective geometric description that satisfies the continuity requirements is provided, the discontinuity problem of the mesh discretization can be avoided. Using the IGA concept developed, it is found that B-Spline or NURBS-based geometry description can improve the smoothness and accuracy of the parametric model. There is no mature and significant technology on the practice of B-spline-based or NURBS-based numerical model description for very complex objects [59]. Therefore, this still needs to be developed, especially in the case of dealing with complex parametric structures characterized by more irregularities and deformations.

Since the accurate free-form property of B-splines is so beneficial, this research chooses this method to generate the parametric model for the FEA. As it is shown in the green parts of Figure 4, the parametric model is created here as a mediation tool. The approximation task of point clouds measured from TLS can be accomplished by B-spline fitting. In this case, the advanced parametric model can meet the continuity requirement and, subsequently, be adopted in commercial FE computation.

Attention should be paid to the fact that the choices of measurement methods are diversified. Many excellent and efficient methods have been developed, e.g., TLS, digital photogrammetry, and radar technology [60]. The data measured to generate the parametric FE model in this research is acquired by means of TLS. In other references, the radar, the photogrammetry, and the integration technologies of multi-sensors can also be alternatively developed to measure the geometry surface information [38,39,40,41,42,43]. Analogously, the T-spline method, and NURBS approximation can be an alternative to the B-spline methods in fitting the geometric surface.

There are five main sections in this paper. Section 1 describes the geometric generation and advanced measurement methods regarding structural problems. A contrast analysis is carried out to understand the advantages and disadvantages of different methods comprehensively. Meanwhile, the kernel motivation and the overall framework of this paper are declared. Section 2 discusses the key methodology of the parametric model and applies the novel parametric model in the FEA. The process to acquire the data measured, the parametric function of B-spline approximation, and the model generation steps are listed in detail in this section. Section 3 applies the parametric model regarding structural problems into a FE computation and computes the results in both the static structural and dynamic computations. Section 4 discusses the static and dynamic computation results and analyzes the model quality based on point clouds, the errors between the parametric model and the simplified model, and the result deviations between both models regarding the static and dynamic computation. Section 5 draws the main conclusions of this paper.

2. Methodology of the Parametric Model

The methodology of the parametric computation regarding the structural problems contains four main steps in total: the measurement, the parametric approximation, the parametric modelling, and the computation, as is shown in Figure 5. The common computation method constructs the target object as a simplified form in the CAD modelling software, which only keeps basic features regarding some sizes, while the deformed details are omitted. Measurement experiments with accurate sensors are desired, which aims to reconstruct the deformed information. Measured data is converted to the available parameter set and then extracted into different features. The 3D model can be fitted by the B-spline volume directly if the model shape is not complex [28], which is not common in structural problems. Therefore, volume features are further decomposed into surface patches. The B-spline surface approximation is applied to fit the decomposed surfaces. The final parametric surfaces are only acceptable when the fitting accuracy is satisfactory. The parametric surfaces are converted to the CAD modelling software. The parametric geometric modelling is processed based on the parametric surfaces and transformed into the computational model. Lastly, the computation progress can be accomplished after the mesh generation and boundary condition application. Kernel links are listed in the next sections in details and applied into a composite structural building.

2.1. Feature Acquisition of the Object

In the feature acquisition progress, it is initially necessary to obtain the global feature of the target object through measurement methods by advance sensors. A TLS-based sensor system is adopted in this step. Secondly, the task is to transform the large datasets into the key and useful coordinate values. Thirdly, the patch decomposition is carried out due to the discontinuity characteristics in the right angle parts, which is ubiquitous for complex building and mechanical structures. There are also some windows, entrance and exit passages, and non-removable parts on the surfaces. Finally, as a result, it is necessary to recognize the boundary features that represent windows, passages, and non-removable parts. In the following, the detailed experiment and operations are introduced to obtain the object feature with the application of these steps.

The researched object is used on the basis of a composite structure regarding a warehouse building with a north-wall length of about 16 m, south-wall length of about 16.5 m, width of about 14.2 m, height of 3.6 m, and wall thickness of 0.2 m. The point cloud of this structure is sample data from PointCab GmbH which offers powerful software to make the processing of high-resolution point clouds easy. It also offers open access for users to download sample data regarding point clouds. The aim is to research weak situations and positions in this deformed structure by combining TLS with the FE computation. The TLS equipment (FARO Focus^S 350) and laser tracker equipment (FARO Laser Tracker Xi) were used here to obtain the 3D point clouds [46,61,62,63]. The systematic distance error of this type TLS is$\pm 1 \mathrm{mm}$. The vertical and horizontal resolution of the TLS is$0.009°$. The maximum vertical scanning speed of the TLS is$97\mathrm{Hz}$. The laser tracker, which has an angular accuracy of$18 \mathsf{\mu }\mathrm{m}+3 \mathsf{\mu }\mathrm{m}/\mathrm{m}$ in angle measurement performance, is applied for the purpose of reference and validation.

The scanning positions are separated into three parts, as it is shown in Figure 6 from the top view. The north direction is marked as N here. There is a total of 13 positions for scanning totally. From the ground-based scanning positions outside the building, there are five positions to capture the outside details: 1, 2, 3, 4, and 5. There are also five scanning positions inside the building: 6, 7, 8, 9, and 10. There are normally only one or two scanning positions inside such a building. However, there are four forced pillars inside the building. Consequently, the number of scanning positions is increased to five to capture the details of the forced pillar. The final scanning positions 11, 12, and 13 are based on the building roof.

The point clouds measured from PointCab contain not only the data needed, but also the unexpected noisy points. A smoothing operation is commonly applied when the initial images are dealt with. This is known as blurring. The role of the smoothing operation here is to reduce the noisy points in point clouds from the TLS. The Gaussian Filter applying the OpenCV function “GaussianBlur” is adopted in this research. The size value and application of the function can be referred to in [64]. The detailed structure of this warehouse building can be found in the point clouds of Figure 7 after the smoothing operation. There are coordinate, scalar, RGB color, and normal values in the scanning results. The kernel information in Figure 7 is the coordinate value of every point which is finally extracted as a P matrix by deleting the unexpected values in MATLAB R2018b, see Equation (1).

$$\mathit{P}=\left[\begin{array}{ccc}{x}_1& y_1& z_1\\ \begin{array}{c}{x}_2\\ ⋮\end{array}& \begin{array}{c}{y}_2\\ ⋮\end{array}& \begin{array}{c}{z}_2\\ ⋮\end{array}\\ x_k& y_k& z_k\end{array}\right]$$

(1)

The focus of this research lies in the reconstruction of the complex composite structures with unknown deformations based on the parametric solution. Hence, points on the volume are firstly decomposed into points on the surface patches. The parametric solid volume is assembled again by these patches after accomplishing the appropriate surface approximations. This is the solution to reconstruct the parametric solid model in this research. The decomposed patches are shown in Figure 8.

As shown in Figure 8, the global volume is decomposed into five parts which are marked as east, west, south, north, and pillar. The door, window, and lower window features are represented as D1 to D5, W1 to W20, and L1 to L6, respectively. The following fitting process is based on the continuous B-spline surfaces. Consequently, surfaces connected with right-angle interface layer should be subdivided into smaller patches. The east outside wall in Figure 9, which is implemented in the open access software CloudCompare V2.8, is subdivided into five patches by blue dotted layers. The purpose is to guarantee the continuity of every independent patch. As a result, the subdivision is adopted in other parts in Figure 8, which is followed by the same method in Figure 9. The boundary extractions of the window and door features apply the boundary threshold method [6] using the principle of Delaunay triangulation [64,65] with a threshold value of 10 mm. The aim of this extraction is to recognize the boundary features of different surfaces for the future trimming step in the B-spline surface.

2.2. B-Spline Approximation

2.2.1. B-Spline Basis Function

Given a knot vector$\Xi$, the B-spline basis function${\left\{{\mathrm{N}}_{i,p}\right\}}_{i=1}^n$ starts with a zeroth order basis function ($p=0$) as is shown in Equation (2). It is known as the Cox–deBoor recursion function [66]. The values of${\xi }_i$ are elements of knot vectors. Here, the parameter$\xi$ is the coordinate in the parametric space that is needed to satisfy the relationship in Equation (2). The parameter$\xi$ varies from the minimum value to the maximum value along the calculated curve.

$${\mathrm{N}}_{i,0}\left(\xi \right)=\left\{\begin{array}{c}1 if {\xi }_i\le \xi <{\xi }_{i+1}\\ 0 otherwise\end{array}\right.$$

(2)

The B-spline basis functions ${\left\{{\mathrm{N}}_{i,p}\right\}}_{i=1}^n$ are defined recursively by Equation (3), for$p\ge 1$. $p$ is the degree of the basic function. Here, it is possible for the status during the molecular formula computation in Equation (3) to behave in special situation as$0/0$in the ratio of the numerator and the denominator. The convention of this special situation is to default this as 0 in the recursion work. For detailed explanations, please refer to [12,67].

$${\mathrm{N}}_{i,p}\left(\xi \right)=\frac{\xi -{\xi }_i}{{\xi }_{i+p}-{\xi }_i}{N}_{i,p-1}\left(\xi \right)+\frac{{\xi }_{i+p+1}-\xi }{{\xi }_{i+p+1}-{\xi }_{i+1}}{N}_{i+1,p-1}\left(\xi \right)$$

(3)

2.2.2. B-Spline Curves and Surfaces

The last remaining step is to create the control points$Q_i$ of the B-spline. They form the so-called control polygon in the calculation. Interested readers can refer to [12,13,63] for more details regarding the B-spline parameter calculation.

The B-spline curve$C\left(\xi \right)$ can be constructed by linearizing the B-spline basis functions, as is given by Equation (4). Coefficients of the basis functions are referred to as control points.

$$C\left(\xi \right)={{\displaystyle \sum }}_{i=0}^n{N}_{i,p}\left(\xi \right)Q_i$$

(4)

$${\mathit{S}}^{\mathit{B}}=S\left(\xi ,\eta \right)={\displaystyle \sum }_{i=0}^n{\displaystyle \sum }_{j=0}^m{N}_{i,p}\left(\xi \right)M_{j,q}\left(\eta \right)Q_{i,j}$$

(5)

The B-spline surface${\mathit{S}}^{\mathit{B}}$ is represented as Equation (5), which is a tensor product of B-spline curves. The control points of the B-spline surface form a control net$Q_{i,j}$,$i=0,1,\cdots ,n$,$j=0,1,\cdots ,m.$ $N_{i,p}\left(\xi \right)$ and$M_{j,q}\left(\eta \right)$ are basis functions of B-spline curves which approximate B-spline curves in two directions.

2.2.3. B-Spline Approximation

The first B-spline approximation task is based on curves which are of significance to describe different feature boundaries of the structure, for example, windows and doors. The degree value of basis functions is 3 [6,67]. The values of the control point number of different feature boundaries regarding the east side wall in Figure 9 are listed in

Table 1. Numbers of control points of B-spline curves regarding feature parts.

Region	Side	n+1	Region	Side	n+1
D1	Left	20	D3	Upper	30
	Upper	9		Right	28
	Right	20		-	-
D2	Left	20	D4	Left	28
	Upper	9		Upper	36
	-	-		Right	28

. Four feature regions have to be extracted on the east outside surface. Ground-based sides of each feature are ignored. Other sides obtain reasonable numbers of controlled points as is shown in Table 1, which are referred to in detail in [12,63]. The extracted B-spline feature curves are then applied as the trimming and projection boundaries to the B-spline surfaces.

Aiming at reconstructing the solid parametric model of composite structures, the surface-based approximation is a key step in the reconstruction step. Taking the south side surface as an example, the deformed surface is approximated as a parametric surface with the application of B-splines, as is shown in Figure 10a. The degree value of the basis functions is 3. The degree value can be adjusted according to the approximation demands. The choice in this manuscript is based on previous papers [6,63]. The control points are distributed in two directions in the surface fitting process, as mentioned in Equation (5). In the south surface example, numbers of control points are 36 in the height direction and 165 in the length direction. The accuracy characteristic of B-spline surfaces in [6,63,67] is proved which adopts 20 control points in the approximation progress regarding an arch structure with the length of 2 m, which indicated that there is a control point every 0.1 m. The determination regarding the quantity of control points affects the quality of the approximated surface. Detailed explanations regarding the parameter and model selection problem are given in [63]. The length of this surface is about 16.5 m. The height is about 3.6 m. Hence, control points are added to the model every 0.1 m in this example. The standard deviations are calculated in Section 4.1, which indicates that the quality of the parametric model is more accurate than the simplified model. Consequently, the number of the control points can be considered as reasonable and reliable in this computation.

Parametric surface trimming is a complex activity [68]. The shape of B-spline curves regarding the feature parts are not absolutely compatible with the shape of the B-spline surface due to the uncertainty of the boundary extraction progress. Hence, the surface trimming step is of great significance to obtain accurate boundary recognition regarding the feature parts in the solid model. The final B-spline surface with window features after trimming is shown in Figure 10b.

3. Finite Element Method

3.1. Computations

The final reconstructed parametric model and the designed simplified model are shown in Figure 11. The simplified model here is a similar model that keeps the same length, width, height, and thickness as the measured data. It also ensures the same small feature parts, for example, windows, doors, and holes. The most obvious difference is that the parametric model contains the detailed deformed surfaces while the simplified model is constructed with flat surfaces only. Deformed surfaces regarding this building are obtained by the TLS. The description of these deformed surfaces takes advantage of the B-spline approximation to fit the point clouds. The surface features can be found in the dark parts of Figure 10a,b. The detailed analysis regarding the model quality and accuracy is discussed in Section 4.

It should be noticed that the final reconstructed model only generates the detailed geometric surfaces from the TLS data. The internal information regarding the material should be taken into account in the future.

Here, the pillar parts consist of reinforced concrete, as is shown in Figure 12. Hence, the reinforced concrete characteristic [69] is added into the final composite structures based on the parametric and simplified models. The diameter of the reinforcement is 25 mm. One pillar contains four steel bars. The concrete density applied here is$2400 \mathrm{kg}/{\mathrm{m}}^3$. The Young’s modulus and Poisson’s ratio of the concrete are$30 \mathrm{GPa}$ and $0.18$, respectively. Additionally, the reinforcement is structural steel material, which contains the density of $7850 \mathrm{kg}/{\mathrm{m}}^3$, Young’s modulus of $20 \mathrm{GPa}$, and Poisson’s ratio of $0.3$. The reference temperature of the reinforcement is the same as that of the concrete. The properties of the concrete and the reinforcement are listed in detail in

Table 2. The main properties of the computation.

(Concrete)
Concrete density (kg/m3)	2400
Coefficient of the thermal expansion (1/°C)	0.000014
Reference temperature (°C)	22
Young’s modulus (GPa)	30
Poisson’s ratio	0.18
(Steel)
Steel density (kg/m3)	7850
Young’s modulus (GPa)	20
Poisson’s ratio	0.3

. The reconstructed composite models after compounding the characteristics of the reinforcement and the concrete are finally computed in the ANSYS Workbench software.

There are many future behaviors and loading types of the structural problems. The structural loading and small vibrations are common destructive performances in engineering [48,70]. Hence, the force-loading and vibration progresses are two main reasons leading to the possible future behaviors. This manuscript compute two kinds of mechanical behaviors based on the static and dynamic analysis.

In the static analysis based on the force-loading computation, surfaces in contact with the ground are defined as fixed supports. The roof surface is loaded with a 2 kN force in the negative direction of coordinate z. The overall loading process is based on the static structure computation. The aim of this computation is to find the weak zones in different models and validate the significance of the parametric model in a static structure analysis.

The dynamic analysis in this computation focuses mainly on the vibration computation based on the modal computation. It is combined with the random vibration module. The aim is to evaluate the response characteristics of different model structures in dynamic processes. Here, surfaces in contact with the ground are set as fixed supports in the dynamic computation.

3.2. Static Computation Results

The deformation results regarding both models in the static structure computation are shown in Figure 13. The red color bar in this figure indicates the largest deformation zones in the deformation progress. The maximum deformation of the parametric model is 0.193 mm, which appears in the north-central location. The maximum deformation 1.36 mm of the simplified model appears in the central part, which is seven times the value of that in the parametric model. One of the possible reason is the wall weight difference between the two models. From the perspective of Figure 13, small deformations can be found in the north side wall of the parametric model. It should be noticed that both models contain the similar deformation diffusion and development laws generally.

3.3. Dynamic Computation Results

In the dynamic computation, the focus of this paper is based on the modal analysis regarding the structure which is a basis of the vibration analysis. It indicates natural frequencies and mode shapes of this structure. Figure 14 shows the deformation shape in different modes of two models.

Generally, the deformation contour based on the vibration computation is different with those based on the static computation. The modal shapes in different modal conditions regarding the same model are completely inconsistent. The maximum modal deformations in the parametric model are 159.94 × 10⁻³ mm, 159.95 × 10⁻³ mm, and 146.06 × 10⁻³ mm, while the maximum modal deformations in the simplified model are 445.18 × 10⁻³ mm, 388.08 × 10⁻³ mm, and 402.23 × 10⁻³ mm, respectively. The development laws regarding the maximum modal deformation are inconsistent. This indicates that the dynamic results in the simplified model computation are invalid. Hence, the accurate parametric model is necessary in the composite structural computation and it is beneficial to improve the reliability and accuracy of the behavior prediction in the FEA progress.

4. Discussion and Analysis

This research is based on the model reconstruction of the composite structures. The evaluation of different models is unknown. Consequently, it is necessary to validate the model quality and describe the main effects due to the model deviations in following subsections. Section 4.1 discusses the quality of the reconstructed model itself and the deviation between the measured point cloud data and reconstructed models. Section 4.2 analyzes the deviations between two researched models. Section 4.3 extracts the deformation contour of both models into one figure and analyzes development laws regarding the deformation. Section 4.4 carries out a detailed discussion regarding the stress result. Section 4.5 analyzes the modal features of the dynamic computation.

4.1. Model Quality Based on Measured Data

Model errors based on different surfaces are listed in detail in

Table 3. Approximation deviations of surfaces.

Surfaces	Maximum Deviation (mm)	Standard Deviation (mm)
(Parametric model)
East outside	37	1.07
West outside	34	0.65
South outside	26	1.15
North outside	27	1.09
Roof outside	26	0.94
Pillar part	51	0.82
(Simplified model)
East outside	55	15.13
West outside	51	10.58
South outside	46	21.38
North outside	48	22.56
Roof outside	43	14.59
Pillar part	53	18.23

in order to determine the model qualities. The error analysis$d_i$ is based on the comparison between the point cloud matrix P and the analyzed model surfaces. The standard deviation$\delta$ is calculated by Equation (6). In the parametric model based on B-splines, the standard deviations of all fitting surfaces are about or within 1 mm. The maximum standard deviation is 1.15 mm which is on the south outside surface. The maximum standard deviation in the simplified model is 22.56 mm which is on the north outside surface.

$$\delta =\sqrt{\frac{{{\displaystyle \sum }}_{i=1}^k{d}_i^2}{k}}$$

(6)

The maximum deviations in the parametric model are in the interval of [26,51] mm. These values are much larger than the standard deviations. This is due to the remaining noise points on these surfaces. Hence, the maximum deviation also partly contains outliers. The standard deviations of the east, south, and north surfaces are larger than that of other surfaces from Table 3. It is also obvious from Figure 8 that the boundary regions of door or window features are large. The possible reason for the larger standard deviation is that some outliers exist in large boundary regions in the TLS measurement. Both maximum and standard deviations in the simplified model are obvious due to neglecting of the surface deformation in the reconstruction progress of the simplified model.

Generally speaking, the approximation quality of the parametric model behaves superior to that of the simplified model on the basis of surface deviation comparisons, which is indicated in Table 3. The average approximation accuracy of the parametric model can reach the millimeter level. The maximum deviations of all outside surfaces in the simplified model are dozens of times larger than the corresponding standard deviations in the parametric model. This indicates that the simplified model is far from the true measured data.

The standard deviation of the parametric model can still be reduced by increasing the control points and adjusting the order of basis functions in the approximation progress of the composite structures. However, the calculation cost and the workload of the model reconstruction will increase greatly. Meanwhile, the mesh generation and the FE computation will be more complex. As a result, the global computation of the parametric model will perform extremely inefficiently in this condition. In this manuscript, the object researched is a large scale model. It is important to balance the calculation cost and computation accuracy. The present accuracy of the parametric model is much higher than the accuracy of the common simplified model. The determination of the geometric accuracy differs in different scale simulations. The computational model reconstruction of details in composite structures with multi-scales can be improved by adjusting control points and basis function orders based on multi-scales.

4.2. Deviation Analysis of Two Models

Details about the centroid position and moment of inertia regarding to different models are extracted in

Table 4. Properties of solid models.

Parametric Model		Simplified Model
Mass volume (m3)	76.5	Mass volume (m3)	67.8
Local origin point X (m)	−21.2	Local origin point X (m)	0
Local origin point Y (m)	−6.8	Local origin point Y (m)	10.1
Local origin point Z (m)	0.2	Local origin point Z (m)	−9.9
Centroid X (m)	−13.6	Centroid X (m)	7.2
Centroid Y (m)	−15.4	Centroid Y (m)	1
Centroid Z (m)	2.96	Centroid Z (m)	−7.4
Relative centroid X (m)	7.6	Relative centroid X (m)	7.2
Relative centroid Y (m)	−8.6	Relative centroid Y (m)	−9.1
Relative centroid Z (m)	2.76	Relative centroid Z (m)	2.5
Moment of Inertia Ip1 (kg$·$m2)	5.99	Moment of Inertia Ip1 (kg$·$m2)	5.7
Moment of Inertia Ip2 (kg$·$m2)	4.3	Moment of Inertia Ip2 (kg$·$m2)	4.01
Moment of Inertia Ip3 (kg$·$m2)	9.9	Moment of Inertia Ip3 (kg$·$m2)	9.3

. The mass volume refers to the size of the volume filled with model materials, which differs from the common solid enclosed volume. The mass volume of the simplified model is 67.8 m³, which is 88.6% of the mass volume of the parametric model. This indicates the weights of both models are different. As is generally known, the deformation exists when the force is loaded on the non-ideal rigid body. Therefore, the effects due to the different self-weights cause different predeformation effects. Hence, the computation based on the simplified model brings predeviation before the FE computation.

$$\mathit{R}\mathit{C}=\mathit{C}-\mathit{L}\mathit{O}\mathit{P}$$

(7)

The centroid positions are listed in Table 4. The local origin points of both models are set at the endpoint on the north lower side of the east outside surface. The local origin point coordinate in the global coordinate system of the simplified model is (0, 10.1, −9.9) m. The local origin point coordinate in the global coordinate system of the parametric model is (−21.2, −6.8, 0.2) m. The relative centroids in the global coordinate system can be calculated by Equation (7). In this equation, $\mathit{R}\mathit{C}$ is the coordinate of the relative centroid, $\mathit{C}$ is the coordinate of the centroid in global system, and $\mathit{L}\mathit{O}\mathit{P}$ is the global coordinate of the local origin point of each model. The coordinates of each relative centroid are shown in Table 4, respectively. The coordinate difference$\left(∆x, ∆y,∆z\right)$ between two relative centroids is (0.4, 0.5, 0.26) m, respectively. The deviation rate of the centroid coordinate based on the parametric model centroid is about (5.3%, 5.8%, 9.4%). As a result, the deviation rate of the moment of inertia in different directions is (4.8%, 6.7%, 6.1%). While the moment of inertia is greater, the ability to resist deformation of the object is stronger. Hence, the difference of the centroids and moments of inertia is also a possible reason for the deformation deviation after loading.

Surfaces in the simplified model are all based on flat planes. The surface${\mathit{S}}^s$ of the simplified model constructed in this research is described by Equation (8) which is a plane. Hence, the distance of the B-spline surface point$r_B$ in the parametric model to the simplified flat surface${\mathit{S}}^s$ is obtained by Equation (9). $r_B=\left(x_i,y_i,z_i\right)\in {\mathit{S}}^{\mathit{B}}$. The total amount of approximated B-spline points is$l$. The global average deviation$D_{\mathsf{\Sigma }}$ between the B-spline surface${\mathit{S}}^{\mathit{B}}$ and the simplified surface${\mathit{S}}^s$ can be indicated by Equation (10). Hence, the global average deviation$D_{\mathrm{average}}$ regarding 11 surfaces can be calculated by Equation (11). The final average deviation of the global surfaces is 8.76 mm. This value indicates the obvious shortcomings of the simplified model which is applied commonly. The parametric model describes the deformed details on all surfaces with the advantage of B-splines. However, the simplified model reconstructs the surfaces as idealized and flat ones and ignores the deformation information description. This is the main reason for this obvious deviation between both models.

$${\mathit{S}}^s=\left\{\left(x,y,z\right)|ax+by+cz+d=0\right\}$$

(8)

$$∆=\frac{\left|a{x}_i+b{y}_i+c{z}_i+d\right|}{\sqrt{a^2+b^2+c^2}}$$

(9)

$$D_{\mathsf{\Sigma }}=\Vert ∆\mathit{S}\Vert =\Vert {\mathit{S}}^s-{\mathit{S}}^B\Vert =\frac{1}{l}{\displaystyle \sum }_{i=1}^l\frac{\left|a{x}_i+b{y}_i+c{z}_i+d\right|}{\sqrt{a^2+b^2+c^2}}$$

(10)

$$D_{\mathrm{average}}=\frac{1}{11}{{\displaystyle \sum }}_{i=1}^{11}{D}_{\mathsf{\Sigma } i}$$

(11)

In summary, it is of great significance to compare the deviation performance when the geometric models of structural problems are reconstructed before the FE computation. The deviations between the parametric and the simplified models are obvious. The mass volume of the parametric model is larger than that of the simplified model which results in the deviation of the predeformation due to the self-weight. The relative centroid of the parametric model offsets from the simplified model relative centroid in the global coordinate system. The average deviation between two models is finally very large, which leads to influence the parametric model in the model selection progress regarding the structural problems.

4.3. Deformation Analysis Based on Static Structure

Aiming at researching the effect of different geometries on the FE computation, it is of significance to carry out the deformation analysis. The deformation contour figure regarding different deformation values of the roof surface is extracted in Figure 15. The blue lines indicate deformation results based on the parametric model, while the red dashed lines show the deformation results based on the simplified model. The centroids of both models are marked in this figure.

The deformation contour spacing of the parametric model in the L1 positive direction gradually decreases from the maximum deformation zone, while that in the L2 positive direction decreases initially and then increases. There is the same changing law of the deformation contour in the simplified model. Here, 30% of the maximum deformation is defined as the large deformation zone, which is marked as a red solid line in the simplified model and a blue dashed line in the parametric model. It is obvious that the total areas regarding both large deformations are similar. With the decrease of the deformation, the contours of both models are closer. This indicates that the simplified model can also fit a relatively accurate deformation-changing law outside the large deformation zones.

Obvious differences between both models are the maximum deformation positions and large deformation zones. The maximum deformation position in the simplified model is almost in the center of the roof surface. With the application of the parametric model, this position moves to the positive direction of L1 and L2. However, the deviation of two maximum deformation positions in the L2 direction is not as obvious as that in the L1 direction. With the combination of the centroids in Figure 13, the relative position of two centroids is opposite to the relative position of the maximum deformation points. The reason focuses on the effect of the centroids’ position and values of moments of inertia in different directions. With the movement of the centroid to the left and lower sides, here from the simplified model to the parametric model, the anti-deformation ability in the left and lower zones improves in the parametric model. Consequently, the large deformation zone of the parametric model moves to the right and upper direction when compared to the simplified model. The large deformation position changes more obviously in the L1 direction than in the L2 direction. This is due to the geometric characteristics of this object. From the relative perspective in Figure 13, the left and right sides are, respectively, the south and north surfaces of the building. Their architecture and locations are more symmetrical. However, due to the existence of the obvious concave part in the upper zone, this feature structure can better support the loading force of the object. Therefore, the large deformation zone does not move upwards obviously.

4.4. Stress Analysis Based on Static Process

Equivalent stress is of great significance to monitor the strength of the damage possibility in the loading progress, as it is shown in Figure 16.

The color difference in the deformed zones of the surfaces is obvious. This indicates the stress change. There are two concave strip shapes in the red enlarged area of Figure 16, while the equivalent stress contour of this zone is smooth in the simplified model of Figure 16b. The equivalent stress of the upper concave peak part is 125.24 MPa, while the equivalent stress of the concave valley part is 123.93 MPa. The equivalent stress of the lower concave peak part is 131.99 MPa, while the equivalent stress of the concave valley part is 130.49 MPa. Three line examples are extracted in the same concave change part. The length is about 0.2m which is close to the width of the concave region. The equivalent stress along the line length is listed in this figure. It indicates that the deformed concave shape affects the stress results. By contrast, this detailed feature cannot be found on the flat surface of the simplified model. It indicates that the deformed surfaces have a direct effect on the surface stress behaviors in the loading progress. There is a great difference regarding the maximum equivalent stress between the parametric model and simplified model, which is enlarged as the black dashed frame zones in Figure 16. Two black enlarged zones in the simplified model show that the maximum equivalent stress on the roof surface is close to the maximum value of the global model. However, the maximum stress of the parametric roof surface is far less than the global maximum stress value. The damage behavior in the engineering application will be predicted inaccurately if the parametric model is ignored.

The contour shape regarding the large stress zone of both models is generally similar. There is a similar development law of the equivalent stress contour if the specific value is not considered seriously. Consequently, similar models contain similar stress development laws in the same conditions. Hence, the common simplified model can be reasonably adopted in the computation of structural problems when the research focus is on the overall laws because the simplified model is more efficient.

4.5. Vibration Analysis

Modal analysis is the basis of the dynamic structural computation. With the help of the modal analysis, it is helpful to recognize different responses of the structural problems to different dynamic loadings. The vibration analysis involves a wide range of knowledge. However, this research focuses on the analysis of two reconstructed geometric models. Therefore, only a contrast dynamic analysis between two models is carried out based on the frequency and the vibration deformation.

The relative modal data in

Table 5. Contrast vibration analysis between two models.

Mode	Max Deformation	Error 1	Frequency	Error 2
(Parametric model)
1	0.16 mm	-	18.732 Hz	-
2	0.16 mm	-	21.174 Hz	-
3	0.146 mm	-	24.912 Hz	-
(Simplified model)
1	0.445 mm	0.285 mm	8.341 Hz	10.391 Hz
2	0.388 mm	0.228 mm	10.284 Hz	10.89 Hz
3	0.402 mm	0.256 mm	12.938 Hz	11.974 Hz

are beneficial for comparing the vibration characteristics of the structural problems [71]. Error 1 calculates the relative deviation regarding the maximum vibration deformation between the parametric and simplified models. Error 2 calculates the relative deviation regarding the modal frequency of two models. Two errors are obvious, which implies the inadequacies of the common simplified model. By contrast, the maximum deformations of the parametric model in the 1st, 2nd, and 3rd modal conditions are far less than those of the simplified model, while the frequencies in the parametric model are higher than in the simplified model. This indicates that the parametric model is more capable in the dynamic computation of structural problems.

5. Conclusions

This manuscript offers a generic methodology which focuses on the FEA based on the parametric model by approximating 3D actual feature data. 3D actual feature data can be acquired from many efficient measurement methods, e.g., TLS, digital photogrammetry, and radar technology. The methods in approximating parametric surface models accurately are diverse, e.g., the T-spline method, the B-spline method, and NURBS approximation.

It was applied to detect the deformed surface information by TLS which is an accurate and reliable measurement technology in this research. The B-spline method was applied to approximate the measured point clouds data and generate the parametric 3D model of structural problems. The parametric model can be applied in both CAD modelling and CAE analysis innovatively. The model quality and some deviations were discussed. The static and dynamic computations were carried out to imply the advantages and disadvantages of both models regarding the responses of different loadings. The main conclusions are as follows.

1. The numerical parametric model of structural problems satisfies the continuity characteristic in constructing the 3D model with the advantage of the parametric description method B-splines.

2. The parametric model can reserve deformed features of the structures accurately. When the data measured is compared, both the maximum and standard deviations of the parametric model are far less than the simplified model. The mass volume of the simplified model is smaller than that of the parametric model due to the lack of the description of deformed walls. This is one of the main reasons why there are obvious deformations between the parametric and simplified models in the static and dynamic computations.

3. In the static structural analysis, the simplified model is acceptable and it is more convenient and efficient to analyze and predict the overall development law regarding the deformation and stress outside the large deformation zones. The simplified model has obvious shortcomings and inaccuracies in large deformation zones compared to the parametric model in the static structural computation.

4. The FEA computation based on the parametric model is more reliable which benefits from the parametric description of the actual object, e.g., B-Spline method.

5. The parametric model contains fewer errors in the dynamic computation. It indicates that the parametric model is more reasonable and acceptable in the dynamic computation and analysis, especially in the large deformation and stress zones.

6. There is an obvious effect of predeformed parametric surfaces on the equivalent stress of the composite structures in the loading progress. Therefore, the parametric model is more accurate in predicting the future behavior by analyzing the equivalent stress of composite structures.

Attention should be paid to the fact that the FEA model simplification can bring some efficiency to researchers to a certain extent, e.g., when the research purpose is only to predict some simple regularity problems or development trends. However, numerical features, e.g., the deformation, the stress, and the strain, should be focused on significantly in most cases when it is in the accurate design, optimization, and prediction stages. Therefore, the FEA model simplification is not desirable in this case and the accurate parametric model is significantly desired.

In summary, the FEA based on the parametric model is instructive for understanding the structural progress and predicting the damage behaviors in this paper. The novel parametric model method is efficient and beneficial to improve the reliability and accuracy of the FE computation regarding the composite structural problems. Moreover, it is of great significance to apply the FEA parametric model to health to monitor the future damage behavior of the composite structures in the engineering application.