| 查看: 551 | 回复: 3 | |||
[求助]
我的一个程序运行总是不成功，求大神帮忙 已有1人参与
|
| 如题，求助。附件里面是我的代码，求好心人看看帮小弟解答下不胜感激。。。 |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
求各位大神帮忙解释一下审稿人的评论，小硕不太理解，非常感谢~~~
已经有7人回复
- 急求！！！！一个Suzuki反应，求大神帮忙！！！！！
已经有14人回复
- 求大神帮忙解释截断能测试中的问题！
已经有7人回复
- 求大神帮忙解题啊！！！
已经有4人回复
- 求大神们帮忙看看这个程序哪儿需要改啊
已经有7人回复
- 求大神帮我解答，谢谢，必定重谢！！！
已经有4人回复
- 求大神帮忙拟合一个非线性方程，求出模型参数
已经有15人回复
- 新手求助~~求大神帮忙！！
已经有3人回复
- 求大神帮忙解答！
已经有8人回复
- fluent安装，不能用，求助大神帮忙~~~
已经有17人回复
- 在linux下运行ADF出错
已经有18人回复
- 求助大神帮我看看这段程序是干什么的？（就是定义这个函数算什么）
已经有6人回复
- 求大神帮忙“切尾”。。
已经有4人回复
|
! ! SUBROUTINE FOR EARTHQUAKE RESPONSE SPECTAR ! ! CODED BY Y.OHSAKI ! SUBROUTINE ERES(NH,H,ND1,NT,T,ND2,DT,NN,DDY,ND3,IND,QMAX,RES) ! DIMENSION H(ND1),T(ND2),DDY(ND3),RES(ND2,ND1) DIMENSION EMAX(3),RMAX(3) PARAMETER (P2=6.283185) ! ! MAXIMA OF INPUT MOTION ! EMAX(1)=ABS(DDY(1)) EMAX(2)=0. EMAX(3)=0. DDYF=DDY(1) DYF=0. YF=0. DO 110 M=2,NN DDYM=DDY(M) DY=DYF+(DDYF+DDYM)*DT/2 Y=YF+DYF*DT+(DDYF/3.+DDYM/6.)*DT**2 EMAX(1)=AMAX1(EMAX(1),ABS(DDYM)) EMAX(2)=AMAX1(EMAX(2),ABS(DY)) EMAX(3)=AMAX1(EMAX(3),ABS(Y)) DDYF=DDYM DYF=DY YF=Y 110 CONTINUE QMAX=EMAX(IND) ! ! RESPONSE COMPUTATION ! DO 150 L=1,NH DO 140 K=1,NT IF (T(K).EQ.0.)GO TO 130 W=P2/T(K) W2=W*W HW=H(L)*W WD=W*SQRT(1.-H(L)**2) WDT=WD*DT E=EXP(-HW*DT) CWDT=COS(WDT) SWDT=SIN(WDT) A11=E*(CWDT+HW*SWDT/WD) A12=E*SWDT/WD A21=-E*W2*SWDT/WD A22=E*(CWDT-HW*SWDT/WD) SS=-HW*SWDT-WD*CWDT CC=-HW*CWDT+WD*SWDT S1=(E*SS+WD)/W2 C1=(E*CC+HW)/W2 S2=(E*DT*SS+HW*S1+WD*C1)/W2 C2=(E*DT*CC+HW*C1-WD*S1)/W2 S3=DT*S1-S2 C3=DT*C1-C2 B11=-S2/WDT B12=-S3/WDT B21=(HW*S2-WD*C2)/WDT B22=(HW*S3-WD*C3)/WDT RMAX(1)=2.*HW*ABS(DDY(1))*DT RMAX(2)=ABS(DDY(1))*DT RMAX(3)=0. DXF=-DDY(1)*DT XF=0. DO 120 M=2,NN DDYM=DDY(M) DDYF=DDY(M-1) X=A12*DXF+A11*XF+B12*DDYM+B11*DDYF DX=A22*DXF+A21*XF+B22*DDYM+B21*DDYF DDX=-2.*HW*DX-W2*X RMAX(1)=AMAX1(RMAX(1),ABS(DDX)) RMAX(2)=AMAX1(RMAX(2),ABS(DX)) RMAX(3)=AMAX1(RMAX(3),ABS(X)) DXF=DX XF=X 120 CONTINUE RES(K,L)=RMAX(IND) GO TO 140 130 RES(K,L)=0. IF(IND.EQ.1) RES(K,L)=EMAX(1) 140 CONTINUE 150 CONTINUE RETURN END DIMENSION H(3),T(35),DDY(4),RES(35,3) DATA NH/3/,H/0.,0.05,0.1/,IND/1/,DT/0.02/ DATA NT/35/,T/0.00,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,& 0.45,0.50,0.55,0.60,0.65,0.70,0.75,0.80,0.85,& 0.90,0.95,1.00,1.20,1.40,1.60,1.80,2.00,2.20,& 2.40,2.60,2.80,3.00,3.50,4.00,4.50,5.00/ open(5,file='data.dat',status='old') REWIND(5) ! READ(5,501)(DDY(M),M=2,5) CALL ERES(NH,H,3,NT,T,35,DT,NN,DDY,4,IND,QMAX,RES) STOP ! 501 FORMAT(F10.3) END 不好意思，程序没传上。上面的就是我的代码 |
2楼2014-10-10 18:33:34
weiwei43
至尊木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 18316.1
- 红花: 6
- 帖子: 737
- 在线: 94.7小时
- 虫号: 250282
- 注册: 2006-05-13
- 专业: 机械动力学
3楼2014-10-17 08:28:57
weiwei43
至尊木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 18316.1
- 红花: 6
- 帖子: 737
- 在线: 94.7小时
- 虫号: 250282
- 注册: 2006-05-13
- 专业: 机械动力学
4楼2014-10-17 08:34:42
