Mathematical Problems in Engineering

Volume 2008 (2008), Article ID 181548, 17 pages

http://dx.doi.org/10.1155/2008/181548

## Identification of Stochastic Loads Applied to a Nonlinear Dynamical System Using an Uncertain Computational Model

^{1}Université Paris-Est, Laboratoire Modélisation et Simulation Multi Echelle, MSME FRE3160 CNRS, 5 bd Descartes, 77454 Marne-la-Vallée, France^{2}R&D Department Analyse Mecanique et Acoustique, EDF, 1 avenue du Géneral de Gaulle, 92140 Clamart Cedex, France

Received 6 April 2008; Accepted 4 June 2008

Academic Editor: Jose Balthazar

Copyright © 2008 C. Soize and A. Batou. 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.

#### Abstract

This paper deals with the identification of stochastic loads applied to a nonlinear dynamical system for which a few experimental responses are available using an uncertain computational model. Uncertainties are induced by the use of a simplified computational model to predict the responses of the real system. A nonparametric probabilistic approach of both parameter uncertainties and model uncertainties is implemented in the simplified computational model in order to take into account uncertainties. The level of uncertainties is identified using the maximum likelihood method. The identified stochastic simplified computational model which is obtained is then used to perform the identification of the stochastic loads applied to the real nonlinear dynamical system. A numerical validation of the complete methodology is presented.

#### 1. Introduction

This paper is devoted to the identification of a stochastic load applied to a nonlinear dynamical system for which a few measurements of its responses are available and for which an uncertain simplified computational model is used. In the dynamical system, the uncertainties are taken into account in the context of the probability theory. Consequently, the uncertain simplified computational model is in fact a stochastic simplified computational model for which the input is a stochastic process (stochastic load) and for which the linear operators of the computational model are random. This identification is then performed using the stochastic simplified computational model which allows the responses of the real system to be predicted and, then, the stochastic loads to be identified in minimizing a certain distance between the experimental responses and the random responses calculated with the stochastic simplified computational model. In fact, the methodology presented is developed in the context of the nonlinear dynamical analysis of tube bundles in pressurized water reactor. The stochastic loads applied to the tubes which have to be identified are then induced by the turbulent flow. Since such a nonlinear dynamical system is very complex, the computational model developed cannot exactly represent the complexity of the system. Consequently, the identification is not performed using a computational model which has the capability to accurately predict the experimental responses but is performed using a simplified computational model containing model errors. In order to perform a robust identification of the stochastic loads with respect to model uncertainties in the nonlinear dynamical computational model, a probabilistic model of uncertainties allowing both parameter uncertainties and model uncertaintes to be taken into account is introduced. The responses of the computational model are then random and the randomness is due to the stochastic loads and is due to the stochasticity of the system. In a first step, the probability model of uncertainties in the computational model is identified using the maximum likelihood method. We then deduce a stochastic computational model which allows a robust identification of stochastic loads to be carried out with respect to uncertainties in the nonlinear computational model. The second step is devoted to the stochastic inverse problem consisting in identifying the stochastic loads. From a theoretical and methodological point of view, we then present a complete probabilistic construction and the associated methodology to solve an inverse problem consisting of the identification of a Gaussian stationary stochastic process which is the input of a continuous nonlinear dynamical system with random operators and for which the stochastic output is measured. It should be noted that, if the parametric probabilistic approach is usual to take into account system parameter uncertainties, in the present paper, both the system parameter uncertainties and the model uncertainties are taken into account using a nonparametric probabilistic approach consisting in directly modeling the linear operators of the dynamical system by random operators using the random matrix theory.

Section 2 deals with the construction of the mean computational model. In Section 3, the probabilistic model of the stochastic loads is introduced. Section 4 is devoted to the identification of the stochastic load. The last section presents a numerical validation of the methodology proposed.

#### 2. Mean Computational Model

Let be the domain of the dynamical system having a nonlinear behavior due to the presence of elastic stops located to several points of the part of the boundary of . The domain is decomposed in two bounded open subdomains of : the subdomain and the subdomain . The subdomain is constituted of a three-dimensional linear viscoelastic medium with instantaneous memory and there are elastic stops located at points in the boundary of . In addition, the subsystem occupying the subdomain is fixed on the part of its boundary . The outward unit normal of is denoted by . The subdomain is constituted of a three-dimensional linear viscoelastic medium with instantaneous memory, fixed on the part of its boundary . The outward unit normal of is denoted by . Consequently, each uncoupled subsystem and does not have rigid body displacement. These two subsystems are coupled on the common coupling interface . One then has and . We are interested in constructing the stationary random responses of the nonlinear stochastic dynamical system excited by stationary stochastic processes. Consequently, we will not introduce the initial conditions and we will assume that the time parameter belongs to .

##### 2.1. Mean Boundary Value Problems

###### 2.1.1. Mean Boundary Value Problems for the Linear Subsystem

Let be the cartesian coordinates and be the displacement field for the linear subsystem at time . The external prescribed volumetric and surface forces fields applied to and to its boundary are denoted by and . The stress tensor is written as where is the linearized strain tensor. The fourth-order tensors and verify the usual properties of symmetry and positiveness [1]. Then, the displacement field verifies, for all and for , the mean boundary value problemin which a dot means the partial time derivative and a double dot means the double partial time derivative, where is the forces induced by the subsystem on via the coupling interface . One has used the classical convention for summations over repeated latin indices. The parameter is the mass density for the subsystem .

###### 2.1.2. Mean Boundary Value Problems for the Nonlinear Subsystem

Let be the displacement field for the nonlinear subsystem at time . The external prescribed volumetric and surface forces fields applied to and to its boundary are denoted by and . Since is occupied by a linear viscoelastic material with instantaneous memory, the stress tensor is written as where is the linearized strain tensor. The fourth-order tensors and verify, as above, the usual properties of symmetry and positiveness. Then, the displacement field verifies, for all and for , the mean boundary value problemin which is the forces induced by the subsystem on via the coupling interface . The forces— represent the actions exerted by the elastic stop located at point on the subsystem and is the surface Dirac measure such that, for all continuous function defined on , one has . The parameter is the mass density for the subsystem .

###### 2.1.3. Interface Conditions for the Coupling of with

The coupling conditions on are written as

##### 2.2. Mean Finite Element Model

The mean finite element method [2] is applied to the variational formulation of the boundary value problems defined in Section 2.1.

###### 2.2.1. Mean Finite Element Model for Subsystem

The vector of the DOF of the subsystem is written as , where is the vector of the internal DOF and where is the vector function of the coupling DOF on the interface. From (2.1), it can be deduced that the mean finite element model of subsystem is written asin which , and are, respectively, the positive-definite symmetric real positive mass, damping, and stiffness matrices. The vectors and of the external forces and of the coupling forces are written as and .

###### 2.2.2. Mean Finite Element Model for Subsystem

Similarly to Section 2.2.1, the vector of the DOF of the subsystem is written as , where is the vector of the internal DOF and where is the -valued function of the coupling DOF. From (2.2), it can be deduced that the mean finite element model of subsystem is written asin which , and are, respectively, the positive-definite symmetric real positive mass, damping, and stiffness matrices. The vectors , , and of the external forces, of the coupling forces and of the nonlinear forces, are written as , and .

###### 2.2.3. Interface Conditions for the Coupling of with

The finite element discretization of the interface conditions defined by (2.3) yields

##### 2.3. Reduced Mean Matrix Model

The continuous linear subsystem (linear dynamical subsystem) contains elastic modes in the frequency band of analysis. In addition, the computational model of the continuous linear subsystem is uncertain (presence of both the system parameter uncertainties and the model uncertainties). As we have explained in Section 1, these uncertainties are taken into account using the nonparametric approach of uncertainties which requires a reduced matrix order model (see [3–6]). Since we have to represent the effects of this substructure on the nonlinear substructure through the coupling interface, it is natural to use the Craig Bampton method [7] in order to reduce the finite element model of subsytem . Finally, to reduce the computational cost of the coupled system, subsystem is also reduced with the same technique.

###### 2.3.1. Reduced Mean Matrix Model for Subsystem

The following change of coordinates is introduced:in which is the real matrix whose columns are the first elastic modes for the subsystem with a fixed coupling interface. Those modes are associated with the first eigenvalues such thatwhere and are the internal DOF blocks of the matrices and , where is an matrix, where is the unity matrix and where is an -vector. Let . Then, the vector is a solution of the reduced mean computational modelin which the matrices , and are positive-definite symmetric real matrices.

###### 2.3.2. Reduced Mean Matrix Model for Subsystem

Using the same reduction method and introducing the elastic modes of the linear subsystem with fixed interface and without elastic stops, the -vector verifies the following matrix equation:in which the matrices , and are positive-definite symmetric real matrices.

###### 2.3.3. Transient Dynamical Response of the Reduced Nonlinear Computational Model

Let be the total number of DOF for the nonlinear computational model. The -vector of the mean nonlinear computational model is written asin which the matrix is constructed by the assemblage of and . Let . Then, using the coupling conditions defined by (2.6), the -vector is a solution of the reduced nonlinear dynamical systemwith

#### 3. Stochastic Nonlinear Computational Model Including System Uncertainties and Identification

In this part, firstly the nonparametric probabilistic approach will be used to take into account both data uncertainties and model uncertainties in the reduced mean computational model of the linear subsystem of the computational model. This approach which has recently been introduced consists in replacing the mass, damping, and stiffness matrices of reduced mean computational model by random matrices for which the probability distributions are explicitly given by the theory and for which a generator of independent realizations is known. Such an approach has been validated for many cases. For the details concerning the nonparametric probabilistic approach, one refers the reader, for instance, to [3–6]. In such an approach, the levels of uncertainties for the mass, damping, and stiffness random matrices are defined by the dispersion parameters which are defined below. Secondly, these dispersion parameters will be identified using the maximum likelihood method. Finally, the stochastic nonlinear computational model will be introduced and deduced from Section 2.2.3. It should be noted that (2.1) only the linear subsystem is assumed to be uncertain and (2.2) the mean nonlinear subsystem is representative and consequently, that both data uncertainties and model uncertainties can be neglected. If such an assumption was not verified, then the nonparametric probabilistic approach of uncertainties could always be implemented without any difficulties in this nonlinear subsystem (see for instance [8, 9]).

##### 3.1. Stochastic Linear Subsystem Modeling Uncertainties

Therefore, the matrices , and of the reduced mean computational model are replaced by the random matrices , and defined on a probability space and such thatin which is the set of all the positive-definite symmetric matrices, where is the mathematical expectation and where is the Frobenius norm such that with , is the conjugate of and is the trace for matrices. Let be the random matrix denoting , or . The probability distribution of the random matrix depends on the dispersion parameter related to the coefficient of variation of the random matrix by the equationThe dispersion parameter allows the level of uncertainites of the random matrix to be controlled. It can be found in [3, 5] an algebraic representation of random matrix which allows independent realizations to be explicitly constructed in order to solve the random equations by the Monte Carlo method. For each random matrix, this random generator depends only on the mean value, on the dimension of the matrix, and on the dispersion parameter. Such an approach is used in this paper.

##### 3.2. Identification of the Dispersion Parameters

As explained in Section 3.1, the probability distributions of the random matrices (and then of the random generators) depend on the vector of the dispersion parameters which is identified using the measurements. The observation of the stochastic computational model is defined introducing the random complex dynamic stiffness matrix of the linear subsystem written asThen the random condensed dynamical stiffness matrix of the linear subsystem on the coupling interface is such that . Taking into account the properties of the probabilistic model, it can be shown that, for all fixed in , the random matrix is invertible almost surely and the random variable defined byexists and has a finite mean value. This random variable gives a measure over of the dynamical effects of subsystem on the subsystem at the coupling interface. It should be noted that the random variable depends on because the probability distributions of the random matrices , and depend on . Let be the probabilty density function of the random variable with respect to . For any fixed in and for any value of the vector belonging to the admissible set of the dispersion parameters, the value of the probability density function is estimated by using the above probabilistic model and the Monte Carlo simulation. The corresponding deterministic experimental value of is calculated using experimental data. The method used to identify vector is the maximum likelihood method (see for instance [10]) for the random variable for which is one realization. We then have to solve the following optimization problem:in which is the identified value of .

##### 3.3. Random Transient Dynamical Response of the Stochastic Nonlinear Computational Model

Using the probabilistic model defined in Section 3.1, the deterministic (2.11) to (2.15) give the following stochastic nonlinear computational model:in which, for all fixed , the vector-valued random variable verifiesand where the random matrices , and are written as

#### 4. Identification of Stochastic Loads

The transient load defined by corresponding to the displacement vector is modeled by a stochastic process . Since all the degrees of freedom of the computational model are not excited by this stochastic load, we then introduce the usual projection operator Proj in order to extract the vector of the nonzero random components of the random vector . This equation can easily be inversed and yields .

##### 4.1. Construction of the Stochastic Load

The stochastic load is modeled by an -valued Gaussian stationary centred second-order stochastic process defined on a probability space different from the probability space . In addition, it is assumed that the stochastic process is mean square continuous on , physically realizable (causal) and for which its matrix-valued autocorrelation function is integrable on . This stochastic process is then completely defined by its matrix-valued spectral density function which is a continuous and integrable function on and which is in values in the set of all the positive hermitian matrices. In addition, we will assume that for all in , the matrix is with values in the set of all the positive definite hermitian matrices. Since the stochastic process is assumed to be physically realizable, the matrix valued spectral density function must satisfy the following usual inequality [11, 12]:The numerical simulation of independent realizations for (trajectories) can easily be generated by using adapted algorithms (see, e.g., [13, 14]).

##### 4.2. Stochastic Equation for Simulation of Responses

We have to identify the stochastic process in presence of uncertainties in the linear subsystem . This identification consists in identifying the matrix-valued spectral density function which completely describes the stochastic process. This stochastic inverse problem is formulated as a stochastic optimization problem. Such an identification is performed using the stochastic equation deduced from (3.6) to (3.8) with (2.15) in which the deterministic load is replaced by the stochastic load . We then have to construct the -valued stationary solution (corresponding to ) which is written asin which subindex is relative to the stationary solution and where the stationary stochastic process satisfies the stochastic equationin which and are the mean-square first and second derivative of . For the identification of , for all fixed, we introduce the -valued random variable which represents the observations of the stochastic computational model made up of components of the vector-valued random response . Thus there exists a projection from into such that . For all in , the stationary stochastic process is such that , where the stationary stochastic process is such thatwhere is a stochastic process defined in Section 4.1. Section 4.3 is devoted to the identification of . In order to perform this identification, we need to introduce an observation relative to the stochastic equation and which is useful to construct the cost function. For all , the matrix-valued spectral density function can be estimated. Generating independent realizations of the random matrices , and , the matrix-valued spectral density function is estimated by the Monte Carlo simulation method. For all , one has

##### 4.3. Identification of the Stochastic Loads

The identification is performed in introducing a parametric representation of this function which is rewritten asin which is the admissible set of the parameter with values in where is the number of unknown scalar parameters which have to be identified and where is a given function from into . Therefore, the identification of the stochastic load consists in identifying the -valued vector . Let be the -valued stationary stochastic process which is measured for the manufactured real system and corresponding to the observation stochastic process . The matrix-valued spectral density function of this stochastic process is estimated using the periodogram method. Then, the parameter is estimated in minimizing the distance between the matrix-valued spectral density function calculated with the stochastic computational model and the experimental matrix-valued spectral density function. We then have to solve the following optimization problemin which is the identified value of the vector .

#### 5. Application

In this section, a numerical simulation of a simple example is presented in order to validate the methodology developed in this paper.

##### 5.1. Data for the Experimental Model

The measurements are generated by an experimental model which is made up of one linear subsystem and one nonlinear subsystem. The linear subsystem is made up of four parallel beams fixed at their ends. The nonlinear subsystem is made up of a beam also fixed at its ends, parallel to the other beams and with one transversal symmetric elastic stop (two identical stops, see Figure 1). The five beams are linked by three transversal grids, each grid being modeled by four transversal springs (see Figure 1). Therefore, the coupling interface between the two subsystems is composed of three points located in the neutral fiber of the beam of the nonlinear subsystem. Each beam is modeled by eight Euler beam finite elements of equal lengths for which the DOF of the two nodes at the ends of the beam are locked. The twelve springs defining the three transversal grids are modeled by twelve spring elements. The two elastic stops are modeled by two springs. We are only interested in the -direction displacements of the beam of the nonlinear subsystem (see Figure 1). Consequently, each beam has 14 DOF (-translation and -rotation). The total number of the free DOF for the linear subsystem is then 59 and the total number of the free DOF for the nonlinear subsystem is then 14. The beam of the nonlinear subsystem is exited by seven transversal forces applied following the -direction. The vectors of these seven nonzero components are denoted by . The stochastic process is a second-order centred stationary Gaussian stochastic process for which its matrix-valued spectral density function is such that (2.1) for all in , is a constant on the frequency band of analysis , and (2.2) for all and in , where in which is the distance between the two excited points and the value of is equal to the quarter of the beam length. In the frequency band of analysis , there are 21 eigenfrequencies for the linearized coupled system made up of the linear subsystem coupled with the linear beam of the nonlinear subsystem (nonlinear subsystem without the stops), for which the first three eigenfrequencies are Hz, Hz, and Hz.

##### 5.2. Data for the Mean Computational Model

This part is devoted to the construction of a simplified mean computational model for the nonlinear dynamical system described in Section 5.1. This simplified mean computational model will be used to identify the stochastic loads. It consists in modeling (2.1) the four beams and the three transversal grids of the linear subsystem of the experimental model by a unique equivalent linear Euler beam and by three equivalent springs (see Figure 2) and (2.2) the linear beam with elastic stops of the nonlinear subsystem of the experimental model by a linear beam with two springs for the elastic stops. The section of the equivalent beam for the linear subsystem is arbitrarily chosen and its Young's modulus and its mass density are identified so that the three first eigenfrequencies of this mean computational model are the same that the three first eigenfrequencies of the experimental model. Note that only the three first eigenfrequencies are correctly fitted and consequently, there are model uncertainties in this simplified mean computational model, which are taken into account as explained in Section 3. It should be noted that the objective of this paper is not to construct an accurate mean computational model in order to exactly represent the experimental model, but to test the validity of the use of a simplified mean computational model in order to represent a much more complex system. After identification, the first three eigenfrequencies of the simplified mean computational model are Hz, Hz, and Hz which have to be compared to the first three eigenfrequencies Hz, Hz, and Hz of the experimental model.

##### 5.3. Comparison between the Dynamical Responses of the Experimental Model and of the Mean Computational Model for the Same Given Stochastic Load

In this section, it is assumed that the stochastic load is given and the same for the experimental model and for the simplified mean computational model. Then, for the two models, the stationary stochastic responses are calculated in the time interval s using an explicit Euler integration scheme. Let be the point of the non-linear subsystem located at the impact point of the elastic stops. The power spectral density functions of the stochastic -displacement and of the stochastic -rotation in point (see Figure 3) is estimated using the periodogram method. It can be seen that the prevision given by the mean simplified computational model is good enough in the frequency band . Nevertheless, there are significant differences in the frequency band induced by model uncertainties. This is the reason why the model uncertainties are taken into account in order to extend the domain of validity of the simplified mean computational model in the frequency band in order to perform a robust identification of the stochastic loads.

##### 5.4. System Uncertainties Modeling and Dispersion Parameter Identification

The nonparametric probabilistic approach of model uncertainties introduced in Section 3.1 is used for stiffness part of the linear subsystem of the simplified mean computational model. We then have to identify the dispersion parameter . Note that the identification procedure which is proposed is independent of the stochastic loads. The estimation of the probability density function in (3.5) is carried out with 200 realizations for the Monte Carlo simulation. Figure 4 shows the likelihood function calculated using (3.5) with . The maximum is reached for .

##### 5.5. Case of an Unknown Stochastic Load and Its Identification

In this section, the responses of the experimental model are given (those constructed in Section 5.3) and the stochastic load is assumed to be unknown and has to be identified using the uncertain simplified computational model, that is to say the stochastic simplified computational model for which the dispersion parameter has been identified in Section 5.4. We begin defining a model as simple as possible for the stochastic load introduced in Section 4.1. We have then chosen to model as in which is a -force and is a -moment applied to the middle of the beam of the nonlinear subsystem (see Figure 5). This force and this moment are independent second-order centred stationary Gaussian stochastic processes. So, they are both completely defined by their power spectral density functions and . The matrix-valued spectral density function of the stochastic process is then defined byIt is assumed that the function is constant in the frequency band of analysis and is such that (4.1) is verified. The experimental stochastic process defined in Section 4.3 is composed of stochastic -displacements. Taking into account (4.6), the function which is a constant diagonal hermitian matrix, can then be rewritten for all in asin which the admissible set . This vector is identified using the trial method to solve the optimization problem defined by (4.7). Such a method consists in calculating the cost function for 100 values of the vector . Figure 6 shows the graph of the function which allows the optimal value to be determined. The confidence region associated with a probability level of the reponse of the stochastic simplified computational model on which the identified stochatic load is applied can then be estimated. The comparison between the experimental responses with the responses constructed with the stochastic simplified computational model is given in Figure 7. This figure displays the confidence region of the power spectral density function of the stochastic -displacement and the stochastic -rotation for point .

#### 6. Conclusions

We have presented a methodology and its validation to perform the identification of a stochastic loads applied to a complex nonlinear dynamical system for which a few measurements of its responses are available. To carry out this identification, a simplified computational model of the real system is introduced. Since such a simplified computational model induces model uncertainties, a probabilistic model of these uncertainties is introduced in the simplified computational model. The identification of the stochastic loads is then performed using this stochastic computational model which takes into account model uncertainties and consequently, we have validated a method to perform a robust identification with respect to model uncertainties. It should be noted that the nonlinear dynamical system used for this validation is representative of real industrial systems and then validates the methodology proposed.

#### References

- R. Ohayon and C. Soize,
*Structural Acoustics and Vibration*, Academic Press, San Diego, Calif, USA, 1998. View at Zentralblatt MATH · View at MathSciNet - O. C. Zienkiewicz and R. L. Taylor,
*The Finite Element Method*, McGraw-Hill, New York, NY, USA, 4th edition, 1989. View at Zentralblatt MATH · View at MathSciNet - C. Soize, “A nonparametric model of random uncertainties for reduced matrix models in structural dynamics,”
*Probabilistic Engineering Mechanics*, vol. 15, no. 3, pp. 277–294, 2000. View at Publisher · View at Google Scholar - C. Soize, “Maximum entropy approach for modeling random uncertainties in transient elastodynamics,”
*The Journal of the Acoustical Society of America*, vol. 109, no. 5, pp. 1979–1996, 2001. View at Publisher · View at Google Scholar - C. Soize, “Random matrix theory for modeling uncertainties in computational mechanics,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 194, no. 12–16, pp. 1333–1366, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. Soize and H. Chebli, “Random uncertainties model in dynamic substructuring using a nonparametric probabilistic model,”
*Journal of Engineering Mechanics*, vol. 129, no. 4, pp. 449–457, 2003. View at Publisher · View at Google Scholar - R. R. Craig, Jr. and M. C. C. Bampton, “Coupling of substructures for dynamic analysis,”
*AIAA Journal*, vol. 6, no. 7, pp. 1313–1319, 1968. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. Desceliers, C. Soize, and S. Cambier, “Nonparametric parametric model for random uncertainties in nonlinear structural dynamics: application to earthquake engineering,”
*Earthquake Engineering and Structural Dynamics*, vol. 33, no. 3, pp. 315–327, 2004. View at Publisher · View at Google Scholar - R. Sampaio and C. Soize, “On measures of nonlinearity effects for uncertain dynamical systems-application to a vibro-impact system,”
*Journal of Sound and Vibration*, vol. 303, no. 3–5, pp. 659–674, 2007. View at Publisher · View at Google Scholar - R. J. Serfling,
*Approximation Theorems of Mathematical Statistics*, Wiley Series in Probability and Mathematical Statistic, John Wiley & Sons, New York, NY, USA, 1980. View at Zentralblatt MATH · View at MathSciNet - Yu. A. Rozanov,
*Stationary Random Processes*, Holden-Day, San Francisco, Calif, USA, 1967. View at Zentralblatt MATH · View at MathSciNet - C. Soize,
*The Fokker-Planck Equation for Stochastic Dynamical Systems and Its Explicit Steady State Solutions*, vol. 17 of*Series on Advances in Mathematics for Applied Sciences*, World Scientific, River Edge, NJ, USA, 1994. View at Zentralblatt MATH · View at MathSciNet - M. Shinozuka, “Simulation of multivariate and multidimensional random processes,”
*The Journal of the Acoustical Society of America*, vol. 49, no. 1, pp. 357–367, 1971. View at Publisher · View at Google Scholar - F. Poirion and C. Soize, “Numerical methods and mathematical aspects for simulation of homogeneous and non homogeneous Gaussian vector fields,” in
*Probabilistic Methods in Applied Physics*, pp. 17–53, Springer, Berlin, Germany, 1995. View at Publisher · View at Google Scholar