有限元程序输出文件的问题

shmily1231988

这是网上下载的一个有限元程序  现有初始输入文件（1.dat），此程序无法正确读入，请各位帮忙指导一下，小弟是刚接触Fortran，希望能得到各位耐心的指导，谢谢

楚香饭

我这里运行后是这样：
                    P R O G R A M    F L a g S H y P

          Finite element LArGe Strain HYperelasticity Program

Is the problem starting from scratch (y/n) ?:
y
Enter the data file name   :
1.dat
Enter the results file name:
2.dat
Enter the restart file name:
3.dat
Unknown element type
elinfo: Unknown element type

你呢？

珊瑚虫

根据程序的内容，1.dat 的数据格式与程序的读入数据格式完全不符，尤其是节点坐标数据位置，少了一列
ps ：这个程序写的一般，建议楼主选择其他程序作为有限元入门。
下面是根据猜测调整的输入文件格式，注意这个文件格式也有问题，节点坐标列少1列数字，因为没有具体读程序所以不太清楚该写什么，原来的1.dat里面也没有。此外，在这个里面我假定原来的1.dat中的20行数字是荷载，至于约束信息，面荷载，重力荷载等没看出来对应原1.dat中的哪一个数字，所以我写的0。希望对楼主有用。

plane strain for TRI 0.41
quad4
113  
   1    0.678412413137      0.500000000000   
   2    0.500000000000      0.678412413137   
   3    0.321587586863      0.500000000000   
   4    0.500000000000      0.321587586863   
   5    0.670005307592      0.554121941501   
   6    0.645576306020      0.603143241596   
   7    0.606009913237      0.643502221087   
   8    0.555805069669      0.669460270744   
   9    0.445878058499      0.670005307592   
  10    0.396856758404      0.645576306020   
  11    0.356497778913      0.606009913237   
  12    0.330539729256      0.555805069669   
  13    0.329994692408      0.445878058499   
  14    0.354423693980      0.396856758404   
  15    0.393990086763      0.356497778913   
  16    0.444194930331      0.330539729256   
  17    0.554121941501      0.329994692408   
  18    0.603143241596      0.354423693980   
  19    0.643502221087      0.393990086763   
  20    0.669460270744      0.444194930331
   
  21    0.476725030056      0.543019367480   
  22    0.513790554854      0.521219916119   
  23    0.453359604385      0.511461925367   
  24    0.489858144548      0.484942247136   
  25    0.538867282079      0.455123506627   
  26    0.559325813685      0.548095025781   
  27    0.535597207408      0.565848339208   
  28    0.465384710877      0.592757141494   
  29    0.438456950694      0.561953830320   
  30    0.401067378686      0.507937411315   
  31    0.453476897176      0.457376604743   
  32    0.554336416882      0.507367871062   
  33    0.600762143691      0.518769729388   
  34    0.559360223726      0.592943435486   
  35    0.487803436975      0.627531588946   
  36    0.431082043803      0.611270219314   
  37    0.383819269235      0.542405877584   
  38    0.385191886472      0.467695704391   
  39    0.432465187463      0.400988216053   
  40    0.497701418988      0.382914296090
   
  41    0.591613904625      0.404750128198   
  42    0.609586721790      0.471983784272   
  43    0.587346673599      0.561656353887   
  44    0.524392150154      0.613997358209   
  45    0.455119694805      0.636312262478   
  46    0.405501488131      0.579248595619   
  47    0.374861451911      0.504404437222   
  48    0.407057512507      0.432534358878   
  49    0.457845685577      0.373138707292   
  50    0.545819254379      0.399662651992   
  51    0.635114617340      0.441555171449   
  52    0.624778756084      0.412289672824   
  53    0.500069166892      0.578503866003   
  54    0.411049288406      0.531331050730   
  55    0.415198764344      0.488587372474   
  56    0.482755536212      0.425860603651   
  57     0.00000000000       0.00000000000   
  58     1.00000000000       0.00000000000   
  59     1.00000000000       1.00000000000   
  60     0.00000000000       1.00000000000
   
  61    0.212129992247       0.00000000000   
  62    0.383285904260       0.00000000000   
  63    0.575740015507       0.00000000000   
  64    0.787870007753       0.00000000000   
  65     1.00000000000      0.212129992247   
  66     1.00000000000      0.383285904260   
  67     1.00000000000      0.575740015507   
  68     1.00000000000      0.787870007753   
  69    0.787870007753       1.00000000000   
  70    0.616714095740       1.00000000000   
  71    0.424259984493       1.00000000000   
  72    0.212129992247       1.00000000000   
  73     0.00000000000      0.787870007753   
  74     0.00000000000      0.616714095740   
  75     0.00000000000      0.424259984493   
  76     0.00000000000      0.212129992247   
  77    0.125571160599      0.314124590867   
  78    0.825647614980      0.173373548552   
  79    0.219161802387      0.366385160104   
  80    0.826751637609      0.817904007787
   
  81    0.630098291862      0.847335787503   
  82    0.208140386618      0.825218712790   
  83    0.143677948862      0.647811800867   
  84    0.194821144778      0.126355911515   
  85    0.417293789765      0.130275084140   
  86    0.725623257574      0.333338821405   
  87    0.843391155931      0.374566502743   
  88    0.421280858092      0.875294335300   
  89    0.124574757403      0.479864901533   
  90    0.639751438171      0.157516914797   
  91    0.620401331670      0.680771321871   
  92    0.567835523363      0.736597782176   
  93    0.461555148954      0.754780311019   
  94    0.280359130151      0.586773840539   
  95    0.287186867231      0.530736405085   
  96    0.302491585306      0.497523545978   
  97    0.357451219477      0.344824269526   
  98    0.402798296829      0.311629190989   
  99    0.460336049619      0.284411497270   
100    0.662524086114      0.341728445740
   
101    0.796667726657      0.640430065474   
102    0.681182034971      0.687756097944   
103    0.361641073801      0.721875308426   
104    0.291409750011      0.652101017499   
105    0.284981363614      0.445455648274   
106    0.317438608624      0.388382610606   
107    0.515951744122      0.258658749518   
108    0.602417146410      0.270785237388   
109    0.346354310008      0.218706921237   
110    0.216550550693      0.228800269805   
111    0.333433833722      0.278057947224   
112    0.267204652327      0.309192815458   
113    0.217196172584      0.449224813660
102
   1   2      40    56    39    49
   2   2      25    56    40    50
   3   2      29    28    36    46
   4   2      36    28    35    45
   5   2      23    24    22    21
   6   2      21    29    54    23
   7   2      29    21    53    28
   8   2      27    53    21    22
   9   2      32    26    27    22
  10   2      22    24    25    32
  11   2      54    30    55    23
  12   2      24    23    55    31
  13   2      31    56    25    24
  14   2      31    48    39    56
  15   2      55    38    48    31
  16   2      30    47    38    55
  17   2      54    37    47    30
  18   2      29    46    37    54
  19   2      28    53    44    35
  20   2      27    34    44    53

  21   2      26    43    34    27
  22   2      32    33    43    26
  23   2      50    41    42    25
  24   2      25    42    33    32
  25   2      52    51    42    41
  26   2      41    18    19    52
  27   2      52    19    20    51
  28   2      51    20     1    42
  29   2      50    17    18    41
  30   2      40     4    17    50
  31   2      49    16     4    40
  32   2      39    15    16    49
  33   2      48    14    15    39
  34   2      38    13    14    48
  35   2      47     3    13    38
  36   2      37    12     3    47
  37   2      46    11    12    37
  38   2      36    10    11    46
  39   2      45     9    10    36
  40   2      35     2     9    45

  41   2      44     8     2    35
  42   2      34     7     8    44
  43   2      43     6     7    34
  44   2      33     5     6    43
  45   2      42     1     5    33
  46   1      76    84   110    77
  47   1      77    89    75    76
  48   1      79   113    89    77
  49   1     110   112    79    77
  50   1      87    86    90    78
  51   1      87    78    65    66
  52   1      78    64    58    65
  53   1      78    90    63    64
  54   1      88    82   103    93
  55   1      82    83   104   103
  56   1     107   108    18    17
  57   1     108   100    19    18
  58   1     105   113    79   106
  59   1     106    79   112    97
  60   1     113    95    94    89

  61   1     105    96    95   113
  62   1     112   111    98    97
  63   1     112   110   109   111
  64   1      94   104    83    89
  65   1     109    99    98   111
  66   1      84    85   109   110
  67   1      85   107    99   109
  68   1      90   108   107    85
  69   1      90    86   100   108
  70   1      93    92    81    88
  71   1      92    91   102    81
  72   1     102   101    80    81
  73   1     101    67    68    80
  74   1      87    66    67   101
  75   1      85    62    63    90
  76   1      88    71    72    82
  77   1      81    70    71    88
  78   1      83    74    75    89
  79   1      76    57    61    84
  80   1      68    59    69    80

  81   1      72    60    73    82
  82   1      82    73    74    83
  83   1      80    69    70    81
  84   1      84    61    62    85
  85   1      91     7     6   102
  86   1      92     8     7    91
  87   1      93     2     8    92
  88   1     102     6     5   101
  89   1     103     9     2    93
  90   1     104    10     9   103
  91   1      94    11    10   104
  92   1      95    12    11    94
  93   1      96     3    12    95
  94   1     105    13     3    96
  95   1     106    14    13   105
  96   1      97    15    14   106
  97   1      98    16    15    97
  98   1      99     4    16    98
  99   1     107    17     4    99
100   1      86    20    19   100

101   1      87     1    20    86
102   1     101     5     1    87
2
  1,1
3.45 3.45 0.35  1.278  1.0  0.0  0.0
  2,1
73.1  73.1  0.22  29.959  1.0  0.0  0.0
20    0   0   0  
  57    1    1    0     0
  58    1    1    0     0
  59    1    1    0     0
  60    1    1    0     0
  61    0    1    0     0
  62    0    1    0     0
  63    0    1    0     0
  64    0    1    0     0
  65    1    0    0     0
  66    1    0    0     0
  67    1    0    0     0
  68    1    0    0     0
  69    0    1    0     0
  70    0    1    0     0
  71    0    1    0     0
  72    0    1    0     0
  73    1    0    0     0
  74    1    0    0     0
  75    1    0    0     0

shmily1231988

楚香饭 发表于 2014-11-18 22:34
我这里运行后是这样：
                    P R O G R A M    F L a g S H y P

我运行的也是这样的结果

shmily1231988

谢谢  但是按照你给修改过后的格式  运行还是错误的

shmily1231988

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

plane strain for TRI 0.41
quad4
113
   1   0  0.678412413137      0.500000000000    
   2   0  0.500000000000      0.678412413137    
   3   0  0.321587586863      0.500000000000    
   4   0  0.500000000000      0.321587586863    
   5   0  0.670005307592      0.554121941501    
   6   0  0.645576306020      0.603143241596    
   7   0  0.606009913237      0.643502221087    
   8   0  0.555805069669      0.669460270744    
   9   0  0.445878058499      0.670005307592    
  10   0  0.396856758404      0.645576306020    
  11   0  0.356497778913      0.606009913237    
  12   0  0.330539729256      0.555805069669    
  13   0  0.329994692408      0.445878058499    
  14   0  0.354423693980      0.396856758404    
  15   0  0.393990086763      0.356497778913    
  16   0  0.444194930331      0.330539729256    
  17   0  0.554121941501      0.329994692408    
  18   0  0.603143241596      0.354423693980    
  19   0  0.643502221087      0.393990086763    
  20   0  0.669460270744      0.444194930331    
  21   0  0.476725030056      0.543019367480    
  22   0  0.513790554854      0.521219916119    
  23   0  0.453359604385      0.511461925367    
  24   0  0.489858144548      0.484942247136    
  25   0  0.538867282079      0.455123506627    
  26   0  0.559325813685      0.548095025781    
  27   0  0.535597207408      0.565848339208    
  28   0  0.465384710877      0.592757141494    
  29   0  0.438456950694      0.561953830320    
  30   0  0.401067378686      0.507937411315    
  31   0  0.453476897176      0.457376604743    
  32   0  0.554336416882      0.507367871062    
  33   0  0.600762143691      0.518769729388    
  34   0  0.559360223726      0.592943435486   
  35   0  0.487803436975      0.627531588946    
  36   0  0.431082043803      0.611270219314    
  37   0  0.383819269235      0.542405877584    
  38   0  0.385191886472      0.467695704391    
  39   0  0.432465187463      0.400988216053    
  40   0  0.497701418988      0.382914296090    
  41   0  0.591613904625      0.404750128198    
  42   0  0.609586721790      0.471983784272    
  43   0  0.587346673599      0.561656353887    
  44   0  0.524392150154      0.613997358209    
  45   0  0.455119694805      0.636312262478    
  46   0  0.405501488131      0.579248595619    
  47   0  0.374861451911      0.504404437222    
  48   0  0.407057512507      0.432534358878    
  49   0  0.457845685577      0.373138707292    
  50   0  0.545819254379      0.399662651992    
  51   0  0.635114617340      0.441555171449    
  52   0  0.624778756084      0.412289672824    
  53   0  0.500069166892      0.578503866003    
  54   0  0.411049288406      0.531331050730    
  55   0  0.415198764344      0.488587372474    
  56   0  0.482755536212      0.425860603651    
  57   3   0.00000000000       0.00000000000    
  58   3   1.00000000000       0.00000000000    
  59   3   1.00000000000       1.00000000000    
  60   3   0.00000000000       1.00000000000    
  61   2  0.212129992247       0.00000000000    
  62   2  0.383285904260       0.00000000000    
  63   2  0.575740015507       0.00000000000    
  64   2  0.787870007753       0.00000000000    
  65   1   1.00000000000      0.212129992247    
  66   1   1.00000000000      0.383285904260    
  67   1   1.00000000000      0.575740015507    
  68   1   1.00000000000      0.787870007753    
  69   2  0.787870007753       1.00000000000    
  70   2  0.616714095740       1.00000000000    
  71   2  0.424259984493       1.00000000000    
  72   2  0.212129992247       1.00000000000    
  73   1   0.00000000000      0.787870007753    
  74   1   0.00000000000      0.616714095740    
  75   1   0.00000000000      0.424259984493    
  76   1   0.00000000000      0.212129992247    
  77   0  0.125571160599      0.314124590867    
  78   0  0.825647614980      0.173373548552    
  79   0  0.219161802387      0.366385160104    
  80   0  0.826751637609      0.817904007787    
  81   0  0.630098291862      0.847335787503    
  82   0  0.208140386618      0.825218712790    
  83   0  0.143677948862      0.647811800867    
  84   0  0.194821144778      0.126355911515    
  85   0  0.417293789765      0.130275084140    
  86   0  0.725623257574      0.333338821405    
  87   0  0.843391155931      0.374566502743    
  88   0  0.421280858092      0.875294335300    
  89   0  0.124574757403      0.479864901533    
  90   0  0.639751438171      0.157516914797    
  91   0  0.620401331670      0.680771321871    
  92   0  0.567835523363      0.736597782176    
  93   0  0.461555148954      0.754780311019    
  94   0  0.280359130151      0.586773840539    
  95   0  0.287186867231      0.530736405085    
  96   0  0.302491585306      0.497523545978    
  97   0  0.357451219477      0.344824269526    
  98   0  0.402798296829      0.311629190989    
  99   0  0.460336049619      0.284411497270    
 100   0  0.662524086114      0.341728445740    
 101   0  0.796667726657      0.640430065474    
 102   0  0.681182034971      0.687756097944    
 103   0  0.361641073801      0.721875308426    
 104   0  0.291409750011      0.652101017499    
 105   0  0.284981363614      0.445455648274    
 106   0  0.317438608624      0.388382610606    
 107   0  0.515951744122      0.258658749518    
 108   0  0.602417146410      0.270785237388    
 109   0  0.346354310008      0.218706921237    
 110   0  0.216550550693      0.228800269805    
 111   0  0.333433833722      0.278057947224    
 112   0  0.267204652327      0.309192815458    
 113   0  0.217196172584      0.449224813660
102
  1   2   40    56    39    49
  2   2   25    56    40    50
  3   2   29    28    36    46
  4   2   36    28    35    45
  5   2   23    24    22    21
  6   2   21    29    54    23
  7   2   29    21    53    28
  8   2   27    53    21    22
  9   2   32    26    27    22
 10   2   22    24    25    32
 11   2   54    30    55    23
 12   2   24    23    55    31
 13   2   31    56    25    24
 14   2   31    48    39    56
 15   2   55    38    48    31
 16   2   30    47    38    55
 17   2   54    37    47    30
 18   2   29    46    37    54
 19   2   28    53    44    35
 20   2   27    34    44    53
 21   2   26    43    34    27
 22   2   32    33    43    26
 23   2   50    41    42    25
 24   2   25    42    33    32
 25   2   52    51    42    41
 26   2   41    18    19    52
 27   2   52    19    20    51
 28   2   51    20     1    42
 29   2   50    17    18    41
 30   2   40     4    17    50
 31   2   49    16     4    40
 32   2   39    15    16    49
 33   2   48    14    15    39
 34   2   38    13    14    48
 35   2   47     3    13    38
 36   2   37    12     3    47
 37   2   46    11    12    37
 38   2   36    10    11    46
 39   2   45     9    10    36
 40   2   35     2     9    45
 41   2   44     8     2    35
 42   2   34     7     8    44
 43   2   43     6     7    34
 44   2   33     5     6    43
 45   2   42     1     5    33
 46   1   76    84   110    77
 47   1   77    89    75    76
 48   1   79   113    89    77
 49   1  110   112    79    77
 50   1   87    86    90    78
 51   1   87    78    65    66
 52   1   78    64    58    65
 53   1   78    90    63    64
 54   1   88    82   103    93
 55   1   82    83   104   103
 56   1  107   108    18    17
 57   1  108   100    19    18
 58   1  105   113    79   106
 59   1  106    79   112    97
 60   1  113    95    94    89
 61   1  105    96    95   113
 62   1  112   111    98    97
 63   1  112   110   109   111
 64   1   94   104    83    89
 65   1  109    99    98   111
 66   1   84    85   109   110
 67   1   85   107    99   109
 68   1   90   108   107    85
 69   1   90    86   100   108
 70   1   93    92    81    88
 71   1   92    91   102    81
 72   1  102   101    80    81
 73   1  101    67    68    80
 74   1   87    66    67   101
 75   1   85    62    63    90
 76   1   88    71    72    82
 77   1   81    70    71    88
 78   1   83    74    75    89
 79   1   76    57    61    84
 80   1   68    59    69    80
 81   1   72    60    73    82
 82   1   82    73    74    83
 83   1   80    69    70    81
 84   1   84    61    62    85
 85   1   91     7     6   102
 86   1   92     8     7    91
 87   1   93     2     8    92
 88   1  102     6     5   101
 89   1  103     9     2    93
 90   1  104    10     9   103
 91   1   94    11    10   104
 92   1   95    12    11    94
 93   1   96     3    12    95
 94   1  105    13     3    96
 95   1  106    14    13   105
 96   1   97    15    14   106
 97   1   98    16    15    97
 98   1   99     4    16    98
 99   1  107    17     4    99
100   1   86    20    19   100
101   1   87     1    20    86
102   1  101     5     1    87
2
1,1
3.45  3.45  0.35  1.278  1.0  0.0  0.0
2,1
73.1  73.1  0.22  29.959  1.0  0.0  0.0
20  0  0  0   
57    1    1    0     0
58    1    1    0     0
59    1    1    0     0
60    1    1    0     0
61    0    1    0     0
62    0    1    0     0
63    0    1    0     0
64    0    1    0     0
65    1    0    0     0
66    1    0    0     0
67    1    0    0     0
68    1    0    0     0
69    0    1    0     0
70    0    1    0     0
71    0    1    0     0
72    0    1    0     0
73    1    0    0     0
74    1    0    0     0
75    1    0    0     0
76    1    0    0     0

  

爷且无辜

首先你要搞明白这个程序解决的是什么问题：静力？动力？，然后搞明白是平面问题还是空间问题？平面问题节点坐标可以只有两列，然后看看这个程序用的是什么单元，单元类型搞明白。如果初次学习FEM建议自己编一个小程序，比如很基本的平面三角形三节点单元求静力问题。不建议上来就读别人的程序。

珊瑚虫

shmily1231988 发表于 2014-11-19 16:05
[mw_shl_code=fortran,true]plane strain for TRI 0.41
quad4
113

plane strain for TRI 0.41
quad4
113
   1   0  0.678412413137      0.500000000000   
   2   0  0.500000000000      0.678412413137   
   3   0  0.321587586863      0.500000000000   
   4   0  0.500000000000      0.321587586863   
   5   0  0.670005307592      0.554121941501   
   6   0  0.645576306020      0.603143241596   
   7   0  0.606009913237      0.643502221087   
   8   0  0.555805069669      0.669460270744   
   9   0  0.445878058499      0.670005307592   
  10   0  0.396856758404      0.645576306020   
  11   0  0.356497778913      0.606009913237   
  12   0  0.330539729256      0.555805069669   
  13   0  0.329994692408      0.445878058499   
  14   0  0.354423693980      0.396856758404   
  15   0  0.393990086763      0.356497778913   
  16   0  0.444194930331      0.330539729256   
  17   0  0.554121941501      0.329994692408   
  18   0  0.603143241596      0.354423693980   
  19   0  0.643502221087      0.393990086763   
  20   0  0.669460270744      0.444194930331   
  21   0  0.476725030056      0.543019367480   
  22   0  0.513790554854      0.521219916119   
  23   0  0.453359604385      0.511461925367   
  24   0  0.489858144548      0.484942247136   
  25   0  0.538867282079      0.455123506627   
  26   0  0.559325813685      0.548095025781   
  27   0  0.535597207408      0.565848339208   
  28   0  0.465384710877      0.592757141494   
  29   0  0.438456950694      0.561953830320   
  30   0  0.401067378686      0.507937411315   
  31   0  0.453476897176      0.457376604743   
  32   0  0.554336416882      0.507367871062   
  33   0  0.600762143691      0.518769729388   
  34   0  0.559360223726      0.592943435486   
  35   0  0.487803436975      0.627531588946   
  36   0  0.431082043803      0.611270219314   
  37   0  0.383819269235      0.542405877584   
  38   0  0.385191886472      0.467695704391   
  39   0  0.432465187463      0.400988216053   
  40   0  0.497701418988      0.382914296090   
  41   0  0.591613904625      0.404750128198   
  42   0  0.609586721790      0.471983784272   
  43   0  0.587346673599      0.561656353887   
  44   0  0.524392150154      0.613997358209   
  45   0  0.455119694805      0.636312262478   
  46   0  0.405501488131      0.579248595619   
  47   0  0.374861451911      0.504404437222   
  48   0  0.407057512507      0.432534358878   
  49   0  0.457845685577      0.373138707292   
  50   0  0.545819254379      0.399662651992   
  51   0  0.635114617340      0.441555171449   
  52   0  0.624778756084      0.412289672824   
  53   0  0.500069166892      0.578503866003   
  54   0  0.411049288406      0.531331050730   
  55   0  0.415198764344      0.488587372474   
  56   0  0.482755536212      0.425860603651   
  57   3   0.00000000000       0.00000000000   
  58   3   1.00000000000       0.00000000000   
  59   3   1.00000000000       1.00000000000   
  60   3   0.00000000000       1.00000000000   
  61   2  0.212129992247       0.00000000000   
  62   2  0.383285904260       0.00000000000   
  63   2  0.575740015507       0.00000000000   
  64   2  0.787870007753       0.00000000000   
  65   1   1.00000000000      0.212129992247   
  66   1   1.00000000000      0.383285904260   
  67   1   1.00000000000      0.575740015507   
  68   1   1.00000000000      0.787870007753   
  69   2  0.787870007753       1.00000000000   
  70   2  0.616714095740       1.00000000000   
  71   2  0.424259984493       1.00000000000   
  72   2  0.212129992247       1.00000000000   
  73   1   0.00000000000      0.787870007753   
  74   1   0.00000000000      0.616714095740   
  75   1   0.00000000000      0.424259984493   
  76   1   0.00000000000      0.212129992247   
  77   0  0.125571160599      0.314124590867   
  78   0  0.825647614980      0.173373548552   
  79   0  0.219161802387      0.366385160104   
  80   0  0.826751637609      0.817904007787   
  81   0  0.630098291862      0.847335787503   
  82   0  0.208140386618      0.825218712790   
  83   0  0.143677948862      0.647811800867   
  84   0  0.194821144778      0.126355911515   
  85   0  0.417293789765      0.130275084140   
  86   0  0.725623257574      0.333338821405   
  87   0  0.843391155931      0.374566502743   
  88   0  0.421280858092      0.875294335300   
  89   0  0.124574757403      0.479864901533   
  90   0  0.639751438171      0.157516914797   
  91   0  0.620401331670      0.680771321871   
  92   0  0.567835523363      0.736597782176   
  93   0  0.461555148954      0.754780311019   
  94   0  0.280359130151      0.586773840539   
  95   0  0.287186867231      0.530736405085   
  96   0  0.302491585306      0.497523545978   
  97   0  0.357451219477      0.344824269526   
  98   0  0.402798296829      0.311629190989   
  99   0  0.460336049619      0.284411497270   
100   0  0.662524086114      0.341728445740   
101   0  0.796667726657      0.640430065474   
102   0  0.681182034971      0.687756097944   
103   0  0.361641073801      0.721875308426   
104   0  0.291409750011      0.652101017499   
105   0  0.284981363614      0.445455648274   
106   0  0.317438608624      0.388382610606   
107   0  0.515951744122      0.258658749518   
108   0  0.602417146410      0.270785237388   
109   0  0.346354310008      0.218706921237   
110   0  0.216550550693      0.228800269805   
111   0  0.333433833722      0.278057947224   
112   0  0.267204652327      0.309192815458   
113   0  0.217196172584      0.449224813660
102
  1   2   40    56    39    49
  2   2   25    56    40    50
  3   2   29    28    36    46
  4   2   36    28    35    45
  5   2   23    24    22    21
  6   2   21    29    54    23
  7   2   29    21    53    28
  8   2   27    53    21    22
  9   2   32    26    27    22
10   2   22    24    25    32
11   2   54    30    55    23
12   2   24    23    55    31
13   2   31    56    25    24
14   2   31    48    39    56
15   2   55    38    48    31
16   2   30    47    38    55
17   2   54    37    47    30
18   2   29    46    37    54
19   2   28    53    44    35
20   2   27    34    44    53
21   2   26    43    34    27
22   2   32    33    43    26
23   2   50    41    42    25
24   2   25    42    33    32
25   2   52    51    42    41
26   2   41    18    19    52
27   2   52    19    20    51
28   2   51    20     1    42
29   2   50    17    18    41
30   2   40     4    17    50
31   2   49    16     4    40
32   2   39    15    16    49
33   2   48    14    15    39
34   2   38    13    14    48
35   2   47     3    13    38
36   2   37    12     3    47
37   2   46    11    12    37
38   2   36    10    11    46
39   2   45     9    10    36
40   2   35     2     9    45
41   2   44     8     2    35
42   2   34     7     8    44
43   2   43     6     7    34
44   2   33     5     6    43
45   2   42     1     5    33
46   1   76    84   110    77
47   1   77    89    75    76
48   1   79   113    89    77
49   1  110   112    79    77
50   1   87    86    90    78
51   1   87    78    65    66
52   1   78    64    58    65
53   1   78    90    63    64
54   1   88    82   103    93
55   1   82    83   104   103
56   1  107   108    18    17
57   1  108   100    19    18
58   1  105   113    79   106
59   1  106    79   112    97
60   1  113    95    94    89
61   1  105    96    95   113
62   1  112   111    98    97
63   1  112   110   109   111
64   1   94   104    83    89
65   1  109    99    98   111
66   1   84    85   109   110
67   1   85   107    99   109
68   1   90   108   107    85
69   1   90    86   100   108
70   1   93    92    81    88
71   1   92    91   102    81
72   1  102   101    80    81
73   1  101    67    68    80
74   1   87    66    67   101
75   1   85    62    63    90
76   1   88    71    72    82
77   1   81    70    71    88
78   1   83    74    75    89
79   1   76    57    61    84
80   1   68    59    69    80
81   1   72    60    73    82
82   1   82    73    74    83
83   1   80    69    70    81
84   1   84    61    62    85
85   1   91     7     6   102
86   1   92     8     7    91
87   1   93     2     8    92
88   1  102     6     5   101
89   1  103     9     2    93
90   1  104    10     9   103
91   1   94    11    10   104
92   1   95    12    11    94
93   1   96     3    12    95
94   1  105    13     3    96
95   1  106    14    13   105
96   1   97    15    14   106
97   1   98    16    15    97
98   1   99     4    16    98
99   1  107    17     4    99
100   1   86    20    19   100
101   1   87     1    20    86
102   1  101     5     1    87
2
1,1
3.45  3.45  0.35  1.278  1.0  0.0  0.0
2,1
73.1  73.1  0.22  29.959  1.0  0.0  0.0
20  0  0  0 1  
57    1    1    0     0
58    1    1    0     0
59    1    1    0     0
60    1    1    0     0
61    0    1    0     0
62    0    1    0     0
63    0    1    0     0
64    0    1    0     0
65    1    0    0     0
66    1    0    0     0
67    1    0    0     0
68    1    0    0     0
69    0    1    0     0
70    0    1    0     0
71    0    1    0     0
72    0    1    0     0
73    1    0    0     0
74    1    0    0     0
75    1    0    0     0
76    1    0    0     0
20     2 4 7 2  2  0

我在上述粗体部分添加了一些数据，根据你的程序，这个地方应该有这些数据，因为没有时间读你的程序，所以我不确定这个数据具体的是什么意义（大概一个是自重的方向，最后一行是荷载增量控制参数）现在程序可以走到计算刚度矩阵了，因为报错说雅克比行列式为负数，我查了一下坐标星系的第二列不是约束信息，你这样写可能是导致错误的原因。约束信息是紧紧跟着荷载信息写的， 你需要重新读一下程序看看坐标怎么写，看程序的方法是 查代码中所有的 read (1    需要注意的是你这个程序在incontr 子程序中也读入了1.dat 的内容 就是最后一行的荷载控制参数。具体取多少你自己查。ps  约束信息直接加到节点坐标信息里面似乎不妥，你再查查是不是第二列就是那个意思。
  

珊瑚虫

爷且无辜 发表于 2014-11-19 16:40
首先你要搞明白这个程序解决的是什么问题：静力？动力？，然后搞明白是平面问题还是空间问题？平面问题节点 ...

同意你的观点

shmily1231988

爷且无辜 发表于 2014-11-19 16:40
首先你要搞明白这个程序解决的是什么问题：静力？动力？，然后搞明白是平面问题还是空间问题？平面问题节点 ...

程序解决的是静力  平面问题，程序中设定的输出结果包括变形后的节点坐标，节点上的力，节点位移，模型的平均位移，高斯点上的柯西应力，单元的平均应力，模型的平均应力。模型的各有刚度系数