# A domain decomposition solver for spectral stochastic finite element: an approximate sparse expansion approach

*Dis Prev Res*2024;3:3.

## Abstract

Spectral stochastic finite element (SSFE) has been widely used in the uncertainty quantification of real-life problems. However, the prohibitive computational burden prevents the application of the method in practical engineering systems because an enormous augmented system has to be solved. Although the domain decomposition method has been introduced to SSFE to improve the efficiency for the solution of the augmented system, there still exist significant challenges in solving the extended Schur complement (e-SC) system from domain decomposition method. In this paper, we develop an approximate sparse expansion-based domain decomposition solver to generalize the application of SSFE. An approximate sparse expansion is first presented for the subdomain-level augmented matrix so that the computational cost in each iteration of the preconditioned conjugate gradient is greatly alleviated. Based on the developed sparse expansion, we further establish an approximate sparse preconditioner to accelerate the convergence of the preconditioned conjugate gradient. The developed approximate sparse expansion-based domain decomposition solver is then incorporated in the context of SSFE. Since the difficulties of solving the e-SC system have been overcome, the developed approximate sparse expansion-based solver greatly improves the computational efficiency of the solution of the e-SC system, and thereby, the SSFE is capable of dealing with large-scale engineering systems. Two numerical examples demonstrate that the developed method can significantly enhance the efficiency for the stochastic response analysis of practical engineering systems.

## Keywords

*,*domain decomposition

*,*approximate sparse expansion

*,*preconditioned conjugate gradient

*,*Schur complement

## INTRODUCTION

Common geological hazards, such as debris flows, landslides, and avalanches, are usually devastating, and the timing, location, and scale of disasters are uncertain, posing a significant threat to people’s lives and property^{[1]}. Therefore, researching the ability of structures to resist disasters is particularly important. However, uncertainty is a common occurrence in various aspects of practical engineering, including structure design, manufacturing, operation, and maintenance^{[2,3]}, which can lead to significant deviations in the structural behavior during structural response analysis and directly result in inaccurate analysis results of the ability of the structure to withstand disasters.

So, the realistic design and analysis of these physical systems must consider uncertainties contributed by various sources such as manufacturing variability, insufficient data, unknown physics, and aging^{[4–7]}. The uncertainties may also significantly influence the predictive capabilities of computer simulations^{[8–12]}. The latest advancements in high-performance computing and sensing technology have stimulated computational simulations with extremely high resolution, providing great possibilities for integrating effective uncertainty quantification methods to achieve realistic and reliable numerical predictions^{[13–18]}. Nevertheless, using standard Monte Carlo simulations for uncertainty quantification in the simulation of random responses to large and complex problems may be time-consuming or impractical.

As an alternative method to Monte Carlo simulation, spectral stochastic finite element (SSFE) has been widely applied in the field of stochastic mechanics^{[19–21]}. Different from various methods that rely on sampling such as the so-called collocation, SSFE is intrusive^{[22,23]}. It estimates the random response of the structure by transforming the random stiffness equation into a coupled set of a series of deterministic equations^{[4]}. In the framework of SSFE, if the uncertain properties of the system follow a Gaussian distribution, then they are quantified with the Karhunen-Loeve (KL) expansion; if not, for example, the uncertain properties follow a lognormal distribution, they will be projected into a polynomial chaos (PC) basis^{[24,25]}. The random response of the system is usually projected into the PC basis because it is unknown and almost never a Gaussian distribution; adopting the Galerkin minimization scheme will establish an augmented system of linear algebraic equations, and only solving this augmented system will determine the coefficients of the basis and then obtain the random response of the system^{[19]}.

However, compared to deterministic equations, the scale of the augmented ones that SSFE needs to solve will be several orders of magnitude larger, because in the analysis of complex structures, the PC expansion (PCE) items of the random response usually reach thousands or even tens of thousands. For such large-scale problems, in the past few years, the customized solution algorithm has been developing continuously, and to some extent, it has achieved successful results in solving such high computational requirements. A kind of method is directly solved by the preconditioned conjugate gradient (PCG) method, and the preconditioned matrix is established using the structure of the coefficient matrix of the augmented system, including the block diagonal preconditioner, successive symmetric over-relaxation, and so on^{[4,25]}. Such preconditions can effectively reduce the number of iterations of PCG, transforming what actually takes time to calculate into a series of deterministic problems, thus greatly improving the solution efficiency of the augmented system.

Nevertheless, with the increasing scale of the augmented system, the deterministic problems demanding resolution through the direct PCG method are growing, which has gradually approached the number of samples that need to be solved by the Monte Carlo method. Domain decomposition method (DDM) provides another essential framework for developing fast and efficient solvers for such applications. The DDM is an important strategy for solving large-scale problems in computational mechanics^{[26,27]} and has been applied to the solution of SSFE in recent years. Using this approach, the solution of augmented systems is transformed into a series of sub-problems about each sub-domain and an extended Schur complement (e-SC) system^{[28]}. The solution of the e-SC system takes up most of the calculation time when the number of sub-domains is huge. Hence, the efficiency for the solution of the e-SC system is the key in this kind of method. Till now, all these methods are based on the structure of the coefficient matrix of the e-SC system, and PCG is used to solve the e-SC system, which improves the computational efficiency to a certain extent in the random response analysis of structures^{[29–34]}. However, the above techniques can still not effectively analyze the large-scale structures, and the efficiency is even lower than that of direct conjugate gradient (CG) method. This is because the preconditioner established from the existing methods cannot guarantee the similarity with the e-SC matrix^{[35]}, and as a result, the iteration steps of PCG will be quite huge. In addition, since direct DDM needs to store a large-scale dense matrix in the solution process, the computational cost in each iteration will be quite large, which may further decrease the efficiency of PCG in the stochastic analysis of large-scale structures.

In order to overcome the above deficiencies of the existing methods, we develop an approximate sparse expansion-based domain decomposition solver for stochastic finite element strategy. Firstly, we develop an approximate sparse expansion for the subdomain-level augmented matrix. By utilizing the property of the subdomain-level random stiffness matrix, the subdomain-level augmented matrix can be approximated as a product of a block diagonal matrix and a sparse matrix, which is denoted as approximate sparse expansion in this study. With this expansion, we further establish an algorism of multiplying the e-SC matrix and an arbitrary vector, which may greatly save the computational cost in each iteration in PCG. Secondly, a preconditioner for the PCG solution of the e-SC system is further developed based on the established approximate sparse expansion. Since the preconditioner is constructed as the product of a block diagonal matrix and a sparse matrix, the inverse of the preconditioner can be readily obtained. More importantly, the preconditioner is approximately equal to the e-SC matrix so that the number of iterations can be significantly reduced. As a result, our approach overcomes the difficulties of solving the e-SC system faced by traditional DDM-based SSFE by establishing the approximate sparse preconditioner, and the efficiency for the solution of e-SC system can be greatly improved compared with the traditional DDM-based SSFE. Therefore, the structural stochastic response analysis will become much more efficient.

The remainder of this paper is organized as follows: In Section 2, the primal domain decomposition for SSFE is briefly described. Then, in Section 3, a domain decomposition solver of stochastic finite element method (SFEM) based on the approximate sparse expansion approach is proposed, and the actual calculation steps are summarized. Section 4 contains numerical examples demonstrating the efficiency of the proposed method. Finally, the concluding remarks are given in Section 5.

## DOMAIN DECOMPOSITION-BASED SSFE

### Stochastic finite element: spectral approach

In the response analysis of a structural system with stochasticity in the model parameters, the equilibrium equation is usually expressed as

where

In one particular formulation of the SFEM, according to the KL of the stochastic process ^{[1]}, can be given as

where

The resulting vector

where ^{[25]}:

Substituting Equations (2) and (3) into (1) and forcing the residual to be orthogonal to the approximating space spanned by the PC ^{[36]}:

These equations can be assembled into a matrix of size

where

It is evident that the augmented system is an equation that couples the determinism of finite elements with the randomness of structures. The number of rows and columns of its coefficient matrix are the product of: (1) the spatial degrees of freedom (DOF) after meshing the structure; and (2) the number of expansion terms for random responses using the PC method. When the structural grid is subdivided densely and multiple terms are expanded for random responses, the size of the augmented system can become quite large, which exceeds the computational capacity of current computers. Ultimately, this fact makes it impossible for SSFE to analyze the stochastic response of large-scale structures. Therefore, it is quite necessary to develop efficient algorithms for the augmented system.

### Domain decomposition approach for SSFE

When dealing with large-scale problems, DDM is one of the best candidates due to its proven performance and scalability. In this method, the computational domain ^{[26,27]}.

According to this partition, the DDM numbers, first the subdomains and last the interface, that is, let

where the restriction operators ^{[31]}, expressed as

Rearranging the augmented system Equation (5) based on the order of

where the expansion restriction operators ^{[37]}, determined by

where

Then, using Gaussian elimination technique, Equation (7) becomes a global extended e-SC system, as given in

where the e-SC matrix

By solving the e-SC system Equation (9), the global interface solution coefficients

As given in Equation (9), the key idea of DDM is to transform the original super-large-scale augmented system into a smaller e-SC system so as to reduce the complexity for calculating the stochastic response of structures. However, even if DDM has improved the efficiency of analyzing the stochastic response of structures through the reduction of complexity, since the coefficient matrix is dense and has lost its structure similar to that of the augmented matrix, a suitable preconditioner has been lacking to date. The ultimate result is that the analysis of stochastic response in large-scale structures is not efficient enough.

## APPROXIMATE SPARSE EXPANSION SOLVER FOR DDM-BASED SSFE

As mentioned earlier, the reason why the analysis of stochastic response in large-scale structures is not efficient enough is that, compared with the augmented matrix, the e-SC matrix will become quite dense, and the structure easy to establish preconditioner is lost. As a result, there is no suitable preconditioner for the e-SC system till now. If a precondition suitable for the structure of the e-SC system can be established according to the new structure of the e-SC matrix, its solving efficiency will be greatly improved, and then the efficiency of structural random response analysis will be improved. So, in Section 3.1, we first establish the approximate sparse matrix for the augmented matrix at the subdomain level, and then in Section 3.2, we further develop an approximate sparse preconditioner for the e-SC system. Since the new preconditioner is established, the computational efficiency of the e-SC system will be greatly improved; thus, the analysis of structural stochastic response will be more efficient.

### Approximate sparse expansion of subdomain-level augmented matrix

In this section, we develop an approximate sparse expansion of the subdomain-level augmented matrix

In order to establish the approximate sparse expansion, we first develop an improved mean-based approximation. Considering the original mean-based approximation ^{[37]}, because it is a diagonal matrix, given as

In order to establish the approximate value

where

The value of

This result further indicates that if the matrix

Then, its sub-matrices can also be expressed in this approximate form. The submatri

It means that any matrix

Then, we establish the approximate relation matrix

The submatrix of the relation matrix

where in this equation,

Then, according to the magnitude of these elements, they are sorted, and a certain amount of the big elements are retained. The preserved elements constitute another set

The above steps represent the implementation process of approximate sparsity, and we define the steps as a function, denoted as

Then, the approximate sparse relationship matrix of the sub-domain augmented matrix will be established through the following steps as given in Algorism 1 (in the APPENDIX). Now, all the elements where each node is located are defined as a set

Through Algorism 1, a series of relation matrices of random stiffness matrices at the subdomain level are obtained as

As given in Equation (18), approximate sparse expansion

#### Multiplication of e-SC matrix and vector

As given in Equation (9), the cost of solving the e-SC system by PCG depends on: (1) the computational efficiency of the multiplication of the e-SC matrix and vector; and (2) the preconditioner of the system. So, in this section, in order to improve the efficiency of solving e-SC systems, we develop a new algorism of multiplication of the e-SC matrix and vector based on approximate sparse expansion in Equation (18). The algorism is efficient, avoids the direct involvement of large matrices in operations, and significantly reduces the cost of stochastic response analysis of Structure.

During each iteration step in PCG, the multiplication of the e-SC matrix with a temporary vector,

In this work, the step of directly multiplying the e-SC matrix with the vector is transformed into a series of steps of multiplying the small-scale augmented matrix with the vector, and the coefficient matrix is the solution problem of the small-scale augmented system. By rewriting the algorithm of multiplying these small-scale augmented matrices with vectors and proposing an approximate sparse preconditioned matrix for small-scale augmented systems, the computational efficiency of the multiplication of e-SC matrices and vectors is improved.

Considering the multiplication of an e-SC matrix and a temporary vector

In the actual calculation process, the e-SC matrix at the subdomain level

Then, the multiplication of argument matrix

The solution of a series of small-scale augmented systems in each iteration step of PCG, which are not large in scale but need to be solved many times in the iterative process, will lead to the time-consuming calculation and finally have a great impact on the computational efficiency of the e-SC system. So, in this work, PCG is used to solve the equation, and the matrix-vector multiplication step is conducted by Algorism 2. Then, the approximate sparse matrix in Equation (18) is used as the preconditioner of the solution.

In recent years, the parallelization of DDM is usually calculations involving the interior of sub-domains. In order to calculate larger-scale structures, the strategy of multi-computer parallelization is often adopted. The computations inside each subdomain usually do not interfere with each other. Therefore, in the process of DDM, this portion of the calculations is typically parallel, and because of this, in this work, the steps of multiplying e-SC matrices and temporary vectors are finally transformed into the following parallel steps of Algorism 3 (in the APPENDIX).

As given in Algorism 3, by virtue of the approximate sparse expansion of the subdomain-level augmented matrix in Equation (18), we establish an algorism of multiplying the e-SC matrix and vector for the PCG solver of e-SC system Equation (9). The algorithm can avoid direct operations on large dense matrices, improve the computation efficiency of each iteration step in PCG, and thus improve the efficiency of stochastic response analysis of structures.

### Approximate sparse preconditioner of e-SC system

Once the approximate sparse expansion of the subdomain-level augmented matrix is established, as determined by Equation (18), the approximate sparse preconditioner of e-SC can be established in this section. The preconditioner is expressed in the form of multiplication of a block diagonal matrix with a sparse matrix called the approximate relation matrix. The preconditioner is quite close to the e-SC matrix, which makes the iteration of PCG converge within a relatively small number of iteration steps. Furthermore, because the inverse matrix of the approximation relation matrix can be easily established, the computational complexity of each iteration step is considerably low. These facts make the solving efficiency of the e-SC system extremely high and ultimately substantially improve the efficiency of structural stochastic response analysis.

#### Approximate sparse expansion of e-SC matrix

When solving e-SC systems, the more important factor determining the computational efficiency of PCG is the choice of preconditioned matrix

So, in this section, we derive the approximate form of the e-SC matrix based on the approximate sparse augmented matrix in Equation (18). Firstly, based on the derivation process of the approximately sparse augmented matrix, the approximately sparse forms of the matrices in Equation (10) are further established. Then, the approximation values of each matrix in Equation (10) are substituted into the calculation formula of the e-SC matrix, and the approximate form of the e-SC matrix is obtained through computation.

Consider further adjusting the order of solutions in the augmented system to carry out the overall Schur complement (SC) process, expressed as

Different from Equation (7), in (21), let the solution vector

According to the approximate value of

where

For other submatrices

After dividing the subdomains, we can expand Equation (23) to obtain the block form, expressed as

where

According to the diagonal dominance of

Then, Equations (25) and (22) are substituted into (21) and the SC process is carried out; with some simplification steps, the approximate form of the e-SC matrix can be finally obtained as

where

According to the approximate simplification result in Equation (26), the e-SC matrix will approximate a special structure, which is a multiplication of two matrices: the one is a diagonal matrix composed of m-SC matrix

The relation matrix of the e-SC matrix can be approximated as a sparse matrix with the original matrix, which can be used to establish the precondition for PCG to solve the e-SC system. However, different from establishing the relational matrix of a subdomain augmented matrix, it is usually not feasible to directly multiply the block diagonal matrix

We utilize the similarity between approximate relation matrix in (26) to infer the relationships between various sub-matrices of the actual relationship matrix. Using this relationship, unnecessary calculations can be avoided, thereby improving computational efficiency.

According to Equation (26), if we calculate the theoretical values of each submatrix of the relational matrix based on Equation (26), we can obtain that

To avoid the direct storage of large dense matrices, we derived an approximate sparse expansion form of subdomain-level e-SC through the conclusion of Section 3.1. Each term of the calculation formula

According to Equations (12) and (14), the sub-matrix in (28) can be approximately obtained as

According to Equation (29), by left-multiplying a diagonal matrix

Thus, according to Equations (26), (27) and (30), any sub-matrix of the relation matrix of an e-SC matrix can be denoted as

where the matrix

Due to the approximate diagonal structure given Equation (31), it can be transformed into an approximately sparse form

The matrix

However, another problem is, for the e-SC matrix, when many subdomains are divided, the number of global interface nodes will be quite large, which leads to the large scale of any sub-matrix of the e-SC matrix, along with the m-SC. The fact means that the necessary step of computing the

#### Approximate relation matrix for e-SC matrix

In order to efficiently calculate the approximate values of the relation matrix

(1) *Computation of the inverse of m-SC matrix*

The division of two-level subdomains is based on the division of one-level subdomains, and several adjacent first-level subdomains form a second-level subdomain. The order of e-SC matrix solutions is rearranged by two-level division; then, each sub-matrix has an arrow-shaped block structure according to the element distribution law of its coefficient matrix.

On the basis of finite element mesh, the structure is divided into a certain number

Figure 1. The two-level domain decomposition of a square plate (

After the structure is divided into two levels, the order of Schur complementary system solutions is adjusted according to the results of the division, so that its coefficient matrix becomes a block matrix with an arrow structure element distribution in each submatrix. In order to explain this structure more clearly, firstly, the element distribution of the mean interface stiffness matrix

According to the deterministic domain decomposition, the SC process of Equation (32) will get the following m-SC system, as formulated in

By the two-level domain decomposition, the m-SC matrix

According to the distribution law of zero elements given in Equation (34) and the division of the second-level domain decomposition shown in Figure 1, we reorder the solutions of the interface

where

After rearranging the order of SC matrices according to the new order of interface solution vectors derived in Equation (35), based on the law of element distribution in Equation (34), the matrix will become the following arrow shape expressed as

where

The form of the block matrix given in Equation (36) is quite easy to find the inverse. Because it can be easily expanded into a simpler form by block Cholesky decomposition, denoted as

where

It can be easily observed that the inverse of this matrix can be expressed as

where only a series of internal stiffness matrices for second-level sub-domains

Once the inverse of the m-SC matrix is obtained as Equation (40), the multiplication of the inverse of the m-SC matrix and a vector (or matrix) can be transformed into a more efficient form. Taking multiplication of the matrix and a vector

where only multiplication of the small-scale matrix and vector and inverse of the small-scale matrix exist, and clearly, it is more efficient to multiply than computing the inverse matrix and then multiplying it with the vector directly.

(2) *Determination of relation matrix*

Based on our previous research findings, it has been concluded that the e-SC matrix can be represented in a two-level block form, where each sub-block has the same distribution of elements as the determined SC matrix, obtained as

where each submatrix can be transformed into this arrow-shaped representation akin to Equation (36) through the two-level region decomposition, as derived in

It is noteworthy that substituting matrices of the same number of rows as

Directly approximating sparseness of the relationship matrix would consume a considerable amount of time due to the lowest computational complexity of sorting, which is at least

Since the inverse of second-level SC matrix

where a series of matrix multiplication problems is involved, including the multiplication of Boolean matrix and other matrices. In order to improve the calculation efficiency of this step, we change the original calculation order, Using the properties of Boolean matrix,

Similar to formula a, there is also a Boolean matrix in the calculation of Equation (46), so although what needs to be calculated is the multiplication of

For the last column of the relationship matrix in Equation (44), the interface solution

Bring the calculation results into the sub-formula (3) of Equation (41), the rest of the solution can be obtained as

It can be observed that in Equation (47), the same contents

During the calculations, according to the specified order of computation, Equations (48) and (46) can be written in an approximate sparse form first as

Further approximate sparsity processing is performed on the matrix

In the calculation process, all the computations involved are operations between small matrices and the approximate sparse processing of small matrices. Since the relation matrix of each submatrix is established by calculating the same amount of content separately, in the actual calculation, this part will be calculated in parallel, and the algorithm for finally calculating the relationship matrix of an e-SC matrix is expressed as Algorithm 4 (in the APPENDIX), where we first choose the row and column numbers of the eigen relation matrices by step 5; then, for every eigen relation matrix, which is split into a block matrix akin to Equation (44), the submatrices in the bottom row are first obtained by steps 6 and 8; using the result, the other row can be calculated by steps 9 and 11; steps 12 16 involve creating the approximate value of these submatrices, then use these submatrices to assembly the approximate eigen relation matrix by steps 17 and 18, finally create the approximate relation matrix by step 20.

Once the approximate relation matrix

which denotes the proposed approximate sparse preconditioner that is quite close to the e-SC matrix; therefore, using this matrix as the preconditioner for PCG will lead to convergence within a relatively small number of iterations. On the other hand, the matrix is a product of a diagonal matrix and a sparse matrix with very little condition number, which makes the solution of the preconditioning equation relatively easy. Ultimately, adopting this preconditioner can make PCG more efficient in solving the e-SC system, which greatly improves the solution efficiency of the random response of the structure.

### Approximate sparse approach for DDM-based SSFE

According to the derivation of Sections 3.1 and 3.2, we develop an approximate sparse approach for DDM-based SSFE; the proposed method targets the step of solving e-SC systems with PCG. An efficient algorism for the multiplication of an e-SC matrix and an arbitrary vector, and an approximate sparse preconditioner are established to improve the computing efficiency of PCG, thus making the structural random response more efficient. The following steps show the new computational framework based on the approximate sparse approach for DDM-based SSFE in this work:

In the established algorithm, PCG in step 8 is used to solve the e-SC system; different from the traditional solution, we propose a new calculation scheme based on the approximate sparse stiffness matrix for the two most critical steps in PCG: the e-SC matrix and temporary vector multiplication and the solution of preconditioned equations, to solve the problems that the existing solution methods must save large dense matrices and their calculation efficiency are low.

In step 4, according to the derivation in Section 3.1, the approximate relation matrix of expansion term of the random stiffness matrix in subdomain

On the basis of Section 3.1, after further derivation, the proximate relationship matrix of e-SC matrix

To sum up, by rewriting the step of the e-SC matrix and the vector multiplication based on the approximate sparse random stiffness matrix, we avoid the large dense matrix directly participating in the operation and improve the calculation efficiency. At the same time, based on the approximate sparse random stiffness matrix, the relationship matrix of the e-SC matrix is further established to create a preconditioning matrix. Using this matrix as the preconditioner for solving the e-SC matrix by PCG can greatly improve the calculation efficiency by controlling the iteration steps. Finally, an efficient framework for structural random response analysis is established.

## NUMERICAL EXAMPLES

To assess the computational efficiency of the proposed approximate sparse expansion-based domain decomposition solver for SFEM, we consider two numerical examples. In the two examples, the precision of the results (structural stochastic response) and the distribution of elements in the relation matrix of the e-SC matrix will be verified first, and then the convergence rate, computational time needed, and relative speedup compared with the original DDM-based SSFE will be investigated and manifested. The computational platform used for this investigation is an Intel Core i9-13900 with 24 physical cores at 3.0 GHz, coupled with 128 GB of RAM. The finite element mesh of the two examples is generated using Abaqus; the finite element method (FEM) assembly procedure is adapted from the FEON framework in python.

### Stochastic response of a thin square plate

As a test case, we have considered a stochastic response analysis of a thin square plate. The dimensions of this plate are

where standard deviation is

The physical domain with unit-square geometry is discretized using an unstructured finite element mesh. After partitioning the mesh, the number of nodes is 5, 041, and the count of Quadrilateral elements totals 4, 900; the number of subdomains involves 16, 25, 36, 49, and 64; the mesh and the division of subdomains are shown in Figure 2B illustrating the situation of 49 subdomains.

The random field

Figure 3. KL expansion result of the random field

The response is expressed using sixth order PCE (

Figure 4. The solution process: Coefficients before the first 4 PC basis. (A)

In order to verify the accuracy of the calculated results more clearly, we consider the node No. 2485 (point A in Figure 4A), which has a large average response, and draw the probability density function (PDF) of its displacement in the y direction according to the calculation results of different methods shown as Figure 5.

Figure 5. Stochastic response of y-direction displacement of node No. 2485: (A) the PDF and (B) tail of the PDF.

From the PDF shown in the figure, it can be seen that the random response results in the y-direction of node number 2485 obtained by the proposed method are almost consistent with the traditional SSFE method. Moreover, due to the maximum mean displacement in the y-direction of this node, it can be considered that the random response results of other nodes obtained by the proposed method also meet the accuracy requirements. Therefore, it indicates that the proposed approach has the same accuracy as the traditional SSFE method and can meet the accuracy requirements when the response PCE order is sufficient.

Next, the distribution law of the coefficient matrix elements of the e-SC system is verified. Firstly, the element distribution of the augmented matrix at the subdomain level is verified; the relation matrix

Figure 6. Element distribution of relation matrix of random partial stiffness matrix at the subdomain level: (A)

Figure 6 presents that for the relation matrices of the subdomain-level stochastic stiffness matrices, the non-diagonal elements are almost zero, while the values of the diagonal elements are relatively large compared to the non-diagonal elements, which confirms our findings.

Based on the above verification results, we further verify whether the properties of the relationship matrix of the e-SC matrix are consistent with the conclusion of proof in Section 3.2. Firstly, it is necessary to verify the distribution of the elements in the relation matrix

Figure 7. Element distribution of some relation matrix’s submatrices: (A) the first row, second column

It is evident that the sub-matrices of the relation matrix of the e-SC matrix also exhibit similar patterns of element distribution; that is, the non-diagonal elements are almost zero, while the values of the diagonal elements are relatively large compared to the non-diagonal elements.

Therefore, it can be said that, as demonstrated in this article, the distribution pattern of the elements in the relation matrix has been confirmed. Thus, the approximate sparse preconditioner established in Section 3.2 is approximately equal to the e-SC matrix, thus greatly improving the efficiency of solving the e-SC system.

On this basis, we study the magnitude of improvement in computational efficiency provided by the method proposed in our research. We compared the computational time of the traditional DDM-based SSFE and the proposed method, including the time of the SC process, the direct CG solver and the PCG solver with the proposed preconditioner of the e-SC system, as shown in Table 1.

Statistics of competitional time: DDM-based SSFE and approximate sparse expansion approach

Number of subdomains | Time of SC process | Time of solving e-SC matrix | |

Proposed preconditioner | Traditional method | ||

16 | 1,202.1 | 258.7 | 1,532.7 |

25 | 400.3 | 342.7 | 2,819.3 |

36 | 318.1 | 536.7 | 5,201.9 |

49 | 338.5 | 744.3 | 7,254.2 |

64 | 542.4 | 1,120.6 | 11,533.4 |

As demonstrated in Table 1, it can be obviously seen that compared with the traditional DDM-based SSFE, which uses the direct CG solver to solve the e-SC system, PCG with the proposed approximate sparse preconditioner is highly efficient, and as a result, the efficiency of stochastic response analysis of structures is greatly improved. Especially when there are a large number of subdomains and the solving time of the e-SC system becomes a higher proportion, this improvement in computational efficiency will become more evident from the fact shown in Table 1.

### Stochastic response of a frame structure

The second example features the stochastic response analysis of a frame structure, where we pay greater attention to its computational efficiency under different circumstances, including varying standard deviations and correlation lengths, than in Section 4.1.

As shown in the following Figure 8, the structure we consider in this example is a five-layer frame structure with a layer height of 3 m, the distance between each column is 2 m, the cross-sectional widths of the beams and columns are 0.3 m, thickness is 0.1, and Poisson's ratio is 0.25. The structure is subjected to the self-weight

The KL expansion is used to express the random field; in this example, the number of terms

Figure 9. The coefficient before the first PC basis of solution process: (A) Component in the x-direction (B) Component in the y-direction.

Then, we use PDF of the y-direction displacement at point A and the x-direction displacement at point B to describe the calculated results (points A and B are shown in Figure 9; A is the point with the largest average displacement in the y-direction shown in Figure 9A, while B is the point with the largest average displacement in the x-direction shown in Figure 9B) and compare them with the results of 50000 Monte-Carlo simulations, as shown in Figure 10.

Figure 10. The PDF of (A) x-direction displacement of point A and (B) y-direction displacement of point B.

It is obvious that the solution process calculated by the proposed approach is similar to the MC method in Figure 10, indicating that the proposed one can meet the accuracy requirements when the PCE order of the response is sufficient. Then, we are concerned about the numerical scalability of the improvement of computational efficiency, that is, the computational efficiency improvement of the proposed method under different conditions (the PCGM is considered to have converged when the error satisfies:

Statistics of competitional time of the different standard deviation: DDM-based SSFE and approximate sparse expansion approach

Standard deviation | Time of SC process | Time of solving e-SC matrix | |

Proposed preconditioner | Traditional method | ||

946.4 | 302.1 | 15,893.6 | |

953.1 | 539.4 | 18,106.9 | |

935.4 | 619.8 | 21,011.6 | |

940.2 | 703.5 | 21,933.9 |

Statistics of competitional time of the different correlation lengths: DDM-based SSFE and approximate sparse expansion approach

Correlation length | Time of SC process | Time of solving e-SC matrix | |

Proposed preconditioner | Traditional method | ||

928.6 | 496.4 | 16,555.3 | |

936.3 | 466.1 | 14,879.8 | |

938.2 | 562.9 | 19,194.6 | |

936.4 | 584.5 | 19,527.4 | |

935.4 | 619.8 | 21,011.6 |

From the computational time presented in Table 2 and Table 3, it is evident that the proposed approximate sparse preconditioner exhibits high efficiency; moreover, the degree of improvement does not show significant variations as the standard deviation and the correlation length change, thus indicating its numerical scalability.

## CONCLUSIONS

In this paper, we develop an approximate sparse expansion-based domain decomposition solver in the context of SSFE. By establishing an approximate sparse expansion of the subdomain-level augmented matrix, we transform the multiplication of the e-SC matrix and vector to the operation of a series of subdomain-level augmented matrices and vectors. The transformation can greatly decrease the computational cost in each iteration in PCG. With the approximate sparse expansion, we further develop an approximate sparse preconditioner in the context of PCG for the solution of the e-SC system. The preconditioner is established as a product of a block diagonal matrix and a sparse matrix with low condition number. The convergence of PCG can be remarkably accelerated using the preconditioner. Since the two main difficulties in the existing DDM-based SSFE are overcome, the developed solver addresses the challenge of difficulty solving the e-SC system in traditional methods, and thereby the uncertainty quantification of practical engineering is quite efficient using SSFE. Two numerical examples, including stochastic analysis of a thin square plate and a planar frame structure, have been studied to illustrate the effectiveness of the developed method. In both examples, the developed approximate sparse expansion-based domain decomposition solver enables greatly improving the efficiency of stochastic response analysis, and the improvement of efficiency is not affected by the number of subdomains and the variance of the random field.

Due to our focus on improving the computational efficiency of the proposed algorithm compared to traditional DDM-based SSFE, we only consider the simplest case, assuming that the structure is elastic and only withstands deterministic static loads.

In future work, we will further evaluate the computational performance of the algorithm under complex conditions (such as nonlinear structures subjected to random dynamic loads^{[38]}). Additionally, it is necessary to be aware that the proposed method has certain limitations; for example, the SC process is still inevitably required in establishing the approximate sparse preconditioner. So, in future work, we will also consider developing an approximate SC process to reduce its computational cost while meeting the accuracy requirements.

## DECLARATIONS

### Authors’ contributions

Conceptualization, methodology, analysis and writing: Luo B

Software coding, numerical analysis and validation: Cao W

Methodology, supervision and review: Dai H, Zhou Z

### Availability of data and materials

Not applicable.

### Financial support and sponsorship

This research was supported by a Grant from the National Natural Science Foundation of China (Project 12272109). This support is gratefully acknowledged.

### Conflicts of interest

All authors declared that there are no conflicts of interest

### Ethical approval and consent to participate

Not applicable.

### Consent for publication

Not applicable.

### Copyright

© The Author(s) 2024.

## REFERENCES

1. Cai J, Hu S, Sun F, Tang L, Fan G, Xing H. Exploring the relationship between risk perception and public disaster mitigation behavior in geological hazard emergency management: a research study in Wenchuan county. *Dis Prev Res* 2023;2:21.

2. Du Y, Xu J. Structural reliability analysis with epistemic and aleatory uncertainties via AK-MCS with a new learning function. *Dis Prev Res* 2023;2:15.

3. De Santis Y, Aloisio A, Pasca DP, Fragiacomo M, Dombrowski F. Evaluation of the shear size effect in glued laminated timber using a stochastic FE model calibrated on 17000 glue-line tests. *Constr Build Mater* 2023;399:132488.

4. Ghosh D, Avery P, Farhat C. A FETI-preconditioned conjugate gradient method for large-scale stochastic finite element problems. *Int J Numer Meth Eng* 2009;80:914-31.

5. Zhang R, Dai H. Stochastic analysis of structures under limited observations using kernel density estimation and arbitrary polynomial chaos expansion. *Comput Method Appl Mech Eng* 2023;403:115689.

6. Wang Y, Li J. Probabilistic analysis on stochastic failure modes of reinforced concrete beams under fatigue loading. *J Eng Mech* 2022;148:04022018.

7. Zheng Z, Dai H, Beer M. Efficient structural reliability analysis via a weak-intrusive stochastic finite element method. *Probabilist Eng Mech* 2023;71:103414.

8. Sousedík B, Ghanem RG, Phipps ET. Hierarchical Schur complement preconditioner for the stochastic Galerkin finite element methods. *Numer Linear Algebra Appl* 2014;21:136-51.

9. Peng Y, Zhou T, Li J. Surrogate modeling immersed probability density evolution method for structural reliability analysis in high dimensions. *Mech Syst Signal Proc* 2021;152:107366.

10. Dai H, Zhang R, Beer M. A new method for stochastic analysis of structures under limited observations. *Mech Syst Signal Proc* 2023;185:109730.

11. Pang R, Zhou Y, Chen G, Jing M, Yang D. Stochastic mainshock-aftershock simulation and its applications in dynamic reliability of structural systems via DPIM. *J Eng Mech* 2023;149:04022096.

12. He J, Chen J, Ren X, Li J. Uncertainty quantification of random fields based on spatially sparse data by synthesizing bayesian compressive sensing and stochastic harmonic function. *Mech Syst Signal Proc* 2021;153:107377.

13. Pranesh S, Ghosh D. Addressing the curse of dimensionality in SSFEM using the dependence of eigenvalues in KL expansion on domain size. *Comput Method Appl Mech Eng* 2016;311:457-75.

14. Meng Z, Zhao J, Chen G, Yang D. Hybrid uncertainty propagation and reliability analysis using direct probability integral method and exponential convex model. *Reliab Eng Syst Safe* 2022;228:108803.

15. Kundu A, Diazdelao F, Adhikari S, Friswell M. A hybrid spectral and metamodeling approach for the stochastic finite element analysis of structural dynamic systems. *Comput Method Appl Mech Eng* 2014;270:201-19.

16. Wan Z, Chen J, Li J. Probability density evolution analysis of stochastic seismic response of structures with dependent random parameters. *Probabilist Eng Mech* 2020;59:103032.

17. Mouyeaux A, Carvajal C, Bressolette P, Peyras L, Breul P, Bacconnet C. Probabilistic stability analysis of an earth dam by stochastic finite element method based on field data. *Comput Geotech* 2018;101:34-47.

18. Zheng Z, Dai H. Structural stochastic responses determination via a sample-based stochastic finite element method. *Comput Method Appl Mech Eng* 2021;381:113824.

19. Xiu D. Numerical methods for stochastic computations: a spectral method approach Princeton, NJ: Princeton University Press; 2010.

20. Stefanou G. The stochastic finite element method: past, present and future. *Comput Method Appl Mech Eng* 2008;198:1031-51.

21. Matthies HG, Keese A. Galerkin methods for linear and nonlinear elliptic stochastic partial differential equations. *Comput Method Appl Mech Eng* 2005;194:1295-331.

22. Ghanem RG, Spanos PD. Stochastic finite elements: a spectral approach. New York Berlin Heidelberg: Springer; 1991.

23. Gopalakrishnan S, Chakraborty A, Mahapatra DR. Spectral finite element method. London: Springer; 2012.

24. Dai H, Zheng Z, Ma H. An explicit method for simulating non-Gaussian and non-stationary stochastic processes by Karhunen-Loeve and polynomial chaos expansion. *Mech Syst Signal Proc* 2019;115:1-13.

25. Stavroulakis G, Giovanis DG, Papadopoulos V, Papadrakakis M. A GPU domain decomposition solution for spectral stochastic finite element method. *Comput Method Appl Mech Eng* 2017;327:392-410.

26. Farhat C, Wilson E, Powell G. Solution of finite element systems on concurrent processing computers. *Eng Comput* 1987;2:157-65.

27. Farhat C, Lesoinne M, Pierson K. A scalable dual-primal domain decomposition method. *Numer Linear Algebra Appl* 2000;7:687-714.

28. Sarkar A, Benabbou N, Ghanem R. Domain decomposition of stochastic PDEs: theoretical formulations. *Int J Numer Methods Eng* 2009;77:689-701.

29. Cros JM. A preconditioner for the Schur complement domain decomposition method. 2003. Available from: http://www.ddm.org/DD14/cros.pdf[Last accessed on 21 Mar 2024.

30. Desai A, Khalil M, Pettit C, Poirel D, Sarkar A. Scalable domain decomposition solvers for stochastic PDEs in high performance computing. *Comput Method Appl Mech Eng* 2018;335:194-222.

31. Subber W, Sarkar A. A domain decomposition method of stochastic PDEs: An iterative solution techniques using a two-level scalable preconditioner. *J Comput Phys* 2014;257:298-317.

32. Toselli A, Widlund OB. Domain decomposition methods-algorithms and theory. In Springer series in computational mathematics, Berlin: Springer; 2005.

33. Pranesh S, Ghosh D. A FETI-DP based parallel hybrid stochastic finite element method for large stochastic systems. *Comput Struct* 2018;195:64-73.

34. Mathew TPA. Domain decomposition methods for the numerical solution of partial differential equations. Berlin: Springer; 2008.

36. Babuska I, Tempone R, Zouraris GE. Galerkin finite element approximations of stochastic elliptic partial differential equations. *SIAM J Numer Anal* 2004;42:800-25.

37. Ullmann E. A kronecker product preconditioner for stochastic galerkin finite element discretizations. *SIAM J Sci Comput* 2010;32:923-46.

## Cite This Article

Export citation file: **BibTeX** | **RIS**

**OAE Style**

Luo B, Cao W, Zhou Z, Dai H. A domain decomposition solver for spectral stochastic finite element: an approximate sparse expansion approach. *Dis Prev Res* 2024;3:3. http://dx.doi.org/10.20517/dpr.2023.39

**AMA Style**

Luo B, Cao W, Zhou Z, Dai H. A domain decomposition solver for spectral stochastic finite element: an approximate sparse expansion approach. *Disaster Prevention and Resilience*. 2024; 3(1): 3. http://dx.doi.org/10.20517/dpr.2023.39

**Chicago/Turabian Style**

Bowen Luo, Wen Cao, Zheng Zhou, Hongzhe Dai. 2024. "A domain decomposition solver for spectral stochastic finite element: an approximate sparse expansion approach" *Disaster Prevention and Resilience*. 3, no.1: 3. http://dx.doi.org/10.20517/dpr.2023.39

**ACS Style**

Luo, B.; Cao W.; Zhou Z.; Dai H. A domain decomposition solver for spectral stochastic finite element: an approximate sparse expansion approach. *Dis. Prev. Res.* **2024**, *3*, 3. http://dx.doi.org/10.20517/dpr.2023.39

## About This Article

### Copyright

**Open Access**This article is licensed under a Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, sharing, adaptation, distribution and reproduction in any medium or format, for any purpose, even commercially, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

## Data & Comments

### Data

**Views**

**Downloads**

**Citations**

**Comments**

**2**

### Comments

Comments must be written in English. Spam, offensive content, impersonation, and private information will not be permitted. If any comment is reported and identified as inappropriate content by OAE staff, the comment will be removed without notice. If you have any queries or need any help, please contact us at support@oaepublish.com.

**9**clicks

**2**likes

^{0}