|  | 
 
| 用的代码是课程书上找的 
   !     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
 
 
 
 | 
 |