In this paper, based on the work of L.J.Hayes in 1986, a vectorized algorithm of the global stiffness matrix-vector multiply in finite element structural analysis is developed when we in fact don't form the global stiffness matrix. A numerical example by using the method in iterative solution of finite element equations verify it is very efficient.

This paper presents simiimplicit Runge-kutta type parallel direct integration methods to solve structural dynamic problems using 3-stage-3-order semi-implicit Runge-kutta method and using polynomial preconditioned conjugate gradient method for solving relevant systems of linear algebraic equations. Compared with relevant serical algorithm RK33S on YH-1 computer, parallel algorithm RK33P's speed up is 24~27 when linear systems of equations involved is of order 10³~10⁴.