Метод конечных элементов в расчетах стержневых систем. Дьяков И.Ф - 105 стр.

UptoLike

ВУЗ: 

УлГТУ | Ульяновск

Составители: 

Рубрика: 

105
Примеры подпрограмм
C BЫЧИCЛEHИE BEKTOPA BHУTPEHHИX CИЛOBЫX ФAKTOPOB И НАПРЯЖЕНИЙ
SUBROUTINE VSFST(P)
COMMON /ORT/ CL(3,3),FT(6,6),X,Y,DL
COMMON/MG/FMG(6,6),FMGO(6,6),MAT(20)
COMMON/XAR/IE,IT,T(3),F,ZI,NMAT,WZ,NY
COMMON /RAZM/ KYZO,NS,KELEM,N,LENTA
DIMENSION P(N)
DIMENSION PO(6),FO(6)
REAL MAT
INTEGER T
REWIND 3
WRITE (8,1)
1 FORMAT(/7X,'Bнутренние силовые факторы '/
13X,'Элемент',3X,'Mатрица индексов',4X,
2'Nx',10X,'Qy',10X,'Mz',10X,'G')
DO 5 I=1,KELEM
READ(3)IE,T,NY,NSM,FMGO,FT,F,WZ
NK=NSM+1
NSEM=NY*NSM
NSEO=NY*NS
DO 2 J=1,NY
K1=(T(J)-1)*NS
DO 2 K=1,NS
I1=(J-1)*NS+K
K1=K1+1
2 PO(I1)=P(K1)/1.E5
CALL MPRD(FMGO,PO,FO,NSEO,NSEO,1,6,6,6,1,6,1)
CALL MPRD(FT,FO,PO,NSEM,NSEO,1,6,6,6,1,6,1)
GN=-PO(1)/F
GK=PO(4)/F
IF(WZ)11,11,10
10 GN=GN+PO(3)/WZ
GK=GK-PO(6)/WZ
11 WRITE (8,3) IE,T(1),(PO(J),J=1,NSM),GN
3 FORMAT(I7,I14,F15.2,3F12.2)
WRITE (8,4) T(2),(PO(J),J=NK,NSEM),GK
4 FORMAT(I21,F15.2,3F12.2)
5 CONTINUE
RETURN
END
4.2.6. Ошибки округления вычислений
Как отмечалось, матрица коэффициентов системы уравнений ленточ-
ная, положительно определённая, симметричная. Диагональные коэффи-
циенты матрицы преобладают по значению над побочными, и поэтому
при решении системы уравнений нет необходимости в выборе главного
элемента, т. е. в процессе исключения при делении на элементы, располо-
женные на главной диагонали, элементы образуемой матрицы будут не
слишком велики, что уменьшает ошибки округления вычислений.

Страницы