Fd[3]=F[3]
#Elastic part of F 4.6
#ATENTION F~ is not in the same order of NODAL variable
if (it==0):
Ft=identity(3)
#NOT USE!!! Fd=[[1,0,0,1]]
Fvpt=identity(3)
#print(Ft)
Fvpd[0]=1.
Fvpd[1]=0. #yx
Fvpd[2]=0. #xy
Fvpd[3]=1.
Fvpd[4]=1.
else:
#F is [xx xy yx yy zz]
Fvpd[0]=Fvp[0]
Fvpd[1]=Fvp[2] #yx
Fvpd[2]=Fvp[1] #xy
Fvpd[3]=Fvp[3]