齿轮与润滑油耦合分析出现的问题！

henry_gigi

大家好，采用fem-sph耦合方法做齿轮与润滑油的耦合分析时，出现如下现象（和想象的不一样）：（1）油被带起来后冲击邮箱反弹到主动轮才进入啮合区；（2）被带起一部分油后，后续的油没有想象中的应该在重力的作用下进入搅动区，感觉油的流动性很差；（3）轮齿切入处（图中右下方）的油被带了起来，并与箱壁分离。请大家帮忙分析这些的原因及正确的现象是什么样的？非常感谢！

K文件太大无法上传，除节点和单元外的关键字如下：kg-mm-ms
*DEFINE_COORDINATE_NODES
         1    262008    261925    261973         0         
*CONTROL_SOLUTION
         0
*CONTROL_TERMINATION
    1.0e+2            
*CONTROL_TIMESTEP
         0  0.800000
     0.000         0         0
*CONTROL_THERMAL_SOLVER
         1         1         1     0.000         0       1.0  1.00e-10
*CONTROL_THERMAL_NONLINEAR
         0         0       0.5
*CONTROL_THERMAL_TIMESTEP
         1   1.00000    1.0E-3                        5.0e+2                    
*CONTROL_SPH
         0         0      0.00         3       300   
*CONTROL_CONTACT
     0.100     0.000         0         0         0         0         2         0
         0         0         0         0     0.000         0         0         0
     0.000  0.000000
         1         1
         0         0         0     0.000  1.000000         0  1.000000
         0
*CONTROL_BULK_VISCOSITY
  1.500000  0.060000
*CONTROL_HOURGLASS
         1       0.1
*DATABASE_GLSTAT
       0.1         3
*DATABASE_RCFORC
       0.1         3
*DATABASE_NODOUT
       0.1         3
*DATABASE_TPRINT
       0.1         3
*DATABASE_SLEOUT
       0.1         3
*DATABASE_BINARY_D3PLOT
       0.1
         0
*DATABASE_BINARY_D3THDT
       0.1
*DATABASE_BINARY_INTFOR
       0.1
*DATABASE_EXTENT_BINARY
         0         0         3
         0         0         0         0         0         0         0         0
         0         0     0.000         0         0         0                    
*DATABASE_HISTORY_SPH_SET
         1
*MAT_PLASTIC_KINEMATIC
         1 7.8000E-6 2.1000E+2  0.300000     0.608    21.000     0.000  
         
$*MAT_ELASTIC
$         1 7.8060E-6 2.0400E+2  0.300000     0.000     0.000     0.000         
*MAT_RIGID                                                                     
         2 7.8060E-6 2.0400E+2  0.300000     0.000     0.000     0.000         
  1.000000         7         7
     0.000     0.000     0.000     0.000     0.000     0.000
*MAT_RIGID                                                                     
         3 7.8500E-6 2.0600E+2  0.280000     0.000     0.000     0.000         
   -1.0000         1    111110
     1.000     0.000     0.000     0.000     0.000     0.000
*MAT_RIGID                                                                     
         4 7.8500E-6 2.0600E+2  0.280000     0.000     0.000     0.000         
  1.000000         7         4
     0.000     0.000     0.000     0.000     0.000     0.000
*MAT_NULL
         5 0.9300E-6    0.0000     0.000   0.00000  0.000000 0.0000000
*EOS_GRUNEISEN
         5  1530.000  1.750000                          0.28

*INITIAL_TEMPERATURE_SET
         0  20.00000
*MAT_THERMAL_ISOTROPIC
         6 7.8000E-6
     490.1   47.7e-6   
*PART

         1         3         2         0         0         0         0         6
*SECTION_SHELL
         3
  2.000000  2.000000  2.000000  2.000000
*PART

         2         1         3         0         0         0         0         6
*SECTION_SOLID
         1
*PART

         3         1         1         0         0         0         0         6
*PART

         4         1         4         0         0         0         0         6
*PART

         5         1         1         0         0         0         0         6
*PART

         6         1         1         0         0         0         0         6
*PART

         7         1         1         0         0         0         0         6
*PART

         8         2         5         5         0         0         0         6
*SECTION_SPH
         2  1.050000  0.200000  2.000000     0.0001.0000E+20
*SET_PART_LIST
         1
         2         3         4         5         6         7         
*DEFINE_CURVE
         1         0  1.000000  1.000000
           00.000000           0.0000000
           2.5000000           0.1007000
           5.0000000           0.1861000
           7.5000000           0.2432000
           9.0000000           0.2600300
          10.0000000           0.2633300
          1000.00e+0           0.2640000
*BOUNDARY_PRESCRIBED_MOTION_RIGID_LOCAL
         2         7         0         1  1.000000
*DEFINE_CURVE
         2         0  1.000000  1.000000
$           00.000000           0.0000000
               0.000           3.5000E+2
         1000.000000           3.5000E+2
*LOAD_RIGID_BODY
         4         7         2       1.0         0
$*INITIAL_VELOCITY_RIGID_BODY
$         4         0         0         0         0         0   -0.0463
$*INITIAL_VELOCITY_RIGID_BODY
$         2         0         0         0         0         0    0.2633
$*LOAD_RIGID_BODY
$         2         7         2    0.1755         1
*DEFINE_CURVE
         3
               0.000            9.810e-3
          5000.00000            9.810e-3
*SET_PART_LIST
         2
         8  
*LOAD_BODY_Y
         3  1.000000
*LOAD_BODY_PARTS
         2
*CONTACT_ERODING_NODES_TO_SURFACE
         1         1         4         0         0         0         0         0
       0.0       0.0       0.0       0.0       0.0         1
       1.0       1.0       0.0       0.0       1.0       1.0       1.0       1.0         
         0         1         0
*CONTACT_AUTOMATIC_NODES_TO_SURFACE
         1         1         4         2         0         0         0         0
       0.0       0.0       0.0       0.0       0.0         1
       1.0       1.0       0.0       0.0       1.0       1.0       1.0       1.0         
         0         1         0
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_THERMAL_ID
         1                                                                     
         7         3         3         3         0         0         1         1
     0.002     0.002     0.000    3.0e-1    20.000         0     0.000     0.000
     0.000     0.000     0.000     0.000     0.000     0.000     0.000     0.000
    0.0000    0.0000   4.77e-6  0.100000  1.000000       1.0


         1         1
*BOUNDARY_SPC_SET
         2         0         0         0         1         1         1         1
*BOUNDARY_SPC_SET
         1         0         0         0         1         1         1         1