- 积分
- 0
- 注册时间
- 2013-5-8
- 仿真币
-
- 最后登录
- 1970-1-1
|
程序如下, 问号的地方看不懂,恳请高手用浅显易懂的语言点拨。。
C 2-NODE BEAM C1 ELEMENT
C WITH GENERALIZED SECTION PROPERTIES AND PRESSURE
C
SUBROUTINE UEL(RHS,AMATRX,SVARS,ENERGY,NDOFEL,NRHS,NSVARS,
1 PROPS,NPROPS,COORDS,MCRD,NNODE,U,DU,V,A,JTYPE,TIME,DTIME,
2 KSTEP,KINC,JELEM,PARAMS,NDLOAD,JDLTYP,ADLMAG,PREDF,NPREDF,
3 LFLAGS,MLVARX,DDLMAG,MDLOAD,PNEWDT,JPROPS,NJPROP,PERIOD)
C
INCLUDE ’ABA_PARAM.INC’
C
DIMENSION RHS(MLVARX,*),AMATRX(NDOFEL,NDOFEL),SVARS(*),PROPS(*),
1 ENERGY(7),COORDS(MCRD,NNODE),U(NDOFEL),DU(MLVARX,*),V(NDOFEL),
2 A(NDOFEL),TIME(2),PARAMS(*),JDLTYP(MDLOAD,*),ADLMAG(MDLOAD,*),
3 DDLMAG(MDLOAD,*),PREDEF(2,NPREDF,NNODE),LFLAG(4),JPROPS(*)
C
C DEFINE PARAMETER AND VARIABLE
C
DIMENSION H(1,4)
PARAMETER(ZERO=0.D0,ONE=1.D0,TWO=2.D0,THREE=3.D0,FOUR=4.D0,
1 SIX=6.D0,EIGHT=8.D0,TWELVE=12.D0)
C
C CALCULATE LENGTH AND GET PROPERTIES
C
DL=COORDS(1,2)-COORDS(1,1)
DL2=DL**2
E=PROPS(1)
I=PROPS(2)
P=PROPS(3) !ASSIGIN PRESSURE
YMAX=PROPS(4)
C
C INITIALIZE RHS AND AMATRX
C
NSVINT=NSVARS/2
DO K1=1,4
RHS(K1,1)=ZERO
DO K2=1,4
AMATRX(K1,K2)=ZERO
END DO
END DO
DO K1=1,6
SVARS(K1)=ZERO
END DO
C
C AMATRX
Stress analysis in a beam using the UEL 10
C
AMATRX(1,1)=TWELVE
AMATRX(1,2)=SIX*DL
AMATRX(1,3)=-TWELVE
AMATRX(1,4)=SIX*DL
AMATRX(2,1)=SIX*DL
AMATRX(2,2)=FOUR*DL2
AMATRX(2,3)=-SIX*DL
AMATRX(2,4)=TWO*DL2
AMATRX(3,1)=-TWELVE
AMATRX(3,2)=-SIX*DL
AMATRX(3,3)=TWELVE
AMATRX(3,4)=-SIX*DL
AMATRX(4,1)=SIX*DL
AMATRX(4,2)=TWO*DL2
AMATRX(4,3)=-SIX*DL
AMATRX(4,4)=FOUR*DL2
DO K1=1,4
DO K2=1,4
AMATRX(K1,K2)=AMATRX(K1,K2)*E*I/(DL**3)
END DO
END DO
C
C RHS
C
DO K1=1,4
DO K2=1,4
RHS(K1,1)=RHS(K1,1)-AMATRX(K1,K2)*U(K2) ?U(K2)是啥意思? 整句话是啥意思?
END DO
END DO
RHS(1,1)=RHS(1,1)-(-P*DL/TWO) ?整句话是啥意思?
RHS(2,1)=RHS(2,1)-(-P*DL2/TWELVE)
RHS(3,1)=RHS(3,1)-(-P*DL/TWO)
RHS(4,1)=RHS(4,1)-(P*DL2/TWELVE)
C
C CALCULATE STRESS
C
CALL STRESS(SVARS,PROPS,DL,DL2,AMATRX,U)
C
C OUTPUT
C
IF(ABS(KINC) .LT. ABS(SVARS(5))) THEN
KINC=SVARS(5) ?SVARS(5) 是啥意思?整句话啥意思?
KSTEP=JELEM ?整句话啥意思?
END
IF
IF(ABS(KINC) .LT. ABS(SVARS(6))) THEN
KINC=SVARS(6)
KSTEP=JELEM+1 ?整句话啥意思?
Stress analysis in a beam using the UEL 11
END IF
OPEN(UNIT=16,FILE=’C:\TEMP\STRESS.DAT’,STATUS=’UNKNOWN’,
$ FORM=’FORMATTED’)
IF(JELEM .EQ. 1) THEN
WRITE(16,*) ’-------------------------NEXT STEP
$ -------------------------’
END IF
WRITE(16,*) ’ELEMENT NUMBER’,JELEM
WRITE(16,*) ’NODE’,JELEM,’FORCE’,SVARS(1)
WRITE(16,*) ’NODE’,JELEM,’MOMENT’,SVARS(2)
WRITE(16,*) ’NODE’,JELEM+1,’FORCE’,SVARS(3)
WRITE(16,*) ’NODE’,JELEM+1,’MOMENT’,SVARS(4)
WRITE(16,*) ’NODE’,JELEM,’NORMAL STRESS’,SVARS(5)
WRITE(16,*) ’NODE’,JELEM+1,’NORMAL STRESS’,SVARS(6)
WRITE(16,*) ’MAX NORMAL STRESS’,KINC,’AT NODE’,KSTEP
WRITE(16,*) ’ ’
C
C END
C
RETURN
END
|
|