Department of Computer Science University of North Carolina at Chapel Hill

Large-Scale Heterogeneous Elastic Body Simulation

Objective

Static analysis of the biological tissue deformation based on 3D medical images.

Approaches

Multi-Resolution (Multi-Grid) Finite Element Method
Parallel Computation on Shared Memory Architecture Machines
Dynamic Switching between Incremental Loading and Arc-Length Continuation Methods
Augmented Lagrangian Method for Satisfying Constraints such as Volume Preservation (incompressibility) and Intersection Avoidance
Fast and Versatile Numerical Integration by Adaptive Placements of Quadrature Points Based on Underlying Voxel Data

Preliminary Implementation and Result

A preliminary version of the algorithm was implemented on a SGI Onyx2. The program is written in C++. The C code that computes the derivatives of the elastic energy is automatically generated by a Maple program. Therefore, various energy function can be tested without rewriting the program. The example below is a simulation result of 10x10x10 finite element mesh pushed from the bottom. The program was successfully applied to a 40x40x40 mesh.

MPEG ANIMATION (532 KB) Quick Time

Current Project Members:

For more information, contact Gentaro Hirota <hirota@cs.unc.edu> or Ming C. Lin <lin@cs.unc.edu>
Last Content Update: May 10, 1999
Last Content Review: Feb. 11, 2000

Copyright 1999. Personal use of this material is permitted. However, permissi on to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or redistribution to servers or li sts, or to reuse any copyrighted component of this work in other works must be o btained from the author.

This material is presented to ensure timely dissemination of scholarly and te chnical work. Copyright and all rights therein are retained by authors or by oth er copyright holders. All persons copying this information are expected to adher e to the terms and constraints invoked by each author's copyright. In most cases , these works may not be reposted without the explicit permission of the copyrig ht holder.