Civil Engineering Reference
In-Depth Information
(number of degrees of freedom per node), which is set to 1 for potential problems and to
3 for elasticity problems.
Zero coefficient arrays ['
U
] and ['
T
], Determine L
[
and L
K
Colloc_Points: DO i=1,Number of points P
i
Determine number of triangular sub-elements needed
Traingles: DO i=1,Number of triangles
Determine distance of P
i
to Sub-lement and No. of
Gauss points in [and KDirection
Gauss Points xsi: DO m=1,Number of Gauss in [direction
Gauss Points eta: DO k=1,Number of Gauss in Kdirection
Determine r,dsxr,Jacobian etc. for kernel computation
Node_Points: DO n=1,Number of Element Nodes
Direction_P: DO j=1,2 (direction of force P)
Determine row number for storage
Direction_Q: DO k=1,2 (direction of U,T at Q).
Determine column number for storage
Sum coefficients ['
U
]
IF (n /=
P
i
) sum ['
T
]
Figure 6.18
Structure chart for computation of ['T] and if ['U] if
P
i
is one of the element nodes
Subroutine INTEG3 is divided into two parts. The first part deals with integration
when
P
i
is not one of the nodes of the element over which the integration is made. Gauss
integration in two directions is used here. The integration of
e
ni
e
ni
'
is carried
out concurrently. It should actually be treated separately, because the functions to be
integrated have different degrees of singularity and, therefore, require a different number
'
and
Search WWH ::
Custom Search