On reduced matrix inversion for operator splitting methods

The element-splitting heat-conduction algorithm of Hughes et al. (1983) is modified using the dual-basis approach of Flanagan and Belytschko (1981). The size of the system to be solved is reduced, and the formulation of a stiffness matrix is avoided. The 2D application is demonstrated, and the feasibility of a 3D implementation is indicated.