8node curved shell

[Гость suvorov.as]

Мне нужно сравнить свой оболочечный элемент с чем-нибудь кроме ансиса.

[a_schelyaev 367]

Правильный ответ - неужели есть трудоголики, которые в выходные будут фигней страдать?

;)

[Bonusfrag 104]

неужели есть трудоголики, которые в выходные будут фигней страдать?

Есть)

У каждого свой способ отдыха ))

Но с прочностью ничем не могу помочь.

[a_schelyaev 367]

Копааать огород черепам быстро!

:)

[Bonusfrag 104]

черепам

Я, видимо, не в теме...

[Цветочек 15]

И самое главное, ну не по теме форума вопрос, что означает фраза например 8node curved shell? я знаю 4-5 формулировок таких и 4,8 и 9 узловых элементов, тольк зачем это надо - сравнивать матрицы жесткости, не представляя собственно теории-формулировки элемента. Как я понимаю это продолжение глубокомысленного вопроса на тему - чем отличается в ансис элемент shell A от элемента shell B. Особенно, когда всё кончилось тем что решили, а тем, что нету shear and volumetric locking.

То бишь - каков вопрос, таков и ответ

Прикрепляй к вопросу статью с "твоей" формулировкой, тогда хоть есть над чем копья ломать, а если нету - на то и разговора нету

[Цветочек 15]

Не забудь только, что элементы сплошного тела - что плоские, что трёхмерные - простые элементы, без прикрас(я имею в виду обычные элементы, не исправленные, то есть). Для оболочек и пластин процесс формировования матриц значительно более сложный. Сам понимаешь, там производные 4 порядка, более сложные теории - например, Миндлин, Рейснера. При том, появляется лишняя сложность, уже готовые оболочковые КЭ ведут себя плохо - не проходят patch тесты, а значит использоваться не могут.

Ты предлагаешь разбираться в твоём элементе, и пишешь только числа. Зачем они нам? чтоб сойти с ума?

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

В ансис один и тот же элемент в зависимости от опций может использовать даже разные теории оболочек, а как ты думал.

Какие ссылки? открываешь журныл прикладная механика с 1960 по наши дни, а там статьи, статьи, а потом ещё иностранные журналы до кучи

[Цветочек 15]

Как это хорошо, что вы у нас есть - ИСПА. Что нам ни захочется, в вас всё есть.

[Борман 2 388]

Честно говоря я большего ожидал.

От мирового лидера.

"Бог с ним с ANSYS. Не пытайтесь показать, что какая-либо программа хуже вашей программы.

Покажите лучше, чем ваша программа лучше других программ."

знаете кто сказал ?

[Цветочек 15]

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

Я рекоммендую вам использовать книжку The Finite Element Method Using MATLAB - Kwon and Bang. В этой по истине замечательной книжке есть не только теория пластин, но и её применение к МКЭ. На странице 361 вы найдете то, что в общем вам и было нужно с самого начала. Надеюсь вы будете рады

[Цветочек 15]

Я считаю, что для начала нужно освоить простое. Какой смысл говорить сложности, если не понятно, что такое просто, а что сложно? Зачем столько букв в названии, если не понятно, собственно чем теория миндлина отличается от кирхгофа, и что же всё таки называют технической теорией. Уровень форума низок до крайности, а если так, даже если заинтересованным лицам будет ясно то, что должно быть ясно до того, как вопрос был опубликован, я буду только рада. Может быть человек выйдет на то, чего ж он в самом деле не понимает и что ему нужно. А там уже и до tochnog или marc с фортраном недалеко.

А по поводу толщины оболочки как степени свободы? вам это нужно для чего? сами хотите запрограммировать или эффективность проверить?

Shear locking лишены очень многие оболочечные КЭ.

[Цветочек 15]

Думаю вы несколько бахвалитесь, ибо если так, то вам самое время к книжке the finite element method: linear static and dynamic finite element Hughes, ибо конструирование конечных элементов для балок и оболочек это самое сложное и неблагодарное дело,из всего МКЭ. А значит очень скоро мы будем не нужны вам

[Цветочек 15]

В предвкушении скорой разлуки позвольте мне набросать вам списочек книжек по оболочкам и МКЭ,

Elasto-Plastic and Damage Analysis of Plates and Shells George Z. Voyiadjis, Pawel Woelke

Laminated Composite Plates and Shells: 3D Modelling Jianqiao Ye

Theories of Plates and Shells: Critical Review and New Applications Reinhold Kienzler, Holm Altenbach and Ingrid Ott

Finite Rotation Shells: Basic Equations and Finite Elements for Reissner Kinematics Krzysztof Wizniewski

The Mechanics of Solids and Structures - Hierarchical Modeling and the Finite Element Solution Miguel Luiz Bucalem and Klaus-Jurgen Bathe

Уверена, вам нужно знать что сейчас делают в мире по оболочкам и КЭ

[Гость suvorov.as]

Вот она, вот она, оболочка моей мечты

Результат теста (ux, uy, uz, usum)

1.0e-003 *

0 0 0 0

0 0 0 0

0.0284 0.0127 0.3067 0.3082

0.0022 -0.0091 0.2281 0.2283

0.0031 -0.0017 0.0093 0.0099

0.0064 -0.0002 0.1319 0.1320

0.0091 0.0008 0.2490 0.2492

0.0032 -0.1160 0.1222 0.1685

Ansys shell93

NODE UX UY UZ USUM

1 0.0000 0.0000 0.0000 0.0000

2 0.0000 0.0000 0.0000 0.0000

3 0.28488E-04 0.12695E-04 0.30680E-03 0.30838E-03

4 0.21891E-05-0.92212E-05 0.22791E-03 0.22810E-03

5 0.31021E-05-0.16980E-05 0.93071E-05 0.99563E-05

6 0.64340E-05-0.22637E-06 0.13196E-03 0.13211E-03

7 0.92025E-05 0.77262E-06 0.24895E-03 0.24912E-03

8 0.31714E-05-0.11590E-03 0.12217E-03 0.16843E-03

[Гость suvorov.as]

И в чем была ощибка?

Жду новый код.

Это не ошибка. По статье орты интерполируются ф.12.5-13 и это приводит к неточности при большой кривизне.

Нужно орт х и орт y искать в точках интегрирования непосредственно через производные функции формы и координаты. Затем находить орт нормали, подправлять орт y и нормировать все орты.

Я сейчас обнаружил у себя еще один косяк (не влияющий на приведенный результат), когда исправлю выложу код

[Гость suvorov.as]

Вот скрипт:

....

[Гость suvorov.as]

Вы собираетесь все это писать на языке С?

Собираюсь, но не лично и потом.

[Гость suvorov.as]

Но вы будете делать более серьезные тесты?

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

И еще надо провести сравнение на реальных моделях.

Кто все это будет делать?

Если один кривой элемент совпал, то и все остальные совпадут. Объединить все КЭ в глобальную матрицу не так уж и сложно. Полученные глобальные матрицы опять же можно будет сравнить с ансисовскими до решения СЛАУ

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

[Гость suvorov.as]

Как они совпадут когда вы интеграл берете по 2 точкам.

А нужно брать по 3 точкам. Мембранная составляющая у вас неправильно получается.

Так в чем проблема то, поменять потом на три точки и посмотреть, что из этого выйдет? Нюансы то алгоритма построения матрицы элемента выяснены и от количества точек интегрирования они не зависят.

Пока на уровне химии вы подогнали матрицу жесткости одного КЭ. И уже делаете выводы.

А какаой вывод я сделал неправильно? Что с ансисом совпадает? Так я много геометрических конфигураций элемента проверил и в плоскости кривых и изогнутых. Все с ансисом совпадают.

[Гость suvorov.as]

Но делать выводы по 1 матрице жесткости пока рано.

Это же не оболочка Кирхгоффа. На тонкостенных конструкциях химия вылезет.

Проявит так сказать себя.

Так какие выводы то? Что для любой постановки эксперимент совпадет с расчетом? Не далаю я таких выводов и никогда их не сделаю.

В ИСПА уже 20 лет есть тонкая, криволинейная 8-ми узловая оболочка.

Это же тоже не оболочка Киркгофа

DKT в ИСПА только 3 и 4 узловые оболочки.

Все давно протестировано и проверено.

Ну и ладушки.