Transient analysis of 1D inhomogeneous media by dynamic inhomogeneous finite element method

The dynamic inhomogeneous finite element method is studied for use in the transient analysis of one dimensional inhomogeneous media. The general formula of the inhomogeneous consistent mass matrix is established based on the shape function. In order to research the advantages of this method, it is compared with the general finite element method. A linear bar element is chosen for the discretization tests of material parameters with two fictitious distributions. And, a numerical example is solved to observe the differences in the results between these two methods. Some characteristics of the dynamic inhomogeneous finite element method that demonstrate its advantages are obtained through comparison with the general finite element method. It is found that the method can be used to solve elastic wave motion problems with a large element scale and a large number of iteration steps.     

An improved pseudo-static method for seismic resistant design of underground structures

This paper describes a commonly used pseudo-static method in seismic resistant design of the cross section of underground structures. Based on dynamic theory and the vibration characteristics of underground structures, the sources of errors when using this method are analyzed. The traditional seismic motion loading approach is replaced by a method in which a one-dimensional soil layer response stress is differentiated and then converted into seismic live loads. To validate the improved method, a comparison of analytical results is conducted for internal forces under earthquake shaking of a typical shallow embedded box-shaped subway station structure using four methods： the response displacement method, finite element response acceleration method, the finite element dynamic analysis method and the improved pseudo-static calculation method. It is shown that the improved finite element pseudo-static method proposed in this paper provides an effective tool for the seismic design of underground structures. The evaluation yields results close to those obtained by the finite element dynamic analysis method, and shows that the improved finite element pseudo-static method provides a higher degree of precision.     

Seismic damage analysis of the outlet piers of arch dams using the finite element sub-model method

This study aims to analyze seismic damage of reinforced outlet piers of arch dams by the nonlinear finite element(FE) sub-model method. First, the dam–foundation system is modeled and analyzed, in which the effects of infinite foundation, contraction joints, and nonlinear concrete are taken into account. The detailed structures of the outlet pier are then simulated with a refined FE model in the sub-model analysis. In this way the damage mechanism of the plain(unreinforced) outlet pier is analyzed, and the effects of two reinforcement measures(i.e., post-tensioned anchor cables and reinforcing bar) on the dynamic damage to the outlet pier are investigated comprehensively. Results show that the plain pier is damaged severely by strong earthquakes while implementation of post-tensioned anchor cables strengthens the pier effectively. In addition, radiation damping strongly alleviates seismic damage to the piers.     

Study of Finite Element Modeling of Strong Earthquake Activities and Its Preliminary Application—Taking Southwest China as an Example

Based on research result concerning the preparation and activity of strong earthquakes in groups and using the finite element method, a finite element dynamic model for Southwest China is established in this paper. Using this model, the stress adjustment in the whole of the Southwest China region in response to the stress change due to strong earthquake occurrence is studied. The preliminary result shows that many strong earthquakes occurred in areas where the stress heightened after the last strong earthquake. So, the finite element model set up in this paper is useful for judging the regions where strong earthquakes are likely to occur in future.     

Temporal and spatial distribution of dam failure events in China

Reservoirs play a vital role in economic development and flood control. Nevertheless, both human and natural factors may lead to dam failures with catastrophic consequences. Analyzing the data of dam failure events from 1954 to 2003 and using the method of energy spectrum analysis, this paper studies the periodicity of dam failures. The rate of failure of different dam ages is analyzed. The climate is the main factor affecting the rate of dam failure. Climate diagrams are used to analyze the spatial distribution of dam failure events in China. High rate of dam failure occurs with 25-year and 12.5-year periods. The distribution of the percentage of dam failure shows an L-shape as a function of service age. The first 5 years of operation is known as the ＂infant period＂, during which, the probability of dam failure is much higher than during any other periods. The failure rate in areas near or north to the 400 mm annual isopluvial line is notably higher than other areas. In areas with high temperature difference among seasons have a high annual average dam failure rate.     

Reliability and sensitivity analysis of wedge stability in the abutments of an arch dam using artificial neural network

In this study, the seismic stability of arch dam abutments is investigated within the framework of the probabilistic method. A large concrete arch dam is considered with six wedges for each abutment. The seismic safety of the dam abutments is studied with quasi-static analysis for different hazard levels. The Londe limit equilibrium method is utilized to calculate the stability of the wedges in the abutments. Since the finite element method is time-consuming, the neural network is used as an alt...     

Stability of an explicit time-integration algorithm for hybrid tests,considering stiffness hardening behavior

An explicit unconditionally stable algorithm for hybrid tests,which is developed from the traditional HHT-α algorithm,is proposed.The unconditional stability is first proven by the spectral radius method for a linear system.If the value of α is selected within [-0.5,0],then the algorithm is shown to be unconditionally stable.Next,the root locus method for a discrete dynamic system is applied to analyze the stability of a nonlinear system.The results show that the proposed method is conditionally stable for dynamic systems with stiffness hardening.To improve the stability of the proposed method,the structure stiffness is then identified and updated.Both numerical and pseudo-dynamic tests on a structure with the collision effect prove that the stiffness updating method can effectively improve stability.     

The Effect Analysis of Topography on the Spectrum Properties of Anti-Plane Movement

An explicit finite difference method is performed to evaluate the significance of topography effects on the seismic response of viscoelastic sites. The real accelerograms and bell-shape impulse wavelets, slope angle of cliff and angle of incidence on the spectrum property of the incoming waves are all discussed in detail. The conclusions show that the presence of topographic factors and the direction of the incident waves not only significantly affect the peak ground acceleration of a site, but also affect the spectrum properties of the incoming waves, and the effect on the direction of the incident waves is especially remarkable. The study also finds that it is reasonable to analyse spectrum properties with the input of impulse wavelets in the finite element simulation of wave motion.     

Wave passage and incoherency effects on seismic response of high arch dams

The effects of incoherency and wave-passage on the nonlinear responses of concrete arch dams are investigated in this study. A double curvature arch dam is selected as a numerical example. The reservoir is modeled as a compressible material and the foundation is modeled as a massless medium. Ground motion time-histories are artificially generated using the Monte Carlo simulation approach. Four different finite element models (FEM) are considered: uniform excitation; incoherence effect; wave passage effect; and both incoherence and wave passage effects. It was revealed that modeling multiple-supports excitation could have a significant impact on the structural response of the dam by inducing a pseudo-static effect. Also, it was concluded that the coherency effect overshadows the wave passage effect and the results obtained from non-uniform excitation of FEM, including the wave passage effect, is close to the results of the FEM when it is uniformly excited.     

Study of damage to a high concrete dam subjected to underwater shock waves

To study the effect of a strong underwater shock wave on a concrete dam, this research aims to improve hammer impact methods in model tests. Six 1:200 scale models were designed and tested under distributed impact loads. A device was deployed for a direct measurement of the impact force at the upstream face of the dams. The model dam bases were anchored to prevent displacement. The experimental results indicate that the top part of the concrete dam is a weak zone, and the impact failure initiates with a fracture on the top of the dam. The peak value of impact stress increases when the second crack appears in the concrete dam from the upstream face to the downstream face. And, the level of the second crack in the dam body is lower as the peak value of impact stress increases. In this study, dynamic analysis was conducted by calculating the results to verify the effectiveness of a device to directly measure the impact force. This method may be used to approximately forecast the damage of concrete dam and may also be useful in other engineering applications.     

三种介质耦联系统动力反应分析的显式有限元法及其应用

在笔者已建立的流体饱和多孔介质动力分析的显式有限元法的基础上，提出了可分析任意形状的流体饱和多孔介质-单相弹性固体介质-理想流体介质耦联的复杂系统的动力响应的显式有限元方法．该方法建立的有限元方程列式具有解耦特征，不需求解联立方程组，因而极大地提高了计算效率．将这一方法用于分析了考虑库水、坝、淤泥层和基岩这一复杂系统动力相互作用问题的斜坝面的地震响应，并给出了一些计算结果．      

考虑水-饱和土场地-结构耦合时的沉管隧道地震反应分析

用理想流体介质模拟水层、流体饱和多孔介质模拟饱和土地层,在理想流体介质与流体饱和多孔介质相连接边界的连续条件基础上,结合已有的对理想流体介质、流体饱和多孔介质进行动力反应分析的显式有限元方法,开考虑地层与结构的动力相互作用,建立了进行水与场地、结构耦合动力分析的方法。利用该方法对沉管隧道的地震响应进行了研究,重点分析了水深、地质条件等因素对沉管隧道地震反应的影响,并从中得出了一些可供相关人员进行沉管隧道抗震分析时参考的结论。     

库水模型及库水深度对拱坝动力响应的影响

基于接触非线性有限元模型,以锦屏一级拱坝为例,库水分别采用附加质量模型、可压缩流体有限元模型、不可压缩流体有限元模型计算了正常蓄水位及运行低水位时坝体的动力响应,结果表明:库水模型对拱坝动力响应有较大影响,随库水深度的增大,各模型计算结果差异增大;相比于流体可压缩模型,采用不可压缩流体模型所得动力响应普遍偏大;运行低水位工况,由于静水压力减小导致拱效应减弱,从而降低了拱坝的整体性,因此运行低水位工况各缝开度普遍高于正常蓄水位工况,且其拉应力范围较大,因此,运行低水位工况将对抗震设计起控制作用。     

淤积泥砂对水平地运动作用时刚性坝面动压力的影响

文章基于固体介质，流体介质，液固两相饱和介质的时域显示有限元波动分析方法，研究了可压缩库水条件下淤积泥砂对水平地运动作用时刚性坝面上动压力的影响，研究中分析了作为两相介质处理的淤积泥砂的饱和度、厚度、渗透系数和孔隙率对坝面动压力的影响，还比较了作为流程两相介质，单相固体介质、单相重流体介质处理的不同淤积泥砂层模型间的计算结果差异。     

饱和多孔介质三维时域黏弹性人工边界与动力反应分析的显式有限元法

本文基于Biot的饱和多孔介质本构方程,考察具有辐射阻尼的外行球面波,推导了饱和多孔介质三维黏弹性人工边界的法向和切向边界方程;在已有的饱和多孔介质二维显式有限元数值计算方法基础上,提出该理论的三维方法,并开发了实现该三维方法的有限元程序.算例表明饱和多孔介质三维时域黏弹性人工边界与动力反应分析的显式有限元法具有较好的精度和稳定性.     

无限元法在深覆盖层土石坝动力分析中的应用

以复杂地质条件的西藏高原地区旁多水电站实际工程为背景,探讨了有限元法与无限元法在深覆盖层悬挂式防渗墙土石坝结构非线性动力分析中的差异,分析和比较了无限元法和有限元法在西藏高原旁多水电站土石坝结构三维模型地震动响应的计算结果,验证了引入无限元法模拟西藏高原地区地基中复杂地质条件下的无限域或半无限域问题的可靠性和精确性。结果表明：在西藏高原地区复杂地质条件下,无限元法相比有限元法能量弥散现象较为明显。通过分析得到了西藏高原具有深覆盖层坝体结构地震动响应的规律性分析结果,为无限元法在西藏高原地区此类工程中的应用提供参考。     

地下输液管道动力反应分析

本文利用有限元方法,考虑介质、埋设管线及土壤的相互作用,进行了埋地管线动特性和地震动响应研究,得出了土特性、管内液体流动对系统固有频率及动力响应影响的若干初步结论。     

二维土坝的两相动力分析

提出了一种分析饱和土坝动力反应的方法,考虑了土坝的两相介质特性,在固液耦联动力方程的基础上,选取固相位移,液相位移、孔隙水压作为场变量,采用伽辽金加权残数法进行有限元空间离散化,然后在时域上采用Ｗｉｌｓｏｎ－θ法进行逐步积分。该方法不仅能计算出固相位移和液相位移,而且能直接得到孔隙水压的反应过程。文中以一饱和土坝模型进行算例分析,并与将其作为单相介质时的结果进行了比较。该法可用于分析饱和介质的地震     

考虑坝体材料非线性的混凝土拱坝地震反应分析

本文基于可信概率水准的破坏性强震作用,针对小湾高拱坝进行了考虑坝体材料非线性的拱坝地震反应分析。在分析模型中,同时考虑了无约束域地震能量辐射效应和近域地基材料非均匀性的影响。为了实现非线性条件下的静、动力组合分析,利用显式有限元结合修正的黏弹性人工边界的开放系统时域静、动力统一分析方法进行了求解,对在自重作用下的初始静力解计算采用了动力松弛技术。     

考虑堆石料空间变异性的土石坝坝坡地震稳定性随机有限元分析

本文探讨了筑坝堆石料的空间变异性对土石坝坝坡动力稳定性的影响。以新疆某在建高面板堆石坝为例,在蒙特卡洛法的框架下,采用基于局部平均细分法的随机有限元法模拟考虑筑坝堆石料空间变异性时土石坝的地震响应及坝坡滑移情况,通过对比随机有限元法和常规确定性有限元法的计算结果,提出:在地震动作用下,考虑筑坝材料空间变异性时,坝坡滑动体的数量、规模以及滑移量和滑动历时都有不同程度的增大,因而坝坡整体危险程度显著高于不考虑材料空间变异性的情况。坝坡各项动力安全性指标对筑坝材料空间变异性非常敏感;因而,考虑筑坝材料空间变异性时,各项安全性指标的离散性较大。