还是刚塑性流动及欧拉法的问题！请leftzmz及高手指教！！！

firmwave

我做的还是轧制过程的仿真！

最近我做了很多尝试，希望marc能采用欧拉法进行计算，可是好象都没有成功，一旦marc开始算了，就发现他依然采用的是拉格郎日法！
我最近的方法是改变输入文件，在输入文件中加上“r-p flow,1,”，dat文件的开始几行如下：

title               job1
$....MARC input file produced by MSC.Marc Mentat 2005 (32bit)
$...................................
$....input file using extended precision
extended
$...................................
sizing                           1000000     55548    114273         0
elements                    80
version                     11
processor                    1         1         1
$no list
r-p flow,1,    （这里是定义的稳态分析）
couple                       0
heat                         0         0         1         2         1         1
all points
no echo                      1         2         3
films                     1012
setname                   1189
end
$...................
solver
         0         1         0         0         0         0         0         0         0         0         0         0
optimize                     9

按理由，我采用了刚塑性流动的稳态分析，marc应当采用欧拉法进行计算，可是到了正式算的时候，我发觉他依然不是用的欧拉法，很奇怪的是，在output文件中，他清楚的标明了我采用了刚塑性流动的非稳态分析（见下面的out文件），也就是说，marc依然使用拉格朗日法进行计算。

*************************************************
                    *************************************************

                     program sizing and options requested as follows

                     element type requested*************************        80
                     element type requested*************************        39
                     number of elements in mesh*********************     55548
                     number of nodes in mesh************************    114273
                     max number of elements in any dist load list***         0
                     maximum number of point loads******************         0
                     large displacement analysis flagged************
                     thermal stress analysis flagged****************
                     load correction flagged or set*****************
                     number of lists of distributed loads***********         3
                     values stored at all integration points********
                     tape no.for input of coordinates + connectivity         5
                     no.of different materials    2 max.no of slopes        20
                     heat transfer analysis,  extrapolation flag, **         1
                     latent heats per element changed to ***********         1
                     gradient of scaler field printed***************
                     include convective terms in thermal equilibrium
                     number of elements with film coefficient b.c.**      2024
                     maximum number of distributed flux lists*******         3
                     rigid-plastic non-steady state flow solution******
                                                （怎么变成了非稳态？？？？？？？？）
                     number of points on shell section *************        11
                     geometry updated after each load step**********
                     new style input format will be used************
                     maximum number of set names is*****************      1189
                     coupled thermal-mechanical analysis flagged****
                     number of processors used *********************         1
                     extended precision input is used **************
                     MSC.Marc input version ************************        11
                     suppress echo of list items *******************
                     suppress echo of bc summary *******************
                     suppress echo of nurbs data *******************

                     end of parameters and sizing
                    *************************************************
                    *************************************************

我做过了e7x1c.dat这个例子，成功了，然后我稍微改变模型，定义了接触，这时我发觉它也不再采用欧拉法了，因为他也是rigid-plastic non-steady state flow solution******，晕！！！

我就奇怪了，为什么我明明要求的是刚塑性稳态分析，结果却成了非稳态分析呢？难道用欧拉法计算时不能定义接触吗？如果这样，那marc 的帮助文档volume A的地119页上面的例子不是哄人的吗？

急等高手指教！

leftzmz

按个人的理解在marc中用欧拉法求解塑性变形中一下问题应该不是它的强项，那个例子我看过了，你注意没有，它的变形量是很小的。对于大变形问题可以肯定的是marc中用欧拉法是很难实现的。我也注意到例子中只采用定义边界条件的方法来定义接触条件，但是说用欧拉法计算时不能定义接触似乎是很难让人接受的，deform中不就是采用定义接触的方法吗？可以看出你在这方面是花了一些功夫，期待你对marc的内涵有更深的理解。

firmwave

谢谢leftzmz！！！
这点我也注意到了，并且，这个题目很不完整，我再延长分析时间，加几个载荷步时，他就不能算了！
我这几天还会继续调试，到时有问题或者有什么进展的话，再向你请教！