! ********************************************************************
!
!  Данная программа входит в набор программ SMs1.inp - SMs10.inp 
!                          учебного пособия:
!   Наседкин А.В. "Конечно-элементное моделирование на основе ANSYS.
!   Программы решения статических задач сопротивления материалов
!   c вариантами индивидуальных заданий" //
!   Ростов-на-Дону: УПЛ РГУ, 1998. 44 с. 
!   
!  Автор программ - Наседкин Андрей Викторович,
!   Ростовский госуниверситет, кафедра математического моделирования,
!   nasedkin@ms.math.rsu.ru или nasedkin@ns.math.rsu.ru
!
! ********************************************************************
/batch           ! Пакетный режим /BATCH
/com,
/com,                 Файл SMs8.inp
/com,    ПРОГРАММА РАСЧЕТА ЗАДАЧИ СОПРОТИВЛЕНИЯ МАТЕРИАЛОВ
/com,    ПО ТЕМЕ: "ВНЕЦЕНТРЕННОЕ РАСТЯЖЕНИЕ-СЖАТИЕ"
/com,                  (версия 1.1)
/com,
/com,   Сборник задач по сопротивлению материалов /
/com,   Под ред. В.К. Качурина. М.: Наука, 1970.
/com,            (по задаче 10.26.    С.272.)
/com,
/com,   Текстовые результаты сохраняются в файле SMs8.res
/prep7           ! Входим в препроцессор /PREP7
f_r='SMs8'
antype,static    ! Статический тип анализа
                 ! Все данные в системе СГС
h=18           ! Высота поперечного сечения балки
b=0.73         ! Толщина поперечного сечения балки
lbeam=50       ! Длина половины балки
exbeam=2e6     ! Модуль Юнга
mp,ex,1,exbeam ! MAT. PROP. типа MAT=1: модуль Юнга EX=exbeam
et,1,beam54    ! Элементы с TYPE=1 - , балочные КЭ BEAM54
n,1            ! Создаем узлы
n,2,lbeam
area1=b*h      ! Площадь поперечного сечения
iz1=b*h**3/12  ! Момент инерции поперечного сечения
dy1=4          ! Расстояние по оси Y от узла I до центр. оси
pf=9e3         ! Значение силы
r,1,area1,iz1,h/2,h/2,area1,iz1 ! Real Const. с 1 по 6 для балки
rmore,h/2,h/2,,dy1,,dy1         ! Real Const. с 7 по 12 для балки
e,1,2          ! Создаем один КЭ
f,2,fx,pf      ! Прилагаем растягивающую силу
save
finish  
/solu            ! Входим в решатель /SOLUTION
d,1,ux,0 $ d,1,rotz,0  ! Условия, моделирующие симметрию задачи
solve             ! Решаем СЛАУ МКЭ
finish  
/post1
 ! Согласно документации по КЭ BEAM54 
 ! определяем указатели на значения макс. и мин.
 ! напряжений на концах элемента
etable,smaxi,nmisc,1  ! Nmisc,1 - макс. напряжение 
                    ! на сечении в узле I
etable,smaxj,nmisc,3  ! Nmisc,3 - аналог. напряжение в узле J
etable,smini,nmisc,2  ! Nmisc,2 - мин. напряжение 
                    ! на сечении в узле I
etable,sminj,nmisc,4  ! Nmisc,4 - аналог. напряжение в узле J

/output,F_R,res   ! Направляем вывод в файл <F_R>.res
pretab,smaxi,smaxj       ! Печатаем макс. напряжения
pretab,smini,sminj       ! Печатаем мин. напряжения
/output
finish 
/exit