! ********************************************************************
!
!  Данная программа входит в набор программ SMs1.inp - SMs10.inp 
!                          учебного пособия:
!   Наседкин А.В. "Конечно-элементное моделирование на основе ANSYS.
!   Программы решения статических задач сопротивления материалов
!   c вариантами индивидуальных заданий" //
!   Ростов-на-Дону: УПЛ РГУ, 1998. 44 с. 
!   
!  Автор программ - Наседкин Андрей Викторович,
!   Ростовский госуниверситет, кафедра математического моделирования,
!   nasedkin@ms.math.rsu.ru или nasedkin@ns.math.rsu.ru
!
! ********************************************************************
/batch           ! Пакетный режим /BATCH
/com,
/com,                 Файл SMs4.inp
/com,    ПРОГРАММА РЕШЕНИЯ ЗАДАЧИ СОПРОТИВЛЕНИЯ МАТЕРИАЛОВ
/com,               ПО ТЕМЕ: "КРУЧЕНИЕ"
/com,                  (версия 1.1)
/com,
/com,   Сборник задач по сопротивлению материалов /
/com,   Под ред. В.К. Качурина. М.: Наука, 1970.
/com,            (по задаче 4.33.    С.94.)
/com,
/com,       Два варианта конечных элементов - BEAM4 и PIPE16
/com,
/com,   Итоговые результаты записываются 
/com,   в стандартный файл вывода *.out
/com,
/title, TORSION OF A CIRCULAR ROD 
/prep7           ! Входим в препроцессор /PREP7
                 ! Все данные в системе СГС
pi=3.1416
dr=10         ! dr - диаметр стержня
ar=60         ! ar, b2r - геометрические параметры
b2r=40
jpol=pi*dr**4/32 ! константа для полярного момента инерции
et,1,beam4       ! балочные КЭ BEAM4
r,1,pi*dr**2/4,jpol/2   ! REAL CONST :
                 !   площадь поперечного сечения AREA=pi*dr**2/4 ,
                 !   момент инерции IZZ=jpol/2 
rmore,,jpol      !  полярный момент инерции 
                 !  (8-ая REAL-константа для BEAM4)
mp,ex,1,2e6      ! модуль Юнга EX=2e6
mp,gxy,1,8e5     ! модуль сдвига GXY=8e5
n,1              ! Определяем узлы по координатам
n,2,ar
n,3,ar+b2r
n,4,ar+b2r*2
n,5,(ar+b2r)*2
*do,i,1,4        ! Определяем элементы BEAM4
  e,i,i+1
*enddo
finish  
/solu            ! Входим в решатель /SOLUTION
antype,static    ! Статический тип анализа
d,1,all,0,,5,4   ! Задаем условия закрепления в узлах
f,2,mx,8e4       ! Задаем крутящие моменты в узлах
f,4,mx,8e4           
outpr,basic,1    ! Вывод основных результатов в стандартный файл
solve            ! Решаем СЛАУ МКЭ
finish 
/prep7           ! Вариант с КЭ PIPE16
et,1,pipe16      ! КЭ - трубки PIPE16
r,1,dr,dr/2      ! REAL CONST :
                 !   внешний диаметр сечения OD=dr,
                 !   толщина стенок трубы TKWALL=dr/2 
                 !   (сплошная труба)
finish
/solu
solve
finish
/exit