跪问这个分析怎么能够让其收敛啊?

浏览:89083 回答:7
被老板逼疯了,可怜我学习这东西没有几天!!!
三维电弧模拟,走到第二步用计算磁场时就怎么也不会收敛了,用的电流密度是第一步得到的,请各位高手闲暇之余,帮帮小弟,多谢了!
因为转换时清楚物理环境,所以材料特性写了两遍。第二步不收敛,第一步没有问题的。三维网格模型是建立在二维模型旋转的基础上的,所以刚开始用了二维的PLANE67,实际上到三维时这个是不用的。
第一步用的是SOLID69,第二步用的是SOLID97。
Finish
/clear
/PREP7
ET,1,PLANE67
KEYOPT,1,3,1
MPTEMP,1,273,2500,5000,7500,10000,12500
MPTEMP,7,15000,17500,20000
MPDATA,DENS,1,1,1.67,0.2,0.1,0.07,0.04,0.03
MPDATA,DENS,1,7,0.02,0.015,0.014
MPDATA,KXX,1,1,0,0.1,0.2,0.3,0.6,1.6
MPDATA,KXX,1,7,2.4,2.0,2.2
MPDATA,RSVX,1,1,1000,100,1,1.3e-3,3.3e-4,1.6e-4
MPDATA,RSVX,1,7,1.2e-4,9.8e-5,8.5e-5
MPDATA,C,1,1,521,521,521,600,1250,6000
MPDATA,C,1,7,10000,4000,2500
MPDATA,VISC,1,1,2e-5,7e-5,1e-4,2.2e-4,2.7e-4,2e-4
MPDATA,VISC,1,7,6e-5,3e-5,2.6e-5
MPDATA,DENS,2,,7200
MPDATA,KXX,2,,20
MPDATA,RSVX,2,,1/770000
MPDATA,DENS,3,,19350
MPDATA,RSVX,3,,0.0000000528
MPDATA,KXX,3,,173
MPDATA,MURX,1,,1
MPDATA,MURX,3,,1
/REPLOT,RESIZE
TOFFST,0
K,1,0,0,0,
k,2,0.01,0,0,
K,3,0.0002,0,0,
K,4,0.0016,0,0,
K,5,0,0.005,0,
K,6,0.0002,0.005,0,
K,7,0,0.005+0.0014*sqrt(3),0,
K,8,0.0016,0.005+0.0014*sqrt(3),0,
K,9,0.01,0.005+0.0014*sqrt(3),0,
K,10,0.0016,0.005,0,
K,11,0.01,0.005,0,
A,5,6,8,7
A,1,3,6,5
A,3,4,10,6
A,6,10,8
A,8,9,11,10
A,10,11,2,4
LGLUE,ALL
AGLUE,ALL
ASEL,S,AREA,,2,6
AATT,1,,1,0,
ASEL,S,AREA,,1
AATT,3,,1,0,
ALLSEL,ALL
ESIZE,0.0002,0,
MSHKEY,1
AMESH,ALL
/PNUM,KP,1
/PNUM,AREA,1
ET,2,SOLID69
ESIZE,0,20,
TYPE,2
MAT,1
VROTAT,4,5,2,3,6,,1,7,90,,
TYPE,2
MAT,3
VROTAT,1,,,,,,1,7,90,,
ACLEAR,1,6,
/VIEW,1,-1,1,1
APLOT
FINISH
/SOL
TUNIF,10000,
DA,11,TEMP,1000,0
ASEL,S,AREA,,27,28
ASEL,A,AREA,,1
ASEL,A,AREA,,17
ASEL,A,AREA,,9
DA,ALL,TEMP,3000,0
ASEL,S,AREA,,15
ASEL,A,AREA,,19
ASEL,A,AREA,,23
DA,ALL,TEMP,5000,0
DA,ALL,VOLT,0,0
ASEL,S,AREA,,27
NSLA,S,1
F,ALL,AMPS,-160
ALLSEL,ALL
SOLVE
SAVE
FINISH
PHYSICS,CLEAR
/PREP7
ET,2,SOLID97
MPTEMP,1,273,2500,5000,7500,10000,12500
MPTEMP,7,15000,17500,20000
MPDATA,DENS,1,1,1.67,0.2,0.1,0.07,0.04,0.03
MPDATA,DENS,1,7,0.02,0.015,0.014
MPDATA,KXX,1,1,0,0.1,0.2,0.3,0.6,1.6
MPDATA,KXX,1,7,2.4,2.0,2.2
MPDATA,RSVX,1,1,1000,100,1,1.3e-3,3.3e-4,1.6e-4
MPDATA,RSVX,1,7,1.2e-4,9.8e-5,8.5e-5
MPDATA,C,1,1,521,521,521,600,1250,6000
MPDATA,C,1,7,10000,4000,2500
MPDATA,VISC,1,1,2e-5,7e-5,1e-4,2.2e-4,2.7e-4,2e-4
MPDATA,VISC,1,7,6e-5,3e-5,2.6e-5
MPDATA,DENS,2,,7200
MPDATA,KXX,2,,20
MPDATA,RSVX,2,,1/770000
MPDATA,DENS,3,,19350
MPDATA,RSVX,3,,0.0000000528
MPDATA,KXX,3,,173
MPDATA,MURX,1,,1
MPDATA,MURX,3,,1
/REPLOT,RESIZE
TOFFST,0
FINISH
/SOL
TUNIF,10000,
DA,11,TEMP,1000,0
ASEL,S,AREA,,27,28
ASEL,A,AREA,,1
ASEL,A,AREA,,17
ASEL,A,AREA,,9
DA,ALL,TEMP,3000,0
ASEL,S,AREA,,15
ASEL,A,AREA,,19
ASEL,A,AREA,,23
DA,ALL,TEMP,5000,0
ALLSEL,ALL
LDREAD,JS,,,,,,RST
SOLVE
FINISH
SAVE
邀请回答 我来回答

全部回答

(7)
默认 最新
morningsun613

extra displacements options
这一命令转换成命令流格式应该如何表达呢?应该不是SOLCONTROL这个命令吧?
2012年7月30日
评论 点赞
解决非线性分析不收敛的技巧
大家都提到了收敛困难的问题为加速收敛应该注意一下几个问题 17
1 收敛容差ANSYS缺省的收敛准则会根据单元的不同而检查不同的收敛力素和容差例如当采用solid65和link8时,缺省的要检查F和DISP两个力素其容差也是缺省的(Help中有)对于钢筋混凝土结构一般而言其位移比较小仅使用F力素收敛即可但其容差也同时放松一般采用5%即可(缺省是5)命令:cnvtol,f,,0.05,2
2 其它选项的设置
自动时间步打开此选择可以让程序决定子步间荷载增量的大小及其是增加或是减小收敛速度较快(命令autots,1)打开后似乎定义的子步数不起控制作用了
打开线性搜索可以帮助收敛的速度(命令:lnsrch,1)
打开预测器可以帮助收敛的速度(命令red,on)
平衡迭代次数在每一子步中的迭代次数缺省是25,将其增加例如改为50(命令: neqit,50)
NSUBST此值不宜过小否则计算过程中老是调整影响计算速度
当然对于比较简单的算例或是分布模型可能不需要如此多的选项但对于复杂的模型是需要的各位可以试试
影响非线性收敛稳定性及其速度的因素很多:
1、模型——主要是结构刚度的大小。对于某些结构,从概念的角度看,可以认为它是几何不变的稳定体系。但如果结构相近的几个主要构件刚度相差悬殊,在数值计算中就可能导致数值计算的较大误差,严重的可能会导致结构的几何可变性——忽略小刚度构件的刚度贡献。 如出现上述的结构,要分析它,就得降低刚度很大的构件单元的刚度,可以加细网格划分,或着改用高阶单元(BEAM->SHELL,SHELL->SOLID)。构件的连接形式(刚接或铰接)等也可能影响到结构的刚度。
2、线性算法(求解器)。ANSYS中的非线性算法主要有:稀疏矩阵法(SPARSE DIRECT SOLVER)、预共轭梯度法(PCG SOLVER)和波前法(FRONT DIRECT SLOVER)。稀疏矩阵法是性能很强大的算法,一般默认即为稀疏矩阵法(除了子结构计算默认波前法外)。预共轭梯度法对于3-D实体结构而言是最优的算法,但当结构刚度呈现病态时,迭代不易收敛。为此推荐以下算法:
1)、BEAM单元结构,SHELL单元结构,或以此为主的含3-D SOLID的结构,用稀疏矩阵法;
2)、3-D SOLID的结构,用预共轭梯度法;
3)、当你的结构可能出现病态时,用稀疏矩阵法;
4)、当你不知道用什么时,可用稀疏矩阵法。
3、非线性逼近技术。在ANSYS里还是牛顿-拉普森法和弧长法。牛顿-拉普森法是常用的方法,收敛速度较快,但也和结构特点和步长有关。弧长法常被某些人推崇备至,它能算出力加载和位移加载下的响应峰值和下降响应曲线。但也发现:在峰值点,弧长法仍可能失效,甚至在非线性计算的线性阶段,它也可能会无法收敛。
为此,尽量不要从开始即激活弧长法,还是让程序自己激活为好(否则出现莫名其妙的问题)。子步(时间步)的步长还是应适当,自动时间步长也是很有必要的。
4、加快计算速度
在大规模结构计算中,计算速度是一个非常重要的问题。下面就如何提高计算速度作一些建议:
充分利用ANSYS MAP分网和SWEEP分网技术,尽可能获得六面体网格,这一方面减小解题规模,另一方面提高计算精度。
在生成四面体网格时,用四面体单元而不要用退化的四面体单元。比如95号单元有20节点,可以退化为10节点四面体单元,而92号单元为10节点单元,在此情况下用92号单元将优于95号单元。
选择正确的求解器。对大规模问题,建议采用PCG法。此法比波前法计算速度要快10倍以上(前提是您的计算机内存较大)。对于工程问题,可将ANSYS缺省的求解精度从1E-8改为1E-4或1E-5即可。
5、荷载步的设置直接影响到收敛。应该注意以下几点:
1、设置足够大的荷载步(将MAXMIUM SUBSTEP=1000000),可以更容易收敛,避免发散的出现(nsub,nsbstp,nsbmx,nsbmn);
2、设置足够大的平衡迭代步数,默认为25,可以放大到很大(100)(eqit,eqit);
3、将收敛准则调整,以位移控制时调整为0.05,以力控制为0.01(CNVTOL,lab,value,toler,norm,minref)。
4、对于线性单元和无中间节点的单元(SOLID65和SOLID45),关闭EXTRA DISPLACEMENTS OPTIONS(在OPTIONS中)。
5、对于CONCRETE材料,可以关闭压碎功能,将CONCRETE中的单轴抗压强度设置为-1(tadata,mat,shrcf-op,shrcf-cl,UntensSt,UnCompSt(-1))。
2009年5月19日
评论 点赞
不是,是设置非线性运算
2009年5月19日
评论 点赞
hdzw
原帖由 cards101 于 2009-5-18 17:03:00 发表
cnvtol,f,,0.05,2,0.01 !定义收敛条件,使用缺省的VALUE
outres,all,all !输出每一子步的结果
autots,on !打开自动时间步控制
lnsrch,on !打开线性搜索
arclen,off !关闭弧长选项
ncnv,2 !如果不收敛时结束......

什么意思》使用上面的命令检测吗?
2009年5月18日
评论 点赞
cnvtol,f,,0.05,2,0.01 !定义收敛条件,使用缺省的VALUE
outres,all,all !输出每一子步的结果
autots,on !打开自动时间步控制
lnsrch,on !打开线性搜索
arclen,off !关闭弧长选项
ncnv,2 !如果不收敛时结束而不退出
neqit,20 !每一子步中方程的迭代次数限值
pred,on !打开预测器
2009年5月18日
评论 点赞
hdzw
能说的具体点吗?
我试过:
CNVTOL,AMPS,,,,-1,
结果是什么也没有了。
其他的条件改过,要么出错,要么还是不行啊!!
2009年5月18日
评论 点赞
hdzw
怎么没有人可以出手相助啊?!
2009年5月17日
评论 点赞

没解决?试试专家一对一服务

换一批
    App下载
    技术邻APP
    工程师必备
    • 项目客服
    • 培训客服
    • 平台客服

    TOP