切削熱流固偶合計算中斷

ken731101

請問各位牛大 最近再作一個切削加工的題目
目前正在做有關溫度的地方
但每次加入熱指令後...計算只能跑到2秒就停止了  而且溫度上升很快
但是工件有溫度可是刀具卻沒有溫度的上昇
我所使用的是solid 164   *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE
                                      *CONTACT_ERODING_SURFACE_TO_SURFACE
然後使用J-C方程式進行模擬
以下附上我的.K檔  請各位牛大幫幫我哪裡有錯!  感恩
*KEYWORD
*TITLE
                                                                                
$
*CONTROL_SOLUTION
         2
*CONTROL_THERMAL_SOLVER
        1         1         0         0         0         1       0.9
*CONTROL_THERMAL_TIMESTEP
         0       0.0     0.007     0.007         0         0
*CONTROL_CONTACT
         0         2         0         0         0         2         0
         0         0         0         0         0                 
         0         0         0         0         0
*INITIAL_TEMPERATURE_SET
         0        20                                                      
*DATABASE_FORMAT
         0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                               NODE DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*NODE
     $
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             SECTION DEFINITIONS                              $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*SECTION_SOLID
         1         1
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             MATERIAL DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*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,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,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.0,
  0.00  
$
*EOS_GRUNEISEN
3,0.47,1,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.134E-5,-9.0,3.0,0.0,1.5,
  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
0.234,60E-15
*BOUNDARY_PRESCRIBED_MOTION_SET
4,1,-0.600E-02,1,2.0,0,0,0
*DEFINE_CURVE
         1  
  0.000000000000E+00 -0.600E-02
  1.084000000000E+00 -0.600E-02
  2.084500000000E+00 -0.600E-02  
  3.084500000000E+00 -0.600E-02
  4.084500000000E+00 -0.600E-02
  5.084500000000E+00 -0.600E-02
  6.084500000000E+00 -0.600E-02
  7.084500000000E+00 -0.600E-02
  8.084500000000E+00 -0.600E-02
  9.084500000000E+00 -0.600E-02
  1.084500000000E+01 -0.600E-02
  2.084500000000E+01 -0.600E-02
  3.084500000000E+01 -0.600E-02
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              PARTS DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
*PART
Part          1 for Mat         1 and Elem Type         1
         1         1         1         1         0         0         0
$
*PART
Part          2 for Mat         2 and Elem Type         1
         2         1         2         2         0         0         0
$
*PART
Part          3 for Mat         3 and Elem Type         1
         3         1         3         3         0         0         0
$
*PART
Part          4 for Mat         4 and Elem Type         1
         4         1         4         4         0         0         0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             ELEMENT DEFINITIONS                              $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*ELEMENT_SOLID
*BOUNDARY_SPC_SET
         3         0         0         1         0         0         0         0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              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
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              DATABASE OPTIONS                                $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_EXTENT_BINARY
         0         0         3         1         0         0         0         0
         0         0         4         0         0         0
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                         INITIAL VELOCITY DEFINITIONS                         $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*SET_NODE_LIST
         4     0.000     0.000     0.000     0.000
    *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
*END

歡迎大家討論

汪洋1234

我以前也遇到这样的问题，经常出错。有些时候改变时间步或者重划网格会好一些。
实在不行的话，*CONTROL_THERMAL_SOLVER设置slover关键字为1，但是这样的话求解时间会大大增加

ken731101

請問是要再細分一點嗎?
我記得如果slover 設為1的話  好像有些狀態會看不到

ken731101

有人可以幫幫我嗎??

ken731101

有哪位牛大可以幫我看看K檔嗎

ken731101

我努力頂帖  希望有緣人可以幫幫我

ken731101

請問有人可以在提供我一點意見嗎

caat

我考虑温度因素  出现的错误提示是
1、negative jacobian in element *****
2、the solution has a non-positive definite stiffness matrix
谢谢  高人指导

ken731101

8# caat

可以請問你所設定的接觸設定是哪種嗎?

caat

9# ken731101
你好 我做的二维的  自适应
用的接触是  contact 2d surface to surface thermal
另外  我的温度也上升  很快 郁闷  好像  对流边界条件  加不进去   温度散不出去。。。

caat

你的问题 解决了吗 就是温升过快  
还有  我觉得用单元失效  算固热耦合  是不给算的  是不是应为  算温度用隐式算法   单元一失效  就不行了啊  高人指教

ken731101

10# caat 請問可以給我看一下
*contact 2d surf to surf thermal  的卡片參數嗎?
感恩

咯皮

你的问题 解决了吗 就是温升过快  
还有  我觉得用单元失效  算固热耦合  是不给算的  是不是应为  算温度用隐式算法   单元一失效  就不行了啊  高人指教
caat 发表于 2009-3-28 09:56

我也觉得是这原因。

ken731101

那有辦法讓單元不失效嘛?