abaqus帮助文档之地震相应计算分析(11)

0.91600E+01, 0.29060E-01, 0.91700E+01, 0.29690E-01, 0.91800E+01, 0.14230E-01, 0.91900E+01,-0.12200E-02, 0.92000E+01,-0.16680E-01, 0.92100E+01,-0.32130E-01, 0.92200E+01,-0.45700E-01, 0.92300E+01,-0.13870E-01, 0.92400E+01, 0.17960E-01, 0.92500E+01, 0.47150E-01, 0.92600E+01, 0.26170E-01, 0.92700E+01, 0.52000E-02, 0.92800E+01, 0.78900E-02, 0.92900E+01, 0.28450E-01, 0.93000E+01, 0.37120E-01, 0.93100E+01, 0.21660E-01, 0.93200E+01, 0.62100E-02, 0.93300E+01, 0.72900E-02, 0.93400E+01, 0.19390E-01, 0.93500E+01, 0.31480E-01, 0.93600E+01, 0.28400E-01, 0.93700E+01, 0.23050E-01, 0.93800E+01, 0.17710E-01, 0.93900E+01, 0.25890E-01, 0.94000E+01, 0.35920E-01, 0.94100E+01, 0.15840E-01, 0.94200E+01, 0.19940E-01, 0.94300E+01, 0.38850E-01, 0.94400E+01, 0.27140E-01, 0.94500E+01, 0.12400E-01, 0.94600E+01, 0.26080E-01, 0.94700E+01, 0.21520E-01, 0.94800E+01, 0.73800E-02, 0.94900E+01,-0.67600E-02, 0.95000E+01, 0.29100E-02, 0.95100E+01, 0.22790E-01, 0.95200E+01, 0.42660E-01, 0.95300E+01, 0.43300E-01, 0.95400E+01, 0.36470E-01, 0.95500E+01, 0.29640E-01, 0.95600E+01, 0.22810E-01, 0.95700E+01, 0.23260E-01, 0.95800E+01, 0.31310E-01, 0.95900E+01, 0.34750E-01, 0.96000E+01, 0.25110E-01, 0.96100E+01, 0.15460E-01, 0.96200E+01, 0.58100E-02, 0.96300E+01, 0.10960E-01, 0.96400E+01, 0.20280E-01, 0.96500E+01, 0.25400E-01, 0.96600E+01, 0.20260E-01, 0.96700E+01, 0.15110E-01, 0.96800E+01, 0.19570E-01, 0.96900E+01, 0.26280E-01, 0.97000E+01, 0.33000E-01, 0.97100E+01, 0.28510E-01, 0.97200E+01, 0.21720E-01, 0.97300E+01, 0.14940E-01, 0.97400E+01, 0.81500E-02, 0.97500E+01, 0.27500E-02, 0.97600E+01, 0.98100E-02, 0.97700E+01, 0.16870E-01, 0.97800E+01, 0.23930E-01, 0.97900E+01, 0.30990E-01, 0.98000E+01, 0.31740E-01, 0.98100E+01, 0.22200E-01, 0.98200E+01, 0.12670E-01, 0.98300E+01, 0.31300E-02, 0.98400E+01, 0.34500E-02, 0.98500E+01, 0.88400E-02, 0.98600E+01, 0.14230E-01, 0.98700E+01, 0.19620E-01, 0.98800E+01, 0.16040E-01, 0.98900E+01, 0.11790E-01, 0.99000E+01, 0.75300E-02, 0.99100E+01, 0.22980E-01, 0.99200E+01, 0.15160E-01, 0.99300E+01, 0.60700E-02, 0.99400E+01, 0.46700E-02, 0.99500E+01, 0.11950E-01, 0.99600E+01, 0.19220E-01, 0.99700E+01, 0.26500E-01, 0.99800E+01, 0.33780E-01, 0.99900E+01, 0.37510E-01, 0.10000E+02, 0.30640E-01

6、User subroutine UEL used by koyna_std.inp to model hydrodynamic interactions via the added mass technique. SUBROUTINE UEL(RHS,AMATRX,SVARS,ENERGY,

1 NDOFEL,NRHS,NSVARS,PROPS,NPROPS,COORDS, 2 MCRD,NNODE,U,DU,V,A,JTYPE,TIME,DTIME, 3 KSTEP,KINC,JELEM,PARAMS,NDLOAD,JDLTYP, 4 ADLMAG,PREDEF,NPREDF,LFLAGS,MLVARX, 5 DDLMAG,MDLOAD,PNEWDT,JPROPS,NJPROP, 6 PERIOD) C

INCLUDE 'ABA_PARAM.INC'

PARAMETER ( ZERO = 0.D0, HALF = 0.5D0, 1 ONE = 1.D0, SEVEN=7.0D0, EIGHT=8.0D0 ) C

DIMENSION RHS(MLVARX,*), 1 AMATRX(NDOFEL,NDOFEL),

2 SVARS(NSVARS),ENERGY(8),PROPS(*), 3 COORDS(MCRD,NNODE),U(NDOFEL), 4 DU(MLVARX,*),V(NDOFEL),A(NDOFEL), 5 TIME(2),PARAMS(3),JDLTYP(MDLOAD,*),

6 ADLMAG(MDLOAD,*),DDLMAG(MDLOAD,*), 7 PREDEF(2,NPREDF,NNODE),LFLAGS(*), 8 JPROPS(*) C

C USER ELEMENT TO MODEL THE HYDRODYNAMIC C PRESSURE OF RESERVOIR ON DAM USING C WESTERGAARD ADDED MASS TECHNIQUE. C

C PROPERTIES:

C PROPS(1) = Depth of reservoir C PROPS(2) = Y coordinate of water level C PROPS(3) = Density of water C

RESERVOIR_DEPTH = PROPS(1) YCOORD_WATER = PROPS(2) WATER_DENSITY = PROPS(3) YCOORD_ELCENTROID =

* HALF * (COORDS(2,2) + COORDS(2,1)) ELCENTROID_DEPTH =

* YCOORD_WATER - YCOORD_ELCENTROID RHO = (SEVEN/EIGHT) * WATER_DENSITY *

* SQRT(RESERVOIR_DEPTH * ELCENTROID_DEPTH) C

ALEN = ABS(COORDS(2,2)-COORDS(2,1)) AM = HALF*RHO*ALEN C

DO K1 = 1, NDOFEL DO KRHS = 1, NRHS RHS(K1,KRHS) = ZERO END DO

DO K2 = 1, NDOFEL AMATRX(K2,K1) = ZERO END DO END DO

do k1 = 1, nsvars svars(k1) = zero end do C

IF (LFLAGS(3).EQ.1) THEN C Normal incrementation IF (LFLAGS(1).EQ.11

* .OR.LFLAGS(1).EQ.12) THEN C *DYNAMIC