熱固耦合求解終止

ken731101

**** negative jacobian in element        16081
  called from:           3           8
**** negative jacobian in element        16381
  called from:           3           8
**** negative jacobian in element        16381
  called from:           3           8
**** negative jacobian in element        16481
  called from:           3           8
**** negative jacobian in element        16481
  called from:           3           8
**** negative jacobian in element        16781
  called from:           3           8
**** negative jacobian in element        16781
  called from:           3           8
*** Error - the solution has a non-positive definite
stiffnes matrix
     solution time                4.9000E+00
     time step size               1.0000E-01
suggestion - decrease time step
  105080 t 4.8000E+00 dt 1.00E+06 write d3plot file
E r r o r   t e r m i n a t i o n

這是我模擬求解時所出現的E r r o r
久久無法解決這個問題
貼部份K黨的設定
這是切削模擬  刀具為碳化鎢  工件為4043低碳鋼
*KEYWORD
*TITLE
                                                                                
$
*CONTROL_SOLUTION
         2
*CONTROL_THERMAL_SOLVER
        1         1         4         0         0         1       0.9
*CONTROL_THERMAL_TIMESTEP
         0       0.0     0.006     0.006         0         0
*CONTROL_CONTACT
         0         0         2         0         0         2         0
         0         0         0         0         0                 
         0         0         0         0         0
         1
*INITIAL_TEMPERATURE_SET
         0        20
*DATABASE_FORMAT
         0
---------------------------------------
*MAT_JOHNSON_COOK
1,7.80,0.77
7.9E-3,5.1E-3,0.27,0.015,1.05,1.527E+3,20,1.1E-6
0.47E-5,-9.0,3.0,0.0,1.4,
  0.00   
$
*EOS_GRUNEISEN
1,0.47,1.49,0.00,0.00,2.15,0.455,0.0
  0.00
*MAT_THERMAL_ISOTROPIC
         1      7.80        0         0
    465E-8    50E-15
*MAT_JOHNSON_COOK
2,7.80,0.77
7.9E-3,5.1E-3,0.27,0.015,1.05,1.527E+3,20,1.1E-6
0.47E-5,-9.0,3.0,0.0,0.8,
  0.00
$
*EOS_GRUNEISEN
2,0.47,1.49,0.00,0.00,2.15,0.455,0.0
  0.00
*MAT_THERMAL_ISOTROPIC
         2      7.80        0         0
    465E-8    50E-15
*MAT_JOHNSON_COOK
3,7.80,0.77
7.9E-3,5.1E-3,0.27,0.015,1.05,1.527E+3,20,1.1E-6
0.47E-5,-9.0,3.0,0.0,1.4,
  0.00
$
*EOS_GRUNEISEN
3,0.47,1.49,0.00,0.00,2.15,0.455,0.0
  0.00
*MAT_THERMAL_ISOTROPIC
         3      7.80        0         0
    465E-8    50E-15
*MAT_JOHNSON_COOK
4,15.8,2.7
15.0E-2,1.7E-3,0.12,0.016,1.0,2.78E+3,20,1.1E-6
0.235E-5,-9.0,3.0,0.0,2.4,
  0.00
$
*EOS_GRUNEISEN
4,0.47,1.49,0.00,0.00,2.15,0.455,0.0
  0.00
*MAT_THERMAL_ISOTROPIC
4,15.8,0
245E-8,50E-15
*BOUNDARY_PRESCRIBED_MOTION_SET
4,1,-0.300E-02,1,1,0,0,0
*DEFINE_CURVE
         1
  0.000000000000E+00 -0.300E-02
  1.084000000000E+00 -0.300E-02
  2.084500000000E+00 -0.300E-02  
  3.084500000000E+00 -0.300E-02
  4.084500000000E+00 -0.300E-02
  6.084500000000E+00 -0.300E-02
  7.084500000000E+00 -0.300E-02
  8.084500000000E+00 -0.300E-02
  1.084500000000E+01 -0.300E-02
  2.084500000000E+01 -0.300E-02
  3.084500000000E+01 -0.300E-02   
--------------------------------------------------
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              CONTACT DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*CONTACT_ERODING_SURFACE_TO_SURFACE
         1         4         3         3         0         0         0         0
0.1500    0.1000     0.000     0.000     0.000             0 0.000    0.1000E+08
1.000     1.000     0.000     0.000     1.000     1.000     1.000     1.000   
         0         1         0
*CONTACT_ERODING_SURFACE_TO_SURFACE
         2         4         3         3         0         0         0         0
0.1500    0.1000     0.000     0.000     0.000             0 0.000    0.1000E+08
1.000     1.000     0.000     0.000     1.000     1.000     1.000     1.000   
         0         1         0
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE
         3         4         3         3         0         0         0         0
0.1500    0.1000     0.000     0.000     0.000             0 0.000    0.1000E+08
1.000     1.000     0.000     0.000     1.000     1.000     1.000     1.000   
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              CONTROL OPTIONS                                 $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*CONTROL_ENERGY
         2         2         2         2
*CONTROL_SHELL
  20.0             1        -1         1         2         2         1
*CONTROL_TIMESTEP
    0.0000    0.9000         0  0.00      0.00   
*CONTROL_TERMINATION
  30.0             0   0.00000   0.00000   0.00000
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                                 TIME HISTORY                                 $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_RCFORC
0.2000E-01
*DATABASE_BINARY_D3PLOT
0.4000E-01
*DATABASE_BINARY_D3THDT
0.2000E-01
------------------------------------
*INITIAL_VELOCITY_GENERATION
         4         3  0.00    -0.600E-02  0.00      0.00   
  0.00      0.00      0.00      1.00      1.00      1.00             0

ken731101

頂一個  請幫幫我

ken731101

有沒有版主或其他牛人可以幫幫我
給我一點點提示也行

ken731101

有沒有哪位大大可以幫幫我
感恩不盡

sccxhrq

你这个没有看见在接触的地方加接触啊……这样也可以么，求交流，最近也在做切削的温度耦合，遇到的问题很多啊