Решение задач динамики и модальный анализ

[Satellit 0]

Кстати, Satellit решите кубик на 1 000 000 степеней свободы. Ваш компьютер легко потянет. Интересно же время рещения сравнить.

150 частот нашел за 2:55:42

Число степеней свободы 1058610

[Fedor 1 597]

"Вы же профессионал, а на такие мелочи внимание обращаете... Я думаю торг здесь неуместен" - поэтому и обращаю :)

В самом деле, ризоме постмодерна вроде меня это не к лицу, шизофренией попахивает, а у меня и без нее забот хватает :)

[Vova 419]

Nastran TETRA10

2400 MHz

Time 64min 5s

NUMBER OF GRID POINTS = 343711

NUMBER OF CTETRA ELEMENTS = 215862

1 18.498

2 18.658

3 27.557

4 73.042

5 74.641

6 74.981

7 75.446

8 76.005

9 76.531

10 76.791

.........................

Внимание тут исправил!

140 1073.5

141 1080.6

142 1082.1

143 1102.4

144 1164.7

145 1167.4

146 1168.6

147 1168.8

148 1181.6

149 1183.0

150 1188.7

Внимание тут исправил!

Пока тетра считалась, сделал гексу:

Nastran HEXA

2400 MHz

Time 3min 30s

NUMBER OF GRID POINTS = 50048

NUMBER OF CHEXA ELEMENTS = 35951

1 18.339

2 18.698

3 27.466

4 71.993

5 73.588

6 73.946

7 74.380

8 74.937

9 75.909

10 76.427

..............................

140 1075.6

141 1082.2

142 1083.7

143 1104.3

144 1160.8

145 1163.2

146 1164.3

147 1164.6

148 1182.4

149 1189.4

150 1193.7

********************************************************

ABAQUS: TETRA10

Time 36min 5s

1 18.052

2 18.209

3 26.895

4 71.282

5 72.842

6 73.174

7 73.628

8 74.173

9 74.688

10 74.941

............................

140 1074.2

141 1081.3

142 1082.8

143 1103.1

144 1165.7

145 1168.3

146 1169.5

147 1169.8

148 1182.4

149 1183.6

150 1189.5

ABAQUS: HEXA

Time 7min 56s

1 18.357

2 18.750

3 27.517

4 72.092

5 73.712

6 74.071

7 74.508

8 75.069

9 76.011

10 76.613

..................

140 1076.7

141 1083.4

142 1084.8

143 1105.5

144 1162.5

145 1164.9

146 1166.0

147 1166.3

148 1184.0

149 1191.2

150 1195.1

Как будем искать правильный расчет?

Может проголосуем?

[Satellit 0]

Очень хороший результат. Сколько ядер при этом было задействовано?

Все четыре

[Vova 419]

NASTRAN

100 Mode

786,6Hz Hexa,

772,7Hz TETRA10

[Борман 2 375]

Вот, ИСПА, специально для вас порешал задачу (куб 10х10х10, Материал - 2e11,0.3,7850, фиксировал по плоскости, искал первые 10 частот и форм).

На графике зависимость времени решения (y, сек), от кубического корня от степеней свободы (х) .

Кто задумает повторить мой подвиг на серьезной тачке, вот APDL. Менять число в строке N=20, ответ - после решения N задач в массиве TIMES (3 столбик - кол-во узлов, 4 ст. - время решения) 

CUBE.txt

[Борман 2 375]

1 000 000

Я слабо копенгаген, но моя задача умещалась в ОП, и с определенного момента (ххххх DOF) она распухает, и привлекает другие ресурсы компа. Так что не так просто...

[Борман 2 375]

Вы полагаете, что вряд ли дождетесь результата?

Ну да... Шанцев практически нет.

[Fedor 1 597]

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

[Vova 419]

Вова, чтобы понять скоростные качества программы не поленись реши кубик. 20 первых форм будет достаточно.

Полка тетра10 имеет больше миллиона

DEGREES OF FREEDOM 1031061

[Fedor 1 597]

"ИСПА решает данную задачу за 15 минут" - это же нелинейная алгебраическая задача, почему-то упускается из рассмотрения вопрос точности...

[Борман 2 375]

Борман, я уверен, что у вас много свободного времени длинными осенними вечерами.

У меня к вам просьба. Все тот же кубик. Разбейте каждую сторону на 8 элементов. Получится 1944 степени свободы. И решите полную проблему собственных значений. Найдите 1944 собственных значений. Временные параметры не важны. Важно принципиальное решение данной задачи. Но скрывать не буду ИСПА решает данную задачу за 15 минут.

В Ансисе решается за 6мин. 7сек. На картинке - последняя форма, там же частота.

<noindex></noindex>

В аттаче - список частот.

SET.TXT

[Vova 419]

Вова, так совпадение по формам тирочка в тирочку. Остается только радоваться.

Ну тогда поздравляю, что совпадает.

Кошек дома нет, поэтому тренеруюсь на кубиках. Не частоты правда, но всё равно динамика.

Для анимации надо кликнуть на картинку.

Кубик резиновый:

Кубик алюминевый:

Кубик алюминевый с разрушением:

Повторить кто сможет?

[Fedor 1 597]

Хорошие картинки, но последняя неправдоподобна. Должна была образоваться трещина

под цилиндром раньше чем произошел бы срез ...

[Vova 419]

Хорошие картинки, но последняя неправдоподобна. Должна была образоваться трещина

под цилиндром раньше чем произошол бы срез ...

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

[Fedor 1 597]

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

Чем ерундой заниматься, лучше бы исследовали на сколько процентов снижается прочность таких соединений

<noindex>http://www.psk-holding.ru/?action=catalog&id=339</noindex> при растяжении и на сколько при сжатии при различных сталях . И для трудового народа польза и денег сэкономили бы строителям цивилизации :)

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

Можете запатентовать, денег срубить, не забудьте только 10% на деревню дедушке отстегнуть

[Борман 2 375]

Борман, вы в очередной раз подтолкнули меня к исследованиям ИСПА. Так вот, докладываю. Собственно сам расчет идет 8 мин, а все остальное время уходит на сообщения между программой и оболочкой. Вот к чему приводит лишняя болтливость.

"Вещь в себе", однако...

[Борман 2 375]

ISPA

я вот что-то не пойму, когда "динамика" начнется...

[Vova 419]

Борман, я сделал данную модель из стержней и оболочек. Минут 30 ручного труда. Никакая 3D модель не понадобилась. Время расчета – секунды. Так вот, для динамики очень даже приемлемые результаты. И не нужно часами молотить на сверхмощных компьютерах.

Совпадение в пределах 5 процентов.

Так что, вы правильно говорите про ВЕЩЬ.

Для сравнения частот лучше пользоваться абсолютными числами. Осбенно в частотах до 100Гц , а они как правило и интерeсуют больше всего.

Пример 40Гц и 42Гц это большая разница, по крайнее мере в расчетах динамики кузова автомобиля.

Вопрос точности расчетов остается открытым.

ИСПА у вас результатов эксперементов на эту кострукцию случайно нет?

[Борман 2 375]

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