Шаг 11. Получение зависимости прижимающей силы от перемещения
При решении этой задачи фактически было решено еще девять — с промежуточными значениями перемещения u верхней грани. По окончании расчета автоматически загружается решение самой последней из этих задач. Результаты решения промежуточных задач можно просмотреть отдельно (команды General Postproc à Read Results à First Set, Next Set и т. д.).
На каждом шаге представляют интерес перемещение u и прижимающая сила P. Их можно найти, последовательно перебирая решения всех задач, записанных в файл результатов. Однако это долгий и неудобный путь. Гораздо быстрее получить требуемые величины, воспользовавшись скриптом — последовательностью команд ANSYS. Скрипт записывается в отдельный файл и затем запускается командой Utility Menu à File à Read input from.
Получение значения перемещения u на всех шагах решения
Сначала необходимо узнать номер какого-нибудь (любого) узла на верхней грани шара:
1. Показать номера всех узлов: Utility Menu à PlotCtrls à Numbering à Node numbers (On) à OK
2. увеличить верхнюю грань, чтобы увидеть номер какого-нибудь узла на ней (обозначим этот номер N — его надо запомнить): Utility Menu à PlotCtrls à Pan Zoom Rotate
3. Убрать изображение номеров узлов (Utility Menu à PlotCtrls à Numbering à Node numbers (Off) à OK)
4. Изобразить всю модель: 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 надо переименовать, чтобы они не были случайно стерты при анализе результатов второй задачи (а именно, при запуске скриптов).
Дата добавления: 2015-12-20; просмотров: 19; Мы поможем в написании вашей работы! |
Мы поможем в написании ваших работ!