Body Forces - The Boundary Element Method with Programming - page 369

SigK= SigmaK (dxr,r,E,ny,Cdim)
Direction_P: DO ii=1,Cdim
iD= Cdim*(i-1) + ii ! Position in F array
Direction_Q: DO jj=1,Cdim
Node_points: DO n=1,NodelC
F(iD)= F(iD) + Ni(n)*SigK(ii,jj)*Weit*Eps0(jj)
END DO Node_points
END DO Direction_Q
END DO Direction_P
END DO Gauss_points_eta
END DO Gauss_points_xsi
END DO Cells
END DO Colloc_points
DEALLOCATE (xPC)
DEALLOCATE (InciC)
DEALLOCATE (Ni,Elcor)
RETURN
END Subroutine Body_force

In Program General_purpose_BEM we have to insert a call to Body_force as shown

….
END DO &
Elements_1
CALL Body_force(F,CDim,xP,Nodes,Isym,E,ny)
!------------------------------------------------------------
! Add azimuthal integral for infinite regions
!------------------------------------------------------------
…

13.8.1 Input data specification for Body_force

The input data for the cell mesh and the specification of the initial strain is supplied via
file INPUT. If no body force effects are present SUBROUTINE Body_force returns
immediately. The geometry of the cell and the sequence of the input of node numbers is
shown in Figure 13.7.

INPUT DATA SPECIFICATION FOR Body_force

1.0 Node specification
Nodes

Number of cell nodes

2.0 Element specification
Ncells Number of cells

3.0 Initial strain specification
Eps0(1:2) Initial strains

H ,

H

0

0