1.7 Визуализация результатов

В качестве критерия разрушения необходимо взять Seqv. На экран выводятся поля напряженийMain Menu> General Postproc> Plot Results> Nodal Solu…

Рисунок 27 - Поля напряжений на диске

Для корректного определения коэффициента запаса диска необходимо учитывать пределы длительной прочности. Формируется макрос расчета коэффициентов запаса. Для этого задается массив типа Table, позволяющий получать значения между заданными, -Utility Menu> Parameters> Array Parameters> Define.

Рисунок 28 - Опции задания массива

После создания массива, нажимается кнопка EDIT и он заполняется - первым столбцом идут температуры, а вторым соответствующие им пределы прочности. Затем сохраняются результаты заполнения –Apply.

Рисунок 29 - Точки массива

Задается массив для операций макроса Menu>Parameters>Array Parameters>Define-Edit

Рисунок 30 - Опции задания массива

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

*do, I, 1, n_count

_tan = (T_obod - T_stup) / (r3 / 1000-r1 / 1000)

Oper(I, 1) = T_stup + (Nx(i) - r1 / 1000) * _tan

*enddo

Заполняется второй столбец пределами длительной прочности, соответствующимитемпературам для каждого узла.

*do, i, 1, n_count

Oper(i,2) = Dlit (Oper(i,1))

*enddo

Далее следует переход в постпроцессор и заполнение третьего столбца действующими напряжениями.

FINISH

/POST1 ! Переход в постпроцессор

*do, i, 1, n_count

*get, res_, node, i, s, eqv !Задание параметру res_ значения в I узле

Oper(i,3) = res_ ! Заполнение ячейки I 3 го столбца результатом

res_= !Обнуление параметра

*enddo

Производится расчет коэффициентов запаса, равных отношению предела длительной прочности и действующего напряжения в узле.

*do, i, 1, n_count

Oper(i,4) = oper(i,2) / oper(i,3) ! Расчет коэффициента запаса

*enddo

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

/prep7

*do,i,1,n_count

BF,i,TEMP,Oper(i,4)

*enddo

После выполнения макроса отображаем результаты расчета.

Рисунок 31 - Поля коэффициентов запаса