您所在位置：主题：求助：各位高手帮我看一下下面这程序，调试没有错误，运行不出结果，错误提示是在程序里面出现了不相匹配
jiangyjian110
[专家分：0] 发布于 2012-04-03 09:48:00
program test
implicit none
integer :: m,n,i,j,t,s
integer,parameter::length=1
integer,parameter::size=4
real ::ex=3.25E+10,ai=0.054543325
real ::a1,a2,a3
real ::k(size,size)
real ::u(4)
real :: MSE
real :: tra(4)
! 对单元刚度赋值
a1=ex*ai/length**3
a2=ex*ai/length**2
a3=ex*ai/length
k(1,1)=12*a1
k(2,1)=6*a2
k(2,2)=4*a3
k(3,1)=-12*a1
k(3,2)=-6*a2
k(3,3)=12*a1
k(4,1)=6*a2
k(4,2)=2*a3
k(4,3)=-6*a2
k(4,4)=4*a3
do i=1,size
do j=i,size
k(i,j)=k(j,i)
enddo
enddo
!读入模态振型数据
read(*,*) u(1),u(2),u(3),u(4)
!计算模态应变能
do m=1,4
do t=1,4
tra(:)=u(m)*k(t,m)
enddo
enddo
do s=1,4
MSE=tra(s)*u(s)/2
enddo
write(*,*)MSE
end
implicit none
integer :: m,n,i,j,t,s
integer,parameter::length=1
integer,parameter::size=4
real ::ex=3.25E+10,ai=0.054543325
real ::a1,a2,a3
real ::k(size,size)
real ::u(4)
real :: MSE
real :: tra(4)
! 对单元刚度赋值
a1=ex*ai/length**3
a2=ex*ai/length**2
a3=ex*ai/length
k(1,1)=12*a1
k(2,1)=6*a2
k(2,2)=4*a3
k(3,1)=-12*a1
k(3,2)=-6*a2
k(3,3)=12*a1
k(4,1)=6*a2
k(4,2)=2*a3
k(4,3)=-6*a2
k(4,4)=4*a3
do i=1,size
do j=i,size
k(i,j)=k(j,i)
enddo
enddo
!读入模态振型数据
read(*,*) u(1),u(2),u(3),u(4)
!计算模态应变能
do m=1,4
do t=1,4
tra(:)=u(m)*k(t,m)
enddo
enddo
do s=1,4
MSE=tra(s)*u(s)/2
enddo
write(*,*)MSE
end