众说纷纭的错误One or more elements have become highly distorted

xiaaiyin

One or more elements have become highly distorted.  Excessive distortion of elements is usually a symptom indicating the need for corrective action elsewhere.  Try incrementing the load more slowly (increase the number of substeps or decrease the time step size).  You may need to improve your mesh to obtain elements with better aspect ratios.  Also consider the behavior of materials, contact pairs, and/or constraint equations.  If this message appears in the first iteration of first substep, be sure to perform element shape checking.

这个错误很多人都在网上发过贴求助，似乎没有什么满意的答案。
基本上是让重新划分网格，子步设置大些，多添些约束，换材料等。
很多人都是在进行了种种尝试后，还是不行，或是需要极大的改变条件与实际不符合才会不弹出这个错误，但那已经没有意义了。
网上的一些帖子，网友问解决了没，我也遇到了这样的问题，可惜作者工作之后就很少上论坛或是不关注这些了。

我的这个课题全是solid45和solid65单元建立的，也是算到time13就出现单元极度扭曲的错误，但是我看图形单元也没出现什么特别变形的地方。
网格全是规整的六面体，问题肯定也不是网格方面。
想起solid65单元支持大位移，我却给关了，打开单元选型后，仍然不解决问题。

最后就聚焦在子步的问题上了。
为了节省机时，采用重启动设置了几个还原点
例如：
rescontrol,,10,last
rescontrol,,11,last
rescontrol,,12,last

然后重启动就直接
finish
/clear
/filename,keti
/solu
antype,,rest,12,last


TIME,13
D,jiazaidian,uz,-30
kbc=0
NSUBST,100
alls
LSWRITE,13
solve
save

TIME,14
D,jiazaidian,uz,33.5
kbc=0
NSUBST,100
alls
LSWRITE,14
solve
save


TIME,15
D,2972,uz,-33.5
kbc=0
NSUBST,100
alls
LSWRITE,15
solve
save
。。。。
然后就可以计算过去了，也没出现这个单元扭曲警告
我的重启动设置就上面一句话，会不会漏掉什么边界条件，而不是在原基础上分析呢，各位大侠看看！

我把子步设置大些，今天晚上算一晚上，连续计算，看看明天的结果和重启动是否一样，不管结果如何，都将如实发帖进展

xiaaiyin

经过昨天晚上的计算，增加子步也无法实现连续计算！
不知道我的那个重启动（设置很简单），能否实现前面的轴力、边界条件都不变而在其基础上继续计算的功能呢？
请高手们帮帮忙，整了一个月了，还是回到了以前的问题上

htbbzzg

不了解你的具体情况，很难发表什么看法。
不过，从你给出的几句命令来看，每一个新的 Time，你都对约束条件进行了修改 (D 命令)，而按照 ANSYS 的要求，如果想在前一次计算的基础上继续计算的话，约束条件是不能修改的。如果修改了约束条件，ANSYS 将作为新的分析，不再考虑前面的计算结果了，这和你的想法是不一致的。

xiaaiyin

加载的是位移荷载，这个也算成改变约束吗？

xiaaiyin

htbbzzg 发表于 2012-5-16 09:03
不了解你的具体情况，很难发表什么看法。
不过，从你给出的几句命令来看，每一个新的 Time，你都对约束条件 ...

我的试验是一个钢管砼柱和两个钢梁通过高强螺栓连接起来的，具体试验过程，希望大侠能移步稍微看下http://forum.simwe.com/thread-1035096-1-1.html这个介绍，D,jiazaidian,uz,-30。这里的jiazaidian是柱顶的耦合点，通过这个点施加低周反复位移荷载，再计算出反力，从而得到滞回曲线。位移边界条件都没改变，这个D是相当于外荷载施加的，这个也算是约束条件改变吗？

xiaaiyin

htbbzzg 发表于 2012-5-16 09:03
不了解你的具体情况，很难发表什么看法。
不过，从你给出的几句命令来看，每一个新的 Time，你都对约束条件 ...

前辈您好,我建立了一个简单的模型，进行对比，发现这个D命令对结果是没有影响的。
将文件设置为test
finish
/clear
/title,connection
/prep7
/view,1,1,1,1                !设置所有窗口的视图方式
/nerr,0
!采用mm-MPA-N单位制度  !在MPA单位下，力为N，长度为mm，弹模为MPA，质量为t（吨），重力加速度!为9815
et,1,solid65     !定义钢材的单元为solid45
!钢梁H型钢腹板（材料1)
MP,EX,1,2.06e5
MP,PRXY,1,0.3
r,45
tb,bkin,1 !采用双线性随动强化，材料1
TBTEMP,0
tbdata,,270,2.06e4 !定义屈服强度和切向模量
!BLOCK, X1, X2, Y1, Y2, Z1, Z2
block,20,-20,50,-50,20,-20,
esizie,10        
mshape,0,3d
mshkey,1
vatt,1,45,1 !材料7，单元类型2
vmesh,all

!!!!!!!!!!!求解!!!!!!!!!!!!!!!!
/solu
allsel
!lsclear,all !清空所有载荷
antype,static !静力分析
nlgeom,on              !打开大变形选项
!lnsrch,on               !打开线性搜索,与自适应下降互斥
neqit,10000              !每一子步中方程的迭代次数限值            
cnvtol,f,,0.05,2,10    !力收敛准则设定
sstif,off   
eqslv,pcg,1.0e-8,2.0    !使用预条件共轭梯度法，内存需求大，硬盘需求小
nropt,full,,on !采用完全的newton-raphson法，采用自适应下降
outpr,basic,all  
outres,basic,all     
autots,off     !打开自动时间步，以利于采用较少资源获得有效解
kbc,0         !荷载设置为斜坡荷载

!将顶面的uz方向耦合在一起(耦合到边节点)
nsel,s,loc,y,50
cp,200,uz,all   !耦合到1点
nsel,s,node,,1
cm,jiazaidian,node

!底面施加约束，形成一个铰
nsel,s,loc,y,-50
d,all,all,0

施加顶面面荷载
/psymb,adir,1 !查看面法线方向，press是指压力，以面的方向为正，
/pbc,all,,1
allsel
TIME,1
kbc=0
NSUBST,5
nsel,s,loc,y,50
sf,all,pres,50
alls
LSWRITE,1
solve
save

TIME,2
D,jiazaidian,uz,5
kbc=0
NSUBST,5
alls
LSWRITE,2
solve
save


TIME,3
D,jiazaidian,uz,-5
kbc=0
NSUBST,5
alls
LSWRITE,3
solve
save

TIME,4
D,jiazaidian,uz,10
kbc=0
NSUBST,5
alls
LSWRITE,4
solve
save

TIME,5
D,jiazaidian,uz,-10
kbc=0
NSUBST,5
alls
LSWRITE,5
solve
save

清楚刚才的结果，保存一些图片进行对比，仍然文件名为test
finish
/clear
/title,connection
/prep7
/view,1,1,1,1                !设置所有窗口的视图方式
/nerr,0
!采用mm-MPA-N单位制度  !在MPA单位下，力为N，长度为mm，弹模为MPA，质量为t（吨），重力加速度!为9815
et,1,solid65     !定义钢材的单元为solid45
!钢梁H型钢腹板（材料1)
MP,EX,1,2.06e5
MP,PRXY,1,0.3
r,45
tb,bkin,1 !采用双线性随动强化，材料1
TBTEMP,0
tbdata,,270,2.06e4 !定义屈服强度和切向模量
!BLOCK, X1, X2, Y1, Y2, Z1, Z2
block,20,-20,50,-50,20,-20,
esizie,10        
mshape,0,3d
mshkey,1
vatt,1,45,1 !材料7，单元类型2
vmesh,all

!!!!!!!!!!!求解!!!!!!!!!!!!!!!!
/solu
allsel
!lsclear,all !清空所有载荷
antype,static !静力分析
nlgeom,on              !打开大变形选项
!lnsrch,on               !打开线性搜索,与自适应下降互斥
neqit,10000              !每一子步中方程的迭代次数限值            
cnvtol,f,,0.05,2,10    !力收敛准则设定
sstif,off   
eqslv,pcg,1.0e-8,2.0    !使用预条件共轭梯度法，内存需求大，硬盘需求小
nropt,full,,on !采用完全的newton-raphson法，采用自适应下降
outpr,basic,all  
outres,basic,all     
autots,off     !打开自动时间步，以利于采用较少资源获得有效解
kbc,0         !荷载设置为斜坡荷载
rescontrol,,3,last

!将顶面的uz方向耦合在一起(耦合到边节点)
nsel,s,loc,y,50
cp,200,uz,all   !耦合到1点
nsel,s,node,,1
cm,jiazaidian,node

!底面施加约束，形成一个铰
nsel,s,loc,y,-50
d,all,all,0

施加顶面面荷载
/psymb,adir,1 !查看面法线方向，press是指压力，以面的方向为正，
/pbc,all,,1
allsel
TIME,1
kbc=0
NSUBST,5
nsel,s,loc,y,50
sf,all,pres,50
alls
LSWRITE,1
solve
save

TIME,2
D,jiazaidian,uz,5
kbc=0
NSUBST,5
alls
LSWRITE,2
solve
save


TIME,3
D,jiazaidian,uz,-5
kbc=0
NSUBST,5
alls
LSWRITE,3
solve
save
finish

/clear
/filename,test
/solu
antype,,rest,3,last

TIME,4
D,jiazaidian,uz,10
kbc=0
NSUBST,5
alls
LSWRITE,4
solve
save

TIME,5
D,jiazaidian,uz,-10
kbc=0
NSUBST,5
alls
LSWRITE,5
solve
save

htbbzzg

1  D 命令就是定义位移约束的，无论是零位移还是非零位移，都是对位移施加了限制。
   这个问题以前有网友问过我，我建议他做了一个验证，证明了修改位移约束条件后，ANSYS 确实是作为新的分析，与上一步无关了。
2  你的命令流中：
    autots,off     !打开自动时间步，以利于采用较少资源获得有效解
   是关闭自动时间步，与你的注释相反。