`International Journal of Stochastic AnalysisVolume 2012 (2012), Article ID 249201, 14 pageshttp://dx.doi.org/10.1155/2012/249201`
Research Article

## Application of Stochastic Sensitivity Analysis to Integrated Force Method

1Mechanical Engineering Department, The University of Akron, OH 44325, USA
2Ohio Aerospace Institute, Cleveland, OH 44135, USA

Received 18 February 2012; Accepted 18 April 2012

Copyright © 2012 X. F. Wei and S. N. Patnaik. 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

As a new formulation in structural analysis, Integrated Force Method has been successfully applied to many structures for civil, mechanical, and aerospace engineering due to the accurate estimate of forces in computation. Right now, it is being further extended to the probabilistic domain. For the assessment of uncertainty effect in system optimization and identification, the probabilistic sensitivity analysis of IFM was further investigated in this study. A set of stochastic sensitivity analysis formulation of Integrated Force Method was developed using the perturbation method. Numerical examples are presented to illustrate its application. Its efficiency and accuracy were also substantiated with direct Monte Carlo simulations and the reliability-based sensitivity method. The numerical algorithm was shown to be readily adaptable to the existing program since the models of stochastic finite element and stochastic design sensitivity are almost identical.

#### 1. Introduction

As an alternative of the classical stiffness method, the force method [13] was also popular in structural analysis of civil, mechanical, and aerospace engineering because of its accurate estimates of forces in structural analysis. A new formulation in the force method, termed the Integrated Force Method [47] (IFM), was proposed by Patnaik for the analysis of discrete and continuous systems. IFM is a force method [8], which integrates both the system equilibrium equations and the global compatibility conditions together. This method has been successfully applied to many structures for their deterministic models with well-defined parameters. Realizing its potential in structural analysis, it is being further extended to consider investigations on the probabilistic analysis. Meanwhile, the probabilistic sensitivity formulation for IFM is also desirable to be investigated for assessing uncertainty effect in system optimization and identification.

In general, sensitivity analysis of structural systems involves the computation of the partial derivatives of some response function with respect to the design parameters. The sensitivity analysis of structural systems to variations in their parameters is one of the ways to evaluate the performance of structures. It is very important for system optimization, parameter identification, reliability assessment, and so forth in engineering analysis. However, the conventional sensitivity analysis is based on the assumptions of complete determinacy of structural parameters. In reality, the occurrence of uncertainty associated with the system parameters is inevitable and intrinsic. Hence, there is a necessity to estimate the effect of uncertainty in stochastic sensitivity derivatives with respect to random design variables on structural responses.

There has been a great interest in developing various methods for computing response sensitivity of structure [911]. Hien and Kleiber [12] formulated the stochastic design sensitivity problems of structural static in an effective way by using the perturbation approach and adjoint variables method. Numerical algorithms are developed and turn out to be readily adaptable to existing finite element codes. The structural design sensitivity and stochastic finite elements are similar in terms of the methodology and computer implementation, which greatly facilitate the combined analysis. Based on the first-order perturbation method, Ghosh et al. [13] illustrated that the stochastic structural response sensitivity is quite satisfactory compared to Monte Carlo simulation results under small variation of input random parameters.

It was noted that literatures about the stochastic sensitivity analysis are limited and almost all of them focused on the common stiffness method. Therefore, this paper was devoted to develop a set of formulations of the probabilistic sensitivity analysis for Integrated Force Method using the perturbation method. Under the assumption of the spatial homogonous random fields, two approximate formulas of the stochastic sensitivity analysis were implemented for IFM. Another stochastic sensitivity analysis from the reliability-based method was also evaluated for comparisons according to the inverse cumulative distribution function of design variables. Numerical examples are presented in this paper to illustrate the application and differences of these methods. The efficiency and accuracy of these analytical formulas were also compared with results from the direct Monte Carlo simulations and the reliability-based sensitivity method. The numerical algorithm was shown to be readily adaptable to the existing program since the models of stochastic finite element and stochastic design sensitivity are almost identical.

#### 2. Integrated Force Method and Deterministic Sensitivity Formulas

The governing equation of IFM [4] for a continuum is discretized by a finite element model with force and displacement unknowns. The equilibrium equations and the ) compatibility conditions are coupled to obtain the IFM governing equation in static analysis as where is the equilibrium matrix, is the compatibility matrix, is the concatenated flexibility matrix, and is the components load vector. is the components vector of combination of and . The r-component initial deformation vector is . Here, is the n-component initial deformation vector and is the governing matrix.

The solution of (2.1) yields the forces . The displacements are obtained from the force by back-substitution where is the deterministic deformation coefficient matrix defined as rows of . It must be stressed here that (2.1) and (2.2) in IFM represent two key relationships in structural analysis.

Based on the above governing equation (2.1), the differentiating equation (2.1) with respect to the “th” design variable results in

that is, where is the internal force sensitivity vector with respect to , that is, . is the load derivative vector with respect to , that is, . is the governing derivative matrix with respect to , that is, . k is the number of design variables, .

Similarly, the deterministic sensitivity formula of displacement with respect to the “kth” design variable in IFM can be expressed by where is the displacement sensitivity vector with respect to , that is, . is the initial deformation derivative vector with respect to , that is, . is the deformation coefficient derivative matrix with respect to , that is, . Note that the deterministic sensitivity formulas are more complicated than the deterministic analysis ones.

#### 3. Stochastic Sensitivity Analysis of IFM

The stochastic structural parameters are assumed to vary spatially as a homogeneous stochastic field. For example, the Young’s modulus corresponding to the th element with a mean value may be written as where is the fluctuating component of the Young’s modulus, called the normalized primitive random variable of the Young’s modulus. Note that the defined random Young’s modulus has zero mean and standard deviation . Similarly, other random structural parameters can also be defined by the normalized primitive random variables.

It is well known that the perturbation technique had been applied to the stochastic analysis of finite element method for structures. The uncertain parameters in structures can cause the response to be uncertain as well as make the response sensitivity uncertain. Assuming the random variables have small fluctuations, a second-order Taylor’s expansion in relation to the random variables is employed for stochastic response sensitivity calculation in IFM. For example, the Taylor series expansion of the internal force sensitivity vector has the following form: where denote the corresponding deterministic part of the vector at and and are partial derivatives of defined as follows:

Similarly, the Taylor series expansions of , , , and in the deterministic sensitivity (2.4) can be also obtained. Substituting these expressions into (2.4), retaining terms up to the second-order terms, and by equating terms of equal orders, the zeroth-, first- and second-order equations on both sides corresponding to (2.4) can be written by the following set of recursive equations:

Therefore, the first-order approximation for the internal force sensitivity is obtained by truncating the right-hand side of (3.2) after the second term as and the second-order approximation where and are the mean values of the first- and second-order approximate formulations of the internal force sensitivity vector with respect to , respectively. and are the covariance matrices of the first- and second-order approximate formulations of the internal force sensitivity matrices with respect to , respectively. In the same way, the displacement sensitivity formulas of IFM can be also obtained by.

The recursive equations:

The first-order approximation:

The second-order approximation: where and are the mean values of the first- and second-order approximate formulations of the displacement sensitivity vector with respect to , respectively. and are the covariance matrices of the first- and second-order approximate formulations of the displacement sensitivity matrices with respect to , respectively.

It is noted that the stochastic response sensitivity formulas of IFM have similar features of the stochastic response formulas [14]. The stochastic response sensitivity analysis formulas consist of the deterministic and random parts. The zeroth-order sensitivity equations are identical to the deterministic sensitivity equations evaluated at the mean value. The first- and second-order sensitivity equations consist of the recursive items. Moreover, covariance matrices in the first- and second-order approximate formulas have the same forms in expressions with the assumption of and . However, to obtain the acceptable accuracy, the variance of uncertainties must be small enough because only finite terms in the above series were used. In this study, the algorithms of stochastic sensitivity analysis for IFM had been programmed in a closed form to evaluate the mean values and covariance matrices using the symbolic software tools, Maple V.

#### 4. Reliability Sensitivity Analysis of IFM

Since the first two probabilistic moments of responses have been obtained in the stochastic response formulation of IFM, the response value at probability of occurrence can be calculated by the inverse of the cumulative distribution function. By using the reliability-based method, stochastic response sensitivity formulation may be expressed by another form.

Let , and denote the response value at probability of occurrence, the mean value, and standard deviation, respectively. By means of the transformation in the probability theory, the inverse of the cumulative distribution function can be written as where is the cumulative distribution function for the standard normal distribution.

The sensitivity of with respect to the mean of a primitive random variable, , can be given by

It should be noted that the response sensitivity is the sum of the derivatives of the mean response and the standard deviation that is prorated by the inverse standard normal cumulative distribution function . Obviously, the sensitivity formula is a complicated expression, even for the simple structure.

#### 5. Illustration Examples

To clearly demonstrate the methodology application, an example, the six-bar truss with mechanical and thermal loads, is evaluated by the following two cases.(1)Only one random variable, Young’s modulus , is considered in the structure; other variables take their mean values as deterministic variables in the computation.(2)Multirandom variables are considered in the computation.

The six-bar truss shown in Figure 1 has a Young’s modulus with a mean of 68.947 GPa and standard deviation of 3.776 GPa and a coefficient of thermal expansion α with a mean of /°C and standard deviation of /°C. It is subjected to a mechanical load with a mean of 4.448 kN and standard deviation of 0.222 kN at node 1 along the y-direction. The temperature on member 3 has an increasing mean of 37.778°C and standard deviation of 13.299°C. The sizing variables, six cross-sectional areas, have the mean, standard deviation, and correlation coefficient matrix as follows:

Figure 1: Six-bar truss.

It should be noted that the problem has ten random variables defined as follows: where , is the normalized primitive random variable with zero mean and standard deviation and the subscript represents any of random variables. Note that the sixth bar was omitted in the results since it is not subjected to any force.

In Case (1), as it is observed from Figures 2 and 3, after 18,000 simulations the fluctuations of mean value and standard deviation of response sensitivities (here, we just calculated force F1 in the first bar and displacement x1 in the horizontal direction at node 1) with respect to the first bar area A1 are small enough and gradually trend to be smooth with increasing simulations. In fact, the number in simulations has been taken as 25,000.

Figure 2: Fluctuation of mean and standard deviation of force sensitivity in the first bar with respect to A1.
Figure 3: Fluctuation of mean and standard deviation of displacement sensitivity, x1, with respect to A1.

The comparisons of the mean value and standard deviation of response sensitivities with varying c.o.v. of between the perturbation methods and direct Monte Carlo Simulation (DMCS) of IFM are shown in Figures 4 and 5, respectively. It is noted in Figure 4 that in the range of small variance of the Young’s modulus the mean force sensitivity curve of DMCS increases linearly, whereas the first- and second-order perturbation curves overlap each other and remain constant. And the standard deviation force sensitivity curves of three methods are close to each other and linearly increase. In Figure 5, the mean displacement sensitivity curves of both DMCS and the second-order perturbation have an accelerated rate of increase and are close to each other prior to c.o.v. of , whereas the first-order perturbation curve remains constant. The standard deviation displacement sensitivity curve of both the first- and second-order perturbation linearly increases with the coefficient of variation and is close to the DMCS curve before c.o.v. of .

Figure 4: Comparison of mean and standard deviation of force sensitivity in the first bar with varying c.o.v. of uncertain variable E.
Figure 5: Comparison of mean and standard deviation of displacement sensitivity, x1, with varying c.o.v. of uncertain variable E.

It is important to note that the perturbation methods underestimate the response variability after the coefficient of variation E = 0.1. It is consistent with the earlier results as demonstrated in some papers [15, 16].

For easy comparisons, assuming the probability of occurrence as 75%, the force sensitivity in the first bar and displacement sensitivity in the horizontal direction at 1 node with respect to bar areas are shown in Figures 6 and 7, respectively. In Figure 6, the force sensitivity values obtained from the reliability method, the perturbation method, and direct Monte Carlo simulation are almost equal to each other. The first bar force is more sensitive to the second bar area A2, not sensitive to bar areas A3 and A5. The deterministic force sensitivity values are different from results of other methods. In Figure 7, the horizontal displacement sensitivity values from both the perturbation method and Monte Carlo simulation are close to each other, but the reliability method and deterministic sensitivity have different values with results of other methods. The horizontal displacement is very sensitive to the first bar area A1, secondly to A2, and insensitive to other bar areas.

Figure 6: Comparison of force sensitivity in the first bar with one random variable at 75% probability of occurrence.
Figure 7: Comparison of displacement sensitivity at the bar with one random variable at 75% probability of occurrence.

In Case (2), also assuming the probability of occurrence as 75%, the comparisons of the force sensitivity in the first bar and displacement sensitivity in the horizontal direction at 1 node with respect to bar areas are shown in Figures 8 and 9, respectively. In Figure 8, note that there are small differences between the perturbation methods and direct Monte Carlo simulation; moreover, the second-order approximation perturbation method has better results with Monte Carlo simulation. However, the reliability method and deterministic sensitivity are not correlated to DMCS. In Figure 9, the reliability method and deterministic sensitivity have big differences with DMCS, whereas the perturbation methods have a good agreement with DMCS. The first bar force and horizontal displacement in Case (2) have the same sensitive and insensitive design variables in Case (1).

Figure 8: Comparison of force sensitivity in the first bar with Multirandom variables at 75% probability of occurrence.
Figure 9: Comparison of displacement sensitivity at the bar with Multirandom variables at 75% probability of occurrence.

The comparisons of the first bar force and horizontal displacement sensitivities with respect to bar areas between the perturbation methods and direct Monte Carlo simulation are shown in Table 1. The number of simulations in DMCS was taken as 20,000. It was noted that the second-order approximation has a better agreement with DMCS than the first-order approximation on the mean values and standard deviations of response sensitivities.

Table 1: Comparison of stochastic sensitivities between PM and MCS.

#### 6. Conclusions

Using the perturbation method, stochastic sensitivity analysis for Integrated Force Method (IFM) had been formulated with the first- and second-order approximations in this study. Furthermore, the developed algorithm can be easily adopted to fit into the existing deterministic sensitivity analysis program of IFM, only incorporating the random variable part subroutine. Without question, the stochastic analysis should enhance IFM available in the probability field for advanced structure analysis. It was observed that the sensitivity analytical results obtained through the second-order approximation have better accuracy than that through the first-order one compared to Monte Carlo results for either single- or Multirandom variables if the coefficients of variation of input random parameter are small. The reliability-based sensitivity analysis shows a big difference with Monte Carlo results and has a computational complexity, even for the simplest structures.

However, for larger coefficients of variation of random parameters, the perturbation method with the first- and second-order approximations is not effective and the accuracy of their calculations also decreases due to the limitation of truncating error and neglect of the higher-order covariance matrices in the approximation formulas. Thus, the further development of the methodology for providing more accurate calculations in stochastic sensitivity analysis of IFM would be desirable for larger coefficient of variation of randomness.

#### Acknowledgments

The support of Glenn Research Center and help from Dr. James D. Guptill and Mrs. Rula M. Coroneos for this research are gratefully acknowledged.

#### References

1. A. Kaveh, “An efficient flexibility analysis of structures,” Computers & Structures, vol. 22, no. 6, pp. 973–977, 1986.
2. A. Kaveh, “Recent developments in the force method of structural analysis,” Applied Mechanics Reviews, vol. 45, no. 9, pp. 401–418, 1992.
3. A. Kaveh, Structural Mechanics: Graph and Matrix Methods, vol. 6 of Applied and Engineering Mathematics Series, Research Studies Press, Taunton, UK, 1st edition, 1992.
4. S. N. Patnaik and D. A. Hopkins, Strength of Materials: A Unified Theory for the 21st Century, Butterworth-Heinemann, 2004.
5. S. N. Patnaik and D. A. Hopkins, “Integrated force method solution to indeterminate structural mechanics problems,” NASA/TP 1999-207430, 1999.
6. S. N. Patnaik, D. A. Hopkins, and R. M. Coroneos, “Recent advances in the method of forces: integrated force method of structural analysis,” Advances in Engineering Software, vol. 29, no. 3–6, pp. 463–474, 1998.
7. S. N. Patnaik, L. Berke, and R. H. Gallagher, “Compatibility conditions of structural mechanics for finite element analysis,” AIAA Journal, vol. 29, no. 5, pp. 820–829, 1991.
8. S. N. Patnaik, L. Berke, and R. H. Gallagher, “Integrated force method versus displacement method for finite element analysis,” Computers & Structures, vol. 38, no. 4, pp. 377–407, 1991.
9. M. Kleiber, H. Antunz, T. D. Hien, and P. Kowalczyk, Parameter Sensitivity in Nonlinear Mechanics, John Wiley & Sons, New York, NY, USA, 1997.
10. T. D. Hien and M. Kleiber, “Stochastic structural design sensitivity of static response,” Computers & Structures, vol. 38, no. 5-6, pp. 659–667, 1991.
11. P. Kolakowski and J. Holnicki-Szulc, “Sensitivity analysis of truss structures (virtual distortion method approach),” International Journal for Numerical Methods in Engineering, vol. 43, no. 6, pp. 1085–1108, 1998.
12. T. D. Hien and M. Kleiber, “Stochastic design sensitivity in structural dynamics,” International Journal for Numerical Methods in Engineering, vol. 32, no. 6, pp. 1247–1265, 1991.
13. R. Ghosh, S. Chakraborty, and B. Bhattacharyya, “Stochastic sensitivity analysis of structures using first-order perturbation,” Meccanica, vol. 36, no. 3, pp. 291–296, 2001.
14. X. F. Wei, S. N. Patnaik, S. S. Pai, and P. P. Ling, “Extension of integrated force method into stochastic domain,” International Journal of Computational Methods in Engineering Science and Mechanics, vol. 10, no. 3, pp. 197–208, 2009.
15. F. Yamazaki, M. Shinozuka, and G. Dasgupta, “Neumann expansion for stochastic finite element analysis,” Journal of Engineering Mechanics, vol. 114, no. 8, pp. 1335–1354, 1988.
16. W. K. Liu, T. Belytschko, and A. Mani, “Computational method for the determination of the probabilistic distribution of the dynamic response of structures,” in Proceedings of the Pressure Vessels and Piping Conference. Piping, Feedwater Heater Operation, and Pumps (ASME PVP '85), vol. 98, no. 5, pp. 243–248, New Orleans, La, USA, 1985.