This Finite Element Analysis (FEA) stress/strain calculation code uses the finite element method to solve for displacement, stress and strain for two or three dimensional meshes. It allows a range of finite elements types, with differing shape functions and numbers of nodes, as well as a number of analysis modes (e.g. plain strain, axis-symmetric etc.). It uses the popular technique of constructing stiffness matrices for every finite element and incorporating contributions from each of those matrices into a global system matrix. The system matrix is then adjusted to account for boundary conditions with loadings being catered for in a right hand side of the system matrix equation. The solution of displacements to satisfy the system matrix equation is calculated using a conjugate gradient technique. 

Code information: 7000 lines of source and 75 subroutines

Total Parallelization Time using ParaWise : Approximately 2 hours.

User Time: Approximately 20 minutes.

 

Results

The results are for a 3000 element test case . The same code was transported (via ftp) to the different parallel platforms, compiled and linked with CAPLib to create the parallel executable.

bullet

Cray T3D

bullet

IBM SP2

bullet

Transtech Paramid

 

  Cray T3D

Unstructured mesh partition 

(3000 element test case)

Processors

Speed Up
1 -
2 1.78
4 3.39
8 6.55
16 12.34
32 19.47

  IBM SP2

Unstructured mesh partition 

(3000 element test case)

Processors

Speed Up
1 -
2 1.85
4 3.49
8 5.71
16 8.79

  Transtech Paramid

Unstructured mesh partition 

(3000 element test case)

Processors

Speed Up
1 -
2 1.85
4 3.51
8 6.56
16 10.63