LFEMI debug for Out-Of-Plane strain mechanical load
Hand Calculation Verification

WORKING programs:

1. lfemi_oop.dat
2. lfemi_oop.f
3. lfemi_oop.out

IDentification matrix; ID(I,J) where I = Degree Of Freedom (DOF) and J = Global node point number.

ID(I,J) = 0, The ith DOF is defined as "FREE" movement allowed ("active DOF")

ID(I,J) = 1, The ith DOF is defined as "FIXED" movement not allowed

Subroutine EQUNUM numbers Active DOFs as Equation Numbers

LM(I), Local connectivity Matrix connects local DOFs to Global DOFs

Connectivity for K-th Element

DO 20 J=1,3
NND= ID(K,J)
LM(J)=ID(NND,1)
LM(J+3)=ID(NND,2)
LM(J+6)=ID(NND,3)

     20 CONTINUE

Local Element Stiffness Matrix

Element connectivity matrix

Element No. 1, i=1

Element No. 2, i=2

Element No. 3, i=3

Element No. 4, i=4

Global Stiffness Matrix

Global Stiffness Matrix Calculations