Получение значения перемещения u на всех шагах решения
Сначала необходимо узнать номер какого-нибудь (любого) узла на верхней грани шара:
Показать номера всех узлов: Utility Menu PlotCtrls Numbering Node numbers (On) OK
увеличить верхнюю грань, чтобы увидеть номер какого-нибудь узла на ней (обозначим этот номер N — его надо запомнить): Utility Menu PlotCtrls Pan Zoom Rotate
Убрать изображение номеров узлов (Utility Menu PlotCtrls Numbering Node numbers (Off) OK)
Изобразить всю модель: Utility Menu PlotCtrls Pan Zoom Rotate Fit.
Найти перемещение u в каждой из промежуточных задач легче всего при помощи следующего скрипта.
*cfopen,'displ.txt' ! Открываем файл displ.txt
SET,FIRST ! Переходим к первому шагу решения
*do,i,1,10,1 ! Цикл по i от 1 до 10 с шагом 1
qq=UY(N) ! Узнаем перемещение N-го узла
*vwrite,qq ! Записываем значение в файл
(E10.3) ! в этом формате
SET,NEXT ! Переходим к следующему шагу решения
*enddo ! Конец цикла
*CFCLOS ! Закрываем файл displ.txt
Его текст надо набрать в редакторе Notepad (можно без комментариев) и сохранить в каталоге, созданном для задачи Герца. N, указанное в тексте скрипта, необходимо заменить на найденный ранее номер узла.
В результате запуска скрипта (Utility Menu File Read input from, указать имя файла со скриптом) образуется файл displ.txt, в который записываются нужные нам значения u на промежуточных шагах.
Замечание. Файл с этим скриптом не надо называть именем displ.txt (displ.txt — это файл, который создается этим скриптом). Используйте для файла скрипта, например, имя CalcDispl.txt.
Вычисление прижимающей силы p на всех шагах решения
На каждом шаге
надо найти прижимающую силу P
— это неизвестная реакция связи,
соответствующей заданному перемещению u.
P
вычисляется как интеграл от напряжения
y,
действующего на верхней грани полушара.
Его можно вычислить средствами ANSYS
как силу, действующую на выделенные
узлы со стороны выделенных элементов.
Таким образом, сначала необходимо
выделить узлы на верхней грани полушара:
Plot Nodes
Select Entities Nodes By Num/Pick From Full OK, выбрать флажок Box, обвести рамкой узлы, лежащие на верхней грани, нажать OK
Суммарную силу, действующую со стороны элементов на узлы верхней грани (а это и будет P), можно найти, читая подряд все промежуточные решения, и давая каждый раз команду
General Postproc Nodal Calcs Total Force Sum, OK (списать значение FY)
Однако гораздо быстрее с этой же задачей справится второй скрипт:
*cfopen,'force.txt' ! Открываем файл force.txt
SET,FIRST ! Переходим к первому шагу решения
*do,i,1,10,1 ! Цикл по i от 1 до 10 с шагом 1
fsum ! Находим суммарную силу, действующую
! со стороны выделенных элементов
! на выделенные узлы
*get,qq,FSUM,FY ! Записываем Y-компоненту этой силы в
! переменную qq
*vwrite,qq ! Записываем значение в файл
(E10.3) ! в этом формате
SET,NEXT ! Переходим к следующему шагу решения
*enddo ! Конец цикла
*CFCLOS ! Закрываем файл force.txt
В результате запуска скрипта образуется файл force.txt, в который записываются нужные нам значения P на промежуточных шагах.
Замечание. Файл с этим скриптом не надо называть именем force.txt (force.txt — это файл, который создается этим скриптом). Используйте для файла скрипта, например, имя CalcForce.txt.
После вывода необходимых результатов необходимо выделить все объекты (Utility Menu Select Everything.), после чего повторить расчет на другой сетке.
Файлы displ.txt и force.txt надо переименовать, чтобы они не были случайно стерты при анализе результатов второй задачи (а именно, при запуске скриптов).