پروژه فرترن (39526))=

---

شرح پروژه

#پروژه_جدید 

#کدپروژه : 39526

موضوع : پروژه فرترن =

 USE MSIMSL

REAL*8 K(6,6), M(6,6), E, A, L, Iz, ro

REAL*8, Allocatable :: Ksys(:,:), q(:,:), Msys(:,:), EVAL(:), EVEC(:,:)  

ne= 25

nsys= 3*ne

E= 199.9E+9

A= 4.265E-03

L= 6.

L= L/ne

Iz= 6.572E-05

ro= 7849.0474

K=0.

k(1,1)= E*A/L ; K(1,4)= -K(1,1)

K(2,2)= (12.*E*Iz)/(L3) ; K(2,3)= (6.*E*Iz)/(L2) ; K(2,5)= -K(2,2) ; K(2,6)= K(2,3)

K(3,3)= (4.*E*Iz)/L ; K(3,5)= -K(2,3) ; K(3,6)= K(3,3)/2.

K(4,4)= K(1,1)

K(5,5)= K(2,2) ; K(5,6)= -K(2,3)

K(6,6)= K(3,3)

DO i= 2,6

 DO j= 1,i-1

  K(i,j)= K(j,i)

 END DO

END DO

!PRINT* ; PRINT*,'K' ; PRINT "(6E9.2)",K

M=0.

M(1,1)= (ro*A*L)/3. ; M(1,4)= M(1,1)/2.

M(2,2)= ((ro*A*L)/420.)*156. ; M(2,3)= ((ro*A*L)/420.)*22.*L ; M(2,5)= ((ro*A*L)/420.)*54. ; M(2,6)= -((ro*A*L)/420.)*13.*L

M(3,3)= ((ro*A*L)/420.)*4.*L2 ; M(3,5)= -M(2,6); M(3,6)= -((ro*A*L)/420.)*3.*L2

M(4,4)= M(1,1)

M(5,5)= M(2,2) ; M(5,6)= -M(2,3)

M(6,6)= M(3,3)

DO i= 2,6

 DO j= 1,i-1

  M(i,j)= M(j,i)

 END DO

END DO

!PRINT* ; PRINT*,'M' ; PRINT "(6E9.2)",M

Allocate (Ksys(nsys,nsys),q(nsys,1))

Allocate (Msys(nsys,nsys))

Allocate (EVAL(nsys),EVEC(nsys,nsys))

Ksys=0.

Ksys(1:3,1:3)= K(4:6,4:6) ; iL= 1 ; iU=6

DO i= 2,ne

 Ksys(iL:iU,iL:iU)= Ksys(iL:iU,iL:iU)+K(1:6,1:6) ; iL= iL+3 ; iU= iU+3

END DO

!PRINT* ; PRINT*,'Ksys' ; CALL MPRINT(nsys,Ksys)

Msys=0.

Msys(1:3,1:3)= M(4:6,4:6) ; iL= 1 ; iU= 6

DO i= 2,ne

 Msys(iL:iU,iL:iU)= Msys(iL:iU,iL:iU)+M(1:6,1:6) ; iL= iL+3 ; iU= iU+3

END DO

!PRINT* ; PRINT*,'Msys' ; CALL MPRINT(nsys,Msys)

CALL DGVCSP(nsys,Ksys,nsys,Msys,nsys,EVAL,EVEC,nsys)

PRINT* ; PRINT*,'Frequencies' ; PRINT "(/,10F10.2,\)",(SQRT(EVAL(i)),i= nsys,1,-1)

! PRINT* ; PRINT "('Mode shapes',/)" ; CALL MPRINTi(nsys,EVEC)

END 

SUBROUTINE MPRINT(N,M)

REAL*8 M(n,n)

DO i= 1,n

 DO j= 1,n-1

  PRINT "(E9.2,\)",M(i,j) 

 END DO

 PRINT "(E9.2 )",M(i,n)

END DO

END

SUBROUTINE MPRINTi(N,M)

REAL*8 M(n,n)

DO i= 1,n

 DO j= n,2,-1

  PRINT "(E9.2,\)",M(i,j) 

 END DO

 PRINT "(E9.2 )",M(i,1)

END DO

END

: منظور از گام افزایشی ۵۰ تایی  در بند ۱، برای پارامتر ne میباشد

: لطفا هرجا مشخص نبود اعلام کنید من ویس توضیحات بدم که موضوع مشخص بشه

این پروژه شامل 3 فایل مهم است، لطفا قبل از ارسال پیشنهاد حتما نسبت به بررسی این فایل اقدام فرمایید.

مهارت ها و تخصص های مورد نیاز

---

---