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

*ENERGY OUTPUT ETOTAL, ALLVD

****************************** ** OUTPUT FOR QA PURPOSES ******************************

*OUTPUT, HISTORY, TIME INTERVAL=0.05, FILTER=SMOOTH *ELEMENT OUTPUT, ELSET=QA_TEST SP1

*NODE OUTPUT, NSET=QA_TEST U1,U2

****************************** *restart, write *END STEP

9、User subroutine VDLOAD used by koyna_xpl.inp to specify hydrostatic pressure. C

C User subroutine VDLOAD subroutine vdload (

C Read only (unmodifiable) variables -

* nblock, ndim, stepTime, totalTime, * amplitude, curCoords, velocity, dircos, * jltyp, sname, C Write only (modifiable) variable - * value ) C

include 'vaba_param.inc'

parameter ( ywater_level = 91.75d0, * water_density = 1000.0d0, * gravity = 9.81d0) C

dimension curCoords(nblock,ndim), * velocity(nblock,ndim), * dircos(nblock,ndim,ndim), * value(nblock) character*80 sname * C

C USER SUBROUTINE TO APPLY THE WATER

C HYDROSTATIC PRESSURE ON KOYNA CONCRETE DAM. C

do k = 1, nblock

depth = ywater_level - curCoords(k,2) value(k) = amplitude *

* depth * gravity * water_density end do * return end

10、Analysis of the post-seismic state of the Koyna Dam; requires import of the results from koyna_xpl.inp. *HEADING

KOYNA DAM: Static Analysis Units - n, m, sec

************************************************* ** Step 4: Static ** **

** Requires results from koyna_xpl

************************************************* *IMPORT, STEP=3, STATE=YES, UPDATE=NO DAM

*BOUNDARY, OP=MOD NBASE, 1,2 ** **

*STEP, NLGEOM, UNSYMM=YES STEP 4 - Static *STATIC

1.0E-10, 1.0E-10 *DLOAD

DAM, GRAV, 9.81, 0, -1

WDAM, HP4, 900067.6, 91.75, 0

***CONTROLS, ANALYSIS=DISCONTINUOUS *OUTPUT, FIELD, VAR=PRESELECT, FREQ=1 *ELEMENT OUTPUT

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

*ELEMENT OUTPUT, ELSET=EOUT SP1

*NODE OUTPUT, NSET=NOUT U1,U2

*ENERGY OUTPUT ETOTAL, ALLVD

****************************** ** OUTPUT FOR QA PURPOSES

****************************** *OUTPUT, HISTORY

*ELEMENT OUTPUT, ELSET=QA_TEST SP1

*NODE OUTPUT, NSET=QA_TEST U1,U2

****************************** *restart, write *END STEP

References

1、Bhattacharjee, S. S., and P. Léger, “Seismic Cracking and Energy Dissipation in Concrete Gravity Dams,” Earthquake Engineering and Structural Dynamics, vol. 22, pp. 991–1007, 1993. 2、Cervera, M., J. Oliver, and O. Manzoli, “A Rate-Dependent Isotropic Damage Model for the Seismic Analysis of Concrete Dams,” Earthquake Engineering and Structural Dynamics, vol. 25, pp. 987–1010, 1996.

3、Chopra, A. K., and P. Chakrabarti, “The Koyna Earthquake and the Damage to Koyna Dam,” Bulletin of the Seismological Society of America, vol. 63, no.2, pp. 381–397, 1973.

4、Ghrib, F., and R. Tinawi, “An Application of Damage Mechanics for Seismic Analysis of Concrete Gravity Dams,” Earthquake Engineering and Structural Dynamics, vol. 24, pp. 157–173, 1995.

5、Lee, J., and G. L. Fenves, “A Plastic-Damage Concrete Model for Earthquake Analysis of Dams,” Earthquake Engineering and Structural Dynamics, vol. 27, pp. 937–956, 1998.

6、Westergaard, H. M., “Water Pressures on Dams during Earthquakes,” Transactions of the American Society of Civil Engineers, vol. 98, pp. 418–433, 1933.

Tables

Table 2.1.15–1 Material properties for the Koyna dam concrete. Young's modulus: Poisson's ratio: Density: Dilation angle:

Compressive initial yield stress: Compressive ultimate stress: Tensile failure stress:

E = 31027 MPa = 0.15 = 2643 kg/m3 = 36.31o = 13.0 MPa = 24.1 MPa = 2.9 MPa

Table 2.1.15–2 Natural frequencies of the Koyna dam.

Mode Natural Frequency (rad sec–1) Abaqus Chopra and Chakrabarti (1973) 1 2 3 4 18.86 49.97 68.16 98.27 19.27 51.50 67.56 99.73

Figures

Figure 2.1.15–1 Geometry of the Koyna dam.

Figure 2.1.15–2 Finite element mesh.

Figure 2.1.15–3 Koyna earthquake: (a) transverse and (b) vertical ground accelerations.

Figure 2.1.15–4 Concrete tensile properties: (a) tension stiffening and (b) tension damage.

Figure 2.1.15–5 Horizontal crest displacement (relative to ground displacement).

Figure 2.1.15–6 Evolution of tensile damage (Abaqus/Standard); deformation scale factor = 100.