怪异的现象

独孤

我把物体分为66个材料，16个边界，Pre文件如下：
3dxyz 1 1 8 10
para 1 2 82 u
#
element y
para
aec8 8 rc ek eq wb tb hz
alq4 4 eb
#
。。。。。。
39 para aec8
2703000.0     2.280     0.0     0     37.0     0.0
40 para aec8
3958500.0     0.419     705.6     0     37.0     0.0
41 para aec8
2133500.0     0.161     3.4     0     37.0     0.0
42 para aec8
3770000.0     0.209     1010.0     0     37.0     0.0
43 para aec8
2703000.0     2.280     0.0     0     37.0     0.0
44 para aec8
3958500.0     0.419     705.6     0     37.0     0.0
45 para aec8
2133500.0     0.161     3.4     0     37.0     0.0
46 para aec8
3770000.0     0.209     1010.0     0     37.0     0.0
47 para aec8
2703000.0     2.280     0.0     0     37.0     0.0
48 para aec8
3958500.0     0.419     705.6     0     37.0     0.0
49 para aec8
2133500.0     0.161     3.4     0     37.0     0.0
50 para aec8
3770000.0     0.209     1010.0     0     37.0     0.0
51 para aec8
2703000.0     2.280     0.0     0     37.0     0.0
52 para aec8
3958500.0     0.419     705.6     0     37.0     0.0
53 para aec8
2133500.0     0.161     3.4     0     37.0     0.0
54 para aec8
3770000.0     0.209     1010.0     0     37.0     0.0
55 para aec8
2703000.0     2.280     0.0     0     37.0     0.0
56 para aec8
3958500.0     0.419     705.6     0     37.0     0.0
57 para aec8
2133500.0     0.161     3.4     0     37.0     0.0
58 para aec8
3770000.0     0.209     1010.0     0     37.0     0.0
59 para aec8
2703000.0     2.280     0.0     0     37.0     0.0
60 para aec8
3958500.0     0.419     705.6     0     37.0     0.0
。。。。。。。。
但是在Fgidpre文件中41-56号的数组有问题
。。。。。
      rmdata(39,1) = 2703000.0
      rmdata(39,2) =     2.280
      rmdata(39,3) =     0.0
      rmdata(39,4) =     0
      rmdata(39,5) =     37.0
      rmdata(39,6) =     0.0
      rmdata(40,1) = 3958500.0
      rmdata(40,2) =     0.419
      rmdata(40,3) =     705.6
      rmdata(40,4) =     0
      rmdata(40,5) =     37.0
      rmdata(40,6) =     0.0
      rmdata(41,1) = 34
      rmdata(42,1) = 34
      rmdata(43,1) = 34
      rmdata(44,1) = 34
      rmdata(45,1) = 34
      rmdata(46,1) = 34
      rmdata(47,1) = 34
      rmdata(48,1) = 34
      rmdata(49,1) = 34
      rmdata(50,1) = 34
      rmdata(51,1) = 34
      rmdata(52,1) = 34
      rmdata(53,1) = 34
      rmdata(54,1) = 34
      rmdata(55,1) = 34
      rmdata(56,1) = 34
      rmdata(57,1) = 2133500.0
      rmdata(57,2) =     0.161
      rmdata(57,3) =     3.4
      rmdata(57,4) =     0
      rmdata(57,5) =     37.0
      rmdata(57,6) =     0.0
。。。。。。
其中34是边界的热流量
82 para alq4
34;
不知道为什么出现这样的问题？

独孤

问题在于41-56号的数组有问题。这是为什么呢？
如果我手动修改fgidpre中的数组，是否能够解决问题？

zhoulon

你的pre文件的第二行贴写错误了，不是para 1 2 82 u ，而是para 1 2 6 u，你的单元最大材料参数是6

独孤

我已经修改过来了，不过出现的问题还是一模一样。我把pre文件上传了上来，大家帮忙看下？

我手工修改了fgidpre里面错误的那些数组，并且编译，能够正确运行。

zhoulon

我现在没有时间来看，你可以具体看一下fgidpre里面的数组的开辟空间大小，自己仔细读一下程序应该就可以解决了。

独孤

好的，谢谢提供思路！