找回密码
 注册
Simdroid-非首页
查看: 243|回复: 2

求解选项的含义,高手出出招吧

[复制链接]
发表于 2004-6-24 10:29:57 | 显示全部楼层 |阅读模式 来自 黑龙江哈尔滨
下面主要分为4个部分,关于收敛求解选项设置的方法,高手能出出招,挨个作一下解释,并将对应的控制命令作一下简要的介绍吗,新手们也大家一起动手,多多补充。
  
!静态分析输出 S T E P       1     S T A T I C   A N A L Y S I S
  
                                                                                          
  
1。!时间步选项(初始步长,最大最小,户长比例)
          AUTOMATIC TIME CONTROL WITH -
          A SUGGESTED INITIAL TIME INCREMENT OF                 1.00     
          AND A TOTAL TIME PERIOD OF                            1.00     
          THE MINIMUM TIME INCREMENT ALLOWED IS                1.000E-05
          THE MAXIMUM TIME INCREMENT ALLOWED IS                 1.00     
  
2。!自动稳定求解控制
         AUTOMATIC STABILIZATION WITH DISSIPATION FACTOR =   2.000E-03
  
3。!力控制收敛  
CONVERGENCE TOLERANCE PARAMETERS FOR FORCE     
    !力残差系数 CRITERION FOR RESIDUAL FORCE     FOR A NONLINEAR PROBLEM          5.000E-03
     !位移残差系数CRITERION FOR DISP.    CORRECTION IN A NONLINEAR PROBLEM          1.000E-02
     INITIAL VALUE OF TIME AVERAGE FORCE                               1.000E-02
     AVERAGE FORCE     IS TIME AVERAGE FORCE     
     ALTERNATE CRIT. FOR RESIDUAL FORCE     FOR A NONLINEAR PROBLEM    2.000E-02
     CRITERION FOR ZERO FORCE     RELATIVE TO TIME AVRG. FORCE         1.000E-05
     CRITERION FOR RESIDUAL FORCE     WHEN THERE IS ZERO FLUX          1.000E-05
     CRITERION FOR DISP.    CORRECTION WHEN THERE IS ZERO FLUX         1.000E-03
     CRITERION FOR RESIDUAL FORCE     FOR A LINEAR INCREMENT           1.000E-08
     FIELD CONVERSION RATIO                                             1.00     
     CRITERION FOR ZERO FORCE     REL. TO TIME AVRG. MAX. FORCE        1.000E-05
     CRITERION FOR ZERO DISP.    RELATIVE TO CHARACTERISTIC LENGTH     1.000E-08
  
  4。力矩收敛控制
     CONVERGENCE TOLERANCE PARAMETERS FOR MOMENT   
     CRITERION FOR RESIDUAL MOMENT    FOR A NONLINEAR PROBLEM          5.000E-03
     CRITERION FOR ROTATION CORRECTION IN A NONLINEAR PROBLEM          1.000E-02
     INITIAL VALUE OF TIME AVERAGE MOMENT                              1.000E-02
     AVERAGE MOMENT    IS TIME AVERAGE MOMENT   
     ALTERNATE CRIT. FOR RESIDUAL MOMENT    FOR A NONLINEAR PROBLEM    2.000E-02
     CRITERION FOR ZERO MOMENT    RELATIVE TO TIME AVRG. MOMENT        1.000E-05
     CRITERION FOR RESIDUAL MOMENT    WHEN THERE IS ZERO FLUX          1.000E-05
     CRITERION FOR ROTATION CORRECTION WHEN THERE IS ZERO FLUX         1.000E-03
     CRITERION FOR RESIDUAL MOMENT    FOR A LINEAR INCREMENT           1.000E-08
     FIELD CONVERSION RATIO                                             1.00     
     CRITERION FOR ZERO MOMENT    REL. TO TIME AVRG. MAX. MOMENT       1.000E-05
  
     VOLUMETRIC STRAIN COMPATIBILITY TOLERANCE FOR HYBRID SOLIDS       1.000E-05
     AXIAL STRAIN COMPATIBILITY TOLERANCE FOR HYBRID BEAMS             1.000E-05
     TRANS. SHEAR STRAIN COMPATIBILITY TOLERANCE FOR HYBRID BEAMS      1.000E-05
     SOFT CONTACT CONSTRAINT COMPATIBILITY TOLERANCE FOR P&gt0          5.000E-03
     SOFT CONTACT CONSTRAINT COMPATIBILITY TOLERANCE FOR P=0.0         0.100     
     DISPLACEMENT COMPATIBILITY TOLERANCE FOR DCOUP ELEMENTS           1.000E-05
     ROTATION COMPATIBILITY TOLERANCE FOR DCOUP ELEMENTS               1.000E-05
  
  6.时间步控制参数
   TIME INCREMENTATION CONTROL PARAMETERS:
  *** FIRST EQUILIBRIUM ITERATION FOR CONSECUTIVE DIVERGENCE CHECK              8
  *** EQUILIBRIUM ITERATION AT WHICH LOG. CONVERGENCE RATE CHECK BEGINS        10
     EQUILIBRIUM ITERATION AFTER WHICH ALTERNATE RESIDUAL IS USED              9
     MAXIMUM EQUILIBRIUM ITERATIONS ALLOWED                                   16
     EQUILIBRIUM ITERATION COUNT FOR CUT-BACK IN NEXT INCREMENT               10
     MAXIMUM EQUILIB. ITERS IN TWO INCREMENTS FOR TIME INCREMENT INCREASE      4
     MAXIMUM ITERATIONS FOR SEVERE DISCONTINUITIES                            12
     MAXIMUM CUT-BACKS ALLOWED IN AN INCREMENT                                 5
     MAXIMUM DISCON. ITERS IN TWO INCREMENTS FOR TIME INCREMENT INCREASE       6
     CUT-BACK FACTOR AFTER DIVERGENCE                                 0.2500     
     CUT-BACK FACTOR FOR TOO SLOW CONVERGENCE                         0.5000     
     CUT-BACK FACTOR AFTER TOO MANY EQUILIBRIUM ITERATIONS            0.7500     
     CUT-BACK FACTOR AFTER TOO MANY SEVERE DISCONTINUITY ITERATIONS   0.2500     
     CUT-BACK FACTOR AFTER PROBLEMS IN ELEMENT ASSEMBLY               0.2500     
     INCREASE FACTOR AFTER TWO INCREMENTS THAT CONVERGE QUICKLY        1.500     
     MAX. TIME INCREMENT INCREASE FACTOR ALLOWED                       1.500     
     MAX. TIME INCREMENT INCREASE FACTOR ALLOWED (DYNAMICS)            1.250     
     MAX. TIME INCREMENT INCREASE FACTOR ALLOWED (DIFFUSION)           2.000     
     MINIMUM TIME INCREMENT RATIO FOR EXTRAPOLATION TO OCCUR          0.1000     
     MAX. RATIO OF TIME INCREMENT TO STABILITY LIMIT                   1.000     
     FRACTION OF STABILITY LIMIT FOR NEW TIME INCREMENT               0.9500     
  
  *** INDICATES USE OF NON-DEFAULT TIME INCREMENTATION CONTROL PARAMETER
     SELECTED WITH *CONTROLS
  
     RESTART FILE WILL BE WRITTEN EVERY     1  INCREMENTS
  
     ONLY THE LATEST INCREMENT WILL BE RETAINED
  
     THE MAXIMUM NUMBER OF INCREMENTS IN THIS STEP IS                     100
  
          LARGE DISPLACEMENT THEORY WILL BE USED
  
     LINEAR EXTRAPOLATION WILL BE USED
  
     CHARACTERISTIC ELEMENT LENGTH      2.00     
  
     PRINT OF INCREMENT NUMBER, TIME, ETC., TO THE MESSAGE FILE EVERY     1  INCREMENTS
 楼主| 发表于 2004-6-24 10:36:16 | 显示全部楼层 来自 黑龙江哈尔滨

回复: 求解选项的含义,高手出出招吧

Simdroid开发平台
   各位论坛大哥,history data是*.inp的最重要组成部分,特别是必须的一些求解选现的设置,新手对于这些都有很大困惑。
   列位abaqus后生,都盼望斑竹cae519,还有 chenhu等多多高手,给大家在这方面作一刊讲座,好好提高一下.
 楼主| 发表于 2004-6-25 11:13:54 | 显示全部楼层 来自 黑龙江哈尔滨

回复: 求解选项的含义,高手出出招吧

自顶一下
  
帮助上,这些对于.dat 中关于求解选项设置做了很详细的解释.
ABAQUS/Standard output
The controls in effect for an analysis are listed in the data (.dat) and message (.msg) files. Nondefault controls are marked by **. For example, the input  
  
*CONTROLS, PARAMETERS=TIME INCREMENTATION
    7, 10
** ,  
would result in the following output:  
TIME INCREMENTATION CONTROL PARAMETERS:
  *** FIRST EQUIL. ITERATION FOR CONSECUTIVE DIVERGENCE CHECK       7
  *** EQUIL. ITER. AT WHICH LOG. CONVERGENCE RATE CHECK BEGINS     10
     EQUIL. ITER. AFTER WHICH ALTERNATE RESIDUAL IS USED           9
     MAXIMUM EQUILIBRIUM ITERATIONS ALLOWED                       16
     EQUIL. ITERATION COUNT FOR CUT-BACK IN NEXT INCREMENT        10
     MAX EQUIL. ITERS IN TWO INCREMENTS FOR TIME INC. INCREASE     4
     MAXIMUM ITERATIONS FOR SEVERE DISCONTINUITIES                12
     MAXIMUM CUT-BACKS ALLOWED IN AN INCREMENT                     5
     MAX DISCON. ITERS IN TWO INCS FOR TIME INC. INCREASE          6
     CUT-BACK FACTOR AFTER DIVERGENCE                          0.250
     CUT-BACK FACTOR FOR TOO SLOW CONVERGENCE                  0.500
     CUT-BACK FACTOR AFTER TOO MANY EQUILIBRIUM ITERATIONS     0.750
Activating the “line search” algorithm
  
In strongly nonlinear problems the full Newton algorithm used in ABAQUS/Standard may sometimes diverge in the early iterations. In that case it is better to “scale” the solution correction to prevent this behavior. The line search algorithm (discussed in “Convergence criteria for nonlinear problems,” Section 8.3.3) makes it possible to carry out such scaling.
  
By default, the line search algorithm is disabled. Use the PARAMETERS=LINE SEARCH parameter on the *CONTROLS option and set the first data item—the maximum number of zero-finding iterations —to a reasonable value (such as 4) to activate the line search procedure. The default value is .
  
Usage: *CONTROLS, PARAMETERS=LINE SEARCH
  
Defining tolerances for constraint equations
  
Use PARAMETERS=CONSTRAINTS on the *CONTROLS option to set tolerances for constraint equations. This option can be used to set strain compatibility tolerances for hybrid elements, displacement and rotation compatibility tolerances for distributing coupling constraints (specified by the *COUPLING option or by using DCOUP2D/DCOUP3D elements), and compatibility tolerances for softened contact. See “Convergence criteria for nonlinear problems,” Section 8.3.3, for details.
  
Controlling the solution accuracy in direct cyclic analysis
  
Use TYPE=DIRECT CYCLIC on the *CONTROLS option to specify when to impose the periodicity conditions and to set tolerances for stabilized state and plastic ratchetting detections.  
  
Usage: *CONTROLS, TYPE=DIRECT CYCLIC  
, ,  , ,   
  
Imposing the periodicity condition
The iteration number at which the periodicity condition is first imposed, , can be specified on the data line. The default value is  = 1, in which case the periodicity condition is imposed for all iterations from the beginning of an analysis. This data item rarely needs to be reset from its default value.
  
Defining tolerances for stabilized state and plastic ratchetting detections.
The stabilized state detection criteria,  and , can be specified on the data line.  is the maximum allowable ratio of the largest residual coefficient on any terms in the Fourier series to the corresponding average flux norm, and  is the maximum allowable ratio of the largest correction to the displacement coefficient on any terms in the Fourier series to the largest displacement coefficient. The default values are  = 5 × 10–3 and  = 5 × 10–3. The solution converges to a stabilized state if both of these criteria are satisfied.
  
If plastic ratchetting occurs, the shape of the stress-strain curves remains unchanged but the mean value of the plastic strain over a cycle continues to shift from one iteration to the next. In that case it is desirable to use separate tolerances for the constant term in the Fourier series to detect the plastic ratchetting.
  
The plastic ratchetting detection criteria,  and , can also be specified on the data line.  is the maximum allowable ratio of the largest residual coefficient on the constant term in the Fourier series to the corresponding average flux norm, and  is the maximum allowable ratio of the largest correction to the displacement coefficient on the constant term in the Fourier series to the largest displacement coefficient. The default values are  = 5 × 10–3 and  = 5 × 10–3. Plastic ratchetting is expected if the residual coefficients and the corrections to the displacement coefficients on any of the periodic terms are within the tolerances set by  and , respectively, but the maximum residual coefficient on the constant term and the maximum correction to the displacement coefficient on the constant term exceed the tolerances set by  and , respectively.
  
ABAQUS/Standard output
The controls in effect for an analysis are listed in the data (.dat) and message (.msg) files. Non-default controls are marked by **. For example, the input  
  
*CONTROLS, TYPE=DIRECT CYCLIC
   5, 1.0E-4, 1.0E-4, 1.0E-4, 1.E-4   
would result in the following output:  
STABILIZED STATE AND PLASTIC RATCHETTING DETECTION
PARAMETERS FOR FORCE  
** CRIT. FOR RESI. COEFF. ON ANY FOURIER TERMS            1.0E-04
** CRIT. FOR CORR. TO DISP. COEFF. ON ANY FOURIER TERMS   1.0E-04
** CRIT. FOR RESI. COEFF. ON CONSTANT FOURIER TERM        1.0E-04
** CRIT. FOR CORR. TO DISP. COEFF. ON CONST. FOURIER TERM 1.0E-04
  
PERIODICITY CONDITION CONTROL PARAMETER:
** ITERATION NUMBER AT WHICH PERIODICITY CONDITION
** STARTS TO IMPOSE           
您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-10-1 19:22 , Processed in 0.068026 second(s), 12 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表