共查询到20条相似文献，搜索用时 109 毫秒
1.
土层场地对地震动的影响明显,研究成层场地局部复杂地形的波动场特性,对于抗震设防具有重要意义。基于土-结构相互作用理论,含有局部复杂地形的成层场地,可分解成广义结构（即近场复杂地形及其周围土体）和具有规则开挖边界的成层无限地基（即远场）。远场的格林函数可通过精细积分求解频域-波数域内的对偶波动方程获得,利用傅里叶逆变换得到频域-空间域内的柔度矩阵。近场可采用比例边界有限元进行模拟,通过高性能连分式的传递边界求解动刚度。自由场波动可通过传递矩阵法求得,将动刚度和自由场波动位移代入控制方程即可得到散射场的动力响应。数值算例验证了方法的准确性,并利用提出的方法讨论含有软夹层场地局部复杂场地的波动场特性。 相似文献
2.
三维黏弹性介质人工边界研究 总被引：1，自引：0，他引：1
为了得到适于分析非线性结构-地基相互作用问题的人工边界,结合弹性波理论与黏弹性理论,推导出了适于模拟三维黏弹性远场介质辐射阻尼的人工边界,并对其模拟远场地基辐射阻尼的性能及相应地震波动输入方法进行了研究。结果表明：在结构-地基动力相互作用问题中,得到的黏弹性人工边界模拟计算域外远场黏弹性介质的辐射阻尼具有理想精度;不论是简谐波还是非简谐波,地震波动输入都取得了理想结果;在黏弹性介质条件下自由地表处的总场位移幅值n是输入波幅值的 倍。 相似文献
3.
目前隧道及大型地下工程往往在裂隙岩体中开挖,而裂隙与地下空间的距离及裂隙的扩展条件,制约着隧道及地下工程的稳定性。应用能考虑孔洞和裂纹问题的新型边界积分方程与无网格加辽金法结合,建立一种新型的边界无单元法。在该方法中基本的未知量是由边界上的面力和边界上位移密度函数构成的复变量边界函数 。文中应用的边界积分公式和Muskhelishvili的积分公式直接相关。将无网格构造方法引入新型的边界积分方程,建立了新型的边界无单元法。应用该方法详细分析了含隧道和裂纹间相互关系等问题,其数值结果与解析结果吻合很好,说明该方法的正确性和可行性。 相似文献
4.
饱和土半解析边界元法及在双排桩被动隔振中的应用 总被引：1，自引：1，他引：0
针对均质饱和地基中双排桩远场被动隔振问题,首先,基于饱和多孔介质的边界元法,建立了以薄层法（TLM）基本解作为动力Green函数的饱和土半解析边界元法,该方法可有效地分析饱和半空间的土-结构动力相互作用问题。在此基础上,根据双排桩与周围土体的邻接条件,推导了双排桩对Rayleigh波散射的三维边界元方程;运用上述边界元方程,对双排桩的远场被动隔振问题进行了研究。结果表明：双排桩能够有效地降低屏障后的位移振幅,其隔振效果要优于单排桩;排间净距对隔振效果影响不大,而桩间净距则对双排桩隔振效果起控制作用。 相似文献
5.
介绍了联合运用积分方程法与有限元法(简称混合法)来计算三维非均匀介质中电磁场分布的理论方法,并进行了数值模拟实验。混合法的原理是引入一个包围非均匀目标体的虚构边界,在边界内部的场用有限元法模拟,在边界外部(包括边界)的场用积分方程表达,二者在边界上通过场的连续性耦合起来。数值实验结果表明,混合法既能显著地减小网格规模,又能灵活地模拟复杂的介质情况,且计算精度较高。 相似文献
6.
本给出了用边界元法积分方程的算法。通过处例验算表明，这种算法综合了边界元法精度高和用矩最法解积分方程速度快的优点，是一种求解积分方程十分有效的算法。 相似文献
7.
本文给出了用边界元法求积分方程的算法。通过算例验算表明，这种算法综合了边界元法精度高和用矩量法解积分方程速度快的优点，是一种求解积分方程十分有效的算法。 相似文献
8.
土体应变变形是土工测试技术的重要难题之一,加速度积分位移是目前动荷载下原位监测和模型试验中认识土体剪应变响应的一种普遍方法,但缺乏验证与可靠性研究。设计砂质边坡?摇臂?质量块振动台模型试验,实测两地震荷载下多组工况的质量块加速度与位移响应时程,探讨了4种代表性方法求解动态、永久位移的可靠性及偏转影响。结果表明,动态位移求解采用高通滤波是一种可靠途径且受积分方式影响较小,各方法在动态位移工况中积分位移与实测相比,峰值平均误差小于等于25%;双拟法积分位移发生严重偏离,与函数形式和处理过程相关。永久位移求解采用高通滤波导致加速度低频成分缺失而无法得到永久位移;双拟法积分位移与实测相比吻合较好,两地震荷载中误差分别为7%和5%;偏转对动态位移求解可靠性影响较小,但对永久位移影响严重,双拟法呈现发散而失真。试验方法与结论,对认识加速度积分位移求取土体剪应变的可靠性和合理选取数据分析方法具有重要科学价值和指导意义。 相似文献
9.
饱和土体中空沟对移动荷载被动隔振的 2.5D边界元法分析 总被引：1，自引：1，他引：0
根据饱和土Biot理论及移动荷载特性,采用积分变换方法,推导了半无限饱和土体频域-波数内的边界积分方程,讨论了2.5D边界积分方程中Cauchy的奇异性,并分析了空沟对移动荷载引起饱和土体振动的被动隔振效果。通过与已知文献结果比较,验证了算法的正确性。数值结果表明,采用2.5D BEM不仅可以解决具有复杂表面的移动荷载作用下动力响应问题,而且可极大减小计算工作量。对于饱和土体受高速移动荷载作用,空沟更能有效减小土体孔压、竖向位移的振幅;当荷载速度较小时,空沟对土体表面的水平向位移几乎不起衰减作用,甚至有放大效应。另外,随空沟距振源的距离增大,空沟能够更好阻隔土体表面波传播。 相似文献
10.
11.
复合单元法在求解有自由面渗流问题中的应用 总被引：3，自引：0，他引：3
论文指出了在有限元分析过程中地下水自由面上、下的区别,提出可采用复合单元法处理常应变交截单元,并以四面体单元为例,推导出了自由面与单元交截的截面方程,详细说明了复合单元法的原理、系数矩阵的求解方法、推导出体积分和曲面积分的求解公式。 相似文献
12.
In finite element calculations the constitutive model plays a key role. The evaluation of the stress response of the constitutive relation for a given strain increment, which is a time integration in the case of models of the rate type, is a typical sub task in such calculations. Adaptive behaviour of the time integration is essential to assure numerical stability and to control the accuracy of the solution. An adaptive second order semi-implicit method is developed in this paper. Its numerical behaviour is compared with an adaptive second order explicit scheme. The two proposed methods control the local error and guarantee numerical stability of the time integration. We include several numerical geotechnical element tests using hypoplasticity with intergranular strain. The element tests simulate the behaviour of a finite element method based on the displacement formulation. 相似文献
13.
14.
In this paper, The Monte Carlo method is incorporated into the finite element method (FEM) to conduct seepage analysis with a free surface. For the transitional element cut by the free surface, it is used to calculate the composite permeability coefficient, as well as to perform the integration directly using Monte Carlo integration. This new algorithm requires less iteration procedures for convergence. The convergence of the method is also proved for cases where there is a significant difference between the permeability coefficient above and below the free surface. 相似文献
15.
Using the extended finite element method, faults can be introduced into a three-dimensional reservoir-geomechanical model without meshing to assess the potential for fault reactivation associated with industrial activities such as disposal of CO2 or fluid extraction. The residuals to the governing equations include basis functions formed from the enrichment functions for strong and weak discontinuities. The traditional spatial integration scheme is based on an algorithm that partitions each enriched element into a collection of subtetrahedra, which is time and storage intensive. To avoid element partitioning, we adopt the homogeneous numerical integration scheme (referred as HNI hereafter) where the integration of homogeneous monomials over each polyhedron is converted into the integration of the same monomials over the one-dimensional edges of the polyhedron by using Stokes's theorem and Euler's homogeneous function theorem. The integrands for the strong and weak discontinuities are derived to implement the integration scheme in the three-dimensional reservoir-geomechanical model. An example of a jagged fault is presented to reveal the advantages of HNI scheme over the standard element partitioning approach. Several other examples that involve fluid flow, fault sliding, and fault sliding triggered by injection pressure are also presented to demonstrate that accurate and efficient computations are realized by the new integration scheme. 相似文献
16.
Seepage-stress coupling is a key problem in the field of geotechnical engineering, and the finite element method is one of the main methods to study seepage-stress coupling in rock masses. However, the finite element method has issues of poor stability, low efficiency, and low accuracy in the simulation of the seepage-stress coupling problem. In this paper, the homogeneous saturated rock mass is taken as the object to deduce the control equation based on the Biot's theory. Considering the singularity of the coupling matrix, the discrete equation is converted into a precise integral format, and the equation is solved by the precise integration method to avoid instability and low precision. The precise integration method in this paper has good numerical stability, fast convergence speed, and high simulation accuracy, which effectively facilitates the rapid and stable numerical simulation of the seepage-stress coupling problem using the equivalent continuum medium model. The validity and accuracy of the precise integration method for seepage-stress coupling problems are verified by numerical examples. 相似文献
17.
18.
An improved search strategy for the critical slip surface using finite element stress fields 总被引：17，自引：0，他引：17
The finite element method can be used to advantage in slope stability problems. This paper proposes a technique to search for the critical slip surface as well as to define and calculate the factor of safety for the slope, when the finite element method is used to model its formation. First, stresses are estimated at each Gaussian point from the finite element analysis. Then, the global stress smoothing method is applied to get a continuous stress field. Based on this stress field, the factor of safety is calculated for a specified slip surface by a stress integration scheme. An improved search strategy is proposed for a noncircular critical surface which starts with a search method for a circular critical surface. During the search process, points defining a trial slip surface can freely move in the finite element mesh subject to some kinematical constraints. This method can be applied to both the limit equilibrium method and the finite element method. Effects of the slope stress history and soil parameters on the resulting critical surface are investigated. 相似文献
19.
The stability of integration is essential to numerical simulations especially when solving nonlinear problems. In this work, a continuum damage mechanics model proposed by the first author is implemented with an integration method named cutting plane algorithm (CPA) to improve the robustness of the simulation. This integration method is one type of return mapping algorithm that bypasses the need for computing the gradients. We compare the current integration method with the previous direct method, and the result shows that the cutting plane algorithm exhibits excellent performance under large loading rate conditions. To enhance accuracy of the new method, a control procedure is utilized in the implementation of the algorithm based on error analysis. Thereafter, the theory of poromechanics is utilized with the damage model to account for the effects of fluid diffusion. Laboratory tests simulated with finite element method illustrate distinct behaviors of shale with different loading rates and indicate the development of microcrack propagation under triaxial compression. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
20.
扩展有限元法是一种在常规有限元框架内求解强和弱不连续问题的新型数值方法,其原理是在裂尖附近用一些奇异函数和沿裂纹面用阶跃函数加强传统有限元的基,以考虑跨过裂纹的位移场的不连续,该加强策略允许计算网格独立于不连续体几何。讨论了扩展有限元法的一些数值方面,主要包括：水平集法确定界面和加强节点与加强方式、裂尖加强范围的选择、J积分区域的确定和积分方案等。 相似文献