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

Уравнения МКЭ для термического расширения материалов


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

Доброго времени суток.

Я решаю различные задачи в Comsol. Особо глубоко рассматривать вопросы, связанные с математическим описанием различных явлений необходимости не было, т.к. Comsol позволяет моделировать особо об этом не задумываясь.

Однако, последнее время есть необходимость разобраться в некоторых вопросах, связанных с МКЭ.

В хелпе Theory on Solid Mechanics написано что для решения задач напряженно деформированного состояния используется уравнение слабой вариационной формы для "Виртуальной работы". Физический смысл в том, что работа внутренних напряжений (интеграл по объёму от напряжений * деформации) равна работе внешних сил (интеграл от внешних сил * перемещения). (Если я всё правильно понимаю).

А когда мы считаем напряжения в материале, который расширяется под действием температуры, и в нем возникают какие-то напряжения, мы в этом случае используем тот же принцип, или нечто другое?

В частности интересует такой пример: двухслойная пластина нагревается до какой-то температуры, слои имеют разные коэффициенты термического расширения. Возникает коробление пластины. Какой в этом случае используется вариационный принцип?

Где про это почитать, чтобы написано было нормально? а то в Comsol хелп какой-то очень краткий.

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


В хелпе Theory on Solid Mechanics написано что для решения задач напряженно деформированного состояния используется уравнение слабой вариационной формы для "Виртуальной работы". Физический смысл в том, что работа внутренних напряжений (интеграл по объёму от напряжений * деформации) равна работе внешних сил (интеграл от внешних сил * перемещения). (Если я всё правильно понимаю).

Рассматривается функционал полной потенциальной энергии: П(u)=Э(u)-А(u), где Э-энергия деформации, А-работа внешних сил, u - перемещение.

Если u* - реальное смещение, то dП=П(u*+du)-П(u*)=0 для любого du. Внутри программы это уравнение сводится к системе линейных уравнений, и решается, в результате чего находятся u*.

А когда мы считаем напряжения в материале, который расширяется под действием температуры, и в нем возникают какие-то напряжения, мы в этом случае используем тот же принцип, или нечто другое?

Тоже самое. Дело в том, что программа не знает, что нагрузка из себя представляет в реале. Для неё все нагрузки сводятся к сосредоточенным, распределенным по линии (поверхности), распределенным по объему (объемным).

Так вот, температура - это та же объемная нагрузка, но распределенная хитрее чем вес, силы инерции и т.п.

В частности интересует такой пример: двухслойная пластина нагревается до какой-то температуры, слои имеют разные коэффициенты термического расширения. Возникает коробление пластины. Какой в этом случае используется вариационный принцип?

Все одно... Температурной задачи здесь нет. Задавая температуру - вы просто задаете объемную нагрузку.

Ведь задавая собственный вес, вы же не решаете гравитационные уравнения Эйнштейна, ну или кто там их придумал... Я не решаю, во всяком случае.:)

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

Борман, спасибо большое.

Тоже думал, что для всех задач в Solid Mechanics используется один принцип, однако меня смутил в хелпе комсола раздел "Энтропия и термоупругость", в котором вводилась функция свободной энергии термоупругого материала (The free energy for the linear thermoelastic material) как сумма энергии упругой деформации и, видимо, количества сообщенного телу теплоты. Вот я думал, что эта свободная энергия тоже где-то завязана в вариационных принципах.

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

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • Orchestra2603
      Да, в рот мне ноги... Ну, считает же Ансис, считает! Блин.. Как до вас донести простой мэсеж. Мне често, пофигу, делает он там себе факторизацию или еще что-то. Но то, что он нормально считает Ланцошем для случая с вырожденной матрицей жесткости - это факт! А раз он считает, то либо (А) вы хрень несете про факторизацию, либо (Б) он такую факторизацию делат сам без проблем. Третьего тут быть не может, и выходит, что в любом случае ваши тезисы оказываются несостоятельными.   Какая вам нужна факторизация? давайте! заказывайте! На какие вам надо матрицы разложить? с какими свойствами? давайте ваш заказ!   Вы на мой пердак не зарьтесь! Смотрите, чтоб ваш функционировал исправно! А то того и гляди, фекализация случится.  
    • ДОБРЯК
      Вам похоже ничего не нужно, а для алгоритму Ланцоша нужно сделать численную факторизацию. Вам только нужно побольше букв и слов написать. :=) Всё подгорел пердак? :=)
    • Orchestra2603
      При чем здесь это? Речь не идет про эффективные или неэффективнеы алгоритмы. Вы утвержаете, что это просто невозможно.   При чем здесь первые или не первые, все или не все... Какое это имеет значение? У Ансис есть алгоритм Ланцоша, он находит столько, сколько запросите. В Маткаде - он да, находит все. Наверное, можно и не все, но я хз, как это делать. Я не понимаю, какое это имеет отношение к разговору.   Неа Ткните пальцем и объясните, в чем ошибка в конкретном рассуждении   Как я должен решать методом Гаусса задачу на собственные значения? Мне же не решения СЛАУ нужны! Их то, ясен пень, бесконечное множество для системы с вырожденной матрицей. Мне собственные значения нужны! Мне нужно базис найти, чтобы матрица (или матрицы в обобщенном случае) оператора была диагональной.   Сначала скажаите мне, зачем ее дать! В чем идея заключается? Вы все слюной брызжете, но ничего по делу еще не сказали. Какой-то нормальный контраргумент от вас можно толковый услышать?   Как еще изволите вас обслужить? Вы - полегче, уважаемый!
    • ДОБРЯК
      Вам нужно десять раз написать одно и тоже, чтобы вы поняли. Почитайте документацию к программе там же всё написано. Мне не сложно копировать...:=) Те алгоритмы которые вы используете находят все собственные числа и вектора. Для разреженных матриц это приводит к их заполняемости. Смекаете о чем я говорю. Или нужно еще разжевать? Во всех.   Решайте методом Гаусса. Матрица 2х2 в чем проблема сделать численную факторизацию любым методом? :=) Больше дела, меньше слов и букв...
    • mannul
      Скорость резания от 8 до 12 м/мин, подача 0,04 мм/зуб, глубина резания - 1,5 мм.
    • Orchestra2603
      " Стою на асфальте я, в лыжи обутый: То ли лыжи не едут, то ли я ..." (с)   1. Вы утверждаете, что для расчетов на СЗ и СВ в случае вырожденной матрицы жесткости нужно провести ее факторизацию. Без нее - никак! Непонятно, какую и зачем, но надо. Ок, допустим.  2. @Jesse,  я, @Fedor, да и вы сами показываете, что таки СЗ определяются для вырожденной матрицы. Значит, либо ваше утверждение про факторизацию неверно, либо она таки делается тихонько себе за кадром. 3. Еще одно наблюдение... Матлаб говорит, что он делает QZ разложение для матриц, у которых нет положительной определенности, т.е. для которых не работает разложение Холецкого. При этом опять же с собственные числа с одной вырожденной матрицей считаются нормально. Т.е. даже если вдруг и необходима такая факторизация (хотя я сомневаюсь, что вот прям обязательно она нужна), то делается она без особых трудностей.   В каком месте я неверно рассуждаю?   Так и не обязательно их записывать для расчета собственных колебаний. Ну... Можно конечно, но тогде некоторые частоты и формы, которые такой симметрией могут и не обладать, вы потеряете в расчете, поскольку введенные условия симметрии их исключат просто из модели.    
    • gudstartup
    • gudstartup
      у меня рободриллы и чпу и тем не менее picture нет наверно она к роботам относится
    • boomeeeer
      Если компания отечественная, то тоже нет смысла. Роботы, ЧПУ, лазеры и рободрилы
    • gudstartup
      это вряд ли если ось Y угловая а какая у вас редукция у вращающейся оправки? по 2мм нержавейки снимать это нелегкая задача для вашего маленького станочка
×
×
  • Создать...