共查询到20条相似文献，搜索用时 203 毫秒
1.
本文将小波有限元法引入到流体饱和多孔隙介质二维波动方程的正演模拟中,以二维Daubechies小波的尺度函数代替多项式函数作为插值函数,构造二维张量积小波单元.引入一类特征函数解决了Daubechies小波没有显式解析表达式所带来的基函数积分值计算问题,并推导出计算分数节点上Daubechies小波函数值的递推公式,从而构造出由小波系数空间到波场位移空间的快速小波变换.数值模拟结果表明该方法是有效的. 相似文献
2.
3.
采用代数重建技术重建三维物体,提出了一种基于splatting原理的ART权值计算方法,以增强重建图像的空间分辨率.该算法基于三维锥束扫描方式,采用双三次样条局部基函数作为插值核,并利用其足印函数对物体进行插值来计算像素的权值.与常量插值相比,这种方法对物体的灰度值有更精确的近似,在相同的实验条件下可重建出分辨率高的物体.给出了适用于三维锥束重建的splatting算法的具体实现,并对实验数据进行了实物重建及误差分析. 相似文献
4.
5.
针对非均匀介质中热蠕变流动问题,给出了有限单元方法与网格－粒子方法联合求解新技术,即有限单元方法求解欧拉网格节点上的未知量,分布于单元内部作为物质成分标记的粒子反映变形过程.有限元法求解动量方程和连续性方程时引入了速度场和压力场等阶插值的压力场稳定的Petrov Galerkin方法,求解能量方程时采用了流线迎风Petrov Galerkin方法,网格－粒子算法中采用双线性插值与有限单元插值函数对应.有限单元计算与网格－粒子计算相对独立,两种方法计算的数据通过有限单元节点传递.同时,实现了三角形单元的算法和程序,解决了复杂结构条件下不规则网格计算的问题.通过经典方腔热对流问题验证了程序,给出了不规则形态块体沉降算例,并分析了数值解的稳定性. 相似文献
6.
基于开放式有限元系统SiPESC. FEMS,构建了机翼流固耦合分析的程序框架。程序框架通过研发流固耦合界面数据交换模块来调用结构和流体分析功能.在流固耦合分析中采用高阶面元法对流体计算以获得飞行器机翼上各点处的压力系数。耦合界面交互采用径向基函数插值算法进行流固两相的载荷和位移的双向插值。进一步通过调用结构分析模块计算结构响应更新结构构型,实现机翼流固耦合整体流程分析。计算过程中各分析模块之间的数据传递均是基于大规模工程数据库SiPESC. ENGDBS来实现的。通过算例测试验证了在SiPESC平台上进行流固耦合分析的可行性。 相似文献
7.
8.
基于可分离小波理论,由一维Daubechies尺度函数的张量积构造二维Daubechies小波基,并将它作为GPR波动方程求解的插值函数,导出了二维Daubechies小波有限元GPR方程离散格式;通过引入转换矩阵,实现小波系数空间与雷达场值之间转换.引入自由度凝聚技术,有效解决了小波有限元求解中小波单元内部自由度过多的问题,节约了计算量并方便与传统有限元法耦合.然后,详细阐述了Daubechies小波有限元联系系数计算方法,有效解决了小波有限元求解偏微分方程的难点与核心问题.最后,以两个典型GPR模型为例,对比了Daubechies小波有限元与传统有限元的雷达正演剖面图与单道波形图,结果表明:在相同的剖分方式及节点数目条件下,Daubechies小波有限元的紧支性与正交性一定程度上提高了求解效率,它与有限元法求解结果能较好地吻合,验证了Daubechies小波有限元算法的正确性. 相似文献
9.
《地球物理学进展》2017,(5)
二维波动方程谱元法通常采用四边形网格,网格质量和属性建模方式对波场模拟精度和效率有较大影响.常规属性建模需要先建立几何模型再对其进行贴体网格剖分.但存在已知稠密属性控制点却难以建立复杂几何模型的问题.针对此问题,提出对模型进行矩形网格剖分并计算GLL点,利用属性控制点对所有GLL点进行插值以提高属性建模的效率.将属性建模方式归纳为单元属性建模方式和节点属性建模方式.在节点属性建模方式中分析了双线性插值、快速Gauss径向基函数插值法的计算效率,将两种插值算法集成于SPECFEM2D的属性建模程序中.实例表明两种属性建模方式对谱元法波场模拟都是有效的. 相似文献
10.
针对大地电磁正则化反演中正则化因子的选取困难问题提出了自适应正则化反演算法（Adaptive Regularized Inversion Algorithm, ARIA）. 在该算法中， ①提出了一种新的数据方差处理方法：数据方差规范化，使得数据方差的大小只对数据的拟合发生影响，不对数据目标函数和模型约束目标函数的权重产生影响，从而减少了正则化因子取值的影响因素；②提出了粗糙度核矩阵的概念，并给出了由基本结构插值基函数计算粗糙度核矩阵的公式，使得模型目标函数的构建更为简便、直接；③根据数据目标函数、模型约束目标函数和正则化因子之间的关系，提出了两种正则化因子自适应调节方法. 本文详细阐述了最平缓模型约束下的大地电磁一维连续介质反演的ARIA实现，以几个算例的分析比较来说明ARIA的有效性. 相似文献
11.
弹性波矢量波场逆时偏移可以综合利用纵横波场信息,对地下空间进行清晰成像,且对于成像介质没有角度限制,可以对复杂构造进行更清晰的成像.而弹性波逆时偏移中最重要的就是求解波动方程的算法,其直接影响成像的精度以及效率.本文引进电力系统谐波分析中常用的余弦组合窗函数,并通过一种新的优化算法得到了改进的余弦组合窗函数从而得到优化后的有限差分算子.并将此算子应用于解耦方程的矢量波场分离算法从而提高了成像精度.数值测试表明基于新算法的逆时偏移的成像精度和清晰度得到了明显的提高. 相似文献
12.
全局弱式无单元法是在有限单元法基础上发展起来的一种数值模拟方法,它采用局部支持域内的节点信息来构造形函数实现局部精确逼近,摆脱了单元,仅依赖于节点信息,具有预处理简单、模拟精度高、灵活性强的特点,适用于复杂地电条件下直流电阻率正演模拟.本文采用RPIM构造直流电阻率全局弱式无单元法形函数,利用RPIM形函数推导了直流电阻率全局弱式无单元法方程.然后,编制了直流电阻率全局弱式无单元法正演模拟Fortran程序,利用该程序对典型的地电模型进行了正演模拟,并将正演结果与基于线性插值的FEM正演结果及解析解进行对比,结果表明采用RPIM形函数的全局弱式无单元法用于直流电阻率正演模拟的正确性及有效性,且在同等条件下,全局弱式无单元法模拟精度高于矩形剖分的FEM,更有利于指导电法勘探的数据解译;利用该程序对复杂地电模型进行了正演模拟,结果表明全局弱式无单元法对复杂地电模型模拟效果良好,适应性强,灵活性高,可任意加密节点提高模拟精度. 相似文献
13.
对已经提出的四种经常使用的包线函数模型进行了评价,认为对于评价单峰点包线函数,除了要求形式简单以外,其优劣可以用函数对峰点位置和高峰隆起程度的控制能力作为标准。借助结构动力试验中阻尼比的确定方法,提出用平稳度指数来评价单峰点包线函数平缓程度。据此对四种包线函数的平稳度进行了分析,还对某些包线函数提出了一些改进。通过分析指出如果把“金井清谱”和“Penzien谱”等模型,看作是时间域中包线,这两个模型的调节能力都很强,而且比较方便,其中“Penzien谱”模型可以用来表示具有两个峰点的包线。文中的比较分析结果可以用于地震动包线函数选择和模拟。 相似文献
14.
A precise value of the matrix-fracture transfer shape factor is essential for modeling fluid flow in fractured porous media by a dual-porosity approach. The slightly compressible fluid shape factor has been widely investigated in the literature. In a recent study, we have developed a transfer function for flow of a compressible fluid using a constant fracture pressure boundary condition [Ranjbar E, Hassanzadeh H, Matrix-fracture transfer shape factor for modeling flow of a compressible fluid in dual-porosity media. Adv Water Res 2011;34(5):627-39. doi:10.1016/j.advwatres.2011.02.012]. However, for a compressible fluid, the consequence of a pressure depletion boundary condition on the shape factor has not been investigated in the previous studies. The main purpose of this paper is, therefore, to investigate the effect of the fracture pressure depletion regime on the shape factor for single-phase flow of a compressible fluid. In the current study, a model for evaluation of the shape factor is derived using solutions of a nonlinear diffusivity equation subject to different pressure depletion regimes. A combination of the heat integral method, the method of moments and Duhamel’s theorem is used to solve this nonlinear equation. The developed solution is validated by fine-grid numerical simulations. The presented model can recover the shape factor of slightly compressible fluids reported in the literature. This study demonstrates that in the case of a single-phase flow of compressible fluid, the shape factor is a function of the imposed boundary condition in the fracture and its variability with time. It is shown that such dependence can be described by an exponentially declining fracture pressure with different decline exponents. These findings improve our understanding of fluid flow in fractured porous media. 相似文献
15.
波形反演是一种利用全波场信息,通过最小化预测波场和实际波场的残差来揭示地下岩性和构造信息的方法.本文首先简述了常规拟牛顿算法的原理,之后利用一种新的拟牛顿公式对Davidon-Fletcher-Powell(DFP)和Broyden-Fletcher-Goldfarb-Shanno(BFGS)算法进行了修正,改进后的BFGS算法在近似Hessian矩阵逆矩阵时,不仅考虑了梯度和模型信息,还加入了目标函数本身的信息,而且对于每次迭代,基本没有增加计算量.数值试验表明,相对常规拟牛顿方法,修正BFGS算法在保证反演精度的同时,明显提高了反演效率. 相似文献
16.
密度界面反演作为了解地球内部结构的一种重要方法,长期以来都是重力学研究的主要内容.本文结合抛物线密度模型及频率域算法的优点,将抛物线密度函数应用于Parker-Oldenburg算法,经过理论推导得到了抛物线密度模型的频率域公式,从而建立了基于抛物线密度模型的三维密度界面重力异常正反演的算法和流程.理论模型数据试验表明本方法快速、有效,适用于大多数浅部比深部增加更快的实际地壳密度.研究中还利用该方法对川滇地区重力异常进行了反演,获得了该区的莫霍面深度分布,并与接收函数研究结果进行对比分析,进一步验证了本文方法的正确性和有效性. 相似文献
17.
A fiber-section model based Timoshenko beam element is proposed in this study that is founded on the nonlinear analysis of frame elements considering axial,flexural,and shear deformations.This model is achieved using a shear-bending interdependent formulation(SBIF).The shape function of the element is derived from the exact solution of the homogeneous form of the equilibrium equation for the Timoshenko deformation hypothesis.The proposed element is free from shear-locking.The sectional fiber model is constituted with a multi-axial plasticity material model,which is used to simulate the coupled shear-axial nonlinear behavior of each fiber.By imposing deformation compatibility conditions among the fibers,the sectional and elemental resisting forces are calculated.Since the SBIF shape functions are interactive with the shear-corrector factor for different shapes of sections,an iterative procedure is introduced in the nonlinear state determination of the proposed Timoshenko element.In addition,the proposed model tackles the geometric nonlinear problem by adopting a corotational coordinate transformation approach.The derivation procedure of the corotational algorithm of the SBIF Timoshenko element for nonlinear geometrical analysis is presented.Numerical examples confirm that the SBIF Timoshenko element with a fiber-section model has the same accuracy and robustness as the flexibility-based formulation.Finally,the SBIF Timoshenko element is extended and demonstratedin a three-dimensional numerical example. 相似文献
18.
19.
Hongyue Guo Xiaodong Liu 《Stochastic Environmental Research and Risk Assessment (SERRA)》2016,30(7):1875-1887
In this study, we propose a new segmentation algorithm to partition univariate and multivariate time series, where fuzzy clustering is realized for the segments formed in this way. The clustering algorithm involves a new objective function, which incorporates an extra variable related to segmentation, while dynamic time warping (DTW) is applied to determine distances between non-equal-length series. As optimizing the introduced objective function is a challenging task, we put forward an effective approach using dynamic programming (DP) algorithm. When calculating the DTW distance, a DP-based method is developed to reduce the computational complexity. In a series of experiments, both synthetic and real-world time series are used to evaluate the performance of the proposed algorithm. The results demonstrate higher effectiveness and advantages of the constructed algorithm when compared with the existing segmentation approaches. 相似文献
20.
系统地论述了用有限单元法研究复杂地形条件下三维直流电阻率的正演计算技术.首先给出了三维构造中点源电场的边值问题以及相应的变分问题;然后利用有限单元法求解变分问题,采用四面体单元对研究区域进行剖分,在单元中进行三线性函数插值,将变分方程化为线性代数方程组;最后,考虑到节约计算时间,利用对称超松弛顸条件共轭梯度迭代算法求解大型线性方程组,得到了各节点的电位值,进而计算出地表的视电阻率.通过理论模型的计算检验了算法的可行性之后,给出了几种常见纯地形异常的数值模拟结果和一个组合模型的计算结果,其研究工作为研究三维直流电阻率反演奠定了基础. 相似文献