独孤 发表于 2009-3-13 18:06:17

怪异的现象

我把物体分为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;
不知道为什么出现这样的问题？

独孤 发表于 2009-3-13 18:09:31

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

zhoulon 发表于 2009-3-14 09:46:16

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

独孤 发表于 2009-3-14 11:51:28

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

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

zhoulon 发表于 2009-3-16 12:10:15

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

独孤 发表于 2009-3-17 10:22:52

好的，谢谢提供思路！
页: [1]
查看完整版本: 怪异的现象