二维位势多边界问题常量元程序使用说明及举例文档格式.doc

二维位势多边界问题常量元程序使用说明及举例文档格式.doc

《二维位势多边界问题常量元程序使用说明及举例文档格式.doc》由会员分享，可在线阅读，更多相关《二维位势多边界问题常量元程序使用说明及举例文档格式.doc（12页珍藏版）》请在冰豆网上搜索。

输出结果

在GHMAMPC中调用子程序EXTINPC和LOCINPC；

EXTINPC用4个积分点的高斯数值积分法计算H和G矩阵副元素；

LOCINPC计算H和G矩阵的对角元素。

3输入数据

TITLE:

字符型的标题；

N,L,M,NC（k）,k=1,M:

单元总数，内点总数，边界的数量，每个边界的最后结点号；

X（I）,Y（I）,I=1,N:

单元端点坐标；

KODE（I）,FI（I）,I=1,N:

I单元边界条件的类型及相应的已知结点u值或q值，KODE（I）=0表示第I单元的u值已知，KODE（I）=1表示第I单元的q值已知；

CX（I）,CY（I）,I=1,L:

内点的坐标。

4其它标识符

XM,YM:

每个单元的中点（结点）的坐标；

DFI:

方程组的右端向量，解方程后成为解向量（每个单元的未知u值或q值）。

二、例题说明

例题为如下图所示的一圆环，内半径2.5cm，外半径10cm，内边界的已知温度为100℃，外边界的已知温度为0℃，用常量边界元求其稳定温度场，内边界分为32个单元，外边界分为36个单元。

三、例题输入数据

例题输入数据如下：

POTENTIALPROBLEMSWITHMULTIBOUNDARYDOMAINSUSINGCONSTANTELEMENT℃

681023668

10.0000.000

9.8481.736

9.3973.420

8.6605.000

7.6606.428

6.4287.660

5.0008.660

3.4209.397

1.7369.848

0.00010.000

-1.7369.848

-3.4209.397

-5.0008.660

-6.4287.660

-7.6606.428

-8.6605.000

-9.3973.420

-9.8481.736

-10.0000.000

-9.848-1.736

-9.397-3.420

-8.660-5.000

-7.660-6.428

-6.428-7.660

-5.000-8.660

-3.420-9.397

-1.736-9.848

0.000-10.000

1.736-9.848

3.420-9.397

5.000-8.660

6.428-7.660

7.660-6.428

8.660-5.000

9.397-3.420

9.848-1.736

2.452-0.488

2.310-0.957

2.079-1.389

1.768-1.768

1.389-2.079

0.957-2.310

0.488-2.452

0.000-2.500

-0.488-2.452

-0.957-2.310

-1.389-2.079

-1.768-1.768

-2.079-1.389

-2.310-0.957

-2.452-0.488

-2.5000.000

-2.4520.488

-2.3100.957

-2.0791.389

-1.7681.768

-1.3892.079

-0.9572.310

-0.4882.452

0.0002.500

0.4882.452

0.9572.310

1.3892.079

1.7681.768

2.0791.389

2.3100.957

2.4520.488

2.5000.000

00.000

0100.000

2.8750

3.6250

4.3750

5.1250

5.8750

6.6250

7.3750

8.1250

8.8750

9.6250

四、例题输出结果

例题计算的输出结果如下：

*******************************************************************************

POTENTIALPROBLEMSWITHMULTIBOUNDARYDOMAINSUSINGCONSTANTELEMENT

data

numberofboundaryelements=68

numberofinternalpointswherethefunctioniscalculated=10

numberofdiferentboundaries=2

lastnodeofeachboundary=36,68,

coordinatesoftheextermepointsoftheboundaryelements

pointxy

10.10000E+020.00000E+00

20.98480E+010.17360E+01

30.93970E+010.34200E+01

40.86600E+010.50000E+01

50.76600E+010.64280E+01

60.64280E+010.76600E+01

70.50000E+010.86600E+01

80.34200E+010.93970E+01

90.17360E+010.98480E+01

100.00000E+000.10000E+02

11-0.17360E+010.98480E+01

12-0.34200E+010.93970E+01

13-0.50000E+010.86600E+01

14-0.64280E+010.76600E+01

15-0.76600E+010.64280E+01

16-0.86600E+010.50000E+01

17-0.93970E+010.34200E+01

18-0.98480E+010.17360E+01

19-0.10000E+020.00000E+00

20-0.98480E+01-0.17360E+01

21-0.