Гармонический Анализ Одной Детали Из Сборки

[Chilavert 0]

Вобщем в продолжение задачи о нагреве.

Есть корпус (пластиковый), в нем плата из стеклотекстолита, она крепится к корпусу на 3х пластиковых заклепках, плюс в плату еще впаяны контакты из меди, которые влиты в корпус. Т.е. все элементы взаимосвязаны. Все это дело нагревается с 20 до 120 градусов и деформируется из-за разности коэффициентов температурных расширений. Эта задача успешно решена совместно - все работает, причем очень интересная картина деформаций получается (с меня пример).

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

На плату действуют только силы, возникающие при нагреве до 120 градусов. Можно идти двумя, или даже тремя путями - вот я и думаю.

1. В Antype,3 есть вариант приложения гармонического температурного поля, т.е. по идее можно сделать так:

/SOLU

ANTYPE,3 ! Гармонический

HROPT,FULL ! Полный метод

HROUT,ON

LUMPM,0

EQSLV,FRONT,1e-008, ! Прямой решатель !

TUNIF,120, ! Прикладываем температуру в 120 градусов

TREF,20, ! Референсная температура в 20 градусов

HARFRQ,1000,20000, ! Указываем пределы частот в которых исследуем

NSUBST,10, ! Количество подшагов (частот)

KBC,0 ! Линейное изменение частот в пределах интервала

SOLVE

Все бы хорошо, НО, сам Ансис предупреждает - вы действительно хотите приложить гармоническое изменение температуры? И результаты даже на простой задаче получаются неправильными (картина темпераутрных деформаций). И Самое главное - полная модель весит 125 тыс. элементов, а плата - 7 тыс., и даже мой сервак 4хOpteron285 4 Gb RAM, 2x74.5 Raptor RAID затыкается, временные файлы растут, и уже на первой частоте виснем.

2. Вариант два, тоже по классике (по мануалу) - есть вариант с предварительным напряжением.

Учитывая, что: "Files Jobname.EMAT and Jobname.ESAV from the static analysis must be available", записываем файлы матриц и элементов в файлы.

Т.е. в 1м анализе (с температурой) вставляем в код строки:

/PREP7

EWRITE,'file',' ',' ',0 ! Пишем элементы

/SOLU

EMATWRITE,YES ! Даем команду записывать матрицы при решении

PSTRES,ON ! Включаем предварительное напряжение для дальнейшего анализа.

....

SOLVE

Во 2.м анализе код такой:

/SOLU

ANTYPE,3 ! Гармонический

HROPT,FULL ! Полный метод

HROUT,ON

LUMPM,0

EQSLV,FRONT,1e-008, ! Прямой решатель !

HARFRQ,1000,20000, ! Указываем пределы частот в которых исследуем

NSUBST,10, ! Количество подшагов (частот)

KBC,0 ! Линейное изменение частот в пределах интервала

PSTRES,ON ! Включаем предварительное напряжение

SOLVE

Все бы хорошо, но все 125 тыс. элементов решатся неизвестно когда - только матрицы для первой частоты создаются более суток, файл *.tri достигает 50 Гб. А если перед вторым анализом выделить только нужные нам элементы - сразу вываливается с ошибкой - что и понятно, матрицы то и сами элемнты для всех 125 тыс. а есть только 7 тыс.

3. Мой вариант

Вобщем подумал, а если после первого анализа взять перемещения UX, UY, UZ, всех узлов (или только по нижней или верхней части) платы, сохранить в параметры, а потом при гармоническом анализе приложить. получится низкая размерность задачи, чистая гармоничность (что силу прикладывать, что перемещение эквивалентное), и должно сойтись, т.к. нас только плата и интересует. Весь термальный анализ на APDL делался - так что если подумать как это реализвать в коде, должно получиться.

Вопрос в том, КАК не руками сохранить перемещения всех узлов нужных нам элементов в параметры (массив, файл) записать, и потом приложить назад.

Если у кого еще есть какие варианты - с большой благодарностью выслушаю.

Заранее благодарен!