!**************************************************************
!        Файл Buckl4.inp
!   Прохлопывание пологой арки.
!   Нелинейный анализ устойчивости.
! Строится кривая равновесных состояний
! "нагрузка-прогиб" для верхней точки купола арки
!**************************************************************
f_r='Buckl_n'        ! имя файлов для вывода результатов
!vforce=450 ! Значение максимальной силы (половинное значение)
vpres=80	! Значение максимального давления (половинное значение)
nldiv=20          ! Параметр для разбиения линий
! Все входные данные задаются в системе СГС
hl=0.3            ! Высота поперечного сечения арки
bl=1              ! Толщина поперечного сечения арки
hg=3.2            ! Высота купола арки
darc=20           ! Диаметр основания арки

rarc=hg/2+darc*darc/8/hg ! Радиус кривизны арки
pi=4*atan(1)  ! Число PI
alfa=acos((rarc-hg)/rarc)*180/pi ! alfa - угол раствора половины арки

/output,f_r,res	! Сервисный вывод в файл <f_r>.res
*vwrite
(1X,'  Нелинейный анализ')
*vwrite
(1X,'  потери устойчивости пологой арки.')
*vwrite
(1X,'  Nonlinear Buckling Analysis ')
*vwrite
(1X)
*vwrite,darc,hg,hl,bl
(1X,'  DARC=', E12.5,'  HG=', E12.5,'  HL=', E12.5,' BL=', E12.5)
*vwrite,rarc
(1X,'  RARC=', E12.5)
*vwrite
(1X)
/output		! Окончание сервисного вывода

/prep7		
et,1,beam3			! КЭ BEAM3
mp,ex,1,2.02e6		! Модуль Юнга EX=2.02e6
mp,nuxy,1,0.3           ! Коэффициент Пуассона NUXY=0.3
r,1,bl*hl,(bl*hl**3)/12,hl	! R. const для поперечного сечения балки:
					! площадь, момент инерции, толщина
csys,1			! Переход в цилиндрическую систему координат
k,1,rarc,90			! Точка 1 - вершина купола
k,2,rarc,90-alfa		! Точка 2 - основание арки
n,1,rarc,90             ! Узел 1 будет далее узлом приложения силы,
n,2,rarc,90-alfa		! узел 2 - основанием арки
l,1,2,nldiv			! Генерация формы половины арки
! (В цилиндрической системе коорд. линии между точками
!  с одинаковыми радиусами есть дуги окружности)
lmesh,all		! Генерация сетки КЭ на линии
nummrg,node       ! Узлы 1,2 включаются в КЭ модель
finish		
/solu			
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