Jump to content

SPG contact

aleksei_babichev

By aleksei_babichev,
in LS-Dyna

 Share Followers 0

Recommended Posts

aleksei_babichev

aleksei_babichev 0

Posted
Posted

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

3_7_2.zip

Link to post
Share on other sites


UnPinned posts
fzbm

fzbm 44

Posted
Posted

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

image.png

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

Нет 450 км )

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

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

Link to post
Share on other sites
AlexKaz

AlexKaz 1,848

Posted
Posted

Зачем 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 шагов по времени чтобы не ждать неделю.

Link to post
Share on other sites
aleksei_babichev

aleksei_babichev 0

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

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

image.png

размер модели 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 текс у вас прикреплен там что-то исправлено и можно вставить в сам файл *.к и посчитать?

Link to post
Share on other sites
AlexKaz

AlexKaz 1,848

Posted
Posted

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

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

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

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

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

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

  • Нравится 2
Link to post
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.
Note: Your post will require moderator approval before it will be visible.

Guest
Reply to this topic...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

×
 Share
Followers 0
Go to topic listing

  • Recently Browsing   0 members

    No registered users viewing this page.



  • Сообщения

    • A_1
      Круговой массив и подпрограмма. Как?

      By A_1 · Posted

      Станок получит список координат осей отверстий по программе HOLES2. Эта программа позиционирует инструмент к центрам отверстий, а изменённая программа L115 с помощью команд TRANS X=<текущее положение по X> Y=<текущее положение по Y> устанавливает в этих точках начало системы координат для расфрезеровки.   Если не хотите модифицировать подпрограмму L115, то можно сделать "по формулам": def int _i, _n=8; количество отверстий def real _r=125; радиус окружности def real _a=45; угол первого отверстия for _i = 1 to _n trans x=_r * cos(_a) y=_r * sin(_a) L115 _a = _a + 360 / _n endfor trans  
    • vik300001
      Lantek Expert

      By vik300001 · Posted

      согласен, лучше пронест, но ситуация...  Нужно лантеком
    • aleksejkolinko
      Размеры из модели в чертёж. Отсутствует буквенное обозначение допусков

      By aleksejkolinko · Posted

      При переносе размеров из модели в чертёж, почему-то не проставляются буквенные допуски. Пытаюсь освоить этот способ простановки размеров. У кого какой опыт и почему проставляются только числовые значения допусков?
    • Bellka
      Автоматическая разблокировка двери после М0. Heidenhain620

      By Bellka · Posted

      На предыдущей стойке это была команда М6
    • Nod801
      Не создается литейная форма! "Сшивка листов не удалась"

      By Nod801 · Posted

      @timyr8693 по вашей ссылке не получается скачать ... выложите в другое место ваши файлы ... посмотрим что там не так
    • Nikomas
      Круговой массив и подпрограмма. Как?

      By Nikomas · Posted

      Из вашей формулы я не получу список координат.   Очень жаль что для TRANS нельзя использовать полярные координаты (TRANS RP=125 AP=45). Было бы в разы проще. Первое отверстие я нашел с помощью: R1 = SQRT(2) / 2 * 125 А вот второе добавить (+45град на радиусе 125мм) я без понятия как.
    • mortinio
      DMU80p Состояние инструмента неизвестно

      By mortinio · Posted

      Получилось сделать бэкап только после того как притащил стационарный комп с Widows7 и установил последнюю версию TNCremo, с ноутбука старенького с Windows xp ни в какую не проходил пароль. https://disk.yandex.by/d/LvR8ZqTK2YnC2w
    • KsunReh
      Не создается литейная форма! "Сшивка листов не удалась"

      By KsunReh · Posted

      У меня 18 солид так что ваш файл проверить не могу. Но по опыту если вы делаете все действия в правильном порядке, то сшивка листов может не получатся из-за "кривых" (грязно построенных) поверхностей разъема либо из-за плохого стыка этих поверхностей с линией разъема. Это актуально для моделей со сложными формами, надо перестраивать еще раз нормально, искать у себя ошибки.   Могу порекомендовать такую самопроверку. 1) вручную выделяете все поверхности детали относящиеся к матрице и разделяющие поверхности, копируете их (инструмент эквидистанта к поверхности со значением 0) 2) тоже самое для пуансона 3) у вас должно получится 2 поверхности без глюков/разрывов и т.д., проверяете их. Если все ОК, то и автоматическое создание формы должно работать, если не ОК, то исправляете проблемные места. 
    • Andrew_kos
      Полиамид Pa6-gf30

      By Andrew_kos · Posted

      А какова вероятность что характеристики еще можно откопать? Был бы очень благодарен если бы вы написали их.
    • Bellka
      Полное отключение осей на 426

      By Bellka · Posted

      Здравствуйте подскажите пожалуйста как сделать полное программное отключение b и с осей на 426 стоке, 430 и остальные предполагаю будут идентично Ситуация такая, станок на обслуживание частично разбирается, но на включать и тестировать при изменении только  MP 10      : 31         ;Aktive Achsen на 10111 ось С пропадает из системы и не требует разводки, но если физически отключить мотор появляются ошибки какие ещё параметры нужно поменять что бы выпилить полностью?
×
×
  • Create New...