Computational & Technology Resources
an online resource for computational,
engineering & technology publications
Civil-Comp Proceedings
ISSN 1759-3433
CCP: 88
PROCEEDINGS OF THE NINTH INTERNATIONAL CONFERENCE ON COMPUTATIONAL STRUCTURES TECHNOLOGY
Edited by: B.H.V. Topping and M. Papadrakakis
Paper 55

Sensitivity Analysis of Computer Models of Structures with Green's Functions

F. Hartmann and T. Kunow

Civil Engineering, University of Kassel, Germany

Full Bibliographic Reference for this paper
F. Hartmann, T. Kunow, "Sensitivity Analysis of Computer Models of Structures with Green's Functions", in B.H.V. Topping, M. Papadrakakis, (Editors), "Proceedings of the Ninth International Conference on Computational Structures Technology", Civil-Comp Press, Stirlingshire, UK, Paper 55, 2008. doi:10.4203/ccp.88.55
Keywords: Green's functions, functionals, finite elements, discretization error, modeling error, sensitivity.

Summary
A structure consists of infinitely many mechanisms - all bolted and fixed so that the structure can carry the load. But if we release one mechanism, apply a unit rotation or dislocation, then this will induce a movement in the structure and the work done by the load on acting through this movement is equal to M(x) . 1 or V(x) . 1, is equal to the Dirac energy [1]. In finite element (FE) analysis we hinder the movements of the structure and so the mechanism obtains an incorrect measure of how large the Dirac energy really is, and consequently Mh(x) not = M(x) and Vh(x) not = V(x). That is the kinematics of a mesh determines the accuracy of an FE solution
  • mesh = kinematics = accuracy of influence functions = quality of results.
This is motivation to reduce the error G-Gh in the Green's functions. A second source of errors are the model parameters which too have an influence on the shape of the Green's function. So numerical accuracy (verification) and model accuracy (validation) hinge on an accurate approximation of the underlying Green's function.

A change in the bending stiffness EI in one frame element alone implies a global change in the Green's function and normally requires a full new evaluation cycle by integrating over all frame elements of the structure. We propose an alternative formulation where integration is done only over the frame element [x1,x2] that changes to assess the change in u(x) or sigma(x) at an arbitrary point x of the structure.

Normally the Green's functions decay with the distance r = |y-x| between the source point y and the field point x but each Green's function comes with its own domain of influence, its own mode of decay and these patterns - which are specific for each Green's function - encapsulate the design sensitivities of a structure. We suggest that questions such as: 'How accurate is the model and how sensitive is the output to changes in the parameters or how sensitive are the regions with high stress gradients to the model parameters?' should be answered by studying the Green's functions. The plots or charts of these functions are an ideal tool to weight the consequences which modifications of model parameters will have on the stress distribution in a structure.

We provide some numerical examples (bridge structures and plates) to demonstrate the efficiency and accuracy of the proposed method.

References
1
F. Hartmann, C. Katz, "Structural Analysis with Finite Elements", Springer-Verlag Berlin Heidelberg New York, 2nd Ed., 2007. doi:10.1007/s11012-005-6002-5

purchase the full-text of this paper (price £20)

go to the previous paper
go to the next paper
return to the table of contents
return to the book description
purchase this book (price £145 +P&P)