! ********************************************************************
!
!  Данная программа входит в набор программ SMs1.inp - SMs10.inp 
!                          учебного пособия:
!   Наседкин А.В. "Конечно-элементное моделирование на основе ANSYS.
!   Программы решения статических задач сопротивления материалов
!   c вариантами индивидуальных заданий" //
!   Ростов-на-Дону: УПЛ РГУ, 1998. 44 с. 
!   
!  Автор программ - Наседкин Андрей Викторович,
!   Ростовский госуниверситет, кафедра математического моделирования,
!   nasedkin@ms.math.rsu.ru или nasedkin@ns.math.rsu.ru
!
! ********************************************************************
/batch           ! Пакетный режим /BATCH
/com,
/com,                 Файл SMs5.inp
/com,    ПРОГРАММА РЕШЕНИЯ ЗАДАЧИ СОПРОТИВЛЕНИЯ МАТЕРИАЛОВ
/com,               ПО ТЕМЕ: "КРИВЫЕ СТЕРЖНИ"
/com,        (версия 1.1 программы без построения эпюр)
/com,
/com,   Сборник задач по сопротивлению материалов /
/com,   Под ред. В.К. Качурина. М.: Наука, 1970.
/com,            (по задаче 11.18.    С.331.)
/com,
/com,   Итоговые результаты записываются 
/com,   в стандартный файл вывода *.out
/com,
/title, CURVED BAR 
/prep7           ! Входим в препроцессор /PREP7
                 ! Все данные в системе СГС
dr=2          ! dr - диаметр стержня
ar=10         ! ar - радиус кривизны
et,1,pipe18      ! КЭ - искривленные трубки PIPE18
keyopt,1,6,2     ! печать внутренних усилий
r,1,dr,dr/2,ar   ! REAL CONST :
                 !   внешний диаметр сечения OD=dr,
                 !   толщина стенок трубы TKWALL=dr/2 
                 !   (сплошная труба)
                 !   радиус кривизны RADCUR=ar
mp,ex,1,2e6      ! модуль Юнга EX=2e6
mp,gxy,1,8e5     ! модуль сдвига GXY=8e5
n,1              ! Определяем узлы и элементы
n,2,-ar
n,3,0,ar
n,4,ar
e,2,3,1 $ e,3,4,1
finish  
/solu            ! Входим в решатель /SOLUTION
antype,static    ! Статический тип анализа
d,all,uz,0       ! Задаем условия закрепления в узлах
d,2,ux,0 $ d,2,uy,0
d,4,ux,0 $ d,4,uy,0
f,3,fy,-100      ! Задаем силу в узле 2+n_fill
outpr,basic,1    ! Вывод основных результатов в стандартный файл
solve            ! Решаем СЛАУ МКЭ
finish 
/exit