SPG contact

[aleksei_babichev 0]

Прошу помощи, уже не знаю, что и как настроить. Между двумя телами не проходит контакт, вернее сперва вроде увязывается, но потом все равно проникновение идет. Может погоняете у себя что получится....

3_7_2.zip

[fzbm 44]

А что за время окончания у вас такое ?

размер модели 450 метров? Может сначала поиграться на модели с нормальными размерами и временем. Может элементы попробовать 8-и узловые...

Нет 450 км )

Просто обычно ls-dyna для быстрых процессов, скажем 0,1 секунда. А у вас 100 тыс. лет время окончанию расчета. Раньше не встречал такую постановку. Поэтому вопросы )

Вы статейки смотрели? Кто то в ls-dyna решал такие задачки?

[AlexKaz 1 844]

Зачем SPG - не понятно. Здесь отрабатывает обычный section_solid #1.

Скрытый текст

*KEYWORD
*TITLE
$#                                                                         title
LS-DYNA keyword deck by LS-PrePost
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_ID
$#     cid                                                                 title
         7                                                                      
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         3         6         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
      0.01       0.0       0.0       0.0      20.0         0       0.01.00000E20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
       5.0       5.0       0.0       0.0       1.0       1.0       1.0       1.0
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_ID
$#     cid                                                                 title
         8                                                                      
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         3        10         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
      0.01       0.0       0.0       0.0      20.0         0       0.01.00000E20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
       1.0       1.0       0.0       0.0       1.0       1.0       1.0       1.0
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_ID
$#     cid                                                                 title
         9                                                                      
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         6        10         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
      0.01       0.0       0.0       0.0      20.0         0       0.01.00000E20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
       1.0       1.0       0.0       0.0       1.0       1.0       1.0       1.0
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_ID
$#     cid                                                                 title
        10                                                                      
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         3        12         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
      0.01       0.0       0.0       0.0      20.0         0       0.01.00000E20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
       5.0       5.0       0.0       0.0       1.0       1.0       1.0       1.0
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_ID
$#     cid                                                                 title
        11                                                                      
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         6        13         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
      0.01       0.0       0.0       0.0      20.0         0       0.01.00000E20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
       1.0       1.0       0.0       0.0       1.0       1.0       1.0       1.0
*CONSTRAINED_GLOBAL
$#      tc        rc       dir         x         y         z       tol   
         3         7         3       0.0       0.0   45000.0       0.0
         3         7         3       0.0       0.0       0.0       0.0
         3         7         3       0.0       0.0   15000.0       0.0
         3         7         3       0.0       0.0   60000.0       0.0
         2         7         2       0.0       0.0       0.0       0.0
*MAT_PLASTIC_KINEMATIC_TITLE
mat
$#     mid        ro         e        pr      sigy      etan      beta    
         1    3000.01.00000E11      0.251.000000E9       0.0       0.0
$#     src       srp        fs        vp  
       0.0       0.0       0.0       0.0
*MAT_RIGID_TITLE
Ridg
$#     mid        ro         e        pr         n    couple         m     alias
         2    1000.01.00000E12       0.3       0.0       0.0       0.0          
$#     cmo      con1      con2    
       1.0         4         7
$#lco or a1        a2        a3        v1        v2        v3  
       0.0       0.0       0.0       0.0       0.0       0.0
*PART
$#                                                                         title
LeftPart
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         3         7         1         0         1         0         0         0
*SECTION_SOLID
$#   secid    elform       aet    unused    unused    unused    cohoff    unused
         7         1         0                                       0          
*HOURGLASS
$#    hgid       ihq        qm       ibq        q1        q2    qb/vdc        qw
         1         1       0.5         0       1.9       0.6       0.1       0.1
*PART
$#                                                                         title
RightPart
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         6         7         1         0         1         0         0         0
*PART
$#                                                                         title
Bottom
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
        10         7         1         0         1         0         0         0
*PART
$#                                                                         title
Left
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
        12         7         1         0         1         0         0         0
*PART
$#                                                                         title
Right
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
        13         7         1         0         1         0         0         0
*DAMPING_GLOBAL
$#    lcid    valdmp       stx       sty       stz       srx       sry       srz
         02.50000E-6       0.0       0.0       0.0       0.0       0.0       0.0
*DATABASE_BINARY_D3PLOT
$#      dt      lcdt      beam     npltc    psetid      
1.000000E7         0         0         0         0
*CONTROL_CONTACT
$#  slsfac    rwpnal    islchk    shlthk    penopt    thkchg     orien    enmass
       0.1       0.0         2         0         1         1         1         0
$#  usrstr    usrfrc     nsbcs    interm     xpene     ssthk      ecdt   tiedprj
         0         0        25         1       4.0         0         0         0
$#   sfric     dfric       edc       vfc        th     th_sf    pen_sf      
       0.0       0.0       0.0       0.0       0.0       0.0       0.0
$#  ignore    frceng   skiprwg    outseg   spotstp   spotdel   spothin       
         2         0         0         0         0         0       0.0
$#    isym    nserod    rwgaps    rwgdth     rwksf      icov    swradf    ithoff
         0         0         1       0.0       1.0         0       0.0         0
$#  shledg    pstiff    ithcnt    tdcnof     ftall    unused    shltrw    igactc
         0         0         0         0         0                 0.0         0
*CONTROL_ENERGY
$#    hgen      rwen    slnten     rylen     irgen     
         2         2         2         2         2
*CONTROL_OUTPUT
$#   npopt    neecho    nrefup    iaccop     opifs    ipnint    ikedit    iflush
         1         0         0         0       0.0         0       100      5000
$#   iprtf    ierode   tet10s8    msgmax    ipcurv      gmdt   ip1dblt      eocs
         0         0         2        50         0       0.0         0         0
$#   tolev    newleg    frfreq     minfo    solsig    msgflg    cdetol      
         2         0         1         0         0         0      10.0
$# phschng    demden   icrfile   spc2bnd    penout    shlsig   hisnout    engout
         0         0         0         0         0         0         0         0
$#    insf    isolsf      ibsf      issf    mlkbag      
         0         0         0         0         0
*CONTROL_TERMINATION
$#  endtim    endcyc     dtmin    endeng    endmas     nosol     
3.00000E12   10000001.00000E-6       0.0    1000.0         0
*CONTROL_TIMESTEP
$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st
  100000.0      0.45         0       0.0  600000.0         0         0         0
$#  dt2msf   dt2mslc     imscl    unused    unused     rmscl    unused      ihdo
       0.0         0         0                           0.0                   0
*LOAD_GRAVITY_PART
$#     pid       dof        lc     accel      lcdr      stga      stgr    
         3         2         2      5.81         0         0         0
         6         2         2      5.81         0         0         0

 

Расчёт желательно перевести в implicit, вместо постоянного шага по времени dt подставить таблицу *define_curve шагов по времени чтобы не ждать неделю.

[aleksei_babichev 0]

28 минут назад, fzbm сказал:

А что за время окончания у вас такое ?

размер модели 450 метров? Может сначала поиграться на модели с нормальными размерами и временем. Может элементы попробовать 8-и узловые...

Нет 450 км )

Просто обычно ls-dyna для быстрых процессов, скажем 0,1 секунда. А у вас 100 тыс. лет время окончанию расчета. Раньше не встречал такую постановку. Поэтому вопросы )

Вы статейки смотрели? Кто то в ls-dyna решал такие задачки?

Не решали таких задач, не вижу разницы особой в шаге если нет вопросов устойчивости счета, а их пока нет, контакт не работает, а не счет. В этом то и вопрос.... Дальше если контакт будет отрабатываться там уже буду смотреть.

6 минут назад, AlexKaz сказал:

Зачем SPG - не понятно. Здесь отрабатывает обычный section_solid #1.

  Показать содержимое

*KEYWORD
*TITLE
$#                                                                         title
LS-DYNA keyword deck by LS-PrePost
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_ID
$#     cid                                                                 title
         7                                                                      
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         3         6         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
      0.01       0.0       0.0       0.0      20.0         0       0.01.00000E20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
       5.0       5.0       0.0       0.0       1.0       1.0       1.0       1.0
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_ID
$#     cid                                                                 title
         8                                                                      
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         3        10         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
      0.01       0.0       0.0       0.0      20.0         0       0.01.00000E20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
       1.0       1.0       0.0       0.0       1.0       1.0       1.0       1.0
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_ID
$#     cid                                                                 title
         9                                                                      
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         6        10         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
      0.01       0.0       0.0       0.0      20.0         0       0.01.00000E20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
       1.0       1.0       0.0       0.0       1.0       1.0       1.0       1.0
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_ID
$#     cid                                                                 title
        10                                                                      
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         3        12         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
      0.01       0.0       0.0       0.0      20.0         0       0.01.00000E20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
       5.0       5.0       0.0       0.0       1.0       1.0       1.0       1.0
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_ID
$#     cid                                                                 title
        11                                                                      
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         6        13         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
      0.01       0.0       0.0       0.0      20.0         0       0.01.00000E20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
       1.0       1.0       0.0       0.0       1.0       1.0       1.0       1.0
*CONSTRAINED_GLOBAL
$#      tc        rc       dir         x         y         z       tol   
         3         7         3       0.0       0.0   45000.0       0.0
         3         7         3       0.0       0.0       0.0       0.0
         3         7         3       0.0       0.0   15000.0       0.0
         3         7         3       0.0       0.0   60000.0       0.0
         2         7         2       0.0       0.0       0.0       0.0
*MAT_PLASTIC_KINEMATIC_TITLE
mat
$#     mid        ro         e        pr      sigy      etan      beta    
         1    3000.01.00000E11      0.251.000000E9       0.0       0.0
$#     src       srp        fs        vp  
       0.0       0.0       0.0       0.0
*MAT_RIGID_TITLE
Ridg
$#     mid        ro         e        pr         n    couple         m     alias
         2    1000.01.00000E12       0.3       0.0       0.0       0.0          
$#     cmo      con1      con2    
       1.0         4         7
$#lco or a1        a2        a3        v1        v2        v3  
       0.0       0.0       0.0       0.0       0.0       0.0
*PART
$#                                                                         title
LeftPart
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         3         7         1         0         1         0         0         0
*SECTION_SOLID
$#   secid    elform       aet    unused    unused    unused    cohoff    unused
         7         1         0                                       0          
*HOURGLASS
$#    hgid       ihq        qm       ibq        q1        q2    qb/vdc        qw
         1         1       0.5         0       1.9       0.6       0.1       0.1
*PART
$#                                                                         title
RightPart
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         6         7         1         0         1         0         0         0
*PART
$#                                                                         title
Bottom
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
        10         7         1         0         1         0         0         0
*PART
$#                                                                         title
Left
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
        12         7         1         0         1         0         0         0
*PART
$#                                                                         title
Right
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
        13         7         1         0         1         0         0         0
*DAMPING_GLOBAL
$#    lcid    valdmp       stx       sty       stz       srx       sry       srz
         02.50000E-6       0.0       0.0       0.0       0.0       0.0       0.0
*DATABASE_BINARY_D3PLOT
$#      dt      lcdt      beam     npltc    psetid      
1.000000E7         0         0         0         0
*CONTROL_CONTACT
$#  slsfac    rwpnal    islchk    shlthk    penopt    thkchg     orien    enmass
       0.1       0.0         2         0         1         1         1         0
$#  usrstr    usrfrc     nsbcs    interm     xpene     ssthk      ecdt   tiedprj
         0         0        25         1       4.0         0         0         0
$#   sfric     dfric       edc       vfc        th     th_sf    pen_sf      
       0.0       0.0       0.0       0.0       0.0       0.0       0.0
$#  ignore    frceng   skiprwg    outseg   spotstp   spotdel   spothin       
         2         0         0         0         0         0       0.0
$#    isym    nserod    rwgaps    rwgdth     rwksf      icov    swradf    ithoff
         0         0         1       0.0       1.0         0       0.0         0
$#  shledg    pstiff    ithcnt    tdcnof     ftall    unused    shltrw    igactc
         0         0         0         0         0                 0.0         0
*CONTROL_ENERGY
$#    hgen      rwen    slnten     rylen     irgen     
         2         2         2         2         2
*CONTROL_OUTPUT
$#   npopt    neecho    nrefup    iaccop     opifs    ipnint    ikedit    iflush
         1         0         0         0       0.0         0       100      5000
$#   iprtf    ierode   tet10s8    msgmax    ipcurv      gmdt   ip1dblt      eocs
         0         0         2        50         0       0.0         0         0
$#   tolev    newleg    frfreq     minfo    solsig    msgflg    cdetol      
         2         0         1         0         0         0      10.0
$# phschng    demden   icrfile   spc2bnd    penout    shlsig   hisnout    engout
         0         0         0         0         0         0         0         0
$#    insf    isolsf      ibsf      issf    mlkbag      
         0         0         0         0         0
*CONTROL_TERMINATION
$#  endtim    endcyc     dtmin    endeng    endmas     nosol     
3.00000E12   10000001.00000E-6       0.0    1000.0         0
*CONTROL_TIMESTEP
$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st
  100000.0      0.45         0       0.0  600000.0         0         0         0
$#  dt2msf   dt2mslc     imscl    unused    unused     rmscl    unused      ihdo
       0.0         0         0                           0.0                   0
*LOAD_GRAVITY_PART
$#     pid       dof        lc     accel      lcdr      stga      stgr    
         3         2         2      5.81         0         0         0
         6         2         2      5.81         0         0         0

 

Расчёт желательно перевести в implicit, вместо постоянного шага по времени dt подставить таблицу *define_curve шагов по времени чтобы не ждать неделю.

Он не считает неделю, со временем счета все нормально. Конечно потом перейду на неявный решатель да еще и статику, но пока проблема с контактом. 

3 минуты назад, aleksei_babichev сказал:

Не решали таких задач, не вижу разницы особой в шаге если нет вопросов устойчивости счета, а их пока нет, контакт не работает, а не счет. В этом то и вопрос.... Дальше если контакт будет отрабатываться там уже буду смотреть.

Он не считает неделю, со временем счета все нормально. Конечно потом перейду на неявный решатель да еще и статику, но пока проблема с контактом. 

SPG для того что бы два тела когда деформировались не нужна была перестройка, не хотелось пока бы с ней еще разбираться я новичок совсем, а они сильно деформируются на дальних этапах процесса. PS текс у вас прикреплен там что-то исправлено и можно вставить в сам файл *.к и посчитать?

[AlexKaz 1 844]

Сравните карты *SECTION_, *PART, *CONTACT_, *CONTROL_TIMESTEP, *CONTROL_ENERGY.

Отдельно *HOURGLASS_, *DAMPING.

Всё остальное осталось прежним.

24 минуты назад, aleksei_babichev сказал:

они сильно деформируются на дальних этапах процесса

Разваливается, потому что Вы принудительно выставили большой шаг по времени (бОльший, чем критерий Куранта) и не внесли в систему демпфирование, особенно HOURGLASS. Чтобы не разваливалось, или надо переходить в implicite, или корректировать демпфирование.