首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 140 毫秒
1.
n维非线性波动方程对应的线性方程Cauchy问题的解   总被引:1,自引:0,他引:1  
在n维非线性波动方程Cauchy问题的迭代解法中 ,每一个迭代步骤都需要求解一个相应的n维线性波动方程的Cauchy问题。对于n维线性波动方程Cauchy问题而言 ,只要初值适当光滑 ,其解也必适当光滑 ,且在整个半空间t≥ 0上是整体存在的。从较简单的n维线性波动方程Cauchy问题出发 ,利用球面平均法和降维法分别求出n为奇数 (n >1)和n为偶数时其解的表达式。并利用叠加原理求得了一般的n维线性齐次方程Cauchy问题和线性非齐次方程Cauchy问题的解 ,为求解非线性波动方程Cauchy问题奠定了基础  相似文献   

2.
位场向上延拓可归纳为求解函数所满足的Laplace方程,将无网格Galerkin(EFG)方法推广到二维位场延拓的数值计算中,详细论述了EFG方法的基本原理和具体实施过程,建立了无网格Galerkin法求解的离散方程,并进行了数值求解。同时与有限差分(FD)法的数值结果进行了比较,两种方法的求解结果基本吻合。介绍了EFG方法的应用实例,理论模型和实例的数值结果表明,EFG方法在处理二维位场延拓问题时是有效的且具有实现简单的特点。  相似文献   

3.
二维分数阶对流-弥散方程的数值解   总被引:1,自引:0,他引:1  
对二维时间分数阶对流-弥散方程和二维空间分数阶对流-弥散方程分别建立了差分格式,实现了对其的数值求解。针对理想算例进行计算求解,分析了时间和空间分数阶阶数取不同值时的扩散变化规律,验证了各自所描述的时间相关性与空间相关性。同时与传统的二维整数阶对流-弥散方程的求解结果作了对比。当时间和空间分数阶阶数α与γ分别取整数时,二维时间分数阶对流-弥散方程和二维空间分数阶对流-弥散方程都与传统二维整数阶对流-弥散方程的计算结果相同,说明提出的对二维分数阶对流-弥散方程的数值求解方法是可行的。其结果对地下水溶质运移的进一步研究提供了有效的手段。  相似文献   

4.
振荡流底层悬沙运动的数值研究   总被引:1,自引:0,他引:1       下载免费PDF全文
建立了平底振荡流底层立面二维水沙数值模型,利用Smagrionsky(SGS)格子涡模型封闭二维Navier Storkes方程水流运动方程组,控制方程采用SMAC法求解。该模型能较精确地模拟振荡流底层水流流动特性,以及含沙量沿垂线分布和随相位变化的情况,且与水槽实验的实测资料基本吻合。  相似文献   

5.
墩柱结构是一种常见的工程建筑物,正确求解作用在墩柱上的波浪力具有重要的意义。主要考虑大尺度墩柱上的波浪力的计算,应用线性水波绕射理论来计算作用在大尺度墩柱上的波浪力。应用线性小振幅波理论把水波绕射问题转化为一个二维的Helmholtz方程来求解,再应用Nyst?m方法来求解二维的Helmholtz方程,通过求得的数值解与解析解的对比,表明给出的数值模型有很高的精度。  相似文献   

6.
墩柱结构是一种常见的工程建筑物,正确求解作用在墩柱上的波浪力具有重要的意义.主要考虑大尺度墩柱上的波浪力的计算,应用线性水波绕射理论来计算作用在大尺度墩柱上的波浪力.应用线性小振幅波理论把水波绕射问题转化为一个二维的Helmholtz方程来求解,再应用Nystm方法来求解二维的Helmholtz方程,通过求得的数值解与解析解的对比,表明给出的数值模型有很高的精度.  相似文献   

7.
AUSM格式在二维浅水方程求解中的应用   总被引:2,自引:1,他引:1       下载免费PDF全文
引入空气动力学中发展起来的AUSM格式,在三角形网格上采用有限体积法,对二维浅水控制方程进行数值求解;引入过程中,对浅水问题中的源项及露滩问题等采用了相应的特殊处理。通过对二维计算常用的溃坝算例及扎龙湿地实际算例的计算结果分析可知:二维浅水问题中引入该计算格式后,计算结果较好且保持了原格式的优良特性,同时也验证了该格式在二维浅水方程中应用的可行性及可靠性。  相似文献   

8.
沈恺伦  王立忠 《岩土力学》2004,25(4):609-613
以Biot(1941)提出的理想土体的三维固结理论为基础,在水流连续条件中增加了考虑孔隙流体本身的压缩变形,并将这一因素以高饱和度土体的饱和度的变化为表征,写入文中建立的平面应变下的Biot二维固结方程。通过实施Laplace变换和Fourier正余弦变换获得基本方程组,求解得到初始函数表达式;再利用数值反变换得到任意时刻、任意一点的数值解,并侧重分析了高饱和土饱和度的变化对孔隙流体的压缩性能的影响,以及由此而产生的固结全过程的变化,表明了考虑孔隙流体的压缩对分析高饱和度土体的固结行为的意义。该解可以进一步推广到其它类型荷载作用在土体不同位置,以及不同边界条件下多层地基的二维固结分析。  相似文献   

9.
利用有限元法实现了任意方向偶极子源在二维介质中频率域电磁响应的数值模拟,研究了波数取值对模拟结果的影响.通过对构造走向的Fourier变换,将全三维电磁问题,转化为一系列二维问题,并在波数域求解,极大地减小了计算工作量,导出了波数域耦合适用于二维电性介质中任何方向电或磁偶极子响应计算的电磁场方程.针对每个给定的波数,上述耦合电磁场方程用等参有限元方法在x-z平面内求解.采用Fourier逆变换,将波数域解积分,得到空间域电磁场.针对电磁模拟计算中,源点的奇异性,采用具有一定面积的伪δ函数表达源电流分布,使数值解精度得以提高.另外,采用等参有限元,使地下复杂地质体得到准确表达.利用不同波数值对均匀介质与层状介质的模拟结果与解析解的对比,验证了算法的正确性与精度.利用层状介质模型的解析解与数值计算结果的对比,分析了波数的优化取值范围及取值点数对数值模拟结果的影响,考察了算法对非均匀介质的适应性.  相似文献   

10.
一种新的二维电阻率模拟方法   总被引:2,自引:0,他引:2  
目前,对地下介质作电阻率成像需要将研究区域网格化。本文以格林函数为数学工具导出了网格化地电模型二维问题的边界积分方程,该方程适合于电阻率成像。用边界单元法求解该方程可以得出空间电位分布,从而进行视电阻率模拟。计算与实验结果的对比表明,这种方法是正确、可行的。  相似文献   

11.
Elastic lateral dynamic impedance functions are defined as the ratio of the lateral dynamic force/moment to the corresponding lateral displacement/rotation at the top ending of a foundation at very small strains. Elastic lateral dynamic impedance functions have a defining influence on the natural frequencies of offshore wind turbines supported on cylindrical shell type foundations, such as suction caissons, bucket foundations, and monopiles. This paper considers the coupled horizontal and rocking vibration of a cylindrical shell type foundation embedded in a fully saturated poroelastic seabed in contact with a seawater half‐space. The formulation of the coupled seawater–shell–seabed vibration problem is simplified by treating the shell as a rigid one. The rigid shell vibration problem is approached by the integral equation method using ring‐load Green's functions for a layered seawater‐seabed half‐space. By considering the boundary conditions at the shell–soil interface, the shell vibration problem is reduced to Fredholm integral equations. Through an analysis of the corresponding Cauchy singular equations, the intrinsic singular characteristics of the problem are rendered explicit. With the singularities incorporated into the solution representation, an effective numerical method involving Gauss–Chebyshev method is developed for the governing Fredholm equations. Selected numerical results for the dynamic contact load distributions, displacements of the shell, and lateral dynamic impedance functions are examined for different shell length–radius ratio, poroelastic materials, and frequencies of excitation. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

12.
A boundary integral equation method is presented for a rigid cylindrical pipe‐pile of finite length embedded in a transversely isotropic half‐space under lateral loads. In the framework of three‐dimensional elastostatics, the complicated soil‐structure interaction problem is shown to be reducible to three coupled Fredholm integral equations. Through an analysis of the associated Cauchy singular kernels, the intrinsic singular characteristics of the radial, angular, and vertical interfacial load transfers are rendered explicit. By means of a complicated numerical procedure, detailed results on the three‐dimensional load–transfer process are provided for benchmark comparison and practical applications. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

13.
In the framework of elastostatics, a mathematical treatment is presented for the boundary value problem of the interaction of a flexible cylindrical pile embedded in a transversely isotropic half‐space under transverse loadings. Taking the pile region as a stiffened subdomain of an extended half‐space, the formulation of the interaction problem is reduced to a Fredholm integral equation of the second kind. The necessary set of Green's functions for the transversely isotropic half‐space is obtained by means of a method of potentials. The resulting Green's functions are incorporated into a numerical procedure for the solution of the integral equation. The theoretical response of the pile is presented in terms of bending moment, displacement and slope profiles for a variety of transversely isotropic materials so that the effect of different anisotropy parameters can be meaningfully discussed. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

14.
The definitions of complex integrals of Cauchy and Hadamard with the singular point coinciding with the end point of the integration curve are proposed. It is shown that the new integrals satisfy most of the properties of the regular ones, including the change of variables. It is also shown that the Cauchy principal value (CPV) and Hadamard finite-part (HFP) integrals can be considered as a sum of the new type integrals. The application to numerical solution by the boundary element method (BEM) and the complex hypersingular integral equation (CHSIE) for the multiregions of interacting elastic bodies and bodies with cracks and holes is discussed. The different ways to place the collocation points are considered. The numerical results for the problems of circular hole and circular elastic inclusion in infinite plate indicated that the appropriate choice of the approximating functions leads to a high accuracy of the calculation. Applications of the new technique to geomechanics problems are discussed. © 1998 John Wiley & Sons, Ltd.  相似文献   

15.
A micropolar elastoplastic model for soils is formulated and a series of finite element analyses are employed to demonstrate the use of a micropolar continuum in overcoming the numerical difficulties encountered in application of finite element method in standard Cauchy–Boltzmann continuum. Three examples of failure analysis involving a deep excavation, shallow foundation, and a retaining wall are presented. In all these cases, it is observed that the length scale introduced in the polar continuum regularizes the incremental boundary value problem and allows the numerical simulation to be continued until a clear collapse mechanism is achieved. The issue of grain size effect is also discussed. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

16.
陈勋 《探矿工程》2010,37(6):01-03
摆线型大位移井轨道由“直井段—圆弧井段—摆线井段—稳斜井段”构成,与其他类型的大位移井轨道相比,摆线型大位移井具有最小的或很小的摩阻和摩阻力矩。摆线型大位移井轨道设计问题归结为求解一个二元方程组,当未知设计参数为摆线井段初始井斜角或稳斜井段井斜角时,该方程组为非线性方程组,没有解析解。通过数学变换和化简,得到了等价的三角函数方程,根据区间搜索和二分法提出了求该三角函数方程近似解的数值迭代算法。算例表明,该算法具有很好的迭代稳定性,计算精度高、速度快。新算法可用于大位移井轨道设计的计算机软件开发,对于提高软件的适用性和稳定性具有较大的帮助。  相似文献   

17.
液化场地桩-土-桥梁结构地震相互作用简化分析方法   总被引:2,自引:0,他引:2  
液化场地桩-土-桥梁结构地震相互作用分析属于桩基桥梁抗震设计中的一个关键科学问题,而目前尚缺乏合理的简化分析方法。鉴于此,直接针对振动台试验,基于Penzien模型,建立了液化场地桩-土-桥梁结构地震相互作用的数值分析模型与相应的简化分析方法。通过振动台试验验证了数值建模途径与简化计算分析方法的正确性,可用于液化场地桩基桥梁结构地震反应的分析,并且特别考虑砂层中孔压升高引起的砂土承载力衰减效应,推荐了计算参数的合理选取方法;据此进行了桩径、桩土初始模量比、砂土内摩擦角、上部桥梁结构质量等重要参数对液化场地桩-土地震相互作用影响的敏感性分析。研究表明:在液化场地条件下,随桩径和桩土初始模量比的增大,桩的峰值加速度、峰值位移减小,而桩的峰值弯矩则增大;随砂土内摩擦角增大,桩的峰值加速度、峰值弯矩、峰值应力均增大,而桩的峰值位移则减小;随上部结构配重增大,桩的峰值位移、峰值弯矩均增大。  相似文献   

18.
This paper presents a direct boundary element method (BEM), formulated in the Laplace transform space, for plane strain poroelasticity. The paper expands on work by Cheng and Liggett by recasting the theoretical foundation of BEM within the framework of Rice and Cleary's formulation of the Biot theory of poroelasticity. Furthermore, the numerical algorithm is generalized to deal with both interior and exterior domain problems, and a method for indirectly calculating the Cauchy principal value of the singular integrals is presented. Formulae for the stress and flux inside the domain are also derived. Finally, the algorithm is validated by comparing the numerical results with the analytical solution of a borehole subject to a far-field deviatoric stress (exterior domain) and with the solution of Mandel's problem (interior domain). These two examples provide a critical test of the algorithm.  相似文献   

19.
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.  相似文献   

20.
在基于波动方程的有限差分数值模拟中,会不可避免地出现数值频散(也称网格发散)问题。数值频散问题通常会给数值模拟的结果造成严重影响,因此在数值模拟中,应尽量设法消除这种现象。这里在前人的基础上,在基于各向同性介质的弹性波方程的数值模拟过程中,通过引入通量校正方法来解决数值频散问题。由数值模拟结果表明,该方法可以有效地消除数值频散现象,大大改善数值模拟的结果。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号