Перейти к публикации

И снова о потере устойчивости...


Рекомендованные сообщения

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

Ссылка на сообщение
Поделиться на других сайтах


10 часов назад, Борман сказал:
Показать содержимое  

file004.jpg.639864355e3368f666db90849cb6aaac.jpg5c060cf42ced2_.jpg.1b383e935b925e0d500ee3386ee74dab.jpg.6bd48c6ab08055dac3a201a72f5138cc.jpg

 

Hide  

Совпадение очевидно.

И наш РЕАЛЬНЫЙ скачек энергии оказался равен РАБОТЕ ВИРТУАЛЬНОЙ СИЛЫ, КОТОРАЯ НИКОГДА НЕ БЫЛА ПРИЛОЖЕНА, НА ВИРТУАЛЬНЫХ ПЕРЕМЕЩЕНИЯХ, КОТОРЫЕ НИКОГДА НЕ БЫЛИ ПРОЙДЕНЫ.

Приращение потенциальной энергии не обязано равняться работе внешних сил (количество выпитого кофе не всегда переходи в качество).

 

Истинный путь VS все остальные.JPG

Изменено пользователем AlexKaz
Ссылка на сообщение
Поделиться на других сайтах
3 минуты назад, AlexKaz сказал:

Приращение потенциальной энергии не обязано равняться работе внешних сил (количество выпитого кофе не всегда переходи в качество).

1. Как это ?

2. У меня равняется во второй части поста.

Ссылка на сообщение
Поделиться на других сайтах

Тогда ничего не понял.

Допустим, мы зафиксировали стержень-балку в каждом узле по всем степеням свободы. Задали начальные деформации/напряжения. Потенциальная энергия есть? - да. Работа внешних сил есть? - нет.

Допустим № 2, мы катаемся в тележке на американских синусоидальных y=A*sin(w*t) горках. Двигатель тащит нашу тележку и совершает работу P. Когда мы забираемся на очередной y=A приращение нашей потенциальной энергии dE=A*sin(w*t)-A*sin(w*t-2*pi)=0. Приращение потенциальной энергии есть? - нет. Работа внешних сил есть? - да.

Изменено пользователем AlexKaz
Ссылка на сообщение
Поделиться на других сайтах
Цитата

Приращение потенциальной энергии не обязано равняться работе внешних сил (количество выпитого кофе не всегда переходи в качество).

Кофе не энергетический напиток. Тут нужно что-то покрепче для аналогии. Водку или пиво. А лучше то и другое и побольше :)

Речь всегда о приращении энергии в результате работы. А уж какая там она внутри то Эйнштейн знает . Это уже не механика :)  

Изменено пользователем Fedor
Ссылка на сообщение
Поделиться на других сайтах
19 минут назад, AlexKaz сказал:

Задали начальные деформации/напряжения. Потенциальная энергия есть? - да. Работа внешних сил есть? - нет.

Я смотрю не абсолютное значение, а приращение. Приращение вн. энергии равно работе силы на приращении перемещения. Про начальное состояние согласен, но, честно сказать, в это состояние тоже какая-то сила притащила конструкцию.

 

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

Когда мы забираемся на очередной y=A приращение нашей потенциальной энергии dE=A*sin(w*t)-A*sin(w*t-2*pi)=0. Приращение потенциальной энергии есть? - нет. Работа внешних сил есть? - да.

Если бы на горках не было трения, то по синусу можно было бы кататься без двигателя.

Ссылка на сообщение
Поделиться на других сайтах
1 час назад, Fedor сказал:

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

Линейная матрица жесткости не зависит от перемещений и деформаций. А матрица больших перемещений зависит. Но это не говорит о том, что меняется модуль упругости.:biggrin:

Ссылка на сообщение
Поделиться на других сайтах
46 минут назад, Борман сказал:

Про начальное состояние согласен, но, честно сказать, в это состояние тоже какая-то сила притащила конструкцию.

 

Вы на правильном пути. Внешняя сила остается.:biggrin:

Ссылка на сообщение
Поделиться на других сайтах
35 минут назад, ДОБРЯК сказал:

Линейная матрица жесткости не зависит от перемещений и деформаций. А матрица больших перемещений зависит. Но это не говорит о том, что меняется модуль упругости.:biggrin:

От перемещений то не зависит вроде, зависит от произведений их производных судя по формулам. :)

Ссылка на сообщение
Поделиться на других сайтах
4 минуты назад, Fedor сказал:

От перемещений то не зависит вроде, зависит от произведений их производных судя по формулам. :)

Тогда нужно обновить формулы

В 20.11.2018 в 06:54, ДОБРЯК сказал:

Матрица жесткости элемента является важнейшим понятием в приложении МКЭ в механике.

Формула для матрицы жесткости является общей:

image039.gif       (1.1.4.6)

Матрица жесткости image042.gif имеет размерность image044.gif, где image046.gif - число степеней свободы конечного элемента.

Матрица image020.gif - есть матрица деформирования элемента.

Вектор деформации image012.gif, определяемый через перемещения будет имеет следующую  аппроксимацию 

image014.gif (1.1.4.2)

image016.gif

Здесь: image018.gif - вектор степеней свободы элемента (3xn компонент);

image020.gif - матрица деформирования элемента.

----------------------------------------------------------------------

При больших перемещениях деформации нелинейно зависят от перемещений и матрица [В] = image020.gif + [ВL {image018.gif}].

 [ВL {image018.gif}] - зависит от перемещений.

В формуле 1.1.4.6 image020.gif заменяем на [В] и получаем матрицу

[К] = [Ко]  + [Кп]

[Ко] – обычная матрица жесткости

[Кп] – матрица больших перемещений.

Если добавить s] – матрица начальных напряжений, то получим

полную матрицу тангенциальных жесткостей [Кт] = [Ко] + [Кs] + [Кп]

 

Деформации зависят от перемещений (матрица [В]).

А жесткость зависит от деформаций. Матрица упругости [Е] не меняется.:biggrin:

Ссылка на сообщение
Поделиться на других сайтах
1 час назад, Борман сказал:

Приращение вн. энергии равно работе силы на приращении перемещения.

Увидел ошибку в первом случае: работа внешних сил (приложенных напряжений/деформаций) уравновесится с потенциальной энергией.

 

IMHO, надо смотреть уравнение ([K]+lambda*[Kgeom])*{dU}=0 когда вектор приращений не выкидывается (жёсткое нагружение) и решается обобщённая задача. Раскрываем: [K]*{dU} + lambda*[Kgeom]*{dU}=0; {F}+lambda*[Kgeom]*{dU}=0. Здесь {F} = const  - и это внешние силы.

А вот на этом графике даны внутренние силы что-ли?:

12 часа назад, Борман сказал:

file004.jpg.639864355e3368f666db90849cb6aaac.jpg

Кусок lambda*[Kgeom] зависит от {dU} линейно или нелинейно - это абсолютно по-барабану.

 

Можно предположить что произведение lambda*[Kgeom] далеко не всегда будет иметь один и тот же знак когда задано мягкое нагружение. А можно сказать аналогичное о lambda*[Kgeom]*{dU}? Вроде нет, ведь при жёстком {F} = const. Значит вывод - lambda и [Kgeom] меняют знак вместе. В таком случае можно наблюдать появление и отрицательной силы (как на графике выше).

Но почему при решении ([K]+lambda*[Kgeom])*{dU}=0 при мягком нагружении на графике наблюдается только монотонный или не очень но рост силы? Может потому, что решатель отбрасывает в данном случае отрицательные корни?

Тогда, если решатель реально отбрасывает отрицательные корни, будет происходить следующее:

- при мягком нагружении в ([K]+lambda*[Kgeom])*{dU}=0 lambda всегда положительна, часть корней какая-то редиска от нас скрыла.

- при жёстком нагружении [K]*{dU} + lambda*[Kgeom]*{dU}=0 lambda всегда положительна, но при этом [Kgeom] меняет знак, а значит работа внешней силы может падать.

 

IMHO, нада в своём собственном коде составлять все указанные матрицы и смотреть на количество корней. Только не вижу в этом особого смысла, кроме чисто спортивного интереса =)

Изменено пользователем AlexKaz
Ссылка на сообщение
Поделиться на других сайтах

@AlexKaz 

Обязательно прочитаю вдумчиво. Пока запара на работе. :(

Ссылка на сообщение
Поделиться на других сайтах
3 минуты назад, Борман сказал:

@AlexKaz 

Обязательно прочитаю вдумчиво. Пока запара на работе. :(

В вашем случае решаются совсем другие уравнения. В них нет никакого lambda.:bye:

Ссылка на сообщение
Поделиться на других сайтах
Цитата

Деформации зависят от перемещений (матрица [В]).

Вы хоть различайте функции от их производных. Это же как путать расстояние и скорость  :)

 

Цитата

dU

Наверное все-таки просто U без дифференцирования :)

Изменено пользователем Fedor
Ссылка на сообщение
Поделиться на других сайтах

@ДОБРЯК , ну так подскажите, какие именно. А, ну да, вы жеж только за бабло даёте советы.

Зато флудить за бесплатно - это вы умеете.

Ссылка на сообщение
Поделиться на других сайтах
2 минуты назад, AlexKaz сказал:

@ДОБРЯК , ну так подскажите, какие именно. А, ну да, вы жеж только за бабло даёте советы.

Я же много раз писал это уравнение. 

[Кт] = [Ко] + [Кs] + [Кп]

[Кт] * {U} = {F}

Статика. .

Нет тут никакого lambda.

Чем больше грязи и клоунады, тем меньше советов. Это то хотя бы понятно.:bye:

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

- при мягком нагружении в ([K]+lambda*[Kgeom])*{dU}=0 lambda всегда положительна, часть корней какая-то редиска от нас скрыла.

- при жёстком нагружении [K]*{dU} + lambda*[Kgeom]*{dU}=0 lambda всегда положительна, но при этом [Kgeom] меняет знак, а значит работа внешней силы может падать.

И при мягком и при жестком нагружении решается одно и то же уравнение 

([K]+lambda*[Kgeom])*{U}=0

Это если решать начальную потерю устойчивости. 

Ведь можно же и в одной задаче задать и перемещение и силы.:bye:

Ссылка на сообщение
Поделиться на других сайтах
1 час назад, AlexKaz сказал:

надо смотреть уравнение ([K]+lambda*[Kgeom])*{dU}=0 когда вектор приращений не выкидывается (жёсткое нагружение) и решается обобщённая задача. Раскрываем: [K]*{dU} + lambda*[Kgeom]*{dU}=0; {F}+lambda*[Kgeom]*{dU}=0. Здесь {F} = const  - и это внешние силы.

Забавное раскрытие уравнения И еще более забавно, что Здесь {F} = const  - и это внешние силы.:biggrin:

Комментарии излишни.

Ссылка на сообщение
Поделиться на других сайтах
1 час назад, AlexKaz сказал:

надо смотреть уравнение ([K]+lambda*[Kgeom])*{dU}=0 когда вектор приращений не выкидывается (жёсткое нагружение) и решается обобщённая задача

Не очень понимаю контекст. Здесь вы перемешиваете линейный подход и нелинейный. Не совсем корректно.

2 часа назад, AlexKaz сказал:

[K]*{dU} + lambda*[Kgeom]*{dU}=0

Такого не может быть. Вы утверждаете, что есть "траектория бифуркации". Есть только точка бифуркации.

 

 

 

Ссылка на сообщение
Поделиться на других сайтах
5 минут назад, Борман сказал:

Такого не может быть. Вы утверждаете, что есть "траектория бифуркации". Есть только точка бифуркации.

Да он не понимает, что 

"надо смотреть уравнение ([K]+lambda*[Kgeom])*{dU}=0 когда вектор приращений не выкидывается (жёсткое нагружение) и решается обобщённая задача"

вектор приращений перемещений не может не выкидываться. :bye:

Ссылка на сообщение
Поделиться на других сайтах
16 минут назад, Борман сказал:

Такого не может быть. Вы утверждаете, что есть "траектория бифуркации". Есть только точка бифуркации.

Есть упругие перемещения, которые решаются на первом этапе баклинга {dU} - dот их я и имею в виду.

 

Продолжаю рассуждать. Почему при мягком нагружении не существует левой половины графика выше? Потому что конструкция в эти точки никогда не войдёт в статике. А почему не войдёт? Потому что решатель в этих точках при мягком нагружении не сходится. В них не будет решения, потому что появляются скорости и ускорения - конструкция начинает двигаться. Решатель ессно пасует. Но, чуть-чуть пораскинув своими куцыми мозгами, решатель находит такую точку (точки) с такими перемещениями, в которых конструкция устойчива.

Вывод: в своей голове надо разделить задачу баклинга в МКЭ на две подзадачи: 1) статика+проблема сходимости решателя 2) модалка.

Изменено пользователем AlexKaz
Ссылка на сообщение
Поделиться на других сайтах

Присоединяйтесь к обсуждению

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

Гость
Ответить в тему...

×   Вставлено в виде отформатированного текста.   Вставить в виде обычного текста

  Разрешено не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отобразить как ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставить изображения напрямую. Загрузите или вставьте изображения по ссылке.

  • Сейчас на странице   0 пользователей

    Нет пользователей, просматривающих эту страницу.




  • Сообщения

    • gudstartup
      как перевести шпиндель снова первоначальный режим управления скоростью когда ему пофиг число импульсов на оборот? резьбу пока резать не будут а в подпрограмму смены можно что нить добавить или в реализацию m19  чтобы чпу из режима контроля позиционного кодера вышло?
    • gudstartup
    • mactepnew
      Вот я тоже так думал, пока не столкнулся с обратным) Как еще столкнусь, пришлю скриншот. Хаос еще в том, что порядок в дереве построения весьма условный. операция, которая стоит перед групой, имеет бОльший timestamp(пардон, не знаю как это в переводе), чем эта группа, и выполняется позже. При этом все операции в группе имеют timestamp на пару сотен меньше и выполняются раньше, чем две группы, стоящие перед ней. Логика?  Если в эскизе все элементы определены, то можно поставить еще один размер (хочешь переопределить эскиз - да пожалста!), изменить его - и NX самостоятельно выберет, какой размер или ограничение отпустить и весь эскиз летит к чертям. Это ли не хаос?   Назовете хоть один CAD, которым можно пользоваться так? Любой из встреченных мной имеет много неочевидных нюансов, которые конечно не проявляются при выполнении университетского или базового курса. В Creo я на протяжении 5 лет находил всё новые фишки, которые не прописаны в документации. Огромная благодарность Ruslan за пояснения тонкостей. Даже Компас и Solidworks, которые считаются самыми дружественными, оказалось непросто обуздать, когда проходил собеседования и делал тестовые. А ведь университетский курс знал на отлично!   Возвращаясь к NX - вот сегодняшний пример. Утром деталь была оранжевая, единая. Изменил пару уклонов и толщин, стала желтая с элементами оранжевого и частично твердотельной, частично листовой. Я даже не знаю, как сформулировать вопрос в поддержку. 
    • Mikhail_Podgornyj
      Добрый день, кто может поделиться SKF Bearing V1.6 (ANSYS App)? Заранее благодарен!
    • kkk
      Короче, как сделать автоматом не нашел. Раскрасил нужные грани нужными цветами. Затем на вкладке с цветом модели можно выбирать нужный цвет, тогда в модели выбираются грани, покрашенные этим цветом. Можно еще сохранить выбор этих граней, но смысла не вижу, т.к. все легко выбирается по цвету. Затем либо в инструменте Измерить смотрим площадь и как-то ее используем, либо можно макросом (что выкладывал Kir95 ) внести площадь в свойства (на вкладку Конфигурации). Можно еще внести выбранные грани в датчики, но как использовать их в свойствах не нашел.    Макрос прикрепляю.     Площадь выбранных граней в свойства.swp
    • IgnatvKudrin
      Для нашего производства требуются постоянно инструментальные пружины для штампов. Сейчас заказы обрабатываются довольно долго - до 2-х месяцев. А у нас сроки на штампы 1-2 месяца с испытанием, доработкой и т.п. Хотелось бы найти поставщика, у которого уже есть склад с данными пружинами или кто может доставить за 1 месяц.   Пружины стандартные по ISO, но на всякий случай приложил каталог. Вот список пружин, которые нужны очень срочно:   LR 32 x 44 - 1 шт. LR 25 x 51 - 2 шт. LB 63 x 89 - 7 шт. LB 32 x 102 - 10 шт. LB 32 x 89 - 20 шт. LB 32 x 76 - 12 шт. LB 32 x 51 - 6 шт. LB 25 x 127 - 9 шт.   Готов выслушать все предложения. Возможно даже есть те, кто сможет изготовить подобные пружины с такими же параметрами. 20231101.pdf
    • bubblyk
      Знатоки форума, подскажите, как смоделировать простую (на первый взгляд) задачу термомеханики в упругой постановке.   Есть брусок правильный формы ака кирпич, в котором генерируется НЕРАВНОМЕРНОЕ удельное тепловыделение (HGEN). Задняя и передняя поверхности конвективно охлаждаются с известными (но различными) коэффициентами теплоотдачи и температурами жидкости/газа. Свойства кирпича (теплофизика + механика) принимаются зависящими от температуры. Постановка - стационарная Собственно вопрос : как грамотно смоделировать температурное поле в кирпиче с учетом температурного расширения?   Исходный код приложил.     Надеюсь на коллективный разум и опыт! EPTH2D.mac p.s.   Если есть опыт решения этой задачи в Ansys WB, тоже можно сюда приложить командные вставки.
    • OutXL
      Это было давно использовал NC-Link  и данную инструкцию. The Maho.docx
    • Mikhail_Podgornyj
      Добрый день, кто может поделиться SKF Bearing V1.6 (ANSYS App)? Заранее благодарен!
    • alexqwq10
      В данных корпусах ничего лишнего быть не должно. Максимум шпильки заземления, в количестве предусмотренном для конкретного корпуса. Рассматривал очень много подобных корпусов разных производителей и данная проблема актуальна для подавляющего большинства.
×
×
  • Создать...