FE-models for structural solid mechanics analyses can be readily generated from computer images via a 'voxel convesion' method, whereby voxels in a two- or three-dimesional computer image are directly translated to elements in a FE-model. The fact that all elements thus generated are the same creates the possibilities for fast solution algorithm that can compensate for a large number of element. The solving methods described in this paper are based on an iterative solving algorithm in combination with a uniqueelement Element-by-Element (EBE) or with a newly developed Row-by-Row (RBR) matrix-vector multiplication strategy. With these methods it is possible to solve FE-models on the order of 105 3-D brick elements on a workstation and on the order of 106 elements on a Cray computer. The methods are demonstrated for the Boussinesq problem and for FF models that represent a porous trabecular bone structure The results show that the RBR method can be 3.2 times faster than the EBE method. It was concluded that the voxel conversion method in combination with these solving methods not only provides a powerful tool to analyse structures that can not be analysed in another way, but also that this approach can be competitive with traditional meshing and solving techniques.

, , , , , , , , , , ,
International Journal for Numerical Methods in Engineering
Erasmus MC: University Medical Center Rotterdam

van Rietbergen, B., Weinans, H., Huiskes, R., & Polman, B. (1996). Computational strategies for iterative solutions of large fem applications employing voxel data. International Journal for Numerical Methods in Engineering, 2743–2767. Retrieved from http://hdl.handle.net/1765/15304