UEL求解过程中的位移和位移增量的变化规律

蓝牙

这是我编写的一个九节点线弹性的UEL在求解过程中输出的某节点的U和delte(U)
加载形式是在把载荷0.005按照10个时间步加载上去
time 是这样变化的 0  1E-7   0.1   0.2
0调用了一次UEL
1E-7 调用了一次UEL
0.1 以后都是调用2次  

D和delte(U)的规律
0.1时刻以前都是0
0.1第一次D=0.0，delte(U)=0.0
0.1第二次编程了D=0.0005，delte(U)=0.0005

0.2第一次D=0.001，delte(U)=0.0005
0.2第二次编程了D=0.001，delte(U)=0.0005

0.3第一次D=0.0015，delte(U)=0.0005
0.3第二次编程了D=0.0015，delte(U)=0.0005

我总觉得在0.1 第一次uel应该为D=0.0，delte(U)=0.0005才对吧？

帮我解释一下ABAQUS的求解规律









TIME(1),TIME(2),KSTEP,KINC:                                             
  0.000000000000000E+000  0.000000000000000E+000           1           1
DU()                   某节点的Y方向的位移:                             
  0.0000E+00                                                            
U()                    某节点的Y方向的增量位移:                        
  0.0000E+00                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  1.000000000000000E-007  1.000000000000000E-007           1           0
DU()                   某节点的Y方向的位移:                             
  0.0000E+00                                                            
U()                    某节点的Y方向的增量位移:                        
  0.0000E+00                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.100000000000000       0.100000000000000                1           1
DU()                   某节点的Y方向的位移:                             
  0.0000E+00                                                            
U()                    某节点的Y方向的增量位移:                        
  0.0000E+00                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.100000000000000       0.100000000000000                1           1
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.5000E-03                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.200000000000000       0.200000000000000                1           2
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.1000E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.200000000000000       0.200000000000000                1           2
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.1000E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.300000000000000       0.300000000000000                1           3
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.1500E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.300000000000000       0.300000000000000                1           3
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.1500E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.400000000000000       0.400000000000000                1           4
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.2000E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.400000000000000       0.400000000000000                1           4
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.2000E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.500000000000000       0.500000000000000                1           5
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.2500E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.500000000000000       0.500000000000000                1           5
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.2500E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.600000000000000       0.600000000000000                1           6
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.3000E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.600000000000000       0.600000000000000                1           6
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.3000E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.700000000000000       0.700000000000000                1           7
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.3500E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.700000000000000       0.700000000000000                1           7
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.3500E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.800000000000000       0.800000000000000                1           8
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.4000E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.800000000000000       0.800000000000000                1           8
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.4000E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.900000000000000       0.900000000000000                1           9
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.4500E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
  0.900000000000000       0.900000000000000                1           9
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.4500E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
   1.00000000000000        1.00000000000000                1          10
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.5000E-02                                                            
--------------------------------------------------------------------   
TIME(1),TIME(2),KSTEP,KINC:                                             
   1.00000000000000        1.00000000000000                1          10
DU()                   某节点的Y方向的位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的增量位移:                        
  0.5000E-02                                                            
--------------------------------------------------------------------

amani

我觉得1E-7是试计算，所以只调用UEL一次，正式的计算值应该从0.1开始的，调用两次。

蓝牙

2# amani

delta U 和U的变化规律显示



U为当前时间步的位移值，delta（U）只能理解成前一时间不的位移和当前时间步的位移差值，而不能理解为下一时间步位移和当前时间步位移的差值

否则在0.1第一次调用UEl时应该是deltaU=0.0005 而U为0才对


还有两次调用UEl都是干什么的，你有次回帖中提到过，但是我后来一直都没有搜到了，谢谢了
我的子程序就差这部分结果验证了！

fr_zengtao

两次条用UEL 和 UMAT里面的每个STEP的第一个INCREMENT都调用两次UMAT是不是一个意思呢，便于形成初始的刚度矩阵，个人认为UEL的目的是不是就是为了告诉程序你的刚度矩阵是如何形成的呢，

蓝牙

我的重点是U和delta(u)

U应该为为当前时间步的位移值，delta（U）只能理解成前一时间不的位移和当前时间步的位移差值
否则
TIME(1),TIME(2),KSTEP,KINC:                                             
   1.00000000000000        1.00000000000000                1          10
DU()                   某节点的Y方向的增量位移:                             
  0.5000E-03                                                            
U()                    某节点的Y方向的位移:                        
  0.5000E-02           

Delta（U） 应该为零！

turen111

call subroutine这部分可以参看abaqus analysis user manual 15.1.1 User subroutines: overview-User subroutine calls：Subroutines that define material, element, or interface behavior
Most user subroutines that are used to define material, element, or interface behavior are called twice per material point, element, or slave surface node in the first iteration of every increment such that the model's initial stiffness matrix can be formulated appropriately for the step procedure chosen. The subroutines are called only once per material point, element, or slave surface node in each succeeding iteration within the increment.

一定要区分开是在每个增量步（increment）的第一个迭代步（iteration）调用两次uel，umat等，伺候在这个此增量步内，每个迭代步只会调用一次uel，umat等，而不是前面说的每个step，里面的第一个increment

turen111

我赞成du矩阵是上一个收敛增量步位移矩阵u与当前增量步u矩阵之间的差值：
见如下的结果：
KINC
u（1）        un（1）          du（1，1）


           1
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
           1
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
           1
  3.296703296703296E-034  0.000000000000000E+000  3.296703296703296E-034
           1
-2.324086481143608E-036  0.000000000000000E+000 -2.324086481143608E-036
           1
  1.476504708885424E-036  0.000000000000000E+000  1.476504708885424E-036
           1
  1.023026126315315E-037  0.000000000000000E+000  1.023026126315315E-037
           2
  2.046052252630630E-037  1.023026126315315E-037  1.023026126315315E-037
           2
-9.027799165575401E-038  1.023026126315315E-037 -1.925806042872855E-037
           3
-2.828585959430395E-037 -9.027799165575401E-038 -1.925806042872855E-037
           3
  4.088694941963723E-040 -9.027799165575401E-038  9.068686114995038E-038
           4
  9.109573064414676E-038  4.088694941963723E-040  9.068686114995038E-038
           4
-2.914815341588592E-041  4.088694941963723E-040 -4.380176476122582E-040
           5
-4.671658010281442E-040 -2.914815341588592E-041 -4.380176476122582E-040
           5
-6.812626469920116E-044 -2.914815341588592E-041  2.908002715118672E-041
           6
  2.901190088648752E-041 -6.812626469920116E-044  2.908002715118672E-041
           6
-1.047954122312911E-045 -6.812626469920116E-044  6.707831057688825E-044
           7
  6.603035645457534E-044 -1.047954122312911E-045  6.707831057688825E-044
           7
-1.895997707274989E-047 -1.047954122312911E-045  1.028994145240161E-045
           8
  1.010034168167411E-045 -1.895997707274989E-047  1.028994145240161E-045
           8
-3.320570124634503E-049 -1.895997707274989E-047  1.862792006028644E-047
           9
  1.829586304782299E-047 -3.320570124634503E-049  1.862792006028644E-047
           9
  3.692255327314068E-050 -3.320570124634503E-049  3.689795657365910E-049
          10
  4.059021190097317E-049  3.692255327314068E-050  3.689795657365910E-049
          10
  1.796934264712538E-050  3.692255327314068E-050 -1.895321062601529E-050
          11
-9.838679788899094E-052  1.796934264712538E-050 -1.895321062601529E-050
          11
-7.384510654628571E-050  1.796934264712538E-050 -9.181444919341109E-050
          12
-1.656595557396968E-049 -7.384510654628571E-050 -9.181444919341109E-050
          12
  1.846127663657140E-050 -7.384510654628571E-050  9.230638318285711E-050
          13
  1.107676598194285E-049  1.846127663657140E-050  9.230638318285711E-050
          13
-4.919339894460752E-052  1.846127663657140E-050 -1.895321062601748E-050

turen111

感觉abaqus里面你在uel里面是无法得到每一个增量步内不同迭代步之间的位移差值的，