Computational & Technology Resources
an online resource for computational,
engineering & technology publications |
|
Civil-Comp Proceedings
ISSN 1759-3433 CCP: 80
PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON ENGINEERING COMPUTATIONAL TECHNOLOGY Edited by: B.H.V. Topping and C.A. Mota Soares
Paper 68
Line Elements for the Multiphysics Modeling of Axispherical Bodies D.W. Wood, Z.C. Lai, C.S. Long, S. Kok and A.A. Groenwold
Department of Mechanical and Aeronautical Engineering, University of Pretoria, South Africa D.W. Wood, Z.C. Lai, C.S. Long, S. Kok, A.A. Groenwold, "Line Elements for the Multiphysics Modeling of Axispherical Bodies", in B.H.V. Topping, C.A. Mota Soares, (Editors), "Proceedings of the Fourth International Conference on Engineering Computational Technology", Civil-Comp Press, Stirlingshire, UK, Paper 68, 2004. doi:10.4203/ccp.80.68
Keywords: spherical geometry, MEMS, elastostatics, plasticity, spherical harmonics.
Summary
The modeling of axispherical structures is normally carried out using axicylindrical finite element
meshes. When unsymmetrical loads are encountered, the modeling of
spherical bodies becomes far from trivial. It can be accomplished using axicylindrical formulations
that provide for unsymmetrical loading conditions via a Fourier series description of the load.
(Hence, even though the geometry is axicylindrical, all of the displacements of the 3-D continuum are
provided for.) However, for a sphere, axicylindrical meshes are inevitably prone to element distortion,
while the distance to the axis of symmetry varies dramatically for different elements in the mesh.
As a result, convergence rates are low.
Spherically symmetric loads applied to an isotropic sphere result in spherically symmetric displacement fields within the sphere. Since only one non-zero displacement interpolation function, namely , is required to model this condition, a one-dimensional element can be used. In our paper, we present efficient axispherical finite element formulations for the numerical modeling of spherically symmetric structures. The formulation results in an efficiently banded global stiffness matrix . The solution of the equilibrium equations is therefore far less computationally expensive than for axicylindrical or solid models. Formulations are presented, which provide for elastostatic deformation [1], piezoelectric effects [2,3] and material non-linearity in the form of elastic-perfect-plastic behavior [4]. Provision is made for unsymmetrical loading conditions via a suitable series description of the applied loads using spherical harmonics [5]. For the sake of brevity, we focus on the coupling between mechanical, electrostatic and piezoelectric effects during numerical evaluation of the formulations. The advantages and ease of formulation of the proposed formulations as compared to axicylindrical models are emphasized. Said advantages include:
One possible application of our novel element is the modeling of spherical micro-electro-mechanical systems (MEMS). In recent times, MEMS have received significant attention. They are lightweight, ultra small devices, which are widely used for integrated sensor applications, e.g. acceleration sensors for airbags, active actuators for micro air vehicles (MAV's), etc. They have potential applications in micro-surgery, health-monitoring, micro-fabrication, etc. Spherical MEMS devices are perceived to be important in a number of applications. They are unidirectional, an attractive feature when the to-be-sensed quantity is of unknown location in an effectively infinite field. Spherical MEMS are, for example, used on a very large scale in active and passive sonar arrays. While many an axispherical problem can of course be solved analytically, this becomes very difficult, and even infeasible, when multiphysics phenomena and, for example, plasticity effects are considered. The accommodation of asymmetric loads via spherical harmonics further increase the applicability of the one-dimensional elements. References
purchase the full-text of this paper (price £20)
go to the previous paper |
|