二维位势多边界问题常量元程序使用说明及举例Word格式文档下载.doc
《二维位势多边界问题常量元程序使用说明及举例Word格式文档下载.doc》由会员分享,可在线阅读,更多相关《二维位势多边界问题常量元程序使用说明及举例Word格式文档下载.doc(12页珍藏版)》请在冰点文库上搜索。
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.93970E+01-0.34200E+01
22-0.86600E+01-0.50000E+01
23-0.76600E+01-0.64280E+01
24-0.64280E+01-0.76600E+01
25-0.50000E+01-0.86600E+01
26-0.34200E+01-0.93970E+01
27-0.17360E+01-0.98480E+01
280.00000E+00-0.10000E+02
290.17360E+01-0.98480E+01
300.34200E+01-0.93970E+01
310.50000E+01-0.86600E+01
320.64280E+01-0.76600E+01
330.76600E+01-0.64280E+01
340.86600E+01-0.50000E+01
350.93970E+01-0.34200E+01
360.98480E+01-0.17360E+01
370.24520E+01-0.48800E+00
380.23100E+01-0.95700E+00
390.20790E+01-0.13890E+01
400.17680E+01-0.17680E+01
410.13890E+01-0.20790E+01
420.95700E+00-0.23100E+01
430.48800E+00-0.24520E+01
440.00000E+00-0.25000E+01
45-0.48800E+00-0.24520E+01
46-0.95700E+00-0.23100E+01
47-0.13890E+01-0.20790E+01
48-0.17680E+01-0.17680E+01
49-0.20790E+01-0.13890E+01
50-0.23100E+01-0.95700E+00
51-0.24520E+01-0.48800E+00
52-0.25000E+010.00000E+00
53-0.24520E+010.48800E+00
54-0.23100E+010.95700E+00
55-0.20790E+010.13890E+01
56-0.17680E+010.17680E+01
57-0.13890E+010.20790E+01
58-0.95700E+000.23100E+01
59-0.48800E+000.24520E+01
600.00000E+000.25000E+01
610.48800E+000.24520E+01
620.95700E+000.23100E+01
630.13890E+010.20790E+01
640.17680E+010.17680E+01
650.20790E+010.13890E+01
660.23100E+010.95700E+00
670.24520E+010.48800E+00
680.25000E+010.00000E+00
boundadryconditions
nodekodeprescribedvalue
100.00000E+00
200.00000E+00
300.00000E+00
400.00000E+00
500.00000E+00
600.00000E+00
700.00000E+00
800.00000E+00
900.00000E+00
1000.00000E+00
1100.00000E+00
1200.00000E+00
1300.00000E+00
1400.00000E+00
1500.00000E+00
1600.00000E+00
1700.00000E+00
1800.00000E+00
1900.00000E+00
2000.00000E+00
2100.00000E+00
2200.00000E+00
2300.00000E+00
2400.00000E+00
2500.00000E+00
2600.00000E+00
2700.00000E+00
2800.00000E+00
2900.00000E+00
3000.00000E+00
3100.00000E+00
3200.00000E+00
3300.00000E+00
3400.00000E+00
3500.00000E+00
3600.00000E+00
3700.10000E+03
3800.10000E+03
3900.10000E+03
4000.10000E+03
4100.10000E+03
4200.10000E+03
4300.10000E+03
4400.10000E+03
4500.10000E+03
4600.10000E+03
4700.10000E+03
4800.10000E+03
4900.10000E+03
5000.10000E+03
5100.10000E+03
5200.10000E+03
5300.10000E+03
5400.10000E+03
5500.10000E+03
5600.10000E+03
5700.10000E+03
5800.10000E+03
5900.10000E+03
6000.10000E+03
6100.10000E+03
6200.10000E+03
6300.10000E+03
6400.10000E+03
6500.10000E+03
6600.10000E+03
6700.10000E+03
6800.10000E+03
*******************************************************************************
results
boundarynodes
xypotentialpotentialderivative
0.99240E+010.86800E+000.00000E+00-0.72169E+01
0.96225E+010.25780E+010.00000E+00-0.72171E+01
0.90285E+010.42100E+010.00000E+00-0.72170E+01
0.81600E+010.57140