Research Article  Open Access
Ying Liu, Hongguang Li, Yun Li, Huanyu Du, "A ComponentBased Parametric ReducedOrder Modeling Method Combined with Substructural Matrix Interpolation and Automatic Sampling", Shock and Vibration, vol. 2019, Article ID 6407437, 14 pages, 2019. https://doi.org/10.1155/2019/6407437
A ComponentBased Parametric ReducedOrder Modeling Method Combined with Substructural Matrix Interpolation and Automatic Sampling
Abstract
An efficient parametric reducedorder modeling method combined with substructural matrix interpolation and automatic sampling procedure is proposed. This approach is based on the fixedinterface CraigBampton component mode synthesis method (CMS). The novel parametric reducedorder models (PROMs) are developed by interpolating substructural reducedorder matrices. To guarantee the compatibility of the coordinates, we develop a threestep adjustment procedure by reducing the local interface degrees of freedom (DOFs) and performing congruence transformation for the normal modes and interface reduced basis, respectively. In addition, an automatic sampling process is also introduced to dynamically fulfill the predefined error limits. It proceeds by first exploring the parameter space and identifying the sampling points with maximum error indicators for all the parameterdependent substructures. The exact error of the assembled model at the optimal parameter point is subsequently calculated to determine whether the automatic sampling procedure reaches a desired error tolerance. The proposed framework is then applied to the moving coil of electricaldynamic shaker to illustrate the advantage and validity. The results indicate that this new approach can significantly reduce both the offline database construction time and online calculation time. Besides, the automatic procedure can sample the parameter space efficiently and fulfill the stopping criterion dynamically with assurance of the resulting PROM accuracy.
1. Introduction
The finite element method (FEM) has been commonly used in a number of serious, complicated engineering problems. FE models with complex geometry and coupled physical effects often yield highdimensional ordinary differential equations (ODE) to provide the needed accuracy. However, the analysis of fullorder models (FOMs) often requires very extensive computational resources. To alleviate this problem, reducedorder models (ROMs) have been presented and widely used for many largescale engineering applications. At this moment, research studies of reducedorder modeling methods of linear models have matured [1–3]. Also, these techniques have already been applied in some areas like optimization and control [4, 5], uncertainty quantification [6], and parameter identification [7].
However, systems are usually parameterdependent, and the existing reducedorder modeling methods are not applicable for the models with parametric variations, while in some areas like optimization, system models need repeated calculation over a range of parameter values, which restricts the application range of these methods. This has promoted the study of parametric reducedorder modeling technique. This approach is aimed to generate parametric reducedorder models (PROMs) that approximate the original fullorder dynamical system with high fidelity over the range of parameters.
Various parametric reducedorder modeling techniques have been developed over the last two decades [8]. One way to acquire the new reduced system is building the global basis matrices over the parameter space by the singular value decomposition (SVD) or a rankrevealing QR factorization and multiplying the fullorder matrices by the global basis matrices [9]. By contrast, one can also interpolate the local reduced bases at discrete parameter points to obtain the basis matrix for a new parameter value. Direct interpolation may lead to nonorthogonality of the new basis matrix and proves invalid. A better method was proposed by Amsallem and Farhat [10] to interpolate the local bases on a tangent space of the Stiefel manifold. It is verified that the matrices achieved by this method can preserve the desired properties like orthogonality. Another kind of method is to interpolate the local reduced model matrices by using matrix manifold interpolation or direct interpolation [11], but first of all, these methods should perform a congruence transformation of the local reduced matrices in order that they are in the same generalized coordinate system. Then, Panzer et al. [12] used direct Lagrange interpolation and obtained the reduced matrices at the parameter. Different from that, Amsallem et al. [13] projected all the matrices onto a tangent space. After Lagrange interpolation, the new matrices were projected back. Both of the two ways prove effective in achieving the new reduced matrices. Yet another option for acquiring the new reduced system is to interpolate the transfer functions of local reduced models. Baur and Benner [14, 15] took advantage of local reduced transfer functions at the discrete parameter points and obtained the reduced transfer function at the given parameter value by Lagrange interpolation. However, the necessity of explicitly realizing the transfer function makes this method time consuming and not universal.
All the parametric reducedorder modeling methods mentioned above are based on the offline/online decomposition strategy. In the offline stage, ROMs for different parameter values are generated and adjusted, and the database is thereby built. However, as the dimension of parameters increases, the corresponding sampling points increase exponentially and this approach becomes dramatically expensive, especially for largescale FE models. Also, it is not necessary to repeatedly reduce the entire FE model when only partial structures change with parameters. Therefore, it is a computationally cheap way to partition the entire model into several substructures and separate the parts with parametric variations. Among the various substructuring techniques [16–19], the fixedinterface component mode synthesis method (CMS) was developed by Craig and Bampton [20] and improved by Rubin [21]. This method has been widely used due to its numerical stability and accuracy. By combining the parametric reducedorder method and CMS method, we can obtain the ROM database more efficiently in offline stage.
Regarding the parametric reducedorder modeling methods combined with the CBCMS, Park [22] and Zhang [23] developed a Taylor seriesbased method to approximate the componentlevel stiffness and mass matrices. As for the transformation matrices, they collected all the component modes at the sample parameter values and constructed a global basis vector. Hong [24] proposed a new transformation matrix and a novel parameterization technique to capture the elementlevel nonlinearities of bricktype finite elements. Lee and Cho [25] performed congruence transformation for the component modes and interpolated them in a common coordinate system. However, these methods are all aimed at the basis vectors and involve largescale fullorder multiplications between substructural stiffness or mass matrices and collected transformation vectors in the online phase, which is relatively time consuming.
In order to further reduce the online calculation time, a novel componentbased parametric reducedorder modeling method manipulating substructural reducedorder matrices is proposed in this research. A threestep adjustment procedure is developed and implemented for each substructure, and the transformed reducedorder component matrices are deduced, which is coordinate compatible for direct interpolation. Compared with previous approaches, this method is rather efficient since the calculation is conducted at the reducedorder level.
Another drawback of previous methods is that these PROMs might lose their accuracy if the system matrices are polynomially expanded in terms of a parameter. To capture the polynomial dependence, more parameter samples are needed. The most fundamental and critical issue is the sampling strategy of parameter space. The distribution of the parameter sample points should enable the surrogate model to capture the underlying parametric dependence to a maximum extent. Some scholars [26–28] use the greedy approaches to iteratively solve the HDM and determine the configuration at which the current ROM leads to larger errors. However, computing the exact error requires the solution of the HDM, which is time consuming. Instead, error indicators evaluated with ROMs have been developed. To find the optimal sampling points associated with the largest error indicators, one major approach calculates the a posteriori error indicator of large discrete sets of candidate parameter values and locates the points with maximum errors [29, 30]. The other one solves the modelconstrained optimization problem using gradientbased approaches and selects the solution as the sampling point [31]. Both of these two approaches can determine the optimal sampling location in parameter space.
However, Most of the existing automatic sampling approaches are based on the entire structure. Few literatures applied the automatic sampling strategy to componentbased parametric reducedorder modeling framework. The conventional uniform or random sampling procedure is unable to capture the parameter dependence well. Besides, the lack of error analysis provides no guaranty for the PROMs accuracy over the parameter space.
Therefore, in this article, we have also developed an automatic parameter space sampling procedure for the novel componentbased parametric reducedorder modeling method. In order to sample the parameter space automatically, an error indicator for rapid error evaluation is introduced. For each single substructure, we identify the optimal sample point by a direct exhaustive exploration of large discrete candidate parameter sets. All the components at the optimal parameter points are subsequently synthesized, and the exact errors between the directly computed CMS ROM and the interpolated CMS ROM are calculated to monitor the convergence rate and used as the termination criteria.
This paper is organized as follows: Section 2 gives an overview of the CBCMS method. Section 3 describes the componentbased parametric reducedorder modeling method, deduces the transformation matrices, presents the error analysis, and outlines the automatic sampling framework. In Section 4, the novel method is applied to the moving coil of electricaldynamic shaker FE model to illustrate the advantage and validity of the proposed approach.
2. ReducedOrder Modeling with CBCMS Method
For a damped structure, the FE equations of motion can be written in a matrix form aswhere are the mass, damping, and stiffness matrices of the system, respectively, X is the system response in physical coordinate, F is the force vector, and n is the degree of freedom (DOF) of the system.
There are usually more than hundreds of thousands of DOFs for large FE models, which result in large system matrices and require vast amounts of computer memory and process time for solving. In order to reduce computational cost, the fixedinterface CMS method is employed to partition and generate substructures for its simplicity and computational stability.
To apply the CMS method, the complex structure of interest is partitioned into substructures. The generation pass can be simplified if we omit the damping matrix and define master DOFs where loads are applied. The corresponding damping, force, and contact conditions are only included in the assembling stage. Therefore, for a typical undamped component i, the equation of motion can be written aswhere N denotes the number of subdomains; the physical DOFs of the stiffness and mass matrices can be partitioned into interior and interface DOFs aswhere superscript A denotes the interface “active” DOFs including the excitation and response points and O denotes “omitted” DOFs of the interior. Then, the equation of motion for the ith substructure can be expressed as
If the first m loworder modes of the interior DOFs are adopted as the retained component modes, the following coordinate transformation matrix from physical coordinate to modal coordinate can be obtained by
The transformation matrix consists of the loworder normal modes , and the constraint modes can be expressed as
Here the first m loworder normal modes matrix can be obtained by solving the following eigen problem:
And by assuming no external force on the interior DOFs, i.e., , the following constraint modes can be derived through the upper part in equation (4):
Then, we can obtain the reduced stiffness and mass matrices of the ith substructure by multiplying the transformation matrix to the original matrices:
The lowdimensional matrices are then synthesized after all the components are reduced. In general, the final reducedorder matrices of the entire assembled model can be characterized by one common set for the global interface coordinate and one unique set for each component’s internal coordinate. For illustration, we introduce a simple structure that consists of two substructures and one interface. The two components are denoted by using subscript 1 and 2, and the interface is designated by 3. Based on the interface compatibility, the following equation is tenable:
The displacement vectors of the two components and assembled model are then expressed as
The dynamic equation of the entire FE model can thereby be condensed aswhere the structural damping is assumed in the form , the force vector is generated based on the interface DOFs, and the assembled mass and stiffness matrices and are written as
Conclusively, the entire calculation flow of fixedinterface CBCMS method can be summarized as: partition the entire structure, generate the reducedorder substructures, synthesize the components, and calculate the reduced responses.
3. Proposed Parametric ReducedOrder Modeling Method
Global parametric reducedorder models (PROMs) [32] have been developed for fast reanalysis of structures with parametric variability in their properties. However, it is quite computationally expensive to calculate and collect systemlevel eigenvectors at multiple parameter sets. Therefore, Park [22] and Hong [24] developed componentbased PROMs. Although these methods can reduce the offline computational cost to some extent, it is still relatively time consuming in the online phase, especially when the interface DOFs are relatively large scale, since these methods involve largescale fullorder online multiplications between substructural stiffness or mass matrices and collected transformation vectors. Another drawback is that these PROMs might lose their accuracy if the system matrices are polynomially expanded in terms of a parameter. To capture the polynomial dependence, more parameter samples are needed. It is a fundamental problem how the relative error of PROMs evolves as the number of parameter samples grows. These issues are addressed herein by developing a novel PROM using matrix interpolation and automatic sampling procedure. This new method performs the interpolation for the reducedorder substructural matrices after adjusting their coordinates and samples the parameter space dynamically to meet the predefined errors limits. In general, the procedure for this method can be summarized as constructing the ROMs of the substructures at initial parameter values, generating and transforming new components matrices at the next optimal sampling point to establish the compatible offline components database, interpolating the components ROMs for the relevant matrices at a new parameter value in the online phase, and then continuing the rest of the substructuring scheme.
3.1. Adjustment of the Offline ReducedOrder Component Matrices
To further reduce the online computational cost, we conduct the interpolation for the reducedorder substructural matrices. However, direct interpolation may lead to inaccuracy as these matrices are expressed in different modal coordinates and are incompatible. Thus, adjustment of the reducedorder component matrices should be performed at first. In this section, a threestep adjusting method is proposed. We first transform the substructural matrices according to the minimum Frobenius norm criteria between different normal modes. Then, combined with the local interface reduction method, the secondary transformation relation of the constraint DOFs is deduced. The reduced matrices corresponding to the interface DOFs are subsequently adjusted to guarantee that the resulting substructural matrices are compatible to interpolate.
Consider an arbitrary ith parameterdependent substructure. The structural dynamical equation is described in parametric form as follows:
The uniform or random sampling method can be applied to sample the continuous parameter field p and generate a rich set of data with full coverage of the parameter space. An automatic sampling process is proposed in the following section.
Assume now that one sampling method is utilized to obtain a discrete set of parameter samples p_{j} for j = 1,…, P. The largescale FE model is then generated at each corresponding discrete parameter value:
In order to perform the congruence transformation, the previous normal mode matrix should satisfy the orthogonality condition which is guaranteed by a Gram–Schmidt procedure. The reducedorder matrices change accordingly.
Assume that the reduced stiffness and mass matrices can be rearranged as
The elements of the mass matrix can be further expressed aswhere are the normal modes at jth parameter value and are the corresponding constraint modes. The elements of the stiffness matrix can be expressed in the same way.
These different constraint modes are described in physical coordinates such that they are compatible for interpolation, while the normal modes are incompatible since they are expressed in each modal coordinate and represent different spaces. Therefore, the first step of transformation is aimed at adjusting the matrices related to the normal modes .
Assume that one reduced model is chosen as standard, and the corresponding matrices are . In order to calculate the congruence transformation matrix for the jth reduced model, the following minimization problems are introduced:where denotes a set of orthogonal matrices in , the subscript F represents the Frobenius norm, and is a mdimensional matrix.
By expanding equation (19), the equivalent maximization problems is obtained as
According to [33], equation (20) has an analytical solution aswhere and are calculated by the singular value decomposition as follows:
By multiplying , all the normal modes are transformed to coordinate compatible such that
The corresponding projection of the displacement vector can be expressed as
The primary transformation matrix is deduced as follows:
The transformed reduced stiffness and mass matrices can be rewritten as
The elements of the mass matrix in equation (26) can be described as
Therefore, the entire reduced mass matrix is written as
And the reduced stiffness matrix can be expressed in the same way such that
Although the substructural matrices after the firststep transformation are compatible for interpolation, the size of these matrices could be relatively large scale, which is dominated by the constraint DOFs corresponding to the A partition in equations (29) and (30). To address this issue, Castanier et al. [34] suggested the characteristic constraint (CC) modes. However, this technique is performed at system level and thereby not applicable to our case. In the present study, the local interface reduction is implemented for the secondary transformation.
The eigen problem on constraint DOFs of the mass and stiffness matrices of component i at the jth parameter value is constructed by
By selecting the r eigenvectors from the lowest eigenvalue, the reduced basis of the interface DOF is obtained as follows:
Similar to the previous normal mode matrix, these reduced basis vectors should satisfy the orthogonality condition . This is guaranteed by performing a Gram–Schmidt procedure for these vectors.
The displacement vectors are further expressed aswhere is the secondary transformation matrix and can be described as
The reduced stiffness and mass matrices deduced by equations (29) and (30) can be rearranged as
Considering the reduced basis , the elements of the mass matrix can be transformed as follows:
Thus, the entire reduced mass matrix can be rewritten as
It should be noted that the projection has little effect on the results which is proved in [33]. And the reduced stiffness matrix can be expressed as
Closely resembling the first step, the reduced basis vectors represent different local modal spaces and are still incompatible. Therefore, the goal of the third step is to adjust, in a similar way to the first step, the substructural matrices related to these reduced basis vectors.
Assuming that one reduced model is chosen as standard, the following minimization problems is derived:
According to equations (20)–(23), the solution to equation (39) iswhere and are calculated by the singular value decomposition as follows:
By multiplying , all the reduced basis vectors are transformed to coordinate compatible such that
The relation between the displacement vectors in each coordinate system can be expressed as
The tertiary transformation matrix is
The final reduced stiffness and mass matrices can be rearranged as
And the elements of the mass matrix can be deduced as
In the end, the relation between the final reduced stiffness and mass matrices and the original ones is expressed as
After the threestep adjustment of the coordinates, the final reducedorder substructural matrices are compatible for further online interpolation. The online calculation efficiency is thereby significantly improved thanks to the dramatic reduction of computational complexity, compared to the traditional methods.
It should be noted that all the computations mentioned above are conducted at the subdomain level. Moreover, the transformation calculations are performed in the reduced coordinate and thus need little computational cost.
3.2. Online Calculation Using the ReducedOrder Component Matrices
As for the interpolation method, one can use gridded data interpolation for regular parameter samples or scattered data interpolation for irregular parameter values. Since the interpolation is performed for reduced matrices, the process is quite computational efficient; some highorder methods like polynomial or radial basis function interpolation approach are therefore available. The detailed description can be found in [30]. For simplicity, we utilize the Lagrange interpolation method and regular parameter samples in this research.
Consider a reducedorder dataset . The matrices of the reduced system at new parameter values are interpolated by use of Lagrange interpolation function such thatwhere are the Lagrange basis functions satisfying
These new matrices are then assembled in the online stage. The compatibility conditions for the components can be expressed almost identically to those described in equation (10) such that
The relation between the displacement vectors of components and assembly, together with the assembled dynamic equation, can also be expressed in the same way as equations (11)–(14). As the interface DOFs contain all the nodes of interest, we focus on the DOFs denoted by A.
After the reduced response obtained, we project it back to the physical coordinate in order to calculate the actual response of the interface DOFs by
Similar to this, the reduced modal shape vectors of the assembled model can also be projected bywhere is the reduced modal shape vectors of the assembled system and is the corresponding fullorder ones of the interface DOFs in the physical coordinate.
Compared to the traditional methods, the online interpolation of the reducedorder components matrices could reduce the computation complexity significantly, thus providing a better performance with cheaper online computational cost.
3.3. Automatic Sampling Procedure
In the previous sections, it is assumed that the parameter domain is sampled for P discrete parameter values. In order to capture the parameter dependence to the greatest extent, multiple parameter samples are needed. The traditional componentbased parametric reducedorder modeling approach cannot provide the error estimation. Without the error evaluation, there is no guaranty that the parameter space can be sampled effectively and that the PROMs can satisfy the predefined error limits. Therefore, we also propose an automatic sampling procedure for the componentbased parametric reducedorder modeling method.
For the automatic sampling procedure, the definition of errors is the most critical issue. In this research, the exact relative error between the directly computed CMS ROM and the interpolated CMS ROM has to be calculated and applied to monitor the convergence of the algorithm and terminate the procedure automatically. The exact relative errors of eigenvalues, modal shapes, and responses are defined as follows:where are eigenvalues, modal shapes, and responses of the assembled system described in equations (53) and (54).
It should be noted that the subscript p and the in brackets represent the results calculated using the exact model and the PROM at parameter value p, respectively. The symbol denotes the absolute value.
To calculate the exact relative error, we have to generate the corresponding FE model and calculate the directly computed CMS solution. In practice, however, this process is relatively time consuming. Therefore, the exact relative error is only used to monitor whether the current PROM satisfies the predefined error limits at the parameter with maximum errors, but it is still an open question where the point locates. Identification of the parameter points with maximum errors relies on the inexpensive computation of a posteriori error indicator associated with a ROM solution. Therefore, we utilize a computational efficient equivalent error to indicate where to place the sampling point.
The exact errors are defined in global assembled system to evaluate the overall error level, while the error indicators are defined in local component systems and the search process is performed at a subdomain level.
As the frequency can reflect components’ relative error, we use the eigenvalues calculated by spline interpolation as substitute to the exact results and define the relative error between frequencies calculated by spline interpolation and PROMs as error indicator:
Here, is the calculated by spline interpolation for eigenvalues at discrete parameter values of the ith substructure, while is obtained by solving the following orderreduced eigenvalue problem:where and are the reduced interface DOFs partition of the stiffness and mass matrices for the ith substructure described in equations (48) and (49).
Considering spline interpolation is quite an accurate interpolation method, this simplified error indicator can almost represent the exact relative error between interpolated components ROM and the accurate FOM.
The goal of the automatic parameter sampling approach is to estimate, at each iteration N_{iter}, the parameter maximizing the error indicator over the set for each component. These components with maximum error indicators are subsequently assembled to evaluate the exact relative errors.
This approach proceeds in two phases. In the first phase, the error estimator is computed for each component over the set:
Then, we build the following error pool to construct a surrogate model for the error as a function of the parameter values:
At each iteration N_{iter} of the automatic sampling process, the parameter value is chosen as the maximizer:
The parameter values are then used in the next phase to generate corresponding components. After these components are assembled, the overall exact errors are evaluated by equations (55)–(57) for the synthesized system.
To sum up, the automatic sampling procedure can be summarized as evaluating the error indicator over a candidate parameter set for each component, generating new components at the parameters maximizing the posteriori error indicator, assembling the components, and calculating the overall exact errors for the synthesized system.
3.4. Entire Calculation Procedure
The entire automatic calculation procedure is schematically outlined in Figure 1. In what follows, we explain the general idea behind this procedure step by step.
Step 1. Set the maximum number of iterations , maximum number of candidates , and maximum error tolerance .
Step 2. Run the following steps for the parameterdependent substructures to construct the initial offline database. First, choose the start and end points of the parameter domain as initial parameter samples. Then, generate the fullorder substructure matrices and at these two parameter values. The corresponding normal modes and constraint modes are subsequently calculated through the CMS method. In the end, perform the threestep transformation for reducedorder matrices and build the initial compatible ROM database.
Step 3. Use the subroutine to locate the parameter points with maximum error indicators. First, randomly select a set of candidate parameters for each component. Compute and by the Lagrange interpolation method for each candidate parameter. The corresponding are calculated and the parameter points of interest are then calculated through equation (62). We search the parameter candidate sets for the optimal sample points for each substructure. The reader should note that the optimum returned by such a procedure does not need to be very accurate as sampling the FOM in the neighborhood of the parameter of maximum error will likely decrease the error at that parameter. After that, we perform the following automatic sampling and database updating process on the basis of the initial database. At first, we interpolate the previous ROMs database for new substructure matrices and at the parameter value , assemble all the new ROMs, and compute the interpolated . Then, we generate the new model and at and perform the threestep transformation procedure to obtain the new ROM . These new reduced matrices are assembled to calculate the exact . The errors of eigenvalues , modal shape vectors , and responses are determined between the exact and interpolated models. If the final stopping criterion is not fulfilled, update the offline ROMs database and iterate step 3 again until the overall PROMs are satisfactory. If the final stopping criterion is met, then exit.
Compared with the global methods, this componentbased parametric reducedorder modeling approach can significantly reduce the parameter sample points number from to . Besides, in comparison with the classical componentbased PROMs, the proposed one can substantially reduce the online computational complexity. Furthermore, the automatic sampling method can sample the parameter space according to the error indicator and terminate the procedure automatically when the maximum actual error becomes smaller than a specified tolerance. This method can place the sampling points where the error reaches the maximum. Therefore the approach can provide a better performance with smaller relative errors and cheaper offline and online computational cost.
4. Numerical Example
In this section, the novel PROM is illustrated with its application to the moving coil of electricaldynamic shakers. Since the geometric model has perfect cyclic symmetry property, it can be simplified to one over sixteen of the entire model. The simplified moving coil model consists of a drive coil, a stiffener plate, and a table. At the position shown in Figure 2(b), only axial displacement is permitted and the normal displacement at the symmetric boundary is zero. The model is excited by equivalent 1000 N surface force imposed on the interface of the drive coil and the stiffener plate. The measuring points are also marked in Figure 2.
(a)
(b)
For simplicity, we partition the entire model into three substructures, and the subdomains along with the corresponding shape parameters are marked in Figure 3. The material is linear elastic with E = 69 GPa, ν = 0.33, and ρ = 2.7e^{3} Kg/cm^{3}, roughly corresponding to aluminum. The finite element model consists of 15608 SOLID elements and 21246 nodes. The size of the fullorder model is 59487, and the number of interface DOF is 1140, which is 1.9% of the full DOFs.
(a)
(b)
(c)
Regarding the parameter space, each substructure contains one shape parameter. In this application, the thickness of stiffener plates T, changing from 10 mm to 32 mm, the length of drive coil L, changing from 242 mm to 342 mm, and the diameter of table D, changing from 920 mm to 1020 mm, span the threedimensional parameter space. The original parameter configuration is T = 18 mm, L = 292 mm, and D = 970 mm.
It should be noted that even though the shape parameters of the mesh vary, the underlying mesh topology remains the same. The preprocess was conducted in HYPERMESH to parameterize the mesh. To be specific, we first import the simplified geometry model into HYPERMESH and generate the original regular mesh. The HYPERMORPH module is subsequently utilized to parameterize the mesh without changing the underlying topology. Two FE models corresponding to the initial and final parameter sets are generated, respectively. Their node coordinates are then exported for use of generation of new FE model at arbitrary parameter values in ANSYS. The interface and measuring points are set as master nodes. The boundary condition and excitation condition shown in Figure 2(b) are also initialized in HYPERMESH.
The sampling scheme is automatically implemented in MATLAB including calling ANSYS to generate FE model and corresponding substructure matrices, importing the component matrices files, synthesizing the substructures, and conducting the automatic methods outlined in Figure 1. All numerical calculations are performed on a Windows Machine with Intel Core (TM) i54210U CPU at 2.60 GHz and 8 GB of RAM.
In this research, Latin hypercube sampling (LHS) is used to select the random candidates. The number of random candidates is set to N_{c} = 20 for all substructures. The reducedorder basis applied in this approach is built by 10order orthogonal eigenvectors after a Gram–Schmidt process. We set the exact errors of the assembled systems as the stopping criterion. Only components with the maximum error indicator are assembled for evaluation. The performance of this approach is assessed by monitoring the maximum errors of frequency, modal shapes, and responses between the directly computed CBCMS ROM and interpolated CBCMS ROM. The parameter sampling points for each component change from 2 to 7, and this number varies simultaneously for all components.
Figure 4 presents the evolution of exact maximum errors of the assembled system defined in equations (55)–(57). The results indicate that the maximum errors decrease significantly as the parameter sample points increase. It can also be observed that the responses has larger errors, compared with the other two variables. Therefore, it can be used as the stopping criteria of the algorithm.
The convergence performance of the proposed approach is also assessed. There is no guaranty that adding a parameter sample point will necessary lead to a monotonic decrease of maximum errors, but in practice, one can observe that the maximum error will decrease for most additional sample points due to the convergence property of the Lagrange interpolation method exploited in our approach.
In order to illustrate the advantage of the proposed CMSPROMs over the traditional global PROMs developed in [32], we compare the offline PROM database generation time and sampling points’ number.
In the global parametric reducedorder modeling procedure, 7 uniformly distributed sample points are set for each dimension of the threedimensional parameter space as same as the proposed automatic procedure. In ANSYS, it is relatively time consuming to change the nodes coordinate by use of NMODIF command. Therefore, the total time includes FE model generation and model reduction time.
One can observe from Figure 5 that the proposed procedure outperforms the global method thanks to its ability to partition the entire parameter space into separate ones with the help of CMS method. This new method can reduce significantly the sample points’ number from 7^{3} = 343 to 7 ∗ 3 = 21. The offline model generation and reduction time decrease a lot accordingly. The computational benefit of the proposed approach can be apparently observed with a speedup ratio of about 19 times. In general, the proposed approach can improve the calculation efficiency and reduce the offline ROMs database construction time to a great extent.
In addition, this novel CMSPROM surpasses not only the global PROM but also the classical CMSPROM [24] in the aspects of online calculation time and accuracy. To validate the performance, we first compare the online modal analysis time for these two CMSPROMs.
It can be observed from Figure 6 that the online calculation time of classical CMSPROM is almost 120 times larger than its counterpart calculated using the proposed CMSPROM. This is due to the fact that the classical method still involves timeconsuming online multiplications between stiffness or mass matrices and collected transformation vectors. As the number of interface DOFs and parameter samples increases, the calculation time quickly grows. The proposed approach, by contrast, performs the online calculation only for the reducedorder component matrices, thus leading to cheaper online calculation cost.
Furthermore, to assess the accuracy of these CMSPROMs, we subsequently perform accuracy comparisons for first 5order frequencies. The query parameter point is T = 30.31 mm, L = 335.22 mm, and D = 968.19 mm at which the error indicator reaches the maximum.
As shown in Table 1, the proposed approach leads to more accurate PROM than its classical counterpart. The results suggest that the proposed CMSPROM can achieve a higher level of accuracy at about than the classical CMSPROM at about . Moreover, the expected values of the frequency errors confirm the results observed in Figure 4.

Then, we compare the harmonic response for the direct CMS ROM and interpolated CMS ROM. The frequency range is from 1 Hz to 1900 Hz with the frequency interval 0.1 Hz. Rayleigh damping coefficient is applied as alpha = 2.5, beta = 9 ✕ 10^{−5}, and the detecting point locates at the center of top surface as shown in Figure 2(b).
To display the distinct comparison, different scales are adopted for the xaxis. For example, the results between 0 and 1500 Hz are plotted using a base 10 logarithmic scale while for the range of 1500∼1900 Hz, linear scales are used.
From Figure 7, one can observe that the harmonic response predicted by the proposed CMSPROM is almost indistinguishable from its counterpart calculated by direct CMS ROM. In other words, the novel CMSPROM can reproduce the model with almost negligible errors, while the relative errors of the harmonic response obtained using the classical CMSPROM are pretty large within some frequency ranges.
(a)
(b)
To investigate exact errors and assess the accuracy of these CMSPROMs, the first 5order resonance responses are compared in Table 2.

As shown in Table 2, the results associated with the proposed approach show a higher level of accuracy than the other classical approach. The results also reveal that the relative errors of the proposed PROM are below 0.1%, and very good predictions can be obtained using this novel approach.
To sum up, all of the above results illustrate the fact that the novel componentbased parametric reducedorder modeling method can not only sample the parameter space automatically and efficiently but also lead to more accurate ROMs for further calculation.
5. Conclusions
An automatic componentbased parametric reducedorder modeling method has been developed in this paper. On the basis of introducing the CBCMS, we combined it with the congruence transformation method and local interface reduction method to achieve the compatibility between discrete CMS ROMs. The automatic sampling procedure is then described in detail. During this process, we defined the error indicator to locate next optimal sampling points and the exact error to monitor the convergence rate and termination criteria. In addition, the entire analysis procedure was summarized in the schematics. For verification of the presented method, we analyzed the moving coil of electricaldynamic shaker. Numerical experiments illustrate that this novel PROM can provide perfect predictions of frequencies and harmonic responses. The proposed approach can automatically determine the minimum parameter sample points to satisfy the predefined error and combine the benefits of componentbased parametric reducedorder modeling method and those of automatic sampling procedure.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare no conflicts of interest.
Authors’ Contributions
All the authors have contributed equally to the conception and idea of the paper, implementing and analyzing the processing methods, evaluating and discussing the experimental results, and writing and revising this manuscript.
Acknowledgments
This project was supported by the National Natural Science Foundation of China (grant no. 11427801).
References
 W. H. Schilders, H. A. Van der Vorst, and J. Rommes, Model Order Reduction: Theory, Research Aspects and Applications, vol. 13, Springer, Berlin, Germany, 2008.
 W. Gressick, J. T. Wen, and J. Fish, “Order reduction for largescale finite element models: a systems perspective,” International Journal for Multiscale Computational Engineering, vol. 3, no. 3, 2005. View at: Publisher Site  Google Scholar
 A. Quarteroni, in Reduced Order Methods for Modeling and Computational Reduction, vol. 9, Springer, Berlin, Germany, 2014.
 X. Tang, X. Hu, W. Yang, and H. Yu, “Novel torsional vibration modeling and assessment of a powersplit hybrid electric vehicle equipped with a dualmass flywheel,” IEEE Transactions on Vehicular Technology, vol. 67, no. 3, pp. 1990–2000, 2018. View at: Publisher Site  Google Scholar
 T. Zhang and H. G. Li, “Adaptive pole placement control for vibration control of a smart cantilevered beam in thermal environment,” Journal of Vibration and Control, vol. 19, no. 10, pp. 1460–1470, 2013. View at: Publisher Site  Google Scholar
 T. BuiThanh, K. Willcox, and O. Ghattas, “Parametric reducedorder models for probabilistic analysis of unsteady aerodynamic applications,” AIAA Journal, vol. 46, no. 10, pp. 2520–2529, 2008. View at: Publisher Site  Google Scholar
 J. Annoni and P. Seiler, “A method to construct reducedorder parametervarying models,” International Journal of Robust and Nonlinear Control, vol. 27, no. 4, pp. 582–597, 2017. View at: Publisher Site  Google Scholar
 P. Benner, S. Gugercin, and K. Willcox, “A survey of projectionbased model reduction methods for parametric dynamical systems,” SIAM Review, vol. 57, no. 4, pp. 483–531, 2015. View at: Publisher Site  Google Scholar
 A. Bouazzouni, G. Lallement, and S. Cogan, “Selecting a ritz basis for the reanalysis of the frequency response functions of modified structures,” Journal of Sound and Vibration, vol. 199, no. 2, pp. 309–322, 1997. View at: Publisher Site  Google Scholar
 D. Amsallem and C. Farhat, “Interpolation method for adapting reducedorder models and application to aeroelasticity,” AIAA Journal, vol. 46, no. 7, pp. 1803–1813, 2008. View at: Publisher Site  Google Scholar
 J. Degroote, J. Vierendeels, and K. Willcox, “Interpolation among reduced‐order matrices to obtain parameterized models for design, optimization and probabilistic analysis,” International Journal for Numerical Methods in Fluids, vol. 63, no. 2, pp. 207–230, 2010. View at: Google Scholar
 H. Panzer, J. Mohring, R. Eid, and B. Lohmann, “Parametric model order reduction by matrix interpolation,” atAutomatisierungstechnik Methoden und Anwendungen der Steuerungs, Regelungs und Informationstechnik, vol. 58, no. 8, pp. 475–484, 2010. View at: Publisher Site  Google Scholar
 D. Amsallem, J. Cortial, K. Carlberg, and C. Farhat, “A method for interpolating on manifolds structural dynamics reducedorder models,” International Journal for Numerical Methods in Engineering, vol. 80, no. 9, pp. 1241–1258, 2009. View at: Publisher Site  Google Scholar
 U. Baur and P. Benner, “Model reduction for parametric systems using balanced truncation and interpolation,” atAutomatisierungstechnik, vol. 57, no. 8, pp. 411–419, 2009. View at: Publisher Site  Google Scholar
 U. Baur, P. Benner, A. Greiner, J. G. Korvink, J. Lienemann, and C. Moosmann, “Parameter preserving model order reduction for MEMS applications,” Mathematical and Computer Modelling of Dynamical Systems, vol. 17, no. 4, pp. 297–317, 2011. View at: Publisher Site  Google Scholar
 R. J. Guyan, “Reduction of stiffness and mass matrices,” AIAA Journal, vol. 3, no. 2, p. 380, 1965. View at: Publisher Site  Google Scholar
 E. J. Kuhar and C. V. Stahle, “Dynamic transformation method for modal synthesis,” AIAA Journal, vol. 12, no. 5, pp. 672–678, 1974. View at: Publisher Site  Google Scholar
 R. L. Goldman, “Vibration analysis by dynamic partitioning,” AIAA Journal, vol. 7, no. 6, pp. 1152–1154, 1969. View at: Publisher Site  Google Scholar
 R. H. MacNeal, “A hybrid method of component mode synthesis,” Computers & Structures, vol. 1, no. 4, pp. 581–601, 1971. View at: Publisher Site  Google Scholar
 R. Craig and M. Bampton, “Coupling of substructures for dynamic analyses,” AIAA Journal, vol. 6, no. 7, pp. 1313–1319, 1968. View at: Google Scholar
 S. Rubin, “Improved componentmode representation for structural dynamic analysis,” AIAA Journal, vol. 13, no. 8, pp. 995–1006, 1975. View at: Publisher Site  Google Scholar
 K. Park, “Componentbased vibration modeling methods for fast reanalysis and design of complex structures,” University of Michigan, Ann Arbor, MI, USA, 2008, Doctoral dissertation. View at: Google Scholar
 G. Zhang, “Componentbased and parametric reducedorder modeling methods for vibration analysis of complex structures,” University of Michigan, 2005, Doctoral dissertation. View at: Google Scholar
 S. K. Hong, “Parametric reducedorder models of structural dynamics for design, damage detection and structural optimization,” University of Michigan, Ann Arbor, MI, USA, 2012, Doctoral dissertation. View at: Google Scholar
 J. Lee and M. Cho, “An interpolationbased parametric reduced order model combined with component mode synthesis,” Computer Methods in Applied Mechanics and Engineering, vol. 319, pp. 258–286, 2017. View at: Publisher Site  Google Scholar
 P. Binev, A. Cohen, W. Dahmen, R. DeVore, G. Petrova, and P. Wojtaszczyk, “Convergence rates for greedy algorithms in reduced basis methods,” SIAM Journal on Mathematical Analysis, vol. 43, no. 3, pp. 1457–1472, 2011. View at: Publisher Site  Google Scholar
 B. Haasdonk, M. Dihlmann, and M. Ohlberger, “A training set and multiple bases generation approach for parameterized model reduction based on adaptive grids in parameter space,” Mathematical and Computer Modelling of Dynamical Systems, vol. 17, no. 4, pp. 423–442, 2011. View at: Publisher Site  Google Scholar
 K. Veroy and A. T. Patera, “Certified realtime solution of the parametrized steady incompressible NavierStokes equations: rigorous reducedbasisa posteriori error bounds,” International Journal for Numerical Methods in Fluids, vol. 47, no. 89, pp. 773–788, 2005. View at: Publisher Site  Google Scholar
 M. A. Grepl and A. T. Patera, “A posteriorierror bounds for reducedbasis approximations of parametrized parabolic partial differential equations,” ESAIM: Mathematical Modelling and Numerical Analysis, vol. 39, no. 1, pp. 157–181, 2005. View at: Publisher Site  Google Scholar
 M. Geuss, “A Blackbox method for parametric model order reduction based on matrix interpolation with application to simulation and control,” Technische Universität München, Munich, Germany, 2015, Doctoral dissertation. View at: Google Scholar
 T. BuiThanh, “Modelconstrained optimization methods for reduction of parameterized largescale systems,” Massachusetts Institute of Technology, Cambridge, MA, USA, 2007, Doctoral dissertation. View at: Google Scholar
 E. Balmes, F. Ravary, and D. Langlais, “Uncertainty propagation in modal analysis,” in Proceedings of the IMACXXII: A Conference and Exposition on Structural Dynamics, p. 57, Cham, Switzerland, 2004. View at: Google Scholar
 D. Amsallem and C. Farhat, “An online method for interpolating linear parametric reducedorder models,” SIAM Journal on Scientific Computing, vol. 33, no. 5, pp. 2169–2198, 2011. View at: Publisher Site  Google Scholar
 M. P. Castanier, Y.C. Tan, and C. Pierre, “Characteristic constraint modes for component mode synthesis,” AIAA Journal, vol. 39, no. 6, pp. 1182–1187, 2001. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Ying Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.