查看: 650  |  回复: 13
【悬赏金币】回答本帖问题,作者笨人煮酒507将赠送您 10 个金币

笨人煮酒507

新虫 (初入文坛)

[求助] 黄平《弹性流体动压润滑数值计算方法》中线接触等温弹流润滑fortran程序运行的问题

黄平《弹性流体动压润滑数值计算方法》中线接触等温弹流润滑fortran程序运行的问题,求大神指导,非常感谢!这个程序运行不出来,谁知道哪儿有错吗?
PROGRAM GREASELINE
        CHARACTER*1 S,S1,S2
        CHARACTER*16 CDATE,CTIME
        COMMON /COM1/ENDA,A1,A2,A3,Z,C1,C3,CW,LMAX,FF/COM2/EDA0/COM4/X0,XE/COM3/E1,PH,B,U1,U2,R
        DATA PAI,Z,P0/3.14159265,0.68,1.96E8/S1,S2/1HY,1Hy/!预定义数值
        DATA N,X0,XE,W,E1,EDA0,R,Us,CU,C1,FN/129,-4.,1.4,1.E5,2.26E11,0.4,0.0128,0.87,0.67,0.5,0.8/CDATE,CTIME/'The date is','The time is'/!预赋值参数
    OPEN(8,FILE='OUT.DAT',STATUS='UNKNOWN')!建立输出结果文件
1   FORMAT(20X,A12,I2.2,':',I2.2,':',I4.4)
2   FORMAT(20X,A12,I2.2,':',I2.2,':',I2.2,'.',I2.2)
        WRITE(*,*)'参数已经赋值? (Y or N)?'!询问是否在程序中更改参数
        READ(*,'(A)')S
        IF(S.EQ.S1.OR.S.EQ.S2)THEN
        GOTO 10
        ENDIF
        WRITE(*,*)'PH=' !程序运行中更改数值,PH为Hertz接触压力
        READ(*,*)PH
        W=2.*PAI*R*PH*(PH/E1)
        WRITE(*,*)'W=',W
10  CW=N+0.1
    FF=1./FN
        LMAX=ALOG(CW)/ALOG(2.)!将参数量纲一化
        N=2**LMAX+1
        LMIN=(ALOG(CW)-ALOG(SQRT(CW)))/ALOG(2.)
        LMAX=LMIN
        W1=W/(E1*R)
        PH=E1*SQRT(0.5*W1/PAI)
        A1=(ALOG(EDA0)+9.67)
        A2=PH/P0
        A3=0.59/(PH*1.E-9)
        B=4.*R*PH/E1
        ALFA=Z*A1/P0
        G=ALFA*E1
        U=EDA0*US/(2.*E1*R)
        C3=1.6*(R/B)**2*G**0.6*U**0.7*W1**(-0.13)
        ENDA=B**(2.+FF)*(PH/2/EDA0)**FF/R**(1+FF)/US/(2.+FF)
        U1=0.5*(2.+CU)*U
        U2=0.5*(2.-CU)*U
    WRITE(*,*)'B,PH,G,U=',B,PH,G,U
        CW=-1.13*C3
        WRITE(*,*)N,X0,XE,W,E1,EDA0,R,US,PH
        WRITE(8,*)N,W,E1,EDA0,R,US,B,PH,FF
        WRITE(*,40)
40  FORMAT(2X,'Wait  Please',//)
        CALL SUBAK(N)!计算弹性变形系数
        CALL MULTI(N)!计算压力P和H
        STOP
        END
        SUBROUTINE MULTI(N)!计算压力P和H
        REAL*8 X(1100),P(1100),H(1100),RO(1100),POLD(1100),EPS(1100),EDA(1100),R(1100),K(1100),E(1100)
        COMMON /COM1/ENDA,A1,A2,A3,Z,C1,C3,CW,LMAX,FF/COM4/X0,XE/COM3/E1,PH,B,U1,U2,RR
        DATA MK,G0/1,1.570796325/
        NX=N
        DX=(XE-X0)/(N-1.0)
        DO 10 I=1,N
        X(I)=X0+(I-1)*DX
        IF(ABS(X(I)).GE.1.0)P(I)=0.0
        IF(ABS(X(I)).LT.1.0)P(I)=SQRT(1.-X(I)*X(I))
10  CONTINUE
        CALL HREE(N,DX,H00,G0,X,P,H,RO,EPS,EDA)!实现各点的膜厚、密度、粘度、弹性变形的计算
        CALL FZ(N,P,POLD)
14        KK=19
        CALL ITER(N,KK,DX,H00,G0,X,P,H,RO,EPS,EDA,R)!利用雷诺方程进行各点压力的重新计算,并再次调用HREE子程序
        MK=MK+1
        CALL ERROP(N,P,POLD,ERP)!计算迭代前后的压力差
        IF(ERP.GT.1.E-4.AND.MK.LE.800)THEN
        GOTO 14
        ENDIF
        WRITE(*,*)PH,RR,B
105        IF(MK.GE.800)THEN
        WRITE(*,*)'Pressures are not convergent !!!'
        READ(*,*)
        ENDIF
        FM=FRICT(N,DX,X,H,P,EDA)
        H2=1.E3
        P2=0.0
        DO 106 I=1,N
        IF(H(I).LT.H2)H2=H(I)
        IF(P(I).GT.P2)P2=P(I)
106 CONTINUE
    DO 108 I=1,N
    K(I)=P(I)*PH/1.E9
        E(I)=H(I)*B*B*1.E6/RR
108 CONTINUE
        H3=H2*B*B/RR
        P3=P2*PH
110        FORMAT(6(1X,E12.6))
120        CONTINUE
        WRITE(8,*)'P2,H2,P3,H3=',P2,H2,P3,H3
        CALL OUTHP(N,X, K,E)!实现结果的输出功能
        RETURN
        END
        SUBROUTINE OUTHP(N,X, K,E)!实现结果的输出功能
        REAL*8  X(N), K(N),E(N)
        DX=X(2)-X(1)
        DO 10 I=1,N
        WRITE(8,20)X(I),K(I),E(I)
10  CONTINUE
20  FORMAT(1X,6(F20.6,1X))
        RETURN
        END
        SUBROUTINE HREE(N,DX,H00,G0,X,P,H,RO,EPS,EDA)!实现各点的膜厚、密度、粘度、弹性变形的计算
    REAL*8  X(N),P(N),H(N),RO(N),EPS(N),EDA(2200)
        REAL*8  W(2200)
        COMMON /COM1/ENDA,A1,A2,A3,Z,C1,C3,CW,K,FF/COM2/EDA0/COMAK/AK(0:1100)
    DATA KK,NW,PAI1/0,2200,0.318309886/
        IF(KK.NE.0)GOTO 3
        HM0=C3
        H00=0.0
3        W1=0.0
        DO 4 I=1,N
4   W1=W1+P(I)!此压力下的量纲一化载荷W
        C3=(DX*W1)/G0
        DW=1.-C3 !承载力判断值
        CALL DISP(N,NW,K,DX,P,W)!计算各点弹性变形量
        HMIN=1.E3
        DO 30 I=1,N
        H0=0.5*X(I)*X(I)-PAI1*W(I)!
        IF(H0.LT.HMIN)HMIN=H0
        H(I)=H0
30        CONTINUE
    IF(KK.NE.0)GOTO 32
        KK=1
        H00=-HMIN+HM0
32        H0=H00+HMIN
        IF(H0.LE.0.0)GOTO 48
        IF(H0+0.3*CW*DW.GT.0.0)HM0=H0+0.3*CW*DW
        IF(H0+0.3*CW*DW.LE.0.0)HM0=HM0*C3
48        H00=HM0-HMIN
50        DO 60 I=1,N
60        H(I)=H00+H(I)!膜厚计算
        DO 100 I=1,N
        EDA(I)=EXP(A1*(-1.+(1.+A2*P(I))**Z))!各点粘度数值计算
        RO(I)=1.0
        EPS(I)=RO(I)*H(I)**(2+FF)*ENDA/EDA(I)**FF !离散化Reynolds方程中的系数
100        CONTINUE
        RETURN
        END
        SUBROUTINE ITER(N,KK,DX,H00,G0,X,P,H,RO,EPS,EDA,R)!利用雷诺方程进行各点压力的重新计算,并再次调用HREE子程序
        REAL*8 X(N),P(N),H(N),RO(N),EPS(N),EDA(N),R(N)
        COMMON /COM1/ENDA,A1,A2,A3,Z,C1,C3,CW,LMAX,FF/COMAK/AK(0:1100)
        DATA KG1,PAI/0,3.14159265/
        IF(KG1.NE.0)GOTO 5
        KG1=1
        DX1=1./DX
        DX2=DX*DX
        DX3=1./DX2
        DX4=DX1/PAI
        DX5=DX1**(1+FF)
        DXL=DX*ALOG(DX)
        AK0=DX*AK(0)+DXL
        AK1=DX*AK(1)+DXL
5        DO 100 K=1,KK
        D2=0.5*(EPS(1)+EPS(2))
        D3=0.5*(EPS(2)+EPS(3))
        D5=DX1*(RO(2)*H(2)-RO(1)*H(1))
        D7=DX4*(RO(2)*AK0-RO(1)*AK1)
        PP=0.
        DO 70 I=2,N-1
        D1=D2
        D2=D3
        D4=D5
        D6=D7
        IF(I+2.LE.N)D3=0.5*(EPS(I+1)+EPS(I+2))
        D5=DX1*(RO(I+1)*H(I+1)-RO(I)*H(I))
        D7=DX4*(RO(I+1)*AK0-RO(I)*AK1)
        DD=(D1+D2)*DX3
        IF(0.05*DD.LT.ABS(D6))GOTO 10
        RI=-DX5*(D2*SIGN(1.0,(P(I+1)-P(I)))*ABS((P(I+1)-P(I)))**(FF)-D1*SIGN(1.0,(P(I)-P(I-1)))*ABS((P(I)-P(I-1)))**(FF))+D4 !雷诺方程离散形式
        R(I)=RI
        DLDP=-FF*DX5*(D1*ABS((P(I)-P(I-1)))**(FF-1)+D2*ABS((P(I+1)-P(I)))**(FF-1))+D6
        RI=RI/DLDP
        RI=RI/C1
        GOTO 20
10        RI=-DX5*(D2*SIGN(1.0,(P(I+1)-P(I)))*ABS((P(I+1)-P(I)))**(FF)-D1*SIGN(1.0,(P(I)-PP))*ABS((P(I)-PP))**(FF))+D4
        R(I)=RI
        DLDP=-FF*DX5*(2*D1*ABS((P(I)-PP))**(FF-1)+D2*ABS((P(I+1)-P(I)))**(FF-1))+2.*D6
        RI=RI/DLDP

        IF(I.GT.2.AND.P(I-1)-C1*RI.GT.0)P(I-1)=P(I-1)-C1*RI
20  PP=P(I)
        P(I)=P(I)+C1*RI
        IF(P(I).LT.0.0)P(I)=0.0
        IF(P(I).LE.0.0)R(I)=0.0
70        CONTINUE
        CALL HREE(N,DX,H00,G0,X,P,H,RO,EPS,EDA)!再次调用HREE子程序实现各点的膜厚、密度、粘度、弹性变形的计算
100 CONTINUE
        RETURN
        END
        SUBROUTINE DISP(N,NW,KMAX,DX,P1,W) !计算弹性变形
        REAL*8 P1(N),W(NW),P(2200),AK1(0:50),AK2(0:50)
        COMMON /COMAK/AK(0:1100)
        DATA NMAX,KMIN/2200,1/
        N2=N
        M=3+2*ALOG(FLOAT(N))!求和次数 M>3+2ln(n)
        K1=N+KMAX
        DO 10 I=1,N
10  P(K1+I)=P1(I)
        DO 20 KK=KMIN,KMAX-1
        K=KMAX+KMIN-KK !积分系数
        N1=(N2+1)/2
        CALL DOWNP(NMAX,N1,N2,K,P)!下传节点参数到最粗网格
20        N2=N1
        DX1=DX*2**(KMAX-KMIN)
        CALL WI(NMAX,N1,KMIN,KMAX,DX,DX1,P,W) !在粗网格上进行积分
        DO 30 K=KMIN+1,KMAX
        N2=2*N1-1
        DX1=DX1/2.
        CALL AKCO(M+5,KMAX,K,AK1)!将粗网格上积分系数传到细网格节点上
        CALL AKIN(M+6,AK1,AK2)!插值计算上层网格的积分系数
        CALL WCOS(NMAX,N1,N2,K,W)!将在粗网格上得到的积分数值映射到相对应的细网格上
        CALL CORR(NMAX,N2,K,M,1,DX1,P,W,AK1)!利用积分系数插值对映射节点积分值进行修正
        CALL WINT(NMAX,N2,K,W)!插值计算不与粗网格重合的细网格节点的积分数值
        CALL CORR(NMAX,N2,K,M,2,DX1,P,W,AK2)!利用积分系数插值对插值节点积分值进行修正
30  N1=N2
        DO 40 I=1,N
40  W(I)=W(K1+I)
        RETURN
        END
        SUBROUTINE DOWNP(NMAX,N1,N2,K,P)!下传节点参数到最粗网格
        REAL*8 P(NMAX)
        K1=N1+K-1
        K2=N2+K-1
        DO 10 I=3,N1-2
        I2=2*I+K2
10  P(K1+I)=(16.*P(I2)+9.*(P(I2-1)+P(I2+1))-(P(I2-3)+P(I2+3)))/32. !细网格节点压力插值后传到粗网格节点
        P(K1+2)=0.25*(P(K2+3)+P(K2+5))+0.5*P(K2+4)
        P(K1+N1-1)=0.25*(P(K2+N2-2)+P(K2+N2))+0.5*P(K2+N2-1)
        RETURN
        END
        SUBROUTINE WCOS(NMAX,N1,N2,K,W)!将在粗网格上得到的积分数值映射到相对应的细网格上
        REAL*8 W(NMAX)
        K1=N1+K-1
        K2=N2+K
        DO 10 I=1,N1
        II=2*I-1
10  W(K2+II)=W(K1+I)!部分细网格节点上的积分数值未计算,根据已知的粗网格上的数值进行计算
        RETURN
        END
        SUBROUTINE WINT(NMAX,N,K,W)!插值计算不与粗网格重合的细网格节点的积分数值
        REAL*8 W(NMAX)
        K2=N+K
        DO 10 I=4,N-3,2
        II=K2+I
10  W(II)=(9.*(W(II-1)+W(II+1))-(W(II-3)+W(II+3)))/16. !部分细网格节点上的积分数值未计算,根据已知的粗网格上的数值进行计算
        I1=K2+2
        I2=K2+N-1
        W(I1)=0.5*(W(I1-1)+W(I1+1))
        W(I2)=0.5*(W(I2-1)+W(I2+1))
        RETURN
        END
        SUBROUTINE CORR(NMAX,N,K,M,I1,DX,P,W,AK)!利用积分系数插值对映射节点积分值进行修正
        REAL*8 P(NMAX),W(NMAX),AK(0:M)
        K1=N+K
        IF(I1.EQ.2)GOTO 20
        DO 10 I=1,N,2
        II=K1+I
        J1=MAX0(1,I-M)
        J2=MIN0(N,I+M)
        DO 10 J=J1,J2
        IJ=IABS(I-J)
10  W(II)=W(II)+AK(IJ)*DX*P(K1+J)
        RETURN
20  DO 30 I=2,N,2
        II=K1+I
        J1=MAX0(1,I-M)
        J2=MIN0(N,I+M)
        DO 30 J=J1,J2
        IJ=IABS(I-J)
30  W(II)=W(II)+AK(IJ)*DX*P(K1+J)
        RETURN
        END
        SUBROUTINE WI(NMAX,N,KMIN,KMAX,DX,DX1,P,W)!求最粗网格上的数值积分
    REAL*8 P(NMAX),W(NMAX)
        COMMON /COMAK/AK(0:1100)
        K1=N+1
        K=2**(KMAX-KMIN)
        C=ALOG(DX)
        DO 10 I=1,N
        II=K1+I
        W(II)=0.0
        DO 10 J=1,N
        IJ=K*IABS(I-J)
10  W(II)=W(II)+(AK(IJ)+C)*DX1*P(K1+J)
        RETURN
        END
        SUBROUTINE AKCO(KA,KMAX,K,AK1)!将粗网格上积分系数传到细网格节点上
        REAL*8 AK1(0:KA)
        COMMON /COMAK/AK(0:1100)
        J=2**(KMAX-K)
        DO 10 I=0,KA
        II=J*I
10  AK1(I)=AK(II)
        RETURN
        END
        SUBROUTINE AKIN(KA,AK1,AK2)!插值计算上层网格的积分系数
        REAL*8 AK1(KA),AK2(KA)
        DO 10 I=4,KA-3
10  AK2(I)=(9.*(AK1(I-1)+AK1(I+1))-(AK1(I-3)+AK1(I+3)))/16.
        AK2(1)=(9.*AK1(2)-AK1(4))/8.
        AK2(2)=(9.*(AK1(1)+AK1(3))-(AK1(3)+AK1(5)))/16.
        AK2(3)=(9.*(AK1(2)+AK1(4))-(AK1(2)+AK1(6)))/16.
        DO 20 I=1,KA
20  AK2(I)=AK1(I)-AK2(I)
        DO 30 I=1,KA-1,2
        I1=I+1
        AK1(I)=0.0
30  AK1(I1)=AK2(I1)
        RETURN
        END
        SUBROUTINE SUBAK(MM)!计算弹性变形系数
        COMMON /COMAK/AK(0:1100)! 程序所得弹性变形系数AK(I)放在公共区供其它程序使用
        DO 10 I=0,MM
10  AK(I)=(I+0.5)*(ALOG(ABS(I+0.5))-1.)-(I-0.5)*(ALOG(ABS(I-0.5))-1.)!计算公式
        RETURN
        END
        FUNCTION FRICT(N,DX,X,H,P,EDA)
        REAL*8 X(N),H(N),P(N),EDA(N)
        COMMON /COM3/E1,PH,B,U1,U2,R
        DATA TAU0,AT/1.98E7,0.078/
        OPEN (10,FILE='TAU.DAT')
        TP=TAU0/PH
        TE=TAU0/E1
        A=AT/TAU0
        BR=B/R
        FRICT=0.0
        DO 10 I=1,N
        DP=0.0
        A=1.0+AT*P(I)/TAU0
        IF(I.NE.N)DP=(P(I+1)-P(I))/DX
        TAU1=0.5*H(I)*DP*(BR/TP/A)
        TAU2=2.*U1*EDA(I)/(H(I)*BR**2*TE*A)
    FRICT=FRICT+ABS(TAU1)+TAU2
        TAU2=0.05*TAU2
        WRITE(10,5)X(I),TAU1,TAU2,P(I),H(I)
5        FORMAT(5(E12.6,1X))
10  CONTINUE
        FRICT=FRICT*DX*B*TAU0
        RETURN
        END
        SUBROUTINE FZ(N,P,POLD)
        REAL*8 P(N),POLD(N)
        DO 10 I=1,N
10  POLD(I)=P(I)
        RETURN
        END
        SUBROUTINE ERROP(N,P,POLD,ERP)!计算迭代前后的压力差
        REAL*8 P(N),POLD(N)
        SD=0.0
        SUM=0.0
        DO 10 I=1,N
        SD=SD+ABS(P(I)-POLD(I))
        POLD(I)=P(I)
10  SUM=SUM+P(I)
        ERP=SD/SUM
        RETURN
        END

黄平《弹性流体动压润滑数值计算方法》中线接触等温弹流润滑fortran程序运行的问题
QQ图片20160421091007.png
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

笨人煮酒507

新虫 (初入文坛)

没人知道么
2楼2016-04-22 19:09:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
3楼2016-05-08 08:28:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

挖坑埋仪器

新虫 (正式写手)

看编译器和黑窗报错信息啊

发自小木虫Android客户端
4楼2016-05-11 02:50:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

指环王jyj

新虫 (小有名气)

5楼2016-05-11 19:59:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

扬州小天王

新虫 (小有名气)

兄弟,你的问题解决了吗?我也遇到同样的问题了

发自小木虫IOS客户端
6楼2017-03-21 10:15:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

古月山

新虫 (初入文坛)

你的是黄平书上的源程序么?

发自小木虫IOS客户端
7楼2017-04-05 14:11:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

古月山

新虫 (初入文坛)

他的书上主要介绍轴承油膜编程,请问您知道如何对圆盘,圆环进行编程的。

发自小木虫IOS客户端
8楼2017-04-05 14:11:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

竹小min

新虫 (初入文坛)

9楼2017-04-17 11:05:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我的好吗

新虫 (初入文坛)

请问,哪里可以找到这本书的PDF?

发自小木虫Android客户端
10楼2017-12-25 16:53:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 笨人煮酒507 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 投票:  二本博士毕业去向 +6 sunkaier 2018-10-16 7/350 2018-10-19 01:38 by 胡胡的虫
[找工作] 面板行业or电池行业 +13 三堆 2018-10-12 27/1350 2018-10-18 21:42 by 三堆
[职场人生] 袭胸 +3 tang0815 2018-10-18 3/150 2018-10-18 19:39 by caiqi
[硕博家园] 不开心,剑桥的教授回复我!牛津发了一个邮件也没有人回复 +10 liwei_chen 2018-10-18 10/500 2018-10-18 18:28 by longwave
[论文投稿] 一审修回后,状态变为邀请审稿人,然后过去很长一直没变,这是什么情况? 20+3 chgylw 2018-10-12 7/350 2018-10-18 14:58 by ZBBcrazy
[基金申请] 还没毕业和入职之前能申请省基和国基吗? +4 ghvc 2018-10-17 11/550 2018-10-18 14:45 by wudi0729
[硕博家园] 搞不懂 +8 chichi精灵 2018-10-17 12/600 2018-10-18 08:44 by shaoju852
[硕博家园] 后悔!!! +33 刘兜兜。 2018-10-14 34/1700 2018-10-18 01:05 by xxh170925
[外语学习] 雅思资料 +6 Sunshine0327 2018-10-17 11/550 2018-10-18 00:34 by LUJINFENG
[基金申请] 面上项目2年未中,停了;可以申请河南省联合基金吗? +8 fennhuie 2018-10-14 9/450 2018-10-17 22:15 by fennhuie
[硕博家园] 换专业读博可行性多大,有成功的前辈 +13 麦穗睡石头 2018-10-14 18/900 2018-10-17 16:57 by lvchaojia
[职场人生] 贝特瑞和欣旺达选哪个 +8 ayu121 2018-10-14 12/600 2018-10-16 19:07 by xiaomingmin
[分析] 液相测五氟溴苯和十二烷基苯为什么不出峰? +3 chenmin300 2018-10-12 4/200 2018-10-16 08:11 by 繁夏
[硕博家园] 2018.10.13 +5 pengfei_213 2018-10-13 5/250 2018-10-15 09:08 by doctorlan
[教师之家] 重庆大学青年讲师 待遇 +19 cqdx111 2018-10-12 26/1300 2018-10-14 22:53 by wenjie8408
[硕博家园] 求助~与申请港大或南科大与港大联培博士相关 +3 张jiayu 2018-10-14 3/150 2018-10-14 14:01 by 甄士隐207
[硕博家园] 嗯嗯 +25 yeyuhuayue 2018-10-12 51/2550 2018-10-14 08:56 by quan3145
[硕博家园] 读博求助 +4 koh_rong 2018-10-13 4/200 2018-10-13 17:53 by 玉麒麟007
[考博] 【博导见面】求助!明天和导师见面! +3 风追的少年_ 2018-10-12 6/300 2018-10-12 23:48 by Nstudent
[考研] 制药工程专业考研 +4 瓶卍邪 2018-10-12 4/200 2018-10-12 14:01 by 易中华
信息提示
请填处理意见