Generalization of rectangular element stiffness matrix and thermal load vector associated with a(0) + a(1)x + a(2)y + a(3)xy type interpolation rule

The stiffness-matrix formulation for the rectangular finite element described by Melosh (1963) and Weaver and Johnston (1984) is generalized to orthotropic materials with material axes not necessarily coincident with the x and y axes; i.e., the condition d(13) = d(23) = 0 is removed. Also included are explicit expressions for the element load vector associated with nonuniform temperature increase in the element. Applications to the analysis of thermal stresses in thin Si-crystal ribbons subjected to temperature changes with highly nonuniform lengthwise and transverse gradients (Utku et al., 1986) and to the simulation of the thermoviscoelastic behavior of growing Si ribbons (Utku and Ray, 1986) are indicated.