MIDASGTS计算不收敛报错原因和解决方法.doc
《MIDASGTS计算不收敛报错原因和解决方法.doc》由会员分享,可在线阅读,更多相关《MIDASGTS计算不收敛报错原因和解决方法.doc(2页珍藏版)》请在冰点文库上搜索。
1、程序分析时报错可能的原因?
输出结果如下:
>ENTRYLOAD_BLOCK
>ENTRYASSM_STIFF_LOAD(FORMGLOBALEQUILIBRIUMEQUATIONSINBLOCKS)
>Inclass[MultifrontalSolver]:
FactorizationfailedinDSYTRFINCOMPLETEMULTI-FRONTALSOLVER
>IERROR=-1
>不正常的执行中断!
1)问题分析:
出现如上提示时,最常见的原因是模型中出现了不确定的刚体位移。
在分析中必须对模型中所有的节点都定义足够的约束条件,以保证它们在各个平移和转动自由度上都不会出现不确定的刚体位移。
2)类试问题
A-桁架单元不能受弯,而在模型中没有约束它的转动自由度,就出现了不确定位移值。
B-施加接触单元后因为分割节点,导致分割产生的节点没有约束条件,也会出现不确定位移值。
C-单元异常,例如单元过度变形
D-应力应变关系中有负斜率
2、程序计算不收敛可能的原因?
输出结果如下:
X-COORDINATEOFC.G....:
1.49753E+01M
Y-COORDINATEOFC.G....:
1.50000E+01M
Z-COORDINATEOFC.G....:
1.12137E+01M
INCREMENTNO.:
1ITERATIONNO.:
1LOADFACTOR:
0.100E+01DISPL.NORM:
0.478E-02
INCREMENTNO.:
1ITERATIONNO.:
2LOADFACTOR:
0.100E+01DISPL.NORM:
0.516E-02
INCREMENTNO.:
1ITERATIONNO.:
3LOADFACTOR:
0.100E+01DISPL.NORM:
0.249E-01
INCREMENTNO.:
1ITERATIONNO.:
4LOADFACTOR:
0.100E+01DISPL.NORM:
0.165E-01
INCREMENTNO.:
1ITERATIONNO.:
5LOADFACTOR:
0.100E+01DISPL.NORM:
0.953E-02
INCREMENTNO.:
1ITERATIONNO.:
6LOADFACTOR:
0.100E+01DISPL.NORM:
0.318E+02
INCREMENTNO.:
1ITERATIONNO.:
7LOADFACTOR:
0.100E+01DISPL.NORM:
0.300E+00
INCREMENTNO.:
1ITERATIONNO.:
8LOADFACTOR:
0.100E+01DISPL.NORM:
0.353E+00
INCREMENTNO.:
1ITERATIONNO.:
9LOADFACTOR:
0.100E+01DISPL.NORM:
0.847E-01
INCREMENTNO.:
1ITERATIONNO.:
10LOADFACTOR:
0.100E+01DISPL.NORM:
0.133E+00
INCREMENTNO.:
1ITERATIONNO.:
11LOADFACTOR:
0.100E+01DISPL.NORM:
0.715E-01
INCREMENTNO.:
1ITERATIONNO.:
12LOADFACTOR:
0.100E+01DISPL.NORM:
0.964E-01
INCREMENTNO.:
1ITERATIONNO.:
13LOADFACTOR:
0.100E+01DISPL.NORM:
0.629E-01
INCREMENTNO.:
1ITERATIONNO.:
14LOADFACTOR:
0.100E+01DISPL.NORM:
0.324E-01
INCREMENTNO.:
1ITERATIONNO.:
15LOADFACTOR:
0.100E+01DISPL.NORM:
0.454E+01
INCREMENTNO.:
1ITERATIONNO.:
16LOADFACTOR:
0.100E+01DISPL.NORM:
0.385E+00
INCREMENTNO.:
1ITERATIONNO.:
17LOADFACTOR:
0.100E+01DISPL.NORM:
0.164E+00
INCREMENTNO.:
1ITERATIONNO.:
18LOADFACTOR:
0.100E+01DISPL.NORM:
0.111E+00
INCREMENTNO.:
1ITERATIONNO.:
19LOADFACTOR:
0.100E+01DISPL.NORM:
0.512E+02
INCREMENTNO.:
1ITERATIONNO.:
20LOADFACTOR:
0.100E+01DISPL.NORM:
0.269E+00
INCREMENTNO.:
1ITERATIONNO.:
21LOADFACTOR:
0.100E+01DISPL.NORM:
0.151E+00
INCREMENTNO.:
1ITERATIONNO.:
22LOADFACTOR:
0.100E+01DISPL.NORM:
0.896E-01
INCREMENTNO.:
1ITERATIONNO.:
23LOADFACTOR:
0.100E+01DISPL.NORM:
0.339E+01
INCREMENTNO.:
1ITERATIONNO.:
24LOADFACTOR:
0.100E+01DISPL.NORM:
0.120E+00
INCREMENTNO.:
1ITERATIONNO.:
25LOADFACTOR:
0.100E+01DISPL.NORM:
0.129E+00
INCREMENTNO.:
1ITERATIONNO.:
26LOADFACTOR:
0.100E+01DISPL.NORM:
0.711E-01
INCREMENTNO.:
1ITERATIONNO.:
27LOADFACTOR:
0.100E+01DISPL.NORM:
0.271E-01
INCREMENTNO.:
1ITERATIONNO.:
28LOADFACTOR:
0.100E+01DISPL.NORM:
0.588E-01
INCREMENTNO.:
1ITERATIONNO.:
29LOADFACTOR:
0.100E+01DISPL.NORM:
0.108E-01
INCREMENTNO.:
1ITERATIONNO.:
30LOADFACTOR:
0.100E+01DISPL.NORM:
0.700E-02
1)问题分析:
和的关系是在试验基础上由本构模型给定的,在做有限元计算时,它是已知的。
而与的关系是未知的,它决定于和的关系。
程序根据非线性迭代方法多次计算,逐步逼近真实解,最后达到非线性控制里面定义的收敛标准,就收敛了。
而如上结果是不收敛的,结果自然也就不合理。
2)类试问题
A-重复单元
解决方法:
模型-单元-删除-仅删除重复单元
B-节点不耦合
解决方法:
网格-检查网格-自由面(自由边),找到不耦合的地方删除重新划分
C-多余节点
解决方法:
节点-删除-仅删除多余节点
D-局部塑性变形过大
解决方法:
修改土体材料本构或参数