|
用的代码是课程书上找的
! 1、本程序采用3结点三角形单元计算弹性体自振特性问题;
! 2、采用子空间迭代法求解;
! 3、主要输出内容:自振频率和振型。
PROGRAM FEM_TRI3_EIGEN
USE LIBRARY
USE GLOBAL_VARIABLES
USE EIGENMOD
IMPLICIT NONE
! AR:输入文件名(最长15个字符);OK:输出文件名(最长15个字符)
CHARACTER(LEN=15):: AR,OK
! X、Y:结点坐标 AE:材料参数 SK:整体劲度矩阵 DM:集中质量数组 RR:各阶振型数组 RLD:各阶频率数组
REAL(8),ALLOCATABLE:: X(:),Y(:),AE(:,:),SK(:),DM(:),RR(:,:),RLD(:)
! MEO:单元结点信息 JR:结点自由度序号 MA:指示矩阵
INTEGER,ALLOCATABLE::MEO(:,:),JR(:,:),MA(:)
INTEGER:: NS,I
! 屏幕提示:PLEASE INPUT FILE NAME OF DATA 读取输入文件
WRITE(*,"(///A)") 'PLEASE INPUT FILE NAME OF DATA='
READ(*,"(A15)") AR
! 屏幕提示:PLEASE OUTPUT FILE NAME OF DATA 读取输出文件
WRITE(*,"(///A)") 'PLEASE OUTPUT FILE NAME OF DATA='
READ(*,"(A15)") OK
OPEN (5,FILE=AR,STATUS='OLD')
OPEN (7,FILE=OK,STATUS='UNKNOWN') 在这个地方显示断点,显示没有找到文件 file not found unit 5
|
|