共查询到20条相似文献，搜索用时 31 毫秒
1.
时域局部人工边界是以行波概念为基础的,原则上不能用于模拟零频分量（消减运动）。Ｗｏｌｆ提出了能同时模拟零频分量和行波的双渐近－多向透射边界。本文指出,在这一组合中的双渐近边界的高频渐近项是不合理的,还查能降低模拟精度。据此,本文提出一种概念清楚、简便实用的组合方案－多次射透－静力边界,并通过数值试验对此方案进行检验。 相似文献
2.
考虑各向异性的层土-盾构隧道地震反应数值模拟 总被引：3，自引：0，他引：3
在层状各向异性土体-盾构隧道地震反应分析中,引入了横观各向同性弹塑性模型理论,建立了适合于横观各向同性介质的双渐近-多向透射边界条件。针对地铁区间盾构隧道抗震设计的特点,基于横观各向同性弹塑性模型,研制了考虑层状土体各向异性和施工开挖效应,适合于盾构隧道动力计算的各向异性弹塑性动力有限元程序。在程序中对于不同的材料采用了不同的本构关系和不同的单元形式,并采用了关联流动法则和多种屈服准则,可同时进行各向异性土体与地下结构的二维平面应力、平面应变和轴对称问题的静力和动力数值分析。最后利用所研制的程序进行了上海地铁二号线石门一路站附近区间隧道在不同超越概率地震动输入下的隧道反应计算。结果表明,在层状土体-地铁区间隧道的抗震分析中考虑土体各向异性的影响是必要的,所提出的计算模型是可行的。 相似文献
3.
时域高阶双渐近透射边界能够同时模拟层状介质中行波和快衰波的传播,具有很高的计算精度和计算效率.本文将高阶双渐近透射边界推广应用到多层层状地基系统弹性波传播问题的模拟,采用广义特征值分解分析该透射边界的数值稳定性,通过移谱法消除导致数值不稳定的虚假模态.将高阶双渐近透射边界以超单元的形式直接嵌入到近场有限元方程,建立了有限元-高阶双渐近透射边界时域耦合分析模型,并将其应用到重力坝-层状地基动力相互作用分析.数值算例分析结果表明,该时域耦合分析模型具有很高的精度和计算效率,适用于实际重力坝工程的地震响应分析. 相似文献
4.
The incident steady-state SH waves in the neighbourhood of a boundary node point are first expanded in terms of a few plane waves at different incident angles, the motion of the boundary point is then related to those of its adjacent node points to form the multi-directional transmitting boundary conditions, which reduce reflections from the boundary over the entire range of incident angles. This feature of the boundaries is demonstrated by finite element computations of SH motions of a layer fixed at its base, in which the incident wave may travel parallel to the boundary as the frequency of excitation approaches one of the natural frequencies of the system. 相似文献
5.
高阶双渐近时域透射边界能够同时模拟行波和快衰波的传播,并且能够在全频范围内迅速逼近准确解,具有优良的收敛性能和计算效率.本文将动水压力波高阶双渐近透射边界直接嵌入到近场有限元方程中,建立了大坝-库水动力相互作用的直接耦合分析模型.该模型的整体控制方程保留了近场有限元方程系数矩阵对称稀疏的优势,可以方便地利用现有的通用有限元求解器求解.基于有限元开源软件框架体系OpenSees（Open System for Earthquake Engineering Simulation）,编程实现了直接耦合分析模型,并将其应用于二维重力坝、三维拱坝与库水动力相互作用分析.数值算例表明,该直接耦合分析模型具有很高的精度和计算效率. 相似文献
6.
7.
A local time-domain transmitting boundary for simulating cylindrical elastic wave propagation in infinite media 总被引：3，自引：0，他引：3
Finite element simulation of the time-dependent wave propagation in infinite media requires enforcing the transmitting boundary to replace the truncated far-field infinite domain so as to model the effect of the wave radiation towards infinity. This paper proposed a novel local time-domain transmitting boundary for simulating the cylindrical elastic wave radiation problem. This boundary is a mechanical model consisting of the spring, dashpot and mass elements, with the auxiliary degrees of freedom introduced, which is dynamically stable and easily implemented into the commercial finite element codes. Numerical analysis of the cylindrical elastic wave radiation problem indicates that the proposed transmitting boundaries with the order N=3 for cylindrical P and SV waves and with the order N=4 for cylindrical SH wave have very high accuracy, even when the artificial boundary at wave source. The proposed transmitting boundary with order N=0 can be applied approximately to the general two-dimensional infinite elastic wave problems that contain the more complex outgoing wave fields at artificial boundary than the cylindrical waves. The plane-strain Lamb problem is analyzed with the acceptable engineering accuracy achieved. On the other hand, the proposed transmitting boundary with higher order can be a tool to localize the temporal convolution that appears in an exact time-domain transmitting boundary for the general infinite wave problems. This potential applicability is mentioned. 相似文献
8.
A new artificial boundary model based on multi-directional transmitting and viscous-spring artificial boundary theories is proposed to absorb stress waves in a saturated soil foundation in dynamic anal... 相似文献
9.
传统时间域航空电磁全波形正演模拟主要采用间接法(褶积算法)和直接法(时域有限差分方法等),然而褶积算法需要获得精确的电流二阶导数,这给发射电流数据采集工作带来极大挑战;时域有限差分方法受到网格和时间步长的严格限制,缺乏灵活性.为解决这些问题,本文采用时域有限元方法,通过直接改变每个时间道上的瞬时电流强度模拟任意发射波形的电磁响应.由于无需计算电流二阶导数,大大提高了正演结果的精度.利用基于非结构四面体网格的矢量有限元方法和后推欧拉技术对时间域电场扩散方程进行空间和时间离散,实现三维航空电磁时间域全波形的直接正演模拟.由此不仅可以模拟复杂的地电结构,而且基于后推欧拉法的无条件稳定性,可以更加灵活地选取时间步长,提高计算效率.通过与1D数值模拟结果进行对比验证了该方法的准确性.本文对三维柱状体模型上HELITEM MULTIPULSE和VTEM系统实际发射波形电磁响应进行模拟,并与褶积算法的结果进行比较,验证了本文算法模拟实际发射波形电磁响应的优越性.对复杂三维地质体模型上不同发射波形电磁响应进行模拟,验证了时间域有限元算法可有效处理复杂地下地质结构. 相似文献
10.
A local transmitting boundary is presented in a compact form, which can be directly incorporated into finite elements. The accuracy of the boundary is studied thoroughly for a one-dimensional model in order to clarify numerical instabilities introduced by the boundary. Discretization of the model and reflection from the boundary are rigorously considered in the study, and the mechanism of the instability is then illuminated in the frequency domain by the amplification of reflection from the boundary and the multi-reflection of wave motion in a finite computational region. Typical characteristics of the instability in the time domain are illustrated by numerical results of the simple model and explained completely by the mechanism. On the basis of this understanding of the mechanism, a modified transmitting boundary is devised and its stability criterion is given for the one-dimensional model. 相似文献
11.
A three-dimensional finite-element time-domain forward-modelling algorithm is developed to simulate transient electromagnetics excited by grounded-wire sources. The main advantage of this finite-element time-domain algorithm is that full transmitting-current waveforms and complex-shaped sources resulting from topography can be directly dealt with in this algorithm. The models used to test this algorithm include a homogeneous half-space model, a stratified-medium model, the model of a complex conductor at a vertical contact and the Ovoid Zone massive sulfide deposit at Voisey's Bay, Canada. The homogeneous half-space model is used to determine the truncation boundary for a computational domain, and to compare with the electromagnetic responses excited by step-off, step-on and direct current waveforms. For the stratified-medium model, results demonstrate that full transmitting waveforms have strong effects on the observed electromagnetic responses. The model of a complex conductor at a vertical contact is designed for the grounded electrical source airborne transient electromagnetic method and is also used to examine the effectiveness of the broadside and inline configurations for such a vertical, thin plate embedded in the subsurface. The area of the Ovoid Zone massive sulfide deposit possesses non-negligible topography, the effects of which on the shapes of the grounded-wire sources must be taken into account when implementing the finite-element time-domain solver. The results show that both the broadside and inline electromagnetic responses are strongly affected by the massive conductive ore body. 相似文献
12.
从Maxwell方程组出发,推导了探地雷达(GPR)有限元波动方程.阐述了透射边界条件和Sarma边界条件的原理,推导了这两种边界条件的理论公式;通过在衰减层内加入过渡带优化了Sarma边界条件的加载方法,压制了介质区和衰减层交界面处的人为反射.考虑到透射边界条件与Sarma边界条件不同的理论机制,提出了一种结合透射边界条件和Sarma边界条件的混合边界条件,它利用Sarma边界条件对到达边界区域的GPR波能量衰减功能和透射边界对GPR波能量的透射功能,使GPR波经过Sarma边界条件的衰减吸收后,再通过透射边界条件将剩余能量透射出去,集成了二者的优势.并以二维均匀模型中的中心脉冲激励源方式为例,通过Matlab程序实现,以GPR的全波场快照的直观方式,对比了有、无边界条件及不同边界条件对人工截断边界的处理效果,说明了该混合边界条件对到达截断边界处的GPR波的处理优于单一边界条件.最后,以基于混合边界条件的有限单元法对两个典型的GPR地电模型进行了正演模拟,指导了GPR数据处理与工程实践. 相似文献
13.
John L. Tassoulas 《地震工程与结构动力学》2011,40(5):531-550
A half‐space finite element and a transmitting boundary are developed for a water‐saturated layered half‐space using a paraxial boundary condition. The exact dynamic stiffness of a half‐space in plane strain is derived and a second‐order paraxial approximation of the stiffness is obtained. A half‐space finite element and a transmitting boundary are then formulated. The development is verified by comparison of the dynamic stiffness of impermeable and permeable rigid strip foundations with other published results. The advantage of using the paraxial boundary condition in comparison with the rigid boundary condition is examined. It is shown that the paraxial boundary condition offers significant gain and the resulting half‐space finite element and transmitting boundary can represent the effects of a water‐saturated layered half‐space with good accuracy and efficiency. In addition, the numerical method described herein maintains the strengths and advantages of the finite element method and can be easily applied to demanding problems of soil–structure interaction in a water‐saturated layered half‐space. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
14.
Steven D. London 《地球物理与天体物理流体动力学》2013,107(6):616-628
We consider an electrically conducting fluid confined to a thin rotating spherical shell in which the Elsasser and magnetic Reynolds numbers are assumed to be large while the Rossby number is assumed to vanish in an appropriate limit. This may be taken as a simple model for a possible stable layer at the top of the Earth's outer core. It may also be a model for the thin shells which are thought to be a source of the magnetic fields of some planets such as Mercury or Uranus. Linear hydromagnetic waves are studied using a multiple scale asymptotic scheme in which boundary layers and the associated boundary conditions determine the structure of the waves. These waves are assumed to be of the form of an asymptotic series expanded about an ambient magnetic field which vanishes on the equatorial plane and velocity and pressure fields which do not. They take the form of short wave, slowly varying wave trains. The results are compared to the author's previous work on such waves in cylindrical geometry in which the boundary conditions play no role. The approximation obtained is significantly different from that obtained in the previous work in that an essential singularity appears at the equator and nonequatorial wave regions appear. 相似文献
15.
本文在全部入射角范围内系统地讨论了时域局部透射边界真实的反射系数及其与理想稳态和理想暂态反射系数的关系。文中指出，后两者可分别用作前者上，下限的估计；同时，就时域逐步积分计算而言理想暂态较之理想稳态更接近真实情形，简便的理想暂态反射系数可用于在全部入射角范围内度量边界是的误差。 相似文献
16.
Based on the u–p formulation of Biot equation with an assumption of zero permeability coefficient, a high-order transmitting boundary is derived for cylindrical elastic wave propagation in infinite saturated porous media. By this transmitting boundary the total stresses on the truncated boundaries of a numerical model, such as a finite element model, are replaced by a set of spring, dashpot and mass elements, with some additionally introduced auxiliary degrees of freedom. The transmitting boundaries are incorporated into the DIANA SWANDYNE II program and an unconditionally stable implicit time integration algorithm is adopted. Despite the assumption made in the derivation of the transmitting boundary, numerical examples show that it can provide highly accurate results for cylindrical elastic wave propagation problems in infinite saturated porous medium in case the u–p formulation is applicable. Although the direct applications of the proposed transmitting boundary to general two dimensional wave problems in infinite saturated porous media are not highly accurate, acceptable accuracy can still be achieved by placing the transmitting boundary at relatively large distance from the wave source. 相似文献
17.
土—结构相互作用通用程序的改进 总被引：2，自引：0，他引：2
本文提出了用具有工程实用价值的透射边界替换通用程序ＦＬＵＳＨ中一致边界的方法。通过实例证明、二阶透射边界与一致边界具有相当的精度。改进后的程序ＷＭＮＳ－Ｆ１比ＦＬＵＳＨ具有更高的计算效率和实用性。 相似文献
18.
The energy transmitting boundary used in programs such as FLUSH and ALUSH is a very accurate and useful technique for the earthquake response analysis of soil–structure interaction systems. However, it is applicable only to linear analyses or equivalent linear analyses, because it can be calculated only in the frequency domain. The author has proposed methods for transforming frequency-dependent impedance into the time domain. In this paper, an earthquake response analysis method for a soil–structure interaction system, using the energy transmitting boundary in the time domain, is proposed. First, the transform of the transmitting boundary matrices to the time domain using the methods proposed by the author is studied. Then, linear and nonlinear time history earthquake response analyses using the boundary are performed. Through these studies, the validity and efficiency of the proposed methods are confirmed. 相似文献
19.
A three‐dimensional transmitting boundary is formulated in the Cartesian co‐ordinate system. It is developed for the dynamic soil–structure interaction problems of arbitrary shape foundations in laterally heterogeneous strata overlying rigid bedrock. Dynamics of a rectangular rigid surface foundation on a homogeneous stratum is analysed by a hybrid approach in which the finite region including foundation is modelled by the conventional finite element method and the surrounding infinite region by the newly developed transmitting boundary. To demonstrate its strength, the present method is applied to a rectangular foundation in a horizontally heterogeneous ground consisting of two distinct regions divided by and welded along a vertical plane. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
20.
Abstract A spherical α2-dynamo is presented as an expansion in the free decay modes of the magnetic field. In the limit of vanishing viscosity the momentum equation yields various asymptotic expansions for the flow, depending on the precise form of the dissipation and boundary conditions applied. A new form for the dissipation is introduced that greatly simplifies this asymptotic expansion. When these expansions are substituted back into the induction equation, a set of modal amplitude equations is derived, and solved for various distributions of the α-effect. For all choices of α the solutions approach the Taylor state, but the manner in which this occurs can vary, as previously found by Soward and Jones (1983). Furthermore, as hypothesized by Malkus and Proctor (1975), but not previously demonstrated, the post-Taylor equilibration is indeed independent of the viscosity in the asymptotic limit, and depending on the choice of a may be either steady-state or oscillatory. 相似文献