请教蓝博士及其他高手关于EXCITE到ANSYS数据恢复的问题

mingzhangjie

蓝博士,该怎么把EXCITE得到的曲轴数据恢复到ANSYS中啊,
   我在EXCITE计算中建的曲轴模型很简化,没有圆角.所以应力恢复时我另建了半拐曲轴较精确的模型,并用软梁建立其他拐的模型,并保证主自由度点号码与EXCITE计算中建的曲轴模型一致.

  我看到C:\AVL\EXCITE\v6.1\examples\Primer\excite\jh16v240_1.Common\WORK@目录下
曲轴的结果文件中有: [email=13@.INP.NEWM]13@.INP.NEWM[/email]     [email=13@M.INP.TIME]13@M.INP.TIME[/email], 13@MFF.INP.TIME, 13@S.INP.TIME, 13@T.INP.TIME, 13@X2OA.INP.TIME, [email=13@.NAM.DOFT]13@.NAM.DOFT[/email], [email=13@.NAM.GEOM]13@.NAM.GEOM[/email], [email=13@.NAM.MFF]13@.NAM.MFF[/email], [email=13@.NAM.MMAT]13@.NAM.MMAT[/email], [email=13@.NAM.MMAT]13@.NAM.MMAT[/email] (13@是曲轴)

   做曲轴应力恢复时要用到哪个文件啊, 您以前给的那个可执行文件ok_excite2ans2_disp_new.exe怎么用啊,是不是先按照它的提示一步步指定相应文件的地址目录,然后在ANSYS中通过菜单AVL-EXCITE Recovery把上面执行文件生成的结果文件指定给ANSYS就可以啊?

麻烦蓝博士指点一下,谢谢!

[ 本帖最后由 mingzhangjie 于 2007-6-24 10:41 编辑 ]

DrLan

我在EXCITE计算中建的曲轴模型很简化,没有圆角.所以应力恢复时我另建了半拐曲轴较精确的模型,并用软梁建立其他拐的模型,并保证主自由度点号码与EXCITE计算中建的曲轴模型一致.
蓝：很好。

   做曲轴应力恢复时要用到哪个文件啊,
蓝：Utility中的data recovery选择displacement，产生INP4。

您以前给的那个可执行文件ok_excite2ans2_disp_new.exe怎么用啊,是不是先按照它的提示一步步指定相应文件的地址目录,
蓝：将INP4和DOFT文件和该程序拷贝在一个目录下运行，产生ansys可执行命令流。照提示去做。

然后在ANSYS中通过菜单AVL-EXCITE Recovery把上面执行文件生成的结果文件指定给ANSYS就可以啊?
蓝：把网格模型调入后，再调入上面产生的命令流，即可进行瞬态应力计算。

mingzhangjie

蓝博士,在Utility中进行data recovery,出错了,提示如下:


Calling program "EXCITE" and generating EXCITE result file...

                    ==========================================================
                    =                                                        =
                    =                     Excite 6.1 (30104)    Jul 30 2004  =
                    =                                                        =
                    ==========================================================
                    =  OS   Company   :                           Microsoft  =
                    =  OS   Name      :                          Windows XP  =
                    =  OS   Version   :           5.1.2600 (Service Pack 2)  =
                    =  CPU  Model     :              Processor 586 Level 15  =
                    =  Number of CPU`s:                                   1  =
                    =  Host Name      :                                 MZJ  =
                    =  User Name      :                              jackey  =
                    ==========================================================
$ Name of the instruction file:  $ Name of the output file < Return = "00_res_inp4.out" >:
$ EXCITE-6.1 - start
$ General information:
$ Date: 27-Jun-2007  - Time: 10:26:23
$ Available memory -> 700.0000 MB
$ Instruction file -> 00_res_inp4.exc
% RESCTR  TIME &
%  RPM  1000 &
%  STRK  4 &
%  INP 1 &
%  dTre 1 &
%  BEG  0.120000048 &
%  END  0.240000096 &
%  IRES 1
% SETIO BCD
% BODY RES  13@ NAST
%  : DISP DISP ::
% Crankshaft.INP4 |
$ BRES   Data block/file information for NASTRAN-OUT4 output:
$ BRES   Data block -> DISP     : DISP     :: file -> Crankshaft.INP4
$ BRES3  Processing result data for body: 13@
$ BRES3  -> Read/transform data:            [ok]
$ BRES3  -> Interpolation:                  [ok]
$ BRES3  -> Write "DISP    ":               [ok]
% STOP
$ Used CPU time : 2.750000 s
$ EXCITE-6.1 - normal termination!
EXCITE result file >> Crankshaft.INP4 << has been created.
Calling: /cygdrive/c/AVL/EXCITE/v6.1/bin/bin.ia32.p4-unknown-winnt/exciteRecover.exe DataRecover.inp 720 1440
                    ==========================================================
                    =                                                        =
                    =             EXCITE_Recover 1.0 (30107)    Aug  3 2004  =
                    =                                                        =
                    ==========================================================
                    =  OS   Company   :                           Microsoft  =
                    =  OS   Name      :                          Windows XP  =
                    =  OS   Version   :           5.1.2600 (Service Pack 2)  =
                    =  CPU  Model     :              Processor 586 Level 15  =
                    =  Number of CPU`s:                                   1  =
                    =  Host Name      :                                 MZJ  =
                    =  User Name      :                              jackey  =
                    ==========================================================
************************************************************************************
*                           Problem Specification Report                           *
************************************************************************************
    using the following input files
        input file containing EXCITE result matrix:
            Crankshaft.INP4
        input file containing transformation matrix X2OA:
            ../../fem/Crankshaft/Crankshaft_X2OA.OUT4
        input file containing the geometry information (.meg file):
            ../../fem/Crankshaft/Crankshaft_GEOM.meg
    for data recovery to the whole structure in angle domain
    and produce the output files
        output geometry file (to be used as input file for IMPRESS):
            Crankshaft_Disp_Angle.meg
        output result file (.mer file, to be used as input for IMPRESS):
            Crankshaft_Disp_Angle.mer
************************************************************************************
   MESSAGE: start reading all collumns of EXCITE result matrix
            stored in file Crankshaft.INP4 ...
   MESSAGE: finished reading EXCITE result matrix in 1 seconds.
   ERROR: invalid OUPUT4 file ../../fem/Crankshaft/Crankshaft_X2OA.OUT4
          use option -b for binary formatted files and -a for ASCII files
          Program terminated abnormally
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> END OF UTILITY <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<


怎么进行了2个program啊,最后那个错误提示是什么意思啊?是Crankshaft_X2OA.OUT4的问题吗?

另外,create results 的时候又出错了，不知怎么回事，在曲轴１３＠体的结果日志上出现了一个很大的节点号2147483647（日志文件摘录如下），那个节点根本不存在啊，奇怪，我师兄做的模型和我的节点号不一样，竟然也出现这一问题，而且同样也是多了那个节点2147483647！
．．．
$ BRES1  -> Check node 10280 /DoF 3:        [results already available]
$ BRES1  Processing result data for body: Crankshaft
$ BRES1  -> Check node 10280 /DoF 4:        [results already available]
$ BRES1  Processing result data for body: Crankshaft
$ BRES1  -> Check node 10280 /DoF 5:        [results already available]
$ BRES1  Processing result data for body: Crankshaft
$ BRES1  -> Check node 10280 /DoF 6:        [results already available]
% BODY RES 13@ REL 10584 1 2 3 4 5 6
$ BRES1  Processing result data for body: Crankshaft
$ BRES1  -> Check node 10584 /DoF 1:        [results already available]
$ BRES1  Processing result data for body: Crankshaft
$ BRES1  -> Check node 10584 /DoF 2:        [results already available]
$ BRES1  Processing result data for body: Crankshaft
$ BRES1  -> Check node 10584 /DoF 3:        [results already available]
$ BRES1  Processing result data for body: Crankshaft
$ BRES1  -> Check node 10584 /DoF 4:        [results already available]
$ BRES1  Processing result data for body: Crankshaft
$ BRES1  -> Check node 10584 /DoF 5:        [results already available]
$ BRES1  Processing result data for body: Crankshaft
$ BRES1  -> Check node 10584 /DoF 6:        [results already available]
% BODY RES 13@ REL 10585 1 2 3 4 5 6
$ BRES1  Processing result data for body: Crankshaft
$ BRES1  -> Check node 10585 /DoF 1:        [results already available]
$ BRES1  Processing result data for body: Crankshaft
$ BRES1  -> Check node 10585 /DoF 2:        [results already available]
$ BRES1  Processing result data for body: Crankshaft
$ BRES1  -> Check node 10585 /DoF 3:        [results already available]
$ BRES1  Processing result data for body: Crankshaft
$ BRES1  -> Check node 10585 /DoF 4:        [results already available]
$ BRES1  Processing result data for body: Crankshaft
$ BRES1  -> Check node 10585 /DoF 5:        [results already available]
$ BRES1  Processing result data for body: Crankshaft
$ BRES1  -> Check node 10585 /DoF 6:        [results already available]
% BODY RES 13@ ABS 2147483647 1 2 3 4 5 6
$ BRES1  Processing result data for body: Crankshaft

[ 本帖最后由 mingzhangjie 于 2007-6-27 15:28 编辑 ]

tomas

“另外,create results 的时候又出错了，不知怎么回事，在曲轴１３＠体的结果日志上出现了一个很大的节点号2147483647（日志文件摘录如下），那个节点根本不存在啊，奇怪，我师兄做的模型和我的节点号不一样，竟然也出现这一问题，而且同样也是多了那个节点2147483647！”
我再补充两句：我们两个建的FEM模型不同，主节点编号也不同，EXCITE参数也不完全一致，但creat results后都有这个节点2147483647。奇怪！

还有我这几天EXCITE计算的结果中，用show results方式，发电机转子和它相关的REVO，AXBE，ROTX等joint都没有结果输出。我查看了creat results的记录，发现最后两行显示的是
% BODY RES 59@ ABS 2147483647 1 2 3 4 5 6
$ BRES1  Processing result data for body: Crankshaft
59@对应的是曲轴。我想creat results可能没有正常结束，要正常结束的话，会显示，normal termination或end of utility。因为模型中没有2147483647这个节点，创建结果时找不到对应数据，程序就停顿在这了。恰好转子和它相关的REVO，AXBE，ROTX等joint的编号都排在59@的后面，也都没有正确的生成结果，所以我在show results时就看不到这些body和joint的数据。
试着先删去了EXCITE模型中的曲轴，再重新创建曲轴，填好参数，这样曲轴的编号就排在最后了，再creat results，最后两行还是
% BODY RES 68@ ABS 2147483647 1 2 3 4 5 6
$ BRES1  Processing result data for body: Crankshaft
但我在show results时就能够看到转子和相关的REVO，AXBE，ROTX等joint的数据了。
我想这个节点2147483647应该是有问题的。

mingzhangjie

麻烦蓝博士了,呵呵,不知道以前有没有哪位碰到类似问题?

DrLan

上面的一个结果正常结束的。
$ EXCITE-6.1 - normal termination!
EXCITE result file >> Crankshaft.INP4 << has been created.
下面这个结果是用EXCITE进行恢复，跟你没有关系。

% BODY RES 13@ ABS 2147483647 1 2 3 4 5 6
你在13体单元内申请了ABS 2147483647 1 2 3 4 5 6
2147483647是你自己填的，删掉这个申请。

tomas

我试着在曲轴的results中取消Global results和Node results，再creat results、show results，还是在曲轴的absolute、relative目录下有这个2147483647这个节点。呵呵，不知道从哪里出来的。

mingzhangjie

蓝博士,我把产生的命令流文件cs.dat读进ANSYS单拐有限元模型后,进入ANSYS菜单
AVL-EXCITE RECOVERY,指定GID文件,然后进行缩减(读入cs.d后,没必要另外进行瞬态动力分析的相关设置吧?好象读入后就自动设置好了),可是又出错了:
*** ERROR ***                           CP=      4.562   TIME= 14:58:16
File file.esav not found.  Previous solution (required for restarting)  
cannot be recovered


*** ERROR ***                           CP=      4.594   TIME= 14:58:17
Error in restart, RDB file does not exist.                              
  The ANTYPE command is ignored.                                         
/COM,ANSYS RELEASE  8.0    UP20030930       15:01:27    06/28/2007


*** ERROR ***                           CP=   1736.938   TIME= 15:41:20
Present time 0 is less than or equal to the previous time in a         
transient analysis.  


我查了一下，file.esav好象是单元数据文件，ANSYS中瞬态动应力分析怎么还要用到原来单元的数据啊,要是那样的话,我单拐模型单元明显和原曲轴整体模型单元明显不一致,这就没法进行恢复了啊?

GID文件是所有主节点在离散时间点上的各个方向的位移以及载荷等量,可是奇怪的是怎么有的GID文件是0~1440deg,而有的则是720~1440deg啊？

DrLan

恢复模型的节点号和自由度必须与原模型一致。

mingzhangjie

2个模型主自由度节点号和位置都是完全一致的啊,每个主节点都是6个自由度.

通过EXCITE结果得知5#主轴颈载荷最大,所以恢复模型建立的是5#主轴颈以及两边半拐,这样可以吗?
下图是恢复模型图以及EXCITE-RECOVERY对话框.

呵呵,请蓝博士指点一下,谢谢.

[ 本帖最后由 mingzhangjie 于 2007-6-29 11:45 编辑 ]

tomas

关于2147483647这个节点，找到一点线索，这个数等于2的31次方减1，是32位电脑里面最大的整数。不知道这是不是EXCITE程序自动加上的一个节点号，呵呵。

DrLan

通过EXCITE结果得知5#主轴颈载荷最大,所以恢复模型建立的是5#主轴颈以及两边半拐,这样可以吗?
蓝：很好，等你贴出动应力结果啊。有动画最好了。

mingzhangjie

　我用上述模型做了2个步长的应力恢复计算,一个步长用时30多分钟! 这样看来把一个循环算下来,得算几天几夜啊!
  　１．基于以上问题,我进一步缩小模型,只建半拐,并且去掉平衡重,那样一个步长大概9分钟,0~720,共145个载荷步,这样一个晚上多的时间就算完了.可是看到结果不理想就没有保存原结果． ．这样做模型简化有点大，不知是否可行？
  　２．可是得到的结果应力值偏小,vov mises应力值最大才50Mpa，另外，得到的特定节点应力时程曲线不太理想，波动特别大．是什么原因呢？而且，发现在油孔和轴颈蜘蛛腿连接节点附近有较大应力集中．

     ３．我觉得应力集中可能是由于网格太大(18mm加圆角细化),蜘蛛腿数目不够．应力值偏小应该是EXCITE参数设置的问题吧？我算得1~9#主轴颈载荷相差比较大,最大的有80KN,最小的才10KN!
    ４．恢复用的半拐梁模型中软梁参数为E=260Mpa,密度=0,实常数:A=250mm^2，
Ｉ1=5000,Ｉ2=5000,Ｉ3=10000mm^4,其余为０．这个参数是当初建参考点时软梁参数，呵呵．不知道是否合适．
　５．另外,恢复用的半拐梁模型中不要加任何约束吧？

[ 本帖最后由 mingzhangjie 于 2007-7-3 11:08 编辑 ]

DrLan

　我用上述模型做了2个步长的应力恢复计算,一个步长用时30多分钟! 这样看来把一个循环算下来,得算几天几夜啊!
蓝：平衡重等非应力集中区的网格尺寸加大，有望将单元总数将下一半，计算将加快好几倍。

  　１．基于以上问题,我进一步缩小模型,只建半拐,并且去掉平衡重,那样一个步长大概9分钟,0~720,共145个载荷步,这样一个晚上多的时间就算完了.可是看到结果不理想就没有保存原结果． ．这样做模型简化有点大，不知是否可行？
蓝：不要用半拐。另外建议用主轴颈到主轴颈的单拐，不要用曲柄销到曲柄销的单拐（因为这里只有一个spider，边界条件将不够好）。

  　２．可是得到的结果应力值偏小,vov mises应力值最大才50Mpa，另外，得到的特定节点应力时程曲线不太理想，波动特别大．是什么原因呢？而且，发现在油孔和轴颈蜘蛛腿连接节点附近有较大应力集中．
蓝：建议用主轴颈到主轴颈的单拐，主轴颈用整个5层的，建议将最外层的那个spider改大到整个轴颈横截面。两个主轴颈的外层两侧spider都加大。

     ３．我觉得应力集中可能是由于网格太大(18mm加圆角细化),蜘蛛腿数目不够．应力值偏小应该是EXCITE参数设置的问题吧？我算得1~9#主轴颈载荷相差比较大,最大的有80KN,最小的才10KN!
蓝：是分割单拐后造成的，本来拐间还有相互的牵扯关系。这个牵扯关系与spider的定义很有关系。这也是我不建议使用“曲柄销到曲柄销的单拐”的理由。

    ４．恢复用的半拐梁模型中软梁参数为E=260Mpa,密度=0,实常数:A=250mm^2，
Ｉ1=5000,Ｉ2=5000,Ｉ3=10000mm^4,其余为０．这个参数是当初建参考点时软梁参数，呵呵．不知道是否合适．
蓝：把几个I缩小10倍后再试试看是否有明显影响。

　５．另外,恢复用的半拐梁模型中不要加任何约束吧？
蓝：不做半拐。照我上面所述去做。

mingzhangjie

蓝博士,我按您的要求做了主轴颈到主轴颈的单拐.

   我看到有的资料上说恢复计算的曲轴模型要在不重要的地方加3个弹簧约束 ,那么我这个模型还需要加这个约束吗 ?如果要加的话,加的具体位置以及弹簧刚度该怎么把握呢?
   
  麻烦蓝博士指点一下,谢谢!

DrLan

我这里用的是位移边界，来恢复应力。可以不加弹簧。

mingzhangjie

蓝博士,按照您的要求做了单拐-梁模型,分了145个载荷步,每个载荷步间隔5度曲柄转角.
    1.得到的结果如下,感觉应力值还是有点偏小.可能会是什么原因呢？
　　在第14个载荷步(曲柄转角65度)附近时,节点247431有最大应力值,位置在过度圆角处.
　　下图是在第14个载荷步时的应力云图(von mises和 1st principle）

    2.globals参数设置中generate loads中我只选择了一下几项，如图glodals.  这样可以吗？

    3.另外我发现一个问题，在simulation report中列出的大头轴承力以及主轴承力和   impress chart中
　　得到的值不相等，这是怎么回事啊？

[ 本帖最后由 mingzhangjie 于 2007-7-13 21:22 编辑 ]

DrLan

1.得到的结果如下,感觉应力值还是有点偏小.可能会是什么原因呢？
蓝：挺好的啊。你的理由呢？

    2.globals参数设置中generate loads中我只选择了一下几项，如图glodals.  这样可以吗？
蓝：只需产生模型要使用的系统“外”载。

    3.另外我发现一个问题，在simulation report中列出的大头轴承力以及主轴承力和   impress chart中得到的值不相等，这是怎么回事啊？
蓝：一样的。（相对于轴瓦）

mingzhangjie

1。动态应力值一般比静态应力计算得出的结果要大吧，我看到一般机型动态最大应力都有200MPA左右。
     另外，对我这16V240，汽缸爆发压力比较大，得到的动态应力值应该要比一般机型大些吧？
2。哦，那么我只做动态应力计算，使用的外载有作用在连杆小头的汽缸压力，作用在转子上的输出扭矩，产生这几个就可以了吧？
3。应该是相等的吧？我是把impress chart中得到的垂直方向，水平方向2个力用平方和然后开方进行合成，然后再与simulation report中的力进行大小比较的，不知道这样可不可以。

[ 本帖最后由 mingzhangjie 于 2007-7-16 09:10 编辑 ]

DrLan

1。动态应力值一般比静态应力计算得出的结果要大吧，我看到一般机型动态最大应力都有200MPA左右。
蓝：要是你的静态模型是正确的话，是你说的。

     另外，对我这16V240，汽缸爆发压力比较大，得到的动态应力值应该要比一般机型大些吧？
蓝：曲拐的个头还大呢。应该不一定就大。

2。哦，那么我只做动态应力计算，使用的外载有作用在连杆小头的汽缸压力，作用在转子上的输出扭矩，产生这几个就可以了吧？
蓝：是你说的这样。象Primer一样的例子，只使用了加在活塞（连杆小头）的活塞力和曲轴输出扭矩。

3。应该是相等的吧？我是把impress chart中得到的垂直方向，水平方向2个力用平方和然后开方进行合成，然后再与simulation report中的力进行大小比较的，不知道这样可不可以。
蓝：图都是一模一样的。径向的两个方向曲线。你要得到最大力，可以象你说的平分和后开方来寻找最大轴承力。