fig.dat

有限元求解柏松方程。本文采用FORTRAN语言编制程序。程序中大部分变量采用有名公共区存储方式存储

  VARIABLES= "X" "Y" "F" "DIFX" "DIFY" 
 ZONE I=         21J=          21
        .0000       .0000       .0000       .9959       .0000
        .1571       .0000       .1564       .9084      -.1587
        .3142       .0000       .3090       .8780      -.3924
        .4712       .0000       .4540       .8260      -.6160
        .6283       .0000       .5878       .7535      -.8242
        .7854       .0000       .7071       .6625     -1.0121
        .9425       .0000       .8090       .5551     -1.1750
       1.0996       .0000       .8910       .4341     -1.3090
       1.2566       .0000       .9511       .3023     -1.4110
       1.4137       .0000       .9877       .1629     -1.4787
       1.5708       .0000      1.0000       .0195     -1.5106
       1.7279       .0000       .9877      -.1245     -1.5063
       1.8850       .0000       .9511      -.2658     -1.4666
       2.0420       .0000       .8910      -.4008     -1.3934
       2.1991       .0000       .8090      -.5265     -1.2898
       2.3562       .0000       .7071      -.6400     -1.1607
       2.5133       .0000       .5878      -.7391     -1.0125
       2.6704       .0000       .4540      -.8219      -.8542
       2.8274       .0000       .3090      -.8878      -.6982
       2.9845       .0000       .1564      -.9373      -.5624
       3.1416       .0000       .0000      -.9628      -.4871
        .0000       .1571       .0000       .8372      -.0794
        .1571       .1571       .1190       .7573      -.2179
        .3142       .1571       .2353       .7288      -.4205
        .4712       .1571       .3457       .6822      -.6126
        .6283       .1571       .4477       .6186      -.7895
        .7854       .1571       .5387       .5397      -.9468
        .9425       .1571       .6164       .4474     -1.0809
       1.0996       .1571       .6789       .3440     -1.1885
       1.2566       .1571       .7246       .2320     -1.2671
       1.4137       .1571       .7525       .1141     -1.3149
       1.5708       .1571       .7617      -.0070     -1.3310
       1.7279       .1571       .7519      -.1285     -1.3155
       1.8850       .1571       .7234      -.2476     -1.2694
       2.0420       .1571       .6765      -.3618     -1.1949
       2.1991       .1571       .6123      -.4689     -1.0954
       2.3562       .1571       .5320      -.5673      -.9757
       2.5133       .1571       .4368      -.6562      -.8424
       2.6704       .1571       .3282      -.7359      -.7045
       2.8274       .1571       .2074      -.8091      -.5745
       2.9845       .1571       .0747      -.8825      -.4709
       3.1416       .1571      -.0713      -.9099      -.4062
        .0000       .3142       .0000       .6373      -.0602
        .1571       .3142       .0906       .5769      -.1655
        .3142       .3142       .1792       .5554      -.3196
        .4712       .3142       .2634       .5200      -.4658
        .6283       .3142       .3411       .4717      -.6005
        .7854       .3142       .4105       .4116      -.7204
        .9425       .3142       .4698       .3412      -.8225
       1.0996       .3142       .5174       .2622      -.9044
       1.2566       .3142       .5523       .1765      -.9642
       1.4137       .3142       .5734       .0860     -1.0006
       1.5708       .3142       .5802      -.0073     -1.0129
       1.7279       .3142       .5724      -.1014     -1.0011
       1.8850       .3142       .5499      -.1946      -.9662
       2.0420       .3142       .5130      -.2853      -.9097
       2.1991       .3142       .4622      -.3726      -.8343
       2.3562       .3142       .3979      -.4561      -.7437
       2.5133       .3142       .3207      -.5368      -.6430
       2.6704       .3142       .2308      -.6174      -.5386
       2.8274       .3142       .1278      -.7035      -.4391
       2.9845       .3142       .0100      -.8063      -.3547
       3.1416       .3142      -.1267      -.8474      -.2879
        .0000       .4712       .0000       .4856      -.0458
        .1571       .4712       .0691       .4397      -.1257
        .3142       .4712       .1366       .4235      -.2429
        .4712       .4712       .2008       .3966      -.3541
        .6283       .4712       .2601       .3599      -.4566
        .7854       .4712       .3130       .3141      -.5478
        .9425       .4712       .3583       .2604      -.6254
       1.0996       .4712       .3946       .2001      -.6876
       1.2566       .4712       .4212       .1344      -.7329
       1.4137       .4712       .4372       .0648      -.7603
       1.5708       .4712       .4422      -.0072      -.7693
       1.7279       .4712       .4359      -.0805      -.7598
       1.8850       .4712       .4181      -.1539      -.7324
       2.0420       .4712       .3889      -.2266      -.6884
       2.1991       .4712       .3483      -.2986      -.6297
       2.3562       .4712       .2964      -.3705      -.5587
       2.5133       .4712       .2332      -.4443      -.4790
       2.6704       .4712       .1580      -.5239      -.3945
       2.8274       .4712       .0693      -.6161      -.3101
       2.9845       .4712      -.0355      -.7325      -.2304
       3.1416       .4712      -.1615      -.7752      -.1581
        .0000       .6283       .0000       .3703      -.0348
        .1571       .6283       .0527       .3355      -.0955
        .3142       .6283       .1042       .3232      -.1846
        .4712       .6283       .1532       .3028      -.2691
        .6283       .6283       .1985       .2749      -.3469
        .7854       .6283       .2389       .2400      -.4161
        .9425       .6283       .2735       .1991      -.4750
       1.0996       .6283       .3013       .1530      -.5220
       1.2566       .6283       .3216       .1027      -.5561
       1.4137       .6283       .3339       .0492      -.5763
       1.5708       .6283       .3376      -.0065      -.5824
       1.7279       .6283       .3325      -.0636      -.5740
       1.8850       .6283       .3185      -.1216      -.5517
       2.0420       .6283       .2953      -.1802      -.5161
       2.1991       .6283       .2629      -.2398      -.4685
       2.3562       .6283       .2211      -.3016      -.4102
       2.5133       .6283       .1692      -.3681      -.3433
       2.6704       .6283       .1062      -.4437      -.2698
       2.8274       .6283       .0303      -.5352      -.1916
       2.9845       .6283      -.0617      -.6529      -.1098
       3.1416       .6283      -.1750      -.6899      -.0315
        .0000       .7854       .0000       .2827      -.0264
        .1571       .7854       .0403       .2563      -.0725
        .3142       .7854       .0796       .2470      -.1401
        .4712       .7854       .1171       .2316      -.2042
        .6283       .7854       .1517       .2104      -.2633
        .7854       .7854       .1827       .1839      -.3157
        .9425       .7854       .2092       .1528      -.3601
       1.0996       .7854       .2306       .1176      -.3955
       1.2566       .7854       .2462       .0792      -.4208
       1.4137       .7854       .2557       .0382      -.4353
       1.5708       .7854       .2586      -.0048      -.4388
       1.7279       .7854       .2547      -.0492      -.4308
       1.8850       .7854       .2438      -.0948      -.4117
       2.0420       .7854       .2257      -.1418      -.3815
       2.1991       .7854       .2001      -.1907      -.3410
       2.3562       .7854       .1666      -.2430      -.2907
       2.5133       .7854       .1245      -.3013      -.2313
       2.6704       .7854       .0726      -.3698      -.1633
       2.8274       .7854       .0089      -.4546      -.0864
       2.9845       .7854      -.0698      -.5646       .0007
       3.1416       .7854      -.1683      -.5910       .0845
        .0000       .9425       .0000       .2163      -.0200
        .1571       .9425       .0308       .1963      -.0550
        .3142       .9425       .0610       .1892      -.1062
        .4712       .9425       .0897       .1776      -.1548
        .6283       .9425       .1163       .1615      -.1994
        .7854       .9425       .1401       .1415      -.2390
        .9425       .9425       .1605       .1179      -.2724
       1.0996       .9425       .1770       .0913      -.2987
       1.2566       .9425       .1892       .0621      -.3172
       1.4137       .9425       .1967       .0309      -.3272
       1.5708       .9425       .1992      -.0019      -.3283
       1.7279       .9425       .1965      -.0361      -.3203
       1.8850       .9425       .1884      -.0716      -.3030
       2.0420       .9425       .1746      -.1085      -.2763
       2.1991       .9425       .1550      -.1477      -.2402
       2.3562       .9425       .1289      -.1906      -.1947
       2.5133       .9425       .0958      -.2395      -.1395
       2.6704       .9425       .0544      -.2980      -.0736
       2.8274       .9425       .0028      -.3714       .0044
       2.9845       .9425      -.0617      -.4669       .0975
       3.1416       .9425      -.1433      -.4798       .1858
        .0000      1.0996       .0000       .1659      -.0152
        .1571      1.0996       .0237       .1508      -.0416
        .3142      1.0996       .0468       .1455      -.0804
        .4712      1.0996       .0689       .1367      -.1171
        .6283      1.0996       .0894       .1247      -.1508
        .7854      1.0996       .1078       .1097      -.1805
        .9425      1.0996       .1237       .0920      -.2054
       1.0996      1.0996       .1367       .0720      -.2247
       1.2566      1.0996       .1464       .0501      -.2379
       1.4137      1.0996       .1526       .0267      -.2442
       1.5708      1.0996       .1550       .0020      -.2434
       1.7279      1.0996       .1536      -.0237      -.2350
       1.8850      1.0996       .1480      -.0504      -.2187
       2.0420      1.0996       .1382      -.0785      -.1941
       2.1991      1.0996       .1239      -.1085      -.1608
       2.3562      1.0996       .1047      -.1417      -.1181
       2.5133      1.0996       .0800      -.1800      -.0650
       2.6704      1.0996       .0488      -.2263       .0001
       2.8274      1.0996       .0097      -.2846       .0801
       2.9845      1.0996      -.0397      -.3606       .1790
       3.1416      1.0996      -.1027      -.3580       .2702
        .0000      1.2566       .0000       .1278      -.0115
        .1571      1.2566       .0183       .1164      -.0315
        .3142      1.2566       .0362       .1125      -.0607
        .4712      1.2566       .0533       .1060      -.0883
        .6283      1.2566       .0692       .0971      -.1136
        .7854      1.2566       .0836       .0859      -.1358
        .9425      1.2566       .0961       .0729      -.1542
       1.0996      1.2566       .1064       .0582      -.1682
       1.2566      1.2566       .1144       .0421      -.1772
       1.4137      1.2566       .1198       .0250      -.1807
       1.5708      1.2566       .1224       .0071      -.1782
       1.7279      1.2566       .1223      -.0114      -.1694
       1.8850      1.2566       .1192      -.0306      -.1536
       2.0420      1.2566       .1131      -.0505      -.1303
       2.1991      1.2566       .1038      -.0717      -.0987
       2.3562      1.2566       .0912      -.0949      -.0578
       2.5133      1.2566       .0747      -.1216      -.0060
       2.6704      1.2566       .0537      -.1538       .0591
       2.8274      1.2566       .0272      -.1944       .1409
       2.9845      1.2566      -.0064      -.2472       .2443
       3.1416      1.2566      -.0493      -.2284       .3369
        .0000      1.4137       .0000       .0990      -.0087
        .1571      1.4137       .0142       .0905      -.0237
        .3142      1.4137       .0281       .0876      -.0457
        .4712      1.4137       .0415       .0829      -.0664
        .6283      1.4137       .0539       .0764      -.0853
        .7854      1.4137       .0653       .0684      -.1017
        .9425      1.4137       .0753       .0590      -.1151
       1.0996      1.4137       .0838       .0485      -.1250
       1.2566      1.4137       .0906       .0372      -.1308
       1.4137      1.4137       .0956       .0253      -.1321
       1.5708      1.4137       .0987       .0131      -.1283
       1.7279      1.4137       .1000       .0008      -.1190
       1.8850      1.4137       .0993      -.0115      -.1036
       2.0420      1.4137       .0968      -.0238      -.0812
       2.1991      1.4137       .0923      -.0363      -.0509
       2.3562      1.4137       .0859      -.0495      -.0113
       2.5133      1.4137       .0774      -.0639       .0395
       2.6704      1.4137       .0666      -.0809       .1044
       2.8274      1.4137       .0530      -.1018       .1872
       2.9845      1.4137       .0358      -.1290       .2932
       3.1416      1.4137       .0139      -.0936       .3854
        .0000      1.5708       .0000       .0774      -.0065
        .1571      1.5708       .0111       .0710      -.0178
        .3142      1.5708       .0221       .0690      -.0342
        .4712      1.5708       .0326       .0656      -.0497
        .6283      1.5708       .0425       .0611      -.0637
        .7854      1.5708       .0517       .0555      -.0758
        .9425      1.5708       .0599       .0491      -.0854
       1.0996      1.5708       .0671       .0420      -.0921
       1.2566      1.5708       .0732       .0346      -.0955
       1.4137      1.5708       .0781       .0272      -.0950