重启动可以改变求解方式吗

清澈

我现在要做一个对橡胶材料的分析，施加的是位移载荷，想研究下橡胶随着位移载荷变化的相关性能。
目前有个问题是，用普通的打开自动时间步只能算到很小的位移载荷就不收敛了。所以想用弧长法来算。
但是从一开始就打开弧长法的话，运算速度实在是太慢了。。
目前的想法是用两个载荷步来算：第一步采用自动时间步长，算到不收敛处的位移载荷下的行为；然后使用弧长法继续算。

问题是：1、位移载荷貌似不能通过两个载荷步来实现呀。。我做的小例子没有不改变求解方式，第二步也不是在第一步的基础上进行。。
1、师兄建议可以试下重启动。我搜索相关帖子，发现有个说“重启动不能任意改变参数，不是说可以随便根据自己的需要来改变参数之后就可以继续分析，重启动中可以改变边界条件，材料参数，网格密度”。。那改成弧长法能实现吗？

或者该问题大家还有其他更好的办法~~
多谢大家！

清澈

我试了试多载荷步的方法。
发现：如果第一步用弧长，第二步用自动时间步长，可以算；
但是第一步用自动时间步长，第二步用弧长就不行。

例子是一个最简单的单轴拉伸橡胶材料的模型，命令流如下：
/PREP7  
et,1,185
!*  
TB,HYPE,1,1,2,NEO   
TBTEMP,0
TBDATA,,0.1007,0.2,,,,  
BLOCK,0,0.5,0,0.2,0,0.2,
esize,0.1   
vmesh,1
FINISH  
/SOL
ANTYPE,0
NLGEOM,1
NSUBST,20,0,0   
OUTRES,ERASE
OUTRES,ALL,1
TIME,1  
nsel,s,loc,x,0  
d,all,all   
nsel,s,loc,x,0.5
d,all,ux,0.05   
nsel,all
/STATUS,SOLU
SOLVE   
NSUBST,100,0,0  
ARCLEN,1,0,0
AUTOTS,-1.0
TIME,0  
DOFSEL,S,UX
DCUM,ADD,1, ,0,
DOFSEL,ALL  
!*  
nsel,s,loc,x,0.5
d,all,ux,0.05   
nsel,all
/STATUS,SOLU
SOLVE

算完之后发现第二步加的0.05位移载荷根本没有加上去，输出窗口上显示的信息是：


大家知道是什么原因吗？
谢谢大家！

swan321

arclen 是一个很好的办法，只是你需要设置核实的limit以及半径，而且不能使用tabular load 和distributed cpus，
你这种情况，建议设置检测点，检查你的不收敛原因，然后选择合适的算法，
要不你可能还是不知道是否能够成功收敛及结果的正确性。

清澈

swan321 发表于 2013-6-28 16:08
arclen 是一个很好的办法，只是你需要设置核实的limit以及半径，而且不能使用tabular load 和distributed c ...

谢谢你的回复！
能详细说明一下“设置检测点”么？确实很想弄清楚为啥不收敛呀~~
现在就是一点点试着往上加位移载荷，那为什么换成弧长法之后位移载荷加不上，是半径没设置合适呀？

swan321

首先看你的output文件和error文件，对于超弹性材料来说，不收敛主要有以下原因：
1 网格发生畸变，一般此原因是最常见的原因。这个需要你对网格划分的经验，（预判网格变形的趋势，然后划分相应的网格，使得网格在变形之后不至于发生严重的畸变）
2 算法选择。材料的算法选择，对于超弹性材料，ansys有Neo-Hookean, Mooney－Rivlin, Polynomial, Yeoh 以及Ogden算法
每种算法针对的应变变形不同，最好选择合适的算法，

至于监测点。在solu>result tracking 中设置

祝你成功。

清澈

swan321 发表于 2013-7-3 10:09
首先看你的output文件和error文件，对于超弹性材料来说，不收敛主要有以下原因：
1 网格发生畸变，一般此原 ...

谢谢你！
目前算的已经能满足要求啦！
不过对于2，针对不同超弹性材料，每个都有针对的算法？？没看见过呀。。。是在solution选项中设置吗？

swan321

在材料参数输入时选择～

wskypm326

这个是有点恼火