Программа s2.F
* формирование глобальной матрицы жесткости конструкции
dimension cord(100,3),nop(144,4),
*estifm(12,12),rq(300),nbc(50),nfix(50),ndft(1),net(1)
*,ncnt(1),y(51),s(300),a(50,50)*,npr(300),npr1(300)
Продолжение приложения 4
open(16,file='RES2')
open(7,file='NPR',form='unformatted')
open(29,file='NOP',form='unformatted')
READ(29)NOP,CORD,rq,nbc,nfix
close(29)
open(2,file='DT2',form='unformatted')
read(2,*)ne,np,nszf,nb,nband,imax,imx1,nt
close(2)
do 33 i=1,nb
33 write(16,*)'nbc',nbc(i),'nfix',nfix(i)
do 44 i=1,nszf
44 write(16,*)'i=',i,'rq',rq(i)
net(1)=ne
ndft(1)=3
ncnt(1)=4
WRITE(16,8)(NET(I),NDFT(I),NCNT(I),I=1,NT)
8 FORMAT(12I5)
WRITE(16,9)(NBC(I),NFIX(I),I=1,NB)
9 FORMAT(2i8)
ndfg=3
ncnm=4
call FORM(imax,npr,nband,nszf,ne,ndfg,ncnm,nt,
*net,ndft,ncnt,nop,estifm,rq,a,nb,nfix,nbc,ik)
write(7)ik,rq,npr
end file(7)
close(7)
end file(16)
close(16)
stop
end
Программа s3.F
* разложение и решение системы линейных уравнений
dimension nop(144,4),rq(300),y(51),s(300),
*a(50,50),npr(300),npr1(300)
open(1,file='AGLOB',form='unformatted')
open(2,file='F2RES',form='unformatted')
open(16,file='RES3')
open(7,file='NPR',form='unformatted')
read(7)ik,rq,npr
close(7)
open(2,file='DT2',form='unformatted')
read(2,*)ne,np,nszf,nb,nband,imax,imx1,nt
close(2)
ndfg=3
ij=nszf/imax
Продолжение приложения 4
if(ik.ne.imax) ij=ij+1
call RASL(npr,imax,nband,ij,ik,nszf,imx1,y,s,npr1,rq,a)
call RES(nband,imax,ik,a,rq,nszf)
do 7 i=1,np
7 write(16,4)i,(rq((i-1)*ndfg+j),j=1,ndfg)
4 format(i5,3e10.2)
open(3,file='DIS',form='unformatted')
write(3)rq
end file(3)
close(3)
close(2,status='delete')
close(1)
end file(16)
close(16)
stop
end
Программа s4.F
* вычисление напряжений в элементах
dimension nop(144,4),rq(300),force(144,7)
open(16,file='RES4')
open(9,file='NOP',form='unformatted')
read(9)nop
close(9)
open(12,file='B',form='unformatted')
open(3,file='DIS',form='unformatted')
read(3) rq
close(3)
ndfg=3
open(2,file='DT2',form='unformatted')
read(2,*)ne,np,nszf,nb,nband,imax,imx1,nt
close(2)
write(16,*)' УЗЛОВЫЕ ПЕРЕМЕЩЕНИЯ'
write(16,*)' Узел X Y Z '
do 7 i=1,np
7 write(16,4)i,(rq((i-1)*ndfg+j),j=1,ndfg)
4 format(i5,3e10.2)
DO 3 N=1,NE
call STR31(n,rq,FORCE,NSZF,NE,NDFG,nop)
3 CONTINUE
close(9)
close(12)
end file(16)
close(16)
stop
end
Библиографический список
1. Зенкевич О.С. Метод конечных элементов в технике / О.С. Зенкевич – М.: Мир, 1975. 541 с.
2. Галлагер Р. Метод конечных элементов. Основы / Р. Галлагер – М.: Мир, 1984. 428 с.
3. Сегерлинд Л. Применение метода конечных элементов./ Л. Сегерлинд – М.: Мир, 1979. 393 с.
4. Firooz Ghassemit Automatic mesh generation scheme a two- or three-dimensional triangular curved surface // Computers & Structures, No.6,pp.613 – 626, 1982.
5. Cavendish J.C. Automatic triangulation of arbitrary planer domains for the finite element method – Int. J. for Numerical Methods in Engineering, vol.8, (1974).
6. William A. Cook Body oriented (natural) coordinates for generating 3D meshes – Int. J. for Numerical Methods in Engineering, vol.8, pp. 27 – 43, (1974)
dimensional triangular curved surface // Computers & Structures, No.6,pp.613 – 626, 1982.
7. Бате К.Численные методы анализа и метод конечных элементов/ К. Бате, Е. Вилсон – М.: Стройиздат, 1982. 446 с.
8. Pro/MECHANICA. Design Study Reference. – Waltham: Parametric Technology Corporation, 1997. 394 p.
9. Филиппов А.П. Колебания деформируемых систем/ А.П.Филиппов – М.: Машиностроение, 1970. 736 с.
10. Рикардс Р.Б. Метод конечных элементов в теории оболочек и пластин/ Р.Б. Рикардс– Рига.: Зинатне, 1988. 284 с.
11. Ramamurti V., Pattabiraman J. Dynamic behaviour of cilindrical shell with a cutout. – J. Sound and Vibration, 1977, v.52, No.2, – P. 193–200.
ОГЛАВЛЕНИЕ
ВВЕДЕНИЕ 7
ОСНОВНАЯ ИДЕЯ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ 8
ПРЕИМУЩЕСТВА И НЕДОСТАТКИ МКЭ 9
ДИСКРЕТИЗАЦИЯ ОБЛАСТИ 10
ТИПЫ КОНЕЧНЫХ ЭЛЕМЕНТОВ 15
ПРЯМОЙ МЕТОД ЖЕСТКОСТИ 19
УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ 24
АЛГОРИТМЫ ПОСТРОЕНИЯ СЕТОК ДЛЯ РЕШЕНИЯ ЗАДАЧ МЕХАНИКИ ДЕФОРМИРУЕМЫХ ТВЕРДЫХ ТЕЛ 25
СООТНОШЕНИЯ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ В ЗАДАЧАХ ДИНАМИКИ 39
МАТРИЦА ИНЕРТНОСТИ ТРЕУГОЛЬНОГО КОНЕЧНОГО ЭЛЕМЕНТА 41
ОПИСАНИЕ ПРОГРАММЫ РАСЧЕТА ПО МЕТОДУ КОНЕЧНЫХ ЭЛЕМЕНТОВ 46
ПРИМЕР ИСПОЛЬЗОВАНИЯ ПРОГРАММЫ ОПРЕДЕЛЕНИЯ СОБСТВЕННЫХ ЧАСТОТ ТОНКОСТЕННЫХ КОНСТРУКЦИЙ 47
ПРИМЕРЫ ПРАКТИЧЕСКОГО ИСПОЛЬЗОВАНИЯ НЕКОТОРЫХ ТИПОВ КОНЕЧНЫХ ЭЛЕМЕНТОВ ПРИ ИССЛЕДОВАНИИ СТАТИЧЕСКИХ И ДИНАМИЧЕСКИХ СОСТОЯНИЙ КОНСТРУКЦИЙ 53
Пространственные стержневые конструкции 53
Плоская задача теории упругости 57
Построение матрицы жесткости пластинки прямоугольной формы 65
Моделирование оболочечных конструкций 76
Моделирование массивных конструкций 95
ЗАКЛЮЧЕНИЕ 103
Приложение 1 105
Приложение 2 112
Приложение 3 122
Приложение4 126
БИБЛИОГРАФИЧЕСКИЙ СПИСОК 130
Учебное издание
Рукин Юрий Борисович
ОСНОВЫ ПРИМЕНЕНИЯ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ
Издание второе, переработанное и дополненное
В авторской редакции
Компьютерный набор Р.А. Жилина
Подписано к изданию 20.10.2006.
Уч.-изд. л. 7,0.
ГОУВПО «Воронежский государственный технический
университет»