search

luchete80 / EPSol

0 stars 0 forks source link

Boundary conditions application only works for 1 VAR (F), needs to be update to several #90

Closed luchete80 closed 3 years ago

luchete80 commented 3 years ago 
for dy in range(ney+1): 
    inode=dy*dnode
    #print("node",inode)   
    #Deformation gradient F
    for i in range ( var_dim [ 0 ] ):
        idof = var_dim[0] * inode + i
        #print ("idof",idof)
        for j in range(dof):
            Kglob[ idof , j ] = 0
            #Rglob[ j ] -= Kglob[j,idof] * 0 #dU=0, U=1(idof)
            Kglob[ j ,idof ] = 0

        Kglob[idof,idof] = 1
        Rglob[idof  ] = 0           #F INCREMENT (dF) IS NULL!!!!!

    #Sigma is zero, Internal variable s ,      
    if numvars == 2:
        idofs = idof + var_dim[0]
        for i in range(dof):                    
            Kglob[ idofs , i ] = 0
            #Rglob[ i ] -= Kglob[ i, idofs ] * 0 #1 is R(idof)  
            Kglob[ i , idofs ] = 0                

        Kglob[idofs,idofs] = 1
        Rglob[idofs  ] = 0