September 2013
Intermediate to advanced
320 pages
11h 25m
English
Content preview from Finite Element Techniques for Fluid Flow
Become an O’Reilly member and get unlimited access to this title plus top books and audiobooks from O’Reilly and nearly 200 top publishers, thousands of courses curated by job role, 150+ live events each month,
Start your free trial
96 THE FINITE ELEMENT TECHNIQUE
A(IR,ID)=A(IR,ID)+S(IM,IN)
200 CONTINUE
101 C0NTINUE
c
C APPLY BOUNDARY CONDITIONS
C
CALL ZER(B,N,1)
KF = MB-1
DO 300
K
= 1,N4
N5
= IBC(K)
DO 400
K1
= 1,KF
IR = N5-K1
IF(IR.LT.l) GOTO 450
B(IR)=B(IR)-A(IR,(Kl + l))*VIBC(K)
A(IR,(K1 +
1))
=
0
450
IR
=
N5
+
K1
IF(IR.GT.N)G0TO 400
B(IR)
= B(IR)-A(N5,K1 + 1)*VIBC(K)
A(N5,(K1 + 1))=0
400 CONTINUE
A(N5,1)=1
B(N5)
= VIBC(K)
300 CONTINUE
C
C CALL SOLVER
C
CALL BANDSOL(A,B,N,MB).
C
C CALL SUBROUTINE T0 PRINT VALUES OF U
C
CALL 0UTW(B,X,Y,N)
999 STOP
END
REFERENCES
1.
DESAi,
c. s. and
ABEL, J. F.,
Introduction to the
Finite
Element
Method,
Van Nostrand
Reinhold Co (1972 ...