Jump to content

SPG contact


aleksei_babichev

Recommended Posts

aleksei_babichev

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

3_7_2.zip

Link to post
Share on other sites


UnPinned posts

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

image.png

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

Нет 450 км )

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

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

Link to post
Share on other sites
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 шагов по времени чтобы не ждать неделю.

Link to post
Share on other sites
aleksei_babichev
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

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

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

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

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

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

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

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

Create an account or sign in to comment

You need to be a member in order to leave a comment

Create an account

Sign up for a new account in our community. It's easy!

Register a new account

Sign in

Already have an account? Sign in here.

Sign In Now
  • Recently Browsing   0 members

    No registered users viewing this page.

  • Сообщения

    • zerganalizer
    • areg
      и хрен что докажешь!!! - нашли исполнителя. Кончается это тем что сделают как могут. Потом 3 штуки проверяют, если повезет примут, если нет 3 шт из все партии переделают и будут клянчить, чтобы приняли. Лично сам кольцо Ф70 толщиной 0,6 +/- 0,01 точил. Кстати выточил. На мой вопрос а почему нельзя было обернуть калиброванной жестью неполный оборот был ответ - НАДА КОЛЬЦО!!! Майор штаба работающий менеджером приемщиком, он еще на сверлильном фрезеровкой занимался и был этим очень горд (рядом фрезерный с универсальной головкой простаивал, но он просто не знал с какой стороны к нему подойти, а спросить стеснялся). Пардон за флуд. БО НАКИПЕЛО!!!!
    • zerganalizer
      И снова привет, я с извинениями... ввёл в заблуждение общественные массы. После проверки циферок при дальнейшем расчёте я понял, что посчитал НЕ ТУ ТОЧКУ на графике. Ибо 0,13 Маха на хорде 1 метр даёт почти 3 млн. Re. На картинку ранее внёс правку - чего я не там посчитал: Получается путём прикидок по картинке - SU2 показал мой выигрыш более, чем XFLR5 - вместо 200/110=1.8 раза - 2,6 раза!   Обтекание показывает так (я чутка контраст поднял), профиль - GAW (NASA LANGLEY LS(1)-0417 (GA(W)-1) AIRFOIL (ls417-il)).   Свою геометрию не показываю по понятным причинам, но она сильно другая, хотя и ламинизированная по верхней дужки похожим образом.   Я продолжил компудахтать SU2 - выбрал другие 2 точки - Re=2e7 (20 млн) и два угла атаки - 2 градуса, где я в выигрыше и 8 градусов, где XFLR5 слил мой профиль процентов на 30. Сами смотрите: На сей раз я вычленил из всего просчёта на всех Re только нужное. SU2 накомпудахтал в невыгодной позиции в 8 градусов угла атаки прирост (а не проигрыш) 1,39 раза, почти 40%. А в выгодной позиции - на вершине "туннельного эффекта Bartini-Panguver" аж в целых пять (с копейками) раз!!!   Экстраполируя эти данные - я примерно от руки накалялкал ожидаемый туннель по SU2 расчётам.   Мой и оригинальный NASA LANGLEY LS(1)-0417 (GA(W)-1) AIRFOIL (ls417-il) считались на одном конфиге, решатель RANS с вязкостью. Уже по первым ста итерациям сразу стало понятно, насколько мой профиль выигрывает у NASA LANGLEY LS(1)-0417 (GA(W)-1) AIRFOIL (ls417-il), а итераций было 2000.   Так что глупый XFLR5 тем больше занижает мои показатели "туннельного эффекта", чем больше Re. На 3-х млн - раза в полтора, а на 20млн. - в разы.   Теперь вот проверить всё это методом наложения поляр XFLR5 и SU2 на тех же режимах можно будет только после пакетного анализа на SU2, который в прошлый раз не запустился из-за отсутствия питона на компе.   Всем спасибо! Всем приятного флейма! Пошёл за попкорном...
    • Killerchik
      Управляться-то она может по-разному. Самый правильный варик - полноценная подключаемая/отключаемая ось. Но на деле может быть множество других вариантов разной степени колхозности.
    • Snake 60
      @Бестолковый Скинь чертеж с моделькой в личку, посмотрю...
    • Бестолковый
      1. Кривая (ломаная) ОС. 2. Не менее кривое (также ломаное) ПО. 3. Видеокарта. Немного не та что нужна для работы данного ПО  или "дрова" не обновлены.
    • pescarbon
      Здравствуйте. На станке с системой БалтСистем всегда работает импульсная смазка.,хотя должна через промежутки времени. Есть код для включения и отключения смазки. Если кто знает, подскажите пожалуйста что в коде не так
    • Бестолковый
    • gudstartup
      а почему без поддержки поддержка возможно есть только пнр не проведено и поэтому люди не знают как сней быть.
    • Бестолковый
      Есть деталь. Банальная полускоба. Скрины ниже. При создании чертежа из неё "улетают" размерные лини (скрины ниже). И не только в данной детали, а весьма часто в чертежах. SW 2022 SP5.0 Вопрос: почему подобное может происходить?
×
×
  • Create New...