Хитрый Buckling

[Борман 2 376]

Чото не соображу, как сделать LB для такой задачи:

Найди крит. сжимающую силы для трубы, находящейся под подстоянным внутренним давлением ?

 

 

Общий случай: Найти тот самый коэффициент для силы F1 для случая нагружения (F1, F2=const) в разные узлы.

[Vova 419]

Чото не соображу, как сделать LB для такой задачи:

Найди крит. сжимающую силы для трубы, находящейся под подстоянным внутренним давлением ?

 

 

Общий случай: Найти тот самый коэффициент для силы F1 для случая нагружения (F1, F2=const) в разные узлы.

Я бы сделал так:

1. шаг приложил давление

2. шаг приложил перемещение, и смотрел бы силу, при которой труба потеряет устойчивость.

[Fedor 1 597]

У Басова в http://free-books.me/comp/cad_cae/ansys/67558-k-a-basov-ansys-v-primerax-i-zadachax.html   есть пример трубки под наружным давлением. Знак поменять и приложить дополнительную силу. Так то почему нельзя сделать ?  

 

 

       !______________________Решение____________________________________
       allsel
       n=5 ! число мод для проверки устойчивости, если 0 статический анализ **************
       /SOLU   
           ANTYPE,STATIC             ! STATIC ANALYSIS
           *IF,n,GT,0,THEN  $ PSTRES,ON   $   *ENDIF      ! CALCULATE PRESTRESS EFFECTS
           SOLVE
       FINISH

       *IF,n,GT,0,THEN           ! проверка на устойчивость
         ! для строительных конструкций не ниже 1.36 коэффициент запаса по устойчивости (следует из СНИП)    
           /SOLU
              ANTYPE,BUCKLE            
              BUCOPT,LANB,n             ! USE BLOCK LANCZOS EIGENVALUE EXTRACTION METHOD
              MXPAND,n                 
              SOLVE
           FINISH
       *ENDIF

[Борман 2 376]

Мужики, акутально. Кумекаем.

[soklakov 1 475]

Итерациями?

[Борман 2 376]

Итерациями

Чот не очень понял.

[kristeen 355]

Мопед не мой.

[soklakov 1 475]

Чот не очень понял.

Поскольку внутреннее давление постоянное, не увеличивается вместе со сжимающей нагрузкой, то вариант найти коэффициент не  подходит. Поэтому итеративно подбирается такая сжимающая нагрузка, чтобы коэффициент в Buckling'е был равен единице. Это и будет критическая сжимающая нагрузка.

Или вопрос в другом?

[JohnnyLuna82 0]

 

Чот не очень понял.

Поскольку внутреннее давление постоянное, не увеличивается вместе со сжимающей нагрузкой, то вариант найти коэффициент не  подходит. Поэтому итеративно подбирается такая сжимающая нагрузка, чтобы коэффициент в Buckling'е был равен единице. Это и будет критическая сжимающая нагрузка.

Или вопрос в друго

Дело в том что LB учитывает всю приложенную нагрузку в посл шаге-и силы и собств вес и т д-как бы рассплитовать её

[soklakov 1 475]

Дело в том что LB учитывает всю приложенную нагрузку в посл шаге-и силы и собств вес и т д-как бы рассплитовать её

Вы прочитали тему-то? Картинка Вам понятна?

[kol 12]

:)

Такая мысль - задаёмся соотношением давления и силы например у меня 17МПа а сила 100000Н =170 ну или 1/170...

тогда если единичная сила 1Н то давление 1/170 = 0,00588Па (это и прикладываем) и для этого определяем коэффициент. Ща буду пробовать :) ток на чём проверить? в ансисе примера не нашел, а без палочки выручалочки страшно :(

 

близко к истине но пока все : https://cloud.mail.ru/public/2SRLDqrVkDad/buckling.dat

[soklakov 1 475]

Такая мысль - задаёмся соотношением давления и силы например у меня 17МПа а сила 100000Н =170 ну или 1/170... тогда если единичная сила 1Н то давление 1/170 = 0,00588Па (это и прикладываем) и для этого определяем коэффициент. Ща буду пробовать :) ток на чём проверить? в ансисе примера не нашел, а без палочки выручалочки страшно :(

Еще раз смотрим на картинку выше. Добавим к ней пояснений, если она кажется непонятной.

 

Задача 1. Определить критическую силу при давлении 17 МПа.

итерация 1. Прикладываете 17МПа и 1 Н. Находите из Bucklinga коэф-т равный Х.

итерация 2. Прикладываете 17МПа и Х Н. Находите из Bucklinga коэф-т равный Y.

итерация 3. Прикладываете 17МПа и Х*Y Н. Находите из Bucklinga коэф-т равный Z.

И так до тех пор, пока не получите коэффициент 1(единица). Приложенная на этом этапе сила и есть критическая при давлении 17МПа (с учетом всех линейных предположений и прочих упрощений)

 

Задача 2. Определить критическую силу при давлении 8 МПа.

... аналогично.

[kristeen 355]

what's new in 16.0

[Борман 2 376]

Еще раз смотрим на картинку выше. Добавим к ней пояснений, если она кажется непонятной.

Сходится преступно медленно. Берем сжатый стержень с критической силой 100. Получим ее как 50 - постоянная часть и 50 - переменная.

 

 

Хотя... Уверен, есть способ лучше.

[soklakov 1 475]

Хотя... Уверен, есть способ лучше.

 

Может быть как-то так?

Правда, как-то слишком просто.

Проблема, наверное, в том, что когда это будут давление и сила, их сложнее складывать.

[Борман 2 376]

Проблема, наверное, в том, что когда это будут давление и сила, их сложнее складывать.

Да уж. Так как я предлагал (ну в смысле по вашему первоначальному алгоритму) главное угадать первое приближение, тогда нормально сойдется.

 

А угадать его можно, если решать без давления. 

[Борман 2 376]

Правда, как-то слишком просто.

Прикольно мозги запудрил. Как кто то давно пытались в одной теме разобраться почему x^-1 убывает нелинейно. Раз степень отрицательная - значит убывает, раз степень единица - значит линейно. Пол дня думали.... Я его потом чуть не забанил со злости.

[soklakov 1 475]

главное угадать первое приближение, тогда нормально сойдется.

Еще вопрос к желаемой точности. В Вашем варианте справа уже на шестой итерации точность в пределах 1%.

А еще, если итерации выполнять вручную, то можно искусственную подкрутку делать. То есть если видим, что тренд к уменьшению, то уменьшили чуть больше, чем алгоритм предполагает.

[Борман 2 376]

ISPA точно знает, как это проще сделать...

[kristeen 355]

ISPA точно знает

Призовете на свою голову..