Civil Engineering Reference
In-Depth Information
I= JJ(N)
J= KK(N)
B= 0.
IF(I .EQ. J) B= Cdim*C3*DXR(K)
IF(I .EQ. K) B= B + ny*DXR(J)
IF(J .EQ. K) B= B + ny*DXR(I)
B= COSTH *(B - C6*DXR(I)*DXR(J)*DXR(K) )
C= DXR(J)*DXR(K)*ny
IF(J .EQ.K) C= C + C3
C= C*VNORM(I)
B= B+C
C= DXR(I)*DXR(K)*ny
IF(I .EQ. K) C=C + C3
C= C*VNORM(J)
B= B+C
C= DXR(I)*DXR(J)*Cdim*C3
IF(I .EQ. J) C= C - C7
C= C*VNORM(K)
US(K,N)= (B + C)*C5/R/R
IF(Cdim .EQ. 3) US(K,N)= US(K,N)/2./R
END DO &
Stress_components
END DO &
Coor_directions
RETURN
END
The discretised form of equation (9.38) is written as
E
N
E
N
(9.48)
¦
¦
¦¦
e
n
e
n
e
n
e
n
u
P
'
U
t
'
T
u
e
1
n
1
e
1
n
1
where
e
n
³
e
n
³
'
U
U
(
P
,
)
N
dS
(
Q
)
;
'
T
T
(
P
,
)
N
dS
(
Q
)
(9.49)
a
n
a
n
S
S
e
e
The discretised form of equation (9.41) is written as
E
N
E
N
(9.50)
¦
¦
¦¦
e
n
e
n
e
n
e
n
V
P
'
S
t
'
R
u
e
1
n
1
e
1
n
1
where
(9.51)
e
n
³
e
n
³
'
S
S
(
P
,
Q
)
N
dS
(
Q
)
;
'
R
R
(
P
,
Q
)
N
dS
(
Q
)
a
n
a
n
S
S
e
e
Search WWH ::
Custom Search