边界元法程序(计算附加质量).pdfVIP

PROGRAM

BEM_ZJL

USE

IMSL

IMPLICIT

NONE

REAL(8),ALLOCATABLE,DIMENSION(:,:)

::

P0,P1,P2,P3,P4,N1,N2,N3,N01,N02,N03

!面元形心、顶点和面元全局坐标下的局部坐标系的分量和单位向量

REAL(8),ALLOCATABLE,DIMENSION(:,:)

::

AIFLC1,AIFLC2,AIFLC3,AIF1,AIF2,AIF3,AIF,CIF

!局坐标下和总体坐标下面元间的影响系数

REAL(8),ALLOCATABLE,DIMENSION(:,:)

::

VP1,VP2,VP3

!三个方向的速度势分解量

REAL(8),ALLOCATABLE,DIMENSION(:,:,:)

::

P0LC,P1LC,P2LC,P3LC,P4LC

!相对于各坐标系的局坐标

REAL(8),ALLOCATABLE

::

S0(:),NQ1(:),NQ2(:),NQ3(:),N03X(:),N03Y(:),N03Z(:)

!面元面积和法向量的模,N03的三方向表示N03X,N03Y,N03Z

REAL(8),ALLOCATABLE

::

MD1(:),MD2(:),MD3(:)

!三个方向的分布源密度

REAL(8),ALLOCATABLE

::

VP_1(:),VP_2(:),VP_3(:)

!每一个面元的三个方向上的速度势

REAL(8),ALLOCATABLE

::

M_11(:),M_22(:),M_33(:)

!附加质量分解量

REAL(8),ALLOCATABLE

::

ANS(:)

REAL(8)

::

M11,M22,M33

!三个自由度的附加质量

REAL(8)

::

A1,A2,A3,A4,A5,A6,A7,A8,A9,A10,A11,A12

!读入网格数据时使用的中间变量

REAL(8)

::

P4P1,P2P3,P2P1_P4P1,P2P3_P4P3,SIN1,SIN3,S1,S3,P2P1,P3P2,P4P3,P1P4

!计算面元投影面积以及分布源影响系数时使用的中间变量

REAL(8)

::

AIFLC1_1,AIFLC1_2,AIFLC1_3,AIFLC1_4

REAL(8)

::

AIFLC2_1,AIFLC2_2,AIFLC2_3,AIFLC2_4

REAL(8)

::

AIFLC3_1,AIFLC3_2,AIFLC3_3,AIFLC3_4

!!影响系数的四个组成元素（四个边）

REAL(8)

::

C1,C2,C3,C4,H1,H2,H3,H4,M12,M23,M34,M41,R1,R2,R3,R4

!计算方向3的影响系数时使用的中间变量

REAL(8)

::

E_1,E_2,E_3,E_4,F_1,F_2,F_3,F_4,C_1,C_2,C_3,C_4

!计算分布源密度的影响系数Cij时使用的中间变量

INTEGER,PARAMETER

::

Q1=3

REAL(8)

::

A(Q1,Q1)

!解某点的局部坐标方程时的中间数组

REAL(8)

::

S(Q1)

!解某点的局部坐标方程时的中间数组

REAL(8)

::

ANS1(Q1)

!解某点的局部坐标方程时的中间数组

REAL(8)

::

P0IP1,P0IP2,P0IP3,P0IP4

!i面元重心在j面元坐标系下距离界面元四点的距离

INTEGER

::

NEL,I,J,J1,J2,J3,J4

!网格数及循环语句中所使用的整形变量

REAL(8),PARAMETER

::

PI=3.1415926,RO=5.0

!圆周率和球的半径

REAL(8)

::

V

!球的体积

CHARACTER

(LEN=10)

::

DNAME

!定义网格文件

WRITE(*,*)

PLEASE

INPUT

THE

NUMBER

OF

ELEMENTS!

READ(*,*)

NEL

WRITE(*,*)

PLEASE

INPUT

THE

DUCUMENT

MAME.SUCH

AS

XX.TXT!

READ(*,*)

DNAME

ALLOCATE(P0(1:NEL,1:3))

ALLOCATE(P1(1:NEL,1:3))

ALLOCATE(P2(1:NEL,1:3))

ALLOCATE(P3(1:NEL,1: