! ********************************************************************
!
!  Данная программа входит в набор программ SMs1.inp - SMs10.inp 
!                          учебного пособия:
!   Наседкин А.В. "Конечно-элементное моделирование на основе ANSYS.
!   Программы решения статических задач сопротивления материалов
!   c вариантами индивидуальных заданий" //
!   Ростов-на-Дону: УПЛ РГУ, 1998. 44 с. 
!   
!  Автор программ - Наседкин Андрей Викторович,
!   Ростовский госуниверситет, кафедра математического моделирования,
!   nasedkin@ms.math.rsu.ru или nasedkin@ns.math.rsu.ru
!
! ********************************************************************
/batch           ! Пакетный режим /BATCH
/com,
/com,                 Файл SMs6.inp
/com,    ПРОГРАММА РЕШЕНИЯ ЗАДАЧИ СОПРОТИВЛЕНИЯ МАТЕРИАЛОВ
/com,               ПО ТЕМЕ: "КРИВЫЕ СТЕРЖНИ"
/com,         (версия 1.1 программы с построением эпюр)
/com,
/com,   Сборник задач по сопротивлению материалов /
/com,   Под ред. В.К. Качурина. М.: Наука, 1970.
/com,            (по задаче 11.18.    С.331.)
/com,
/title, CURVED BAR 
F_R='SMs6'
/show,F_R,grph   ! Задаем вывод графики в файл <F_R>.grph
/prep7           ! Входим в препроцессор /PREP7
                 ! Все данные в системе СГС
dr=2          ! dr - диаметр стержня
ar=10         ! ar - радиус кривизны
et,1,pipe18      ! КЭ - искривленные трубки PIPE18
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_fill=20        ! число дополнительных узлов на одном КЭ
csys,1           ! Переходим в цилиндрическую систему координат
n,1              ! Определяем узлы по координатам
n,2,ar,180
n,2+n_fill,ar,90
fill,2,2+n_fill,n_fill-1  ! Генерируем n_fill-1 равноотстоящих узлов
                          ! между узлами 2 и 2+n_fill
n,2+n_fill*2,ar
fill,2+n_fill,2+n_fill*2,n_fill-1  
                          ! Генерируем n_fill-1 равноотстоящих узлов
                          ! между узлами 2+n_fill и 2+n_fill*2
csys,0           ! Возвращаемся в декартову систему координат
*do,i,2,1+n_fill*2
  e,i,i+1,1
*enddo
finish  
/solu            ! Входим в решатель /SOLUTION
antype,static    ! Статический тип анализа
d,all,uz,0       
d,2,ux,0 $ d,2,uy,0  ! Задаем условия закрепления в узлах
d,2+n_fill*2,ux,0 $ d,2+n_fill*2,uy,0
f,2+n_fill,fy,-100      ! Задаем силу в узле 2+n_fill
solve            
finish 
/post1
 ! Согласно документации по КЭ PIPE18 
 ! определяем указатели на значения нормальных и поперечных сил и
 ! изгибающих моментов на концах элементов
etable,mforxi,smisc,1  ! Smisc,1 - нормальная сила N в узле I
etable,mforxj,smisc,7  ! Smisc,7 - нормальная сила N в узле J
etable,mforyi,smisc,3  ! Smisc,3 - поперечная сила Q в узле I
etable,mforyj,smisc,9  ! Smisc,9 - поперечная сила Q в узле J
etable,mmomzi,smisc,5  ! Smisc,5 - изгибающий момент M в узле I
etable,mmomzj,smisc,11 ! Smisc,11 - изгибающий момент M в узле J
/output,F_R,res   ! Направляем вывод в файл <F_R>.res
prrsol            ! Печатаем опорные реакции
pretab,mforyi,mforyj       ! Печатаем поперечные силы
pretab,mforxi,mforxj       ! Печатаем нормальные силы
pretab,mmomzi,mmomzj       ! Печатаем изгибающие моменты
/output
/pnum,node,1  ! Нумеровать узлы при графическом показе
eplot         ! Нарисовать элементы
/pbc,u,,1     ! Показывать граничные условия,силы и моменты
/pbc,f,,1
/pbc,m,,1
/psf,pres,norm,1 ! Показывать распределенные нормальные нагрузки
/pnum,elem,1  ! Нумеровать элементы
/pnum,node,0  ! Не нумеровать узлы
/number,1  ! Занумерованные объекты показывать 
           ! только различным цветом без нумерации
/zoom,1,scrn,,,1.7 ! Масштабировать окно для удобства
/pbc,all,,0    ! Не показывать граничные условия, силы и моменты
/psf,pres,norm,0 ! Не показывать распределенные нагрузки
/pnum,elem,0     ! Не нумеровать элементы
/triad,off       ! Не показывать начало и систему координат 
plls,mforyi,mforyj,2  ! Графический показ эпюры поперечных сил Q 
plls,mforxi,mforxj,2  ! Графический показ эпюры нормальных сил N
plls,mmomzi,mmomzj,-2 ! Графический показ эпюры изгибающих моментов M
finish 
/exit