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

Расчет деформаций в трехмерном элементе


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

При работе с деформациями в трехмерном элементе (тетраэдр), наткнулся на одну проблему.

При вычислении деформаций с помощью ansys mechanical apdl 14.5 и вычислении вручную с помощью формул через матрицу B (метод расчета взят из книги О.Зенкевича) результаты не совпадают. Пришел к выводу, что проблема, по моему мнению, в расчете коэффициентов a, b, c, d. При этом, в качестве модели взял простой куб. Также пробовал и на других моделях.

Безымянный.jpg

Написано, что следующие коэффициенты получаются циклической перестановкой индексов  p,I,j,m.  При этом, если просто переставлять индексы, конечный результат не совпадает. Вроде говорят, что при перестановке где-то должен знак меняться.

Может кто-нибудь тоже сталкивался с подобным? Хочу понять сам принцип расчета деформаций в узлах (не элементах), зная перемещения узлов.

P.S. мне это надо для дипломной работы.

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


https://ru.wikipedia.org/wiki/Тензор_деформации   все же просто ... :)

 

"принцип расчета деформаций в узлах (не элементах) " - это невозможно. так как не будет производных. Они будут только в элементах. В узле можно например осреднять от входящих элементов с какими - нибудь весовыми коэффициентами более менее логично отражающими влияние элемента.  Например, чем более углы отклоняются от прямого, тем меньше им доверия :)   Точка в механике сплошной среды понимается как малый объем. Это надо чтобы существовали производные :)  

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

В Ansys не парьтесь, ничего не надо считать, все и так есть и рисуется и выбрать уже осредненные в узлах можно :)

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

Коэффициенты вычисляются в лоб.

Снимок экрана - 10.02.2017 - 12:19:15.png

Снимок экрана - 10.02.2017 - 12:19:43.png

46 минут назад, Fedor сказал:

уже осредненные в узлах можно :)

Наверное на интерполированные @BedaPiton и смотрел...

Если решение сошлось с ANSYS по перемещениям, то дело не в коэффициентах.

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

Для примитивной пирамидки наверное, но кому такая архаика нужна и зачем ? 

И все равно в узле элемента, а в общем случае в узел входит несколько элементов :)

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

Кроме того, в примитивной пирамидке линейные функции перемещений, следовательно при дифференцировании будут константы. и деформации будут постоянными по всему элементу... У Сегерлинда вроде все это расписано :)

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

Позвольте поинтересоваться, а что Вы считаете и зачем вам такой подход? Свой МКЭ пишете?))
Попробуйте для сравнения в параллелепипида посчитать, там вроде и математика и решение проще будут(без поворота площадок и осей) и проще будет сравнить Ансис и аналитику.

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

Да все там однообразно для любых элементов. Базисные функции можно посмотреть тут http://know.sernam.ru/book_mge.php  например .

И выгоднее не вычислять определитель, а разложить матрицу. Быстрее интегрирование, насколько помню :)

 

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

ХЗб я в такие дебри не лез но мне было проще на учебных задачах в универе кубики считать нежели пирамидки) Для понимания так точно.)))

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

Не разобравшись в изопараметрии в мкэ делать нечего. В смысле понимания и программирования. Задачки то можно решать  :)

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

результаты не совпадают

они могут не совпадать по многим причинам..во первых как уже сказал @Fedor деформация и напряжения в элементе постоянны, и расчет в узлах производится либо усреднением по соседним элементам, либо экстрополяцией от центра элемента, и дальнейшем усреднением по узлам. Эстрополяция дело хитрое, по крайней мере в NASTRAN их аж 4 штуки (по умолчанию стоит билинейная), думаю в ANSYS их тоже не один вариант. Помимого всего, при анализе экстраполируемых значений деформации в узлах, стоит внимательно смотреть за качеством элемента. При разбиении на кое как, можно и бред вовсе получить. Поэтому сойтись до знаков после запятой врятли представляется возможным..

Кстати, насколько помню, линейные тетры вообще не экстраполируют в узлы, а осреднение идет только по соседям..

 

 

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

Не разобравшись в изопараметрии в мкэ делать нечего. В смысле понимания и программирования. Задачки то можно решать  :)

Все расчеты МКЭ это всего лишь задачки=)  МКЭ - инструмент-)))

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

Так и человек это животное которое умеет создавать неорганические инструменты для решения задачек.

 

"Интеллект - это способность создавать и применять неорганические инструменты... " - http://society.polbu.ru/reale_westphiloiv/ch142_all.html  :)

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

Так и человек это животное которое умеет создавать неорганические инструменты для решения задачек.

 

"Интеллект - это способность создавать и применять неорганические инструменты... " - http://society.polbu.ru/reale_westphiloiv/ch142_all.html  :)

Никто и не спорит))) Животное с интелектом=)

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

@BedaPiton А вообще, может быть полезена A First Course in the Finite Element Method. Daryl L. Logan стр. 495.

http://kntu.ac.ir/DorsaPax/userfiles/file/Mechanical/OstadFile/dr_nakhodchi/DarylL.LoganAFirstCourse.pdf

 

2017-02-13_21-08-42.png

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

@kolo666 спасибо большое, очень помогло. Теперь все нормально сходится.

Вы там что на питоне пишете собственную решалку*?)
 

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

Вы там что на питоне пишете собственную решалку*?)
 

Тоже самое спросил=)

 

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

Ну почти, только не на питоне, а на яве.

 

На диплом? Зачем это вам))))

 

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

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • vad0000
      Если Вы курсор поставите на OUTPUT 0, и колесом начнёте увеличивать напряжение на аналоговом выходе, то ось Х поедет, или нет? Сфоткайте машинные параметры для оси Х
    • lexa345
      Если фрезерный, то спросите у @andrev  возможно поделиться   
    • Bot
      Семинары «BIM-Практикум 2024» пройдут весной и летом этого года в четырнадцати городах в России и Беларуси. На мероприятии представители строительной отрасли расскажут об опыте применения технологии информационного моделирования в регионе. Российские разработчики представят новые возможности цифровых решений для проектирования и управления строительством. У участников будет возможность обменяться мнением и опытом, проконсультироваться с экспертами отрасли и установить полезные контакты. Зачем Сегодня на развитие строительной отрасли влияют как внутренние факторы, включая небольшой опыт ведения проектов с использованием технологии информационного моделирования, недостаток квалифицированных специалистов, высокую стоимость ошибок, так и внешние обстоятельства: новые государственные требования, риски, связанные с информационной безопасностью, и уход западных вендоров. Участники семинара узнают, как справиться с этими вызовами, в том числе на примере компаний, которые [...] View the full article
    • MOPKOBKA555
      А есть ли возможность в cnckad разделить программу на части? Например сразу обрабатываем 9-ю инструментами, а с новой программы еще недостающими?
    • coolparty
      Привет. Задача - расположить массив отверстий по кривой. Сделал модель из нескольких кривых, не понимаю, как их объединить в одну, пробовал сделать через "объединенная кривая", но при выборе ее в качестве пути, все равно выбирается по сегментно. Всю голову уже сломал.  
    • Bot
      3 и 4 апреля в Екатеринбурге и Челябинске состоялась конференция «Российский программный комплекс T-FLEX PLM: современные методы и подходы к проектированию изделий», организатором которой выступила компания «Топ Системы». Во время мероприятия представители компании «Топ Системы» — генеральный директор Сергей Кураксин и заместитель генерального директора по развитию PLM-технологий Игорь Кочан — представили текущие итоги реализации проектов у заказчиков из различных отраслей машиностроения, поделились перспективами развития программного комплекса T-FLEX PLM, продемонстрировали функциональные возможности перспективных разработок САПР T-FLEX CAD и платформы T-FLEX DOCs по управлению данными и жизненным циклом изделия. В Челябинске вызвал живой интерес доклад Романа Ткачука, начальника отдела САПР, АЗ УРАЛ, во время которого представитель заказчика рассказал о проектном опыте перехода с иностранных систем на комплекс T-FLEX PLM. Цель проекта заключалась в переводе [...] View the full article
    • roiman
      Не одно и тоже, нет?     Не одно и тоже, нет?   Я про это и писал, что вместо того, чтобы решать реальную производственную задачу, мы занимаемся словесной казуистикой. Задача не поставлена, но решение её требуется. А может и не требуется. Начальник сказал попробовать спросить. А может что-то другое. Ценится только своё время, а не чужое...   Ещё раз. Если сложить пазл из данных статьи и вашу картинку, то профиль не может никак участвовать в уравнениях. У вас идёт контакт зуба с роликом цепи. И их положение относительно друг друга не меняется в зависимости от угла поворота "кулачка". И всё натяжение цепи передаётся на пантограф через этот контакт зуба с роликом. Задача профиля - обеспечить натяжение цепи.   Но это, опять же догадки. Потому что непонятно, что это за 2-а других отверстия. Может там не просто поворот этого "кулачка", а его плоское движение.     P.S. Вы зря удалили статью. Я никакого отношение не имею ни к УКВЗ, ни ПК ТС. Опасаться нечего.
    • Александр1979
      На первый взгляд параметры конфигурации осей B, C корректные. У вас станочный ноль осей X, Y находится по центру поворотной оси C?
    • lem_on
    • gudstartup
      какая величина отклонения при развороте плоскости оси b ? кинематику станка проверяли прямолинейность перпендикулярность осей xyz точность поворота осей b и с? 
×
×
  • Создать...