Fortran_平面刚架位移内力求解程序

!

!求解位移

!

program qiujieweiyi

implicit none

real(kind = 8) :: e,a,i,l,b

!e：弹性模量,a：单元面积,i：惯性矩值,l：单元长度,b：整体坐标系到局部坐标系的夹

角

logical :: alive

character(len = 20) :: title

real(kind = 8) :: count

integer :: z,j !循环参数

integer :: error

logical :: notdata = .true. !条件判断参数

integer :: m,n,q,w,g !循环参数

integer :: bec,bej !bec：杆元数目，bej：杆元非固支节点数目

real(kind = 8),allocatable :: zk(:,:),ck(:,:,:),t(:,:,:),f(:),wy(:),d(:,:),fc(:,:)

!zk(6,6,:)：整体刚度矩阵，ck(6,6,:)单元刚度矩阵，t(6,6)：坐标转换矩阵,f(:)：结点受

if(.not. alive)then

write(*,*)trim(filename1),"doesn't exist."

stop

end if

!读取文件，设置参量

open(unit = fileid1,file = filename1)

do m = 1,6

read(fileid1,"(A8,1X,F)",iostat = error)title,count

if(error /= 0)then

write(*,*)"文件出错 1"

exit

select case(title)

e = count

case("单元面积")

l = count

case("节点数目")

bej = count

case default

write(*,*)"数据已读完或不正确的数据"

stop

end select

end do

!设置变量、数组的值

allocate(zk(bej,bej))

allocate(ck(6,6,bec))

allocate(t(6,6,bec))

allocate(x(6,bec))

allocate(f(bej))

allocate(wy(bej))

zk = 0

ck = 0

t = 0

x = 0

d = 0

fc = 0

do q = 1,bec

if(error /= 0)then

write(*,*)"文件出错 2"

notdata = .false.

exit

end if

cosb = cosd(b)

s3 = 6*e*i*sinb*lt*lt

s4 = e*a*sinb*sinb*lt + 12*e*i*cosb*cosb*lt*lt*lt

s5 = 6*e*i*cosb*lt*lt

ck(3,3,q) = s1

ck(1,2,q) = s2

ck(2,1,q) = s2

ck(4,5,q) = s2

ck(5,4,q) = s2

ck(3,4,q) = s3

ck(4,3,q) = s3

ck(4,6,q) = s3

ck(6,4,q) = s3

ck(2,2,q) = s4

ck(5,5,q) = s4

ck(2,3,q) = s5

ck(3,2,q) = s5

ck(2,6,q) = s5

ck(3,3,q) = 2*s6

ck(6,6,q) = 2*s6

ck(1,4,q) = -s1

ck(4,1,q) = -s1

ck(2,4,q) = -s2

ck(4,2,q) = -s2

ck(1,3,q) = -s3

ck(3,1,q) = -s3