【作者】 李湛；

【导师】 栾茂田；

【作者基本信息】 大连理工大学 ， 岩土工程， 2006， 博士

【摘要】 我国西北与西南地区拥有丰富的水利水电资源,同时这些地区的河谷条件、地质条件和环境条件极度复杂,且一般位于高烈度地震活动区。随着我国西北大开发战略计划的实施,越来越多的高坝等大型水利水电工程基础设施正在建设或正处于规划与设计中。土石填筑坝能够有效地就地利用当地筑坝材料,能够适宜于不良地基与河谷条件,具有良好的抗震性能,往往是一种安全经济、具有潜力的坝型。同时在海岸防护与防洪加固中修建了大量的堤坝、河堤与岸堤。但是复杂条件下土石填筑坝的抗震性能与抗震安全性评价一直是设计与建设中的关键技术问题,目前仍尚未得到彻底解决,传统上所采用的拟静力分析方法存在其固有的缺陷,而目前正在发展的有限元动力分析方法尚不成熟而无法在工程实际中推广应用,单凭拟静力抗滑稳定性安全系数无法判断土石填筑坝与堤坝的实际抗震安全度与可靠性,而基于地震永久变形或地震滑移量的堤坝动力稳定性评价方法目前尚不完善而未达到工程实用阶段,因此必须发展理论上合理、工程上实用的土石坝抗震性能分析方法与抗震稳定性评价方法,这不仅对于促进土动力学与岩土地震工程学科的发展具有重要的学术意义,而且对于改进土石填筑坝的抗震设计具有重要的实际应用价值。为此,为了有效地改进这些问题,本文通过理论分析与数值计算等手段对于土石填筑坝地震永久变形分析方法与抗震稳定性评价方法进行了比较深入而系统的探索,主要研究工作包括: 1.基于传统的解析型一维剪切条模型,提出了土石坝地震动力响应分析的离散型剪切条模型。首先将坝体沿坝高离散为一系列具有不同剪切刚度与阻尼比等参数的层状体系,建立了各层的振动控制方程及其初边值条件,进而运用数学物理方程方法进行了数学求解,确定了体系的振动特性,并根据振型分解法和Duhamel积分确定了坝体地震动力响应的线性解答。进而采用等价线性化方法考虑坝料的动力非线性性质,通过对线弹性响应的反复迭代计算使得各层土的剪切模量和阻尼比与其相应的剪应变水平相协调,确定了与非线性坝体系统相等效的线性解答,并以此所得到的地震响应作为非线性地震响应的近似解。最后分别以均质土坝和心墙坝作为算例,进行了具体数值计算,并将所得结果与有限元数值解进行了对比分析,论证了所建议分析方法与计算模型的合理性及适用性。由此表明:离散型剪切条法可以作为一种工程实用分析方法,为土石坝抗震性能简化分析、方案论证阶段的多工况对比计算或工程设计初步阶段的抗震优化设计提供一种简便工具。更多还原

【Abstract】 As the development of water resources in South- and North-Western areas of China where the topographical condition and engineering geological as well as environmental conditions are often rather complex and the seismic activity is sometime high, large-scale major infrastructures such as high dams have been built or to be constructed or under planning or design. As a potential type of dam structures, earth and rock-fill dams (ERFD) or concrete-faced rockfill dams (CFRD) have been most frequently employed in engineering practice. However, the evaluation of earthquake-resistant behavior and seismic safety of these high dams during earthquake is one of the key technical issues in the design and construction of earth or rockfill dams and concrete-faced rockfill dams. It has been not well solved at present. At the same time, embankments are widely used in the river banks and coastal banks as a protective precaution. It has been recognized that the conventional pseudo-static method for evaluating seismic stability of dams or embankments has many shortcomings while the nonlinear dynamic finite element method for analyzing seismic response of dams and embankments subjected to earthquake shaking is somewhat complicated for practicing engineers. At the same time the seismic safety cannot be authentically evaluated by the single factor of safety against sliding during transient overloading. The sliding movement or overall permanent deformation was proposed by Newmark in 1965 to evaluate the overall dynamic stability of dams and embankments. In order to get good understanding and then to offer rational evaluation of earthquake-resistant characteristics and seismic stability of dams or embankments, the numerical methods for evaluating seismic dynamic response and for assessing the earthquake-induced residual deformations or displacements are to be developed in the thesis. The investigation mainly consists of the following portions.1. As a simplified analytical procedure, the shear-beam or shear-slice model is widely employed for evaluating dynamic response of earth dams or embankments. In this thesis, the modification on the equivalent linear technique for analysis of seismic response of dams and embankments which is based on one-dimensional shear-slice model is made. In the proposed method, the discrete shear-slice model is developed to consider the layered distribution of shear moduli and damping ratios of soils along depth. The differential equations for governing the vibration behavior of dam or embankment are solved mathematically under the assumption that the soils of dam display visco-linear-elastic behavior and the closed-form solutions for the natural characteristics and seismic response of the dam are obtained. The resulting linear solutions are incorporated with equivalent linearization technique for seismic response of dam or soil strata. A series of trial-and-error iterations are performed in order to make the moduli and damping ratios of all individual layers be compatible with their strain amplitudes. Then a linear system which is overall equivalent to nonlinear dam system is achieved and the seis-mic response for this equivalent linearization system can be taken as a good approximation of nonlinear seismic response of dam. Finally, example studies are performed for both a homogeneous dam and a clay-core dam and the numerical results are compared with the solutions obtained by the finite element method to illustrate the rationality and applicability of the proposed method.2. A numerical procedure by using elasto-plastic finite element method based on shear strength reduction technique is presented for pseudo-static seismic stability analysis of dams and embankments. In the proposed procedure, two alternatives are employed to determine the seismically-induced inertia loads as following: (1) The inertia forces are computed by using the routine employed in the pseudo-static limit equilibrium slices method proposed in earthquake-resistant design code of earth and rockfill dams; (2) Based on the dynamic response computed by the nonlinear dynamic finite element analysis of dams or embankments, the inertia forces corresponding to the computed nodal acceleration are defined by a given empirical manner. In addition to the seismically-induced inertia loads which is distributed non-uniformly along dam height, the seepage forces or dynamic pore water pressure caused by steady seepage or normal drop of reservoir are taken into account. Under such coupling circumstances, the elasto-plastic finite element analysis is performed to assess the pseudo-static safety factors of the dam by using the shear strength reduction technique. Finally, numerical computations for example dams are conducted to illustrate the rationality and reliability of the proposed analysis method, and the results obtained from two considerations of inertia loads induced by earthquake are compared. It is shown that the finite element procedure for seismic stability analysis of dams or embankments is proposed to prefer to be used in practice.3. In order to evaluate the earthquake-iriduced sliding displacements of a potential sliding mass in dam slope, the average seismic response of the potential sliding mass assessed on the basis of seismic response of dam or embankment subjected to earthquake shaking computed by equivalent dynamic analysis are compared with the critical safety factor to define the condition and state of the transient overloading. Then the rigid-plastic sliding-block model for evaluating sliding movement proposed by Newmark is used estimate the sliding displacement caused by such an overloading. All the sliding displacement caused by all individual overloading are accumulated to obtain the overall movement induced by earthquake. In the proposed method, the Newmark model is modified in order to take account the degradation of strength of soils due to the build-up of pore water pressure caused by earthquake-induced cyclic or transient loading. At the same time, the effect of vertical acceleration of potential sliding mass induced by earthquake shaking on seismic stability of dam which is previously overlooked are taken into account. Two types of dynamic loading including sine-type excitation and earthquake-induced irregular excitation are considered. Finally, numerical analyses for typical cases are made and a number of findings are given through comparison.4. As an alternative mechanism of earthquake effect on dam or embankment, a finiteelement procedure is proposed for evaluating seismically-induced permanent deformation of dam or embankment. In the proposed procedure, both concepts of equivalent nodal force and step-by-step gradually softening modulus are integrated together. The earthquake duration is divided into a certain number of time increments and for each time increments the residual strain and dynamic pore water pressure which is likely induced during previous time increments under undrained condition are estimated on the basis of the stress condition obtained by the finite element dynamic analysis and the empirical patterns of both residual strain and pore water pressure achieved experimentally. Then, the computed accumulative dynamic pore-water pressure at the end of each time increment is used directly to modify the static hyperbolic relationship of stress and strain which is to be used for the next time period, and at the same time, the equivalent nodal forces equivalent to incremental residual strain potential are defined. By using the modified stress-strain relationship, the incremental deformations are computed when the nodal forces equivalent to earthquake effect on dam defined as above are imposed on the dam or embankment. The computed incremental displacements of the dam or embankment for each time increments are accumulated and the accumulative displacements can be regarded as approximation of the residual deformation which is to be initiated by earthquake shaking. In fact, the proposed numerical procedure has taken both the inertia effect and the softening effect of earthquake-induced loading on the permanent deformation of the dam or embankment into consideration.更多还原