Civil Engineering Reference
In-Depth Information
!-------------------------------------------------------------
! Add Diagonal coefficients
!-------------------------------------------------------------
Collocation_points: &
DO
m=1,Ndofs
Nod=0
DO
n=1, Nodes
DO
l=1,Ndof
IF
(m == Ndest(n,l))
THEN
Nod=n
EXIT
END IF
END DO
IF
(Nod /= 0)
EXIT
END DO
DO
k=1,Ndof
DoF=Ndest(Nod,k)
IF
(DoF /= 0)
THEN
IF
(NCode(DoF) == 1)
THEN
Nel=0
Pos=0
DO
i=1,Maxe
DO j=1,Ndofe
IF
(DoF == Ldest(i,j))
THEN
Nel=i
Pos=j
EXIT
END IF
END DO
IF
(Nel /= 0)
EXIT
END DO
F(m) = F(m) - Diag(m,k) * Elres_u(Nel,Pos)
ELSE
Lhs(m,DoF)= Lhs(m,DoF) + Diag(m,k)
END IF
END IF
END DO
END DO &
Collocation_points
!---------------------------------------------------------
! Solve system of equations
!---------------------------------------------------------
CALL
Solve(Lhs,F,u1)
CLOSE
(UNIT=2)
OPEN
(UNIT=2,FILE='BERESULTS',FORM='FORMATTED')
! Gather Element results from global result vector u1
Elements_2:&
DO
nel=1,maxe
D_o_F1: &
DO
nd=1,Ndofe
IF
(Ldest(nel,nd) /= 0)
THEN
Search WWH ::
Custom Search