! ********************************************************************
!
!  Данная программа входит в набор программ SMs1.inp - SMs10.inp 
!                          учебного пособия:
!   Наседкин А.В. "Конечно-элементное моделирование на основе ANSYS.
!   Программы решения статических задач сопротивления материалов
!   c вариантами индивидуальных заданий" //
!   Ростов-на-Дону: УПЛ РГУ, 1998. 44 с. 
!   
!  Автор программ - Наседкин Андрей Викторович,
!   Ростовский госуниверситет, кафедра математического моделирования,
!   nasedkin@ms.math.rsu.ru или nasedkin@ns.math.rsu.ru
!
! ********************************************************************
/batch           ! Пакетный режим /BATCH
/com,
/com,                 Файл SMs3.inp
/com,    ПРОГРАММА РЕШЕНИЯ ЗАДАЧИ СОПРОТИВЛЕНИЯ МАТЕРИАЛОВ
/com,           ПО ТЕМЕ: "РАСЧЕТ ПЛОСКИХ РАМ"
/com,                  (версия 1.1)
/com,
/com,   Сборник задач по сопротивлению материалов /
/com,   Под ред. В.К. Качурина. М.: Наука, 1970.
/com,            (по задаче 9.36.    С.246-247.)
/com,
/com,   Текстовые результаты сохраняются в файле SMs3.res
/com,   Графические результаты выводятся в файл  SMs3.grph
/com,
/title, PLANE FRAME
F_R='SMs3'
/prep7           ! Входим в препроцессор /PREP7
/show,F_R,grph   ! В режиме /BATCH задаем вывод графики
                 ! в файл <F_R>.grph
                 ! Все данные в системе СГС
et,1,beam3       ! Элементы с TYPE=1 - , балочные КЭ BEAM3
r,1,37.5,3800,24 ! REAL CONST. типа REAL=1 : (двутавр N 20а)
                 !   площадь поперечного сечения AREA=28.9 ,
                 !   момент инерции IZZ=2030 ,
                 !   толщина HEIGHT=20
mp,ex,1,2e6      ! MAT. PROP. типа MAT=1: модуль Юнга EX=2e6
lr=500   ! lr, ar, cr - геометрические параметры
ar=400
cr=300
n_fill=16      ! число КЭ, генерируемых 
               ! на участке с распределенной нагрузкой
n,1              ! Определяем узлы по координатам
n,2,,lr
n,3,,lr
n,3+n_fill,ar,lr
n,4+n_fill,ar,lr-cr
fill,3,3+n_fill,n_fill-1  ! Генерируем n_fill-1 равноотстоящих узлов
                          ! между узлами 3 и 3+n_fill
e,1,2            ! Определяем элементы BEAM3
*do,i,3,3+n_fill
  e,i,i+1
*enddo
cp,1,ux,2,3 $ cp,2,uy,2,3   ! Связываем степени свободы UX и UY
                  ! в совпадающих узлах 2 и 3 для задания
                  ! промежуточного шарнира
finish  
/solu            ! Входим в решатель /SOLUTION
antype,static    ! Статический тип анализа
d,1,uy           ! Задаем условия закрепления
d,1,ux
d,1,rotz
d,4+n_fill,uy
d,4+n_fill,ux           
*do,i,2,1+n_fill
 sfbeam,i,,pres,40  ! Задаем распределенную нагрузку
*enddo
     ! Задание сил FY, FX и моментов MZ проводится аналогично
     ! описанному в программе SMs2.inp
solve             ! Решаем СЛАУ МКЭ
save
finish 
/post1
 ! Согласно документации по КЭ BEAM3 при KEYOPT(9)=0
 ! определяем указатели на значения нормальных и поперечных сил и
 ! изгибающих моментов на концах элементов
etable,mforxi,smisc,1  ! Smisc,1 - нормальная сила N в узле I
etable,mforxj,smisc,7  ! Smisc,7 - нормальная сила N в узле J
etable,mforyi,smisc,2  ! Smisc,2 - поперечная сила Q в узле I
etable,mforyj,smisc,8  ! Smisc,8 - поперечная сила Q в узле J
etable,mmomzi,smisc,6  ! Smisc,6 - изгибающий момент M в узле I
etable,mmomzj,smisc,12 ! Smisc,12 - изгибающий момент 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,cp,,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 
  ! Поперечные силы в ANSYS имеют противоположный знак 
  ! по сравнению с принятым в указанном выше задачнике
plls,mforxi,mforxj,2   ! Графический показ эпюры нормальных сил N
plls,mmomzi,mmomzj,2   ! Графический показ эпюры изгибающих моментов M
finish 
/exit