共查询到20条相似文献，搜索用时 37 毫秒
1.
A general viscous‐spring transmitting boundary for dynamic analysis of saturated poroelastic media 
下载免费PDF全文
下载免费PDF全文 A time‐domain viscous‐spring transmitting boundary is presented for transient dynamic analysis of saturated poroelastic media with linear elastic and isotropic properties. The u–U formulation of Biot equation in cylindrical coordinate is adopted in the derivation. By this general viscous‐spring boundary, the effective stress and pore fluid pressure on the truncated boundary of the computational area are replaced by a set of continuously distributed spring and dashpot elements, of which the parameters are defined assuming an infinite permeability and considering the two dilatational waves. Numerical examples demonstrate good absorption of both the two cylindrical dilatational waves by the proposed ‘drained’ boundary. For general two‐dimensional wave propagation problems, acceptable accuracy can still be achieved by setting the proposed boundary relatively far away from the scatter. Numerical comparison shows that the results obtained by using this boundary are more accurate for all permeability values than those by the traditional viscous‐spring or viscous boundaries established for u–U formulation. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
2.
3.
近海风力发电结构是由无限域地基、基础及上部结构组成的体系,有限域透射人工边界和地震动输入模式直接影响结构体系的动力响应。首先,建立透射人工边界-地基-基础-塔筒结构和固定边界-地基-基础-塔筒结构的有限元模型;然后,简要澄清了透射边界地震动输入之所以采用外力而不采取地震动物理量（加速度、速度和位移）的原因所在;最后,对给出的3种作用模式：（1）固定边界地震动输入,（2）只考虑外源输入波作用,（3）同时考虑外源输入波和内源振动的散射波作用,分别进行计算分析和比较。通过数值计算并结合理论分析,（1）澄清了对透射边界采用外力进行地震激励输入的方法的原因,这种输入方法合理可靠,符合实际情况;（2）在地震波等效荷载的生成中,针对圆形横截面地基提出了一种阻尼力、刚度力的便捷性生成方法,大大减小了工作量;（3）从自身振动特点及阻尼设置角度,解释了地震作用下采用透射边界比采用固定边界时风电结构响应减小的原因;（4）指出了地震作用下透射边界地基模型在结构动力响应与倒塌计算中内源振动反射的不容忽视性。 相似文献
4.
We present fully-discrete procedures for computing the impedance functions of rigid massless soil-structure interfaces that are embedded in arbitrarily heterogeneous half-spaces. The finite element method (FEM) is used for obtaining the wave responses of (visco-)elastic half-spaces truncated by Perfectly Matched Layers (PMLs), which provide the wave absorbing boundary conditions. The devised FEM-PML approach is verified in both time and frequency domains by using various benchmark solutions. Requirements on the prescribed input excitations for obtaining accurate impedances in the time domain as well as the relative computational cost of time- and frequency domain solutions are investigated. Accuracy of the implemented PMLs in extracting the impedance functions is also examined in comparison to Lysmer–Kuhlemeyer dashpots; and it was found that this simplified boundary treatment is generally inadequate. The utility of the proposed method is demonstrated by extracting the impedance matrix of rectangular and circular voids embedded in a linearly stiffening half-space. Impedance functions for such complex soil-structure systems are shown to be highly coupled and frequency-dependent due to wave reflections and interference caused by the soil heterogeneity and interface geometry. Fully discrete approaches, such as those proposed herein, are necessary to adequately capture these effects. 相似文献
5.
In modeling of many geomechanics problems such as underground openings, soil-foundation structure interaction problems, and in wave propagation problems through semi-infinite soil medium the soil is represented as a region of either infinite or semi-infinite extent. Numerical modeling of such problems using conventional finite elements involves a truncation of the far field in which the infinite boundary is terminated at a finite distance. In these problems, appropriate boundary conditions are introduced to approximate the solution of the infinite or semi-infinite boundaries as closely as possible. However, the task of positioning the finite boundary in conventional finite element discretization and the definition of the boundary and its conditions is very delicate and depends on the modeller's skill and intuition. Moreover, such a choice is influenced by the size of the domain to be discretized. Consequently, the dimensions of the global matrices and the time required for solution of the problem will increase considerably and also selection of the arbitrary location of truncated boundary may lead to erroneous result. In order to over come these problems, mapped infinite elements have been developed by earlier researchers (Simoni and Schrefier, 1987). In the present work the applicability of infinite element technique is examined for different geomechanics problems. A computer program INFEMEP is developed based on the conventional finite element and mapped infinite element technique. It is then validated using selected problems such as strip footing and circular footing. CPU time taken to obtain solutions using finite element approach and infinite element approach was estimated and presented to show the capability of coupled modeling in improving the computational efficiency. Mesh configurations of different sizes were used to explore the enhancement of both computational economy and solution accuracy achieved by incorporation of infinite elements to solve elastic and elasto-plastic problems in semi-infinite/finite domain as applied to geotechnical engineering. © Rapid Science Ltd. 1998 相似文献
6.
Explicit time integration schemes provide an efficient solution to non-linear dynamic finite element analyses of geotechnical problems especially when high frequency response is important. Such explicit time integration schemes require one of two distinct transmitting boundary formulations to overcome the problem of radiation damping. These are the superposition boundary approach, which involves the cancellation of the reflected waves by combining the solutions of two different boundary conditions and the viscous boundary approach, which involves the absorption of incident wave energy by frequency independent viscous dashpots. The theoretical justification of these two approaches and their means of implementation are reviewed. The solutions obtained using the two different boundary approaches to the problem of a rigid massless circular footing vibrating on an elastic half-space are compared with an independent theoretical solution. The performance of the boundaries for problems involving step loading is also examined and the implications for any loading pattern with a non-zero time average are discussed. 相似文献
7.
殷文 《吉林大学学报(地球科学版)》2008,38(1):144-151
在弹性波频率空间域有限差分数值模拟方面,差分网格及边界条件是影响弹性波模拟成功与否的关键,为了压制数值模拟中的网格频散,采用25点有限差分算子,建立了有限差分矩阵方程,且借鉴匹配层衰减边界条件思想,设计了弹性波频率空间域有限差分数值模拟算法。由于采用高阶有限差分法来提高差分格式的精度,将会导致计算量显著增加,为此,对频率空间域有限差分弹性波数值模拟方法,采用流水线技术与分治策略进行了并行算法研究,提高了计算效率,使得在合理的计算时间内更精确地模拟弹性波在弹性介质中的传播过程。 相似文献
8.
In this paper, a three-dimensional particle-based technique utilizing the discrete element method (DEM) is proposed to study
wave propagation in a dry granular soil column. Computational simulations were conducted to investigate the soil response
to sinusoidal motions with different amplitudes and frequencies. Three types of soil deposits with different void ratios were
employed in these simulations. Different boundary conditions at the base such as rigid bedrock, elastic bedrock, and infinite
medium were also considered. Analysis is done in time domain while taking into account the effects of soil nonlinear behavior.
The computational approach is able to capture a number of essential characteristics of wave propagation including motion amplification
and resonance. Dynamic soil properties were then extracted from conducted simulations and used to predict the response of
the soil using the widely used equivalent linear method program SHAKE and compare its predictions to DEM results. Generally,
there was a good agreement between SHAKE and DEM results except when the exciting frequency was close to the resonance frequency
of the deposit where significant discrepancy in computed shear strains between SHAKE predictions and DEM results was observed. 相似文献
9.
波动方程的系数反演是一种用于识别地下介质物理力学参数的重要方法。过去这方面的研究一般都建立在弹性模型的基础上,本文则提出了粘弹性反演模型,这种模型更能真实地反映波在地下传播的实际情形。在粘弹性反演模型的基础上,我们还提出了一整套行之有效的数值反演方法,并在频域中完成了反演计算。最后的数值计算结果是令人满意的,证明了本文的模型和方法的合理性。 相似文献
10.
结构-地基动力相互作用是结构地震响应分析及安全评估的一个非常重要课题。基于比例边界有限元法,提出了一种新颖的结构-地基动力相互作用的时域模型,即采用比例边界有限元子域模拟近场有限域部分,采用高阶透射边界模拟远场无限域部分。通过采用连分式展开和引入辅助变量,有限域的动力方程采用高阶的静力刚度矩阵和质量矩阵表示。高阶透射边界精确满足无限远处的辐射边界条件,具有全局精确、时间局部和收敛速度快等优点。它是基于改进的连分式法求解无限域动力刚度矩阵而建立的,在时域里表示为一阶常微分方程组。通过联立有限域和无限域的运动方程,建立了结构-地基相互作用的标准动力学方程,采用Newmark法可直接求解。3个算例结果表明,该算法在时域里比黏弹性边界更精确、有效。 相似文献
11.
12.
The method of fundamental solution for 3‐D wave scattering in a fluid‐saturated poroelastic infinite domain 下载免费PDF全文
下载免费PDF全文 Zhongxian Liu Zhikun Wang Alexander H.D. Cheng Jianwen Liang Chuchu Wang 《国际地质力学数值与分析法杂志》2018,42(15):1866-1889
By using a complete set of poroelastodynamic spherical wave potentials (SWPs) representing a fast compressional wave PI, a slow compressional wave PII, and a shear wave S with 3 vectorial potentials (not all are independent), a solution scheme based on the method of fundamental solution (MFS) is devised to solve 3‐D wave scattering and dynamic stress concentration problems due to inhomogeneous inclusions and cavities embedded in an infinite poroelastic domain. The method is verified by comparing the result with the elastic analytical solution, which is a degenerated case, as well as with poroelastic solution obtained using other numerical methods. The accuracy and stability of the SWP‐MFS are also demonstrated. The displacement, hoop stress, and fluid pore pressure around spherical cavity and poroelastic inclusion with permeable and impermeable boundary are investigated for incident plane PI and SV waves. The scattering characteristics are examined for a range of material properties, such as porosity and shear modulus contrast, over a range of frequency. Compared with other boundary‐based numerical strategy, such as the boundary element method and the indirect boundary integral equation method, the current SWP‐MFS is a meshless method that does not need elements to approximate the geometry and is free from the treatment of singularities. The SWP‐MFS is a highly accurate and efficient solution methodology for wave scattering problems of arbitrary geometry, particularly when a part of the domain extends to infinity. 相似文献
13.
TTI介质弹性波频率-空间域有限差分数值模拟 总被引：3，自引：0，他引：3
由周期性薄互层引起的VTI介质是研究比较广泛的一类各向异性介质。当VTI介质对称轴偏离垂向,本构坐标系与观测坐标系不重合时,会形成观测坐标系下的TTI介质。引入25点优化差分算子,推导出二维TTI介质频率域弹性波动方程;为压制边界反射,采用完全匹配层法吸收边界条件,并计算出优化差分系数;最后采用集中力源,模拟了弹性波在TTI介质中的传播过程。从波场快照和地面共炮记录可以看出,笔者采用的数值模拟算法能有效压制数值频散。TTI介质中的波场传播比较复杂,纵波传播相对稳定,横波波前的三分叉现象比较明显,并存在振幅奇异性。当VTI介质的对称轴偏转后,还会增加地面地震记录的复杂性。 相似文献
14.
流体饱和两相多孔介质动力反应计算分析 总被引：2，自引：1，他引：1
基于流体饱和两相多孔介质的弹性波动方程组,运用显式逐步积分格式与局部透射人工边界相结合的时域显式有限元方法对该波动方程组进行求解,对两相介质在输入地震波作用下的弹性动力反应进行计算和分析;对在是否考虑孔隙流体渗流的两种情况下计算得到的两相介质弹性动力反应结果的差异进行对比研究,从而揭示孔隙流体渗流对两相介质动力反应性质的影响。计算结果表明：两相介质弹性动力反应时程的波形与入射地震波的波形相同,且弹性动力反应的峰值出现的时刻对应于入射地震波的峰值出现的时刻;孔隙流体的渗流将对两相介质的弹性动力反应性质产生显著的影响。数值计算同时表明,时域显式有限元方法是进行流体饱和两相多孔介质弹性动力反应计算分析的一种有效的方法。 相似文献
15.
基于复合材料以及连续介质损伤理论,给出了岩石材料的各向异性损伤破坏模型。通过引入与岩石材料单轴加载行为相对应的特征模态构成的四阶对称损伤张量,描述了岩石材料的损伤演化过程,其中对不同主应变方向采用不同的损伤变量,而对同一主应变方向拉压时的损伤则采用不同的损伤变量来描述。在数值模拟岩石破坏过程的程序中,采用了张量分解的方法。将该模型编写用户材料子程序,并嵌入到大型有限元分析程序ABAQUS中,通过ABAQUS/EXPLICIT SOLVER的显式有限元算法求解。利用此程序对岩石材料的单轴压缩进行了数值模拟。 相似文献
16.
研究了二阶声波方程频域PML边界条件和频域变网格并行计算技术。PML边界是一种较为理想的吸收边界方法,多用在求解时域应力速度方程中,但对于频域声波正演,二阶位移方程更常用。从一阶声波方程PML吸收边界条件导出频域二阶位移方程PML边界条件,模拟算例得到的频率切片、时间切片和地震记录对比都说明该边界条件吸收效果很好。频域单炮正演不同频率间是独立的,据此低频部分采用大网格计算,高频采用小网格,实现变网格步长计算技术,这是较时间域正演的一个优势,在保证模拟质量的同时,减少计算量和内存消耗。 相似文献
17.
公路黄土暗穴临界稳定埋深回归的预测方法及其程序实现 总被引：1，自引：0，他引：1
基于动力学基本方程,运用有限单元法（FEM）和NEWMARK隐式积分方法,对动荷载（FWD荷载）作用下路基下伏黄土暗穴动力稳定的临界埋深进行了弹塑性数值分析,考虑到人工截断边界上波反射对计算结果的影响,分别用传递边界和半无限单元来处理侧面和底面边界。在此基础上,对利用有限个数值分析结果进行不同洞径及位置暗穴临界稳定埋深的预测方法和步骤进行了分析,提出了一种基于洞径和中心距的双变量临界稳定埋深计算方法,并用VB6.0编制了计算程序LDMC。为验证程序的有效性,用3 m和5 m洞径的暗穴进行了对比分析,表明计算结果具有较高的精度,完全可满足实际工程的要求,同时指出了使用时应注意的问题。 相似文献
18.
19.
A coupling scheme for boundary and finite elements using joint elements is proposed which includes the consideration of body forces. In this scheme the boundary and joint elements are formulated in a similar way as finite elements (i.e., the equivalent FE procedure). These joint elements are efficiently used to combine different BE regions. For the evaluation of a body forces, two methods are compared on computational efficiency and it is found that the method using Galerkin tensor is more efficient than the method dividing the problem domain into several internal cells. Two main geotechnical problems considering self weight are numerically examined using this coupling procedure. 相似文献
20.
完全匹配层吸收边界（PML）已经被证明是非常有效的边界吸收技术,对体波和面波的吸收都具有非常好的效果,已经被广泛应用于弹性波的数值模拟中。但是在某些情况下传统的PML技术还是存在一定的问题,比如对掠射情况下的体波和窄区域自由表面条件下的面波的吸收等等。在坐标变换中采用复频移拉伸函数的复频移PML可以有效地改善PML边界条件的吸收性能。基于弹性波一阶速度-应力方程,推导了复频移PML的递推卷积实现方法,并采用交错网格高阶有限差分法对其进行了数值模拟,与传统的PML进行了对比。结果表明：传统的PML对掠射情况下的体波和窄区域自由表面条件下的面波吸收不足,会产生虚假反射,影响真实波场;而基于递推卷积的复频移PML算法能够有效地改善困难情况下的吸收效果,并且在实现过程中不用分裂变量,应用更加方便简单。计算卷积时采用递推的形式,推导过程直观易懂,易于编程,而且不会增加计算量,存储量也没有太大的变化。 相似文献