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

BETA = PARAMS(2)

DADU = ONE/(BETA*DTIME**2) C

AMATRX(1,1) = AM*DADU AMATRX(3,3) = AM*DADU RHS(1,1) = RHS(1,1)-AM*A(1) RHS(3,1) = RHS(3,1)-AM*A(3)

ENERGY(1)= HALF*AM*(V(1)*V(1)+V(3)*V(3)) END IF

ELSE IF (LFLAGS(3).EQ.4) THEN C Mass matrix

AMATRX(1,1) = AM AMATRX(3,3) = AM

ELSE IF (LFLAGS(3).EQ.5) THEN C Half-step residual calculation RHS(1,1) = RHS(1,1)-AM*A(1) RHS(3,1) = RHS(3,1)-AM*A(3) ELSE IF (LFLAGS(3).EQ.6) THEN C Initial acceleration calculation AMATRX(1,1) = AM AMATRX(3,3) = AM

ENERGY(1) = HALF*AM*(V(1)*V(1)+V(3)*V(3)) END IF C

RETURN END

7、Analysis of the pre-seismic state of the Koyna dam. These results are imported by koyna_xpl.inp. *HEADING

KOYNA DAM: PRE-SEISMIC ANALYSIS Units - n, m, sec

***************************************** ** Step 1: Gravity load

** Step 2: Hydrostatic pressure load

***************************************** *PREPRINT, MODEL=YES *NODE

1 , 0.00, 0.00 21, 70.00, 0.00 2401, 0.00, 66.50 2421, 19.25, 66.50 3601, 0.00, 91.75

3621, 16.17, 91.75 3801, 0.00, 103.00 3821, 14.80, 103.00 *NGEN, NSET=NBASE 1, 21

*NGEN, NSET=NMID 2401, 2421

*NGEN, NSET=NWL 3601, 3621

*NGEN, NSET=NTOP 3801, 3821

*NFILL, BIAS=1.05

NBASE, NMID, 24, 100 *NFILL, BIAS=0.92 NMID, NWL, 12, 100 *NFILL

NWL, NTOP, 2, 100 *NSET, NSET=NOUT 1, 3801

*ELSET, ELSET=EOUT 2320, 2420

*ELEMENT, TYPE=CPS4R 1, 1, 2, 102, 101

*ELGEN, ELSET=DAM 1, 20, 1, 1, 38, 100, 100

*ELSET, ELSET=WDAM, GENERATE 1, 3501, 100

*SOLID SECTION, ELSET=DAM, MATERIAL=CONCRETE,CONTROLS=HGC 1.0,

*SECTION CONTROLS,NAME=HGC,HOURGLASS=ENHANCED *MATERIAL, NAME=CONCRETE *ELASTIC

3.1027E+10, 0.2 *DENSITY 2643.0

*DAMPING, BETA=0.00323

*CONCRETE DAMAGED PLASTICITY 36.31

*CONCRETE COMPRESSION HARDENING 13.0E+6, 0.0 24.1E+6, 0.001

*CONCRETE TENSION STIFFENING, TYPE=DISPLACEMENT 2.9E+6 ,0

1.94393E+6 ,0.000066185

1.30305E+6 ,0.00012286 0.873463E+6 ,0.000173427 0.5855E+6 ,0.00022019 0.392472E+6 ,0.000264718 0.263082E+6 ,0.000308088 0.176349E+6 ,0.00035105 0.11821E+6 ,0.000394138 0.0792388E+6 ,0.000437744 0.0531154E+6 ,0.000482165

*CONCRETE TENSION DAMAGE, TYPE=DISPLACEMENT 0 ,0

0.381217 ,0.000066185 0.617107 ,0.00012286 0.763072 ,0.000173427 0.853393 ,0.00022019 0.909282 ,0.000264718 0.943865 ,0.000308088 0.965265 ,0.00035105 0.978506 ,0.000394138 0.9867 ,0.000437744 0.99177 ,0.000482165 *BOUNDARY NBASE, 1,2

*STEP, NLGEOM, UNSYMM=YES STEP 1 - GRAVITY LOAD *STATIC

1.0E-10, 1.0E-10 *DLOAD

DAM, GRAV, 9.81, 0, -1

*OUTPUT, FIELD, VARIABLE=PRESELECT, FREQ=10 *ELEMENT OUTPUT

S,PE,LE,PEEQ,PEEQT,DAMAGEC,DAMAGET,SDEG *OUTPUT, HISTORY, VARIABLE=PRESELECT *ELEMENT OUTPUT, ELSET=EOUT SP1

*NODE OUTPUT, NSET=NOUT U

*END STEP **

*STEP, NLGEOM, UNSYMM=YES STEP 2 - HYDROSTATIC LOAD *STATIC

1.0E-10, 1.0E-10 *DLOAD

WDAM, HP4, 900067.6, 91.75, 0 *RESTART, WRITE, FREQUENCY=1 *END STEP

Abaqus/Explicit input files

8、Seismic analysis of the Koyna dam, not including hydrodynamic interactions; requires import of the results from koyna_std_to_xpl.inp. *HEADING

KOYNA DAM: SEISMIC ANALYSIS Units - n, m, sec

*********************************************** ** Step 3: Earthquake

** (without hydrodynamic interactions) **

** Requires results from koyna_std_to_xpl,

** FORTRAN subroutine koyna_hp_vdload, and ** amplitude curves koyna_haccel.inp and ** koyna_vaccel.inp **

** Execution command:

** abaqus job=koyna_xpl oldjob=koyna_std_to_xpl ** user=koyna_hp_vdload double

***********************************************

*IMPORT, STEP=2, INCREMENT=1, STATE=YES, UPDATE=NO DAM

*NSET, NSET=QA_TEST NOUT

*ELSET, ELSET=QA_TEST EOUT

*MATERIAL, NAME=CONCRETE *ELASTIC

3.1027E+10, 0.2 *DENSITY 2643.0

*DAMPING, BETA=0.00323

*CONCRETE DAMAGED PLASTICITY 36.31

*CONCRETE COMPRESSION HARDENING 13.0E+6, 0.000 24.1E+6, 0.001

*CONCRETE TENSION STIFFENING, TYPE=DISPLACEMENT 2.9E+6 ,0

1.94393E+6 ,0.000066185 1.30305E+6 ,0.00012286

0.873463E+6 ,0.000173427 0.5855E+6 ,0.00022019 0.392472E+6 ,0.000264718 0.263082E+6 ,0.000308088 0.176349E+6 ,0.00035105 0.11821E+6 ,0.000394138 0.0792388E+6 ,0.000437744 0.0531154E+6 ,0.000482165

*CONCRETE TENSION DAMAGE, TYPE=DISPLACEMENT 0 ,0

0.381217 ,0.000066185 0.617107 ,0.00012286 0.763072 ,0.000173427 0.853393 ,0.00022019 0.909282 ,0.000264718 0.943865 ,0.000308088 0.965265 ,0.00035105 0.978506 ,0.000394138 0.9867 ,0.000437744 0.99177 ,0.000482165

*AMPLITUDE, NAME=HAMP, INPUT=koyna_haccel.inp *AMPLITUDE, NAME=VAMP, INPUT=koyna_vaccel.inp *FILTER, NAME=SMOOTH, TYPE=BUTTERWORTH 2, 4

*STEP, NLGEOM=YES STEP 3 - EARTHQUAKE *DYNAMIC, EXPLICIT , 10.0 *DLOAD

DAM, GRAV, 9.81, 0, -1 WDAM, P4NU

*BOUNDARY,TYPE=ACCELERATION,AMPLITUDE=HAMP NBASE,1,1,9.81

*BOUNDARY,TYPE=ACCELERATION,AMPLITUDE=VAMP NBASE,2,2,9.81 **

*OUTPUT, FIELD, VAR=PRESELECT, NUMBER INTERVAL=40 *ELEMENT OUTPUT

S,PE,LE,PEEQ,PEEQT,DAMAGEC,DAMAGET,SDEG *OUTPUT,HISTORY

*ELEMENT OUTPUT, ELSET=EOUT SP1 …… 此处隐藏：2760字，全部文档内容请下载后查看。喜欢就下载吧 ……