! ********************************************************************
!
!  Данная программа входит в набор программ SMs1.inp - SMs10.inp 
!                          учебного пособия:
!   Наседкин А.В. "Конечно-элементное моделирование на основе ANSYS.
!   Программы решения статических задач сопротивления материалов
!   c вариантами индивидуальных заданий" //
!   Ростов-на-Дону: УПЛ РГУ, 1998. 44 с. 
!   
!  Автор программ - Наседкин Андрей Викторович,
!   Ростовский госуниверситет, кафедра математического моделирования,
!   nasedkin@ms.math.rsu.ru или nasedkin@ns.math.rsu.ru
!
! ********************************************************************
/batch           ! Пакетный режим /BATCH
/com,
/com,                 Файл SMs10.inp
/com,    ПРОГРАММА РАСЧЕТА ЗАДАЧИ СОПРОТИВЛЕНИЯ МАТЕРИАЛОВ
/com,    ПО ТЕМАМ: "РАСТЯЖЕНИЕ СТЕРЖНЕЙ, ТЕМПЕРАТУРНЫЕ ДЕФОРМАЦИИ"
/com,                  (версия 1.1)
/com,
/com,   Сборник задач по сопротивлению материалов /
/com,   Под ред. В.К. Качурина. М.: Наука, 1970.
/com,            (по задаче 1.89в.  С.40.)
/com,
/com,   Текстовые результаты сохраняются в файле SMs10.res
/prep7           ! Входим в препроцессор /PREP7
f_r='SMs10'
antype,static    ! Статический тип анализа
                 ! Все данные в системе СГС
ax=50             ! Геометрический размер 
ar1=20 $ ar2=10   ! Площади поперечных сечений
mp,ex,1,0.8e6      ! MAT=1: медь
mp,alpx,1,165e-7
mp,ex,2,2e6        ! MAT=2: сталь
mp,alpx,2,125e-7
et,1,link1         ! TYPE=1 : стержневой КЭ LINK1
r,1,ar1 $ r,2,ar2  ! REAL 1 и 2 с разными площадями
n,1,0 $ n,2,ax $ n,3,3*ax ! Создаем узлы
! Создаем КЭ с различными MP и Real Const. 
mat,1 $ real,1 $ e,1,2
mat,2 $ real,2 $ e,2,3
d,all,uy,0    ! Считаем перемещения UY для всех узлов равными 0
d,1,ux,0,,3,2      ! Закрепляем узлы 1 и 3 по оси X
! tref,0  ! Так можно задать начальную температуру TREF=0
          ! (TREF=0 по умолчанию)
bfunif,temp,40     ! Задаем постоянную температуру T: T-TREF=40
! Другие способы задания температуры T (T-TREF=40):
! bf,all,temp,40    !  - для всех узлов 
! bfe,all,temp,1,40 !  - для всех элементов
finish  
/solu            ! Входим в решатель /SOLUTION
solve             ! Решаем СЛАУ
finish  
/post1
 ! Согласно документации по КЭ LINK1 
 ! определяем указатели на значения осевых усилий и
 ! напряжений в стержнях
etable,forax,smisc,1  ! Smisc,1 - осевое усилие 
etable,sigax,ls,1     ! Ls,1 - осевое напряжение
/output,F_R,res       ! Направляем вывод в файл <F_R>.res
pretab,forax,sigax    ! Печатаем усилия и напряжения
/output
finish 
/exit