! ********************************************************************
!
!  Данная программа входит в набор программ SMs1.inp - SMs10.inp 
!                          учебного пособия:
!   Наседкин А.В. "Конечно-элементное моделирование на основе ANSYS.
!   Программы решения статических задач сопротивления материалов
!   c вариантами индивидуальных заданий" //
!   Ростов-на-Дону: УПЛ РГУ, 1998. 44 с. 
!   
!  Автор программ - Наседкин Андрей Викторович,
!   Ростовский госуниверситет, кафедра математического моделирования,
!   nasedkin@ms.math.rsu.ru или nasedkin@ns.math.rsu.ru
!
! ********************************************************************
/batch           ! Пакетный режим /BATCH
/com,
/com,                 Файл SMs9.inp
/com,    ПРОГРАММА РАСЧЕТА ЗАДАЧИ СОПРОТИВЛЕНИЯ МАТЕРИАЛОВ
/com,    ПО ТЕМАМ: "РАСТЯЖЕНИЕ СТЕРЖНЕЙ, НАЧАЛЬНЫЕ ДЕФОРМАЦИИ,
/com,               УРАВНЕНИЯ СВЯЗЕЙ"
/com,                  (версия 1.1)
/com,
/com,   Сборник задач по сопротивлению материалов /
/com,   Под ред. В.К. Качурина. М.: Наука, 1970.
/com,            (по задаче 1.100.  С.44.)
/com,
/com,   Текстовые результаты сохраняются в файле SMs9.res
/prep7           ! Входим в препроцессор /PREP7
f_r='SMs9'
antype,static    ! Статический тип анализа
                 ! Все данные в системе СГС
ax=100 $ ay=100  ! Геометрические размеры 
delt=0.05        ! Начальное удлинение второго стержня
area=20          ! Площади поперечных сечений
mp,ex,1,2e6    ! MAT=1: модуль Юнга EX=2e6
et,1,link1     ! TYPE=1 : стержневой КЭ LINK1
r,1,area           ! REAL=1 : area
r,2,area,delt/(ay-delt) ! REAL=2 : area, 
                        ! ISTRN=delt/(ay-delt) - нач. деформация
n,1,ax $ n,2,2*ax $ n,3,3*ax ! Создаем узлы
n,4,ax,-ay $ n,5,2*ax,-ay $ n,6,3*ax,-ay
! Создаем КЭ с различными Real Const.
real,1 $ e,1,4
real,2 $ e,2,5
real,1 $ e,3,6
d,all,ux,0    ! Считаем перемещения UX для всех узлов равными 0
d,1,uy,0,,3   ! Закрепляем верхние узлы (1,2,3) по оси Y
ce,1,0,5,uy,1,4,uy,-2 ! Задаем 1-ое уравнение связи: UY_5+(-2)*UY_4=0
ce,2,0,6,uy,1,4,uy,-3 ! Задаем 2-ое уравнение связи: UY_6+(-3)*UY_4=0
save
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