memlab (1051146), страница 5
Текст из файла (страница 5)
antype,static
nlgeom,on ! Включение учета больших перемещений
! (геометрической нелинейности)
arcl,on ! Включение Arc-Length метода
nsub,50 ! Задание параметра NSUBST=50 для
! определения радиуса в Arc-Length методе как
! макс. значение нагрузки / NSUBSP
outres,all,all ! Запись всех выходных данных
! на всех шагах решения в базу данных
dk,1,ux,0,,,rotz ! Условия симметрии в вершине купола арки
dk,2,ux,0,,,uy,rotz ! Жесткая заделка опорной точки 2 арки
!f,1,fy,-vforce ! Задание сжимающей силы в верхней точке (1) купола
sfbeam,all,,pres,vpres ! Задание давления
/pbc,all,1 ! Показ всех граничных условий
/psf,pres,norm,1 ! Показ распределенных нагрузок
eplot ! Графический вывод конечно-элементной модели
solve ! Запуск на счет
finish
/post1 ! Обработка результатов построцессором post1
set,last ! Считывание результатов для последнего шага
pldisp,1 ! Изображение деформированного состояния,
! наложенного на недеформированное
finish
/post26 ! Обработка результатов построцессором post26
nsol,2,1,u,y ! Определение переменной номер 2 как
! вертикального перемещения uy в узле 1
abs,2,2 ! Переменная 2 := |переменная 2|
xvar,2 ! Выбор переменной номер 2 в качестве оси x
/axlab,x,Deflection ! Задание подписи Deflection под осью x
!prod,3,1,,,Force,,,2*vforce ! Умножение переменной 1 на 2*vforce
! для получения значений силы
! и присвоение переменной 1 метки Force
!/axlab,y,Force ! Задание подписи Force под осью y
/axlab,y,Pressure ! Задание подписи Force под осью y
prod,3,1,,,Pres,,,2*vpres ! Умножение переменной 1 на 2*vpres
! для получения значений давления
! и присвоение переменной 1 метки Pres
plvar,3 ! Вывод графика переменной 3 от 2
! (кривой "нагрузка-прогиб")
/output,f_r,res,,append
*vwrite
(1X,' Значения точек кривой "нагрузка-прогиб"')
prvar,2,3 ! Печать значений точек предыдущего графика
/output
Результатами выполнения программы Buckl4.inp являются график равновесных состояний «нагрузка-прогиб» (Рис.5.3), где прогиб вычисляется в верхней точке купола, и значения точек этого графика в текстовом файле вывода <f_r>.res. По этим значениям находим, что величина критической нагрузки q, при которой происходит прохлопывание арки есть q=133.557 кГ/см2, что достаточно близко к величине критической нагрузки q1 , полученной по программе Buckl3.inp. (Такая близость результатов вовсе не характерна для задач устойчивости, в чем можно убедиться, проведя вычисления в программах Buckl3.inp и Buckl4.inp для сосредоточенной силы.) Считывая данные для подшагов до прохлопывния арки и после прохлопывания, можно получить на ANSYS графики масштабированных деформируемых состояний до (рис.5.2б) и после (рис.5.2в) хлопка.
Рис. 5.2а. Рис. 5.2б. Рис. 5.3в.
Из сравнения рис.5.2а и 5.2в видно, что форма волнообразования, представленная на рис.5.2а после хлопка не реализуется.
Рис. 5.3.
Используя приведенные программы, аналогичным образом можно исследовать и другие задачи устойчивости деформируемых сред, в частности задачи о прохлопывании мембран и др. При этом необходимо выбрать подходящие КЭ ANSYS [13,14], убедившись что они обладают требуемыми для анализа устойчивости свойствами. Имея рабочие программы, можно также провести различные расчеты при варьировании входных параметров и убедиться в существенной зависимости получаемых результатов от параметров тонкостенности, пологости и др.
Л И Т Е Р А Т У Р А
-
Гроп Д. Методы идентификации систем. М.: Мир, 1979. 294с.
-
Эйкофф П. Основы идентификации систем управления. М.: Мир, 1975. 676с.
-
Ферстер Э., Ренц Б. Методы корреляционного и регрессионного анализа. М.: Финансы и статистика, 1983. 302с.
-
Хампель Ф., Рончетти Э., Рауссеу П. Робастность в статистике. М.: Мир, 1989. 512с.
-
Хьюбер Дж.П. Робастность в статистике. М.: Мир, 1984. 304с.
-
Вольмир А.С. Устойчивость деформируемых систем. М.: Наука, 1967.
-
Валишвилли Н.В. Методы расчета оболочек вращения на ЭЦВМ. М.: Машиностроение, 1976.
-
Григолюк Э.И., Кабанов В.В. Устойчивость оболочек. М.: Наука, 1978. 360с.
-
Статика и динамика сложных структур: Прикладные многоуровневые методы исследований / Вольмир А.С., Куранов Б.А., Турбаивский А.Т. М.: Машиностроение, 1989. 248с.
-
Галлагер Р. Метод конечных элементов. Основы. М.: Мир, 1984. 428с.
-
ANSYS. Basic Analysis Procedures Guide. Rel. 5.3. / ANSYS Inc. Houston, 1994.
-
ANSYS. Commands Reference. Rel. 5.3. / ANSYS Inc. Houston, 1994.
-
ANSYS. Elements Reference. Rel. 5.3. / ANSYS Inc. Houston, 1994.
-
ANSYS. Theory Reference. Rel. 5.3. Ed. P. Kothnke / ANSYS Inc. Houston, 1994.
-
ANSYS. Verification Manual. Rel. 5.3. / ANSYS Inc. Houston, 1994.















