基于LS-DYNA的过盈装配的小案例

以轴承和轴的热装配为例，过盈量为2mm，进行装配

基本思路： 采用*CONTACT_SURFACE_TO_SURFACE_INTERFERENCE_ID进行动态松弛预应力装配

动态松弛计算结果

动态松弛过程.gif

关键计算关键字 卡片已经列下下面，有一定dyna基础的就应该自行实现！

*DEFINE_CURVE_TITLE
stiffness over time dynamic relaxation
$: Contact surface 1 : Factor on dyn rel stiffness vs time
$: X axis : Dynamic rel time           (Units: Time)
$: Y axis : Factor on stiffness        (Units: Scalar, no units)
$: Usage: Stress initialization
$:   label      sidr       sfa       sfo      offa      offo    dattyp
         4         1       0.0       0.0       0.0       0.0         0
$:             xvals               yvals
                 0.0                 0.0
              0.0003          1.0000000
                0.01          1.0000000
$
*DEFINE_CURVE_TITLE
stiffness over time dynamic relaxation
$:   label      sidr       sfa       sfo      offa      offo    dattyp
         5         0       0.0       0.0       0.0       0.0         0
$:             xvals               yvals
                 0.0           1.0000000
               0.010           1.0000000         
$
*CONTACT_SURFACE_TO_SURFACE_INTERFERENCE_ID
$:   label       title
         1contact interference                                                 
$:    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         1        2         3         3         0         0         1         1
$:      fs        fd        dc        vc       vdc    penchk        bt        dt
       0.2       0.0       0.0       0.0      20.0         0       0.0       0.0
$:     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
       0.0       0.0       0.0       0.0       0.0       0.0       0.0       0.0
$:   lcid1     lcid2
         4         5
$

如果实在有困难，请查看如下k文件。