有限差分法及其应用
- 格式:docx
- 大小:330.47 KB
- 文档页数:9
电动力学中的电场分布模拟在电动力学中,电场是一个非常重要的概念,用来描述电荷之间的相互作用。
电场的分布对于理解电磁现象以及解决各种工程问题都具有重要的意义。
为了更好地研究和理解电场分布,科学家们发展了各种电场分布的模拟方法。
本文将介绍几种常见的电场分布模拟方法及其应用。
一、有限元法(Finite Element Method,FEM)有限元法是一种常见的数值计算方法,用于求解偏微分方程和变分问题。
在电场分布模拟中,有限元法可以通过将电场区域划分为有限数量的小元素,然后利用这些小元素的基本信息来近似求解电场分布。
有限元法可以应用于各种复杂的电场问题,并且具有较高的计算精度。
二、有限差分法(Finite Difference Method,FDM)有限差分法是一种基于差分运算的数值计算方法,用于求解偏微分方程。
在电场分布模拟中,有限差分法可以将电场区域划分为离散的网格点,然后利用网格点间的差分运算来逼近求解电场分布。
有限差分法适用于各种简单的电场问题,并且计算速度较快。
三、边界元法(Boundary Element Method,BEM)边界元法是一种基于边界积分方程的数值计算方法,用于求解偏微分方程的边界值问题。
在电场分布模拟中,边界元法可以通过将电场区域划分为有限数量的边界元素,然后利用边界元素上的边界条件来求解电场分布。
边界元法适用于具有无穷远边界条件或者具有局部边界条件的电场问题。
四、有限积分法(Finite Integration Technique,FIT)有限积分法是一种基于积分形式的数值计算方法,用于求解偏微分方程的边界值问题。
在电场分布模拟中,有限积分法可以通过在电场区域中离散采样然后应用积分近似来求解电场分布。
有限积分法可以应用于各种电场问题,并且具有适应性强、计算速度快的特点。
五、快速多极子方法(Fast Multipole Method,FMM)快速多极子方法是一种高效的数值计算方法,用于求解大规模的边界值问题。
偏微分方程的数值求解方法偏微分方程是描述自然现象的重要工具,例如描述热传导、电磁波传播、流体运动等。
然而大多数情况下,这些方程很难通过解析方式求解,因此需要数值求解方法。
本文将介绍偏微分方程的数值求解方法及其应用。
一、有限差分法有限差分法是一种常见的偏微分方程数值求解方法。
它将原本连续的区域离散化,将偏微分方程转化为差分方程。
例如对于一维热传导方程:$$\frac{\partial u}{\partial t} = \alpha\frac{\partial^2 u}{\partial x^2} $$其中 $u(x, t)$ 是温度,$\alpha$ 是热扩散系数。
我们可以选择将空间分成 $N$ 个网格,时间分成 $M$ 个步骤。
则有:$$u_i^{m+1} = u_i^m + \frac{\alpha\Delta t}{\Deltax^2}(u_{i+1}^m - 2u_i^m + u_{i-1}^m)$$其中 $u_i^m$ 表示在位置 $i\Delta x$,时间 $m\Delta t$ 时的温度值。
这是一个显式求解方程,可以直接按照时间步骤迭代计算。
不过由于它的误差可能会增长,因此需要小心选择时间步长和空间步长,以保证误差不会过大。
二、有限元法有限元法是一种更加通用的偏微分方程数值求解方法。
它将连续区域离散化成一些小段,称为单元。
然后针对每个单元,将其上的偏微分方程转化为局部插值函数的方程求解。
例如对于一维波动方程:$$\frac{\partial^2 u}{\partial t^2} = c^2 \frac{\partial^2 u}{\partialx^2}$$我们可以选择将空间分成 $N$ 个网格,用有限元方法将每个网格分成若干个单元。
则对于每个单元 $i$,我们可以得到一个局部插值函数 $u^i(x, t)$ 来近似解该单元上的偏微分方程。
这里不再赘述该函数的形式。
另外,我们还需要满足界面上的连续性和斜率匹配条件,以保证整体解是连续的。
课程论文有限差分法在微分方程中的应用本学期学习了《微分方程数值解》,本书中有限差分法给我留下的印象比较深刻,下边说说自己在方面的一点理解,请老师指正。
1.有限差分法的基本思想:当系统的数学模型建立后,我们面对的主要问题就是微分积分方程的求解。
基本思想是用离散的只含有限个未知量的差分方程组去近似地代替连续变量的微分方程和定解条件,并把差分方程组的解作为微分方程定解问题的近似解。
将原方程及边界条件中的微分用差分来近似,对于方程中的积分用求和或及机械求积公式来近似代替,从而把原微分积分方程和边界条件转化成差分方程组。
2.有限差分法求解偏微分方程的步骤:区域离散,即把所给偏微分方程的求解区域细分成由有限个格点组成的网格,这些离散点称作网格的节点;近似替代,即采用有限差分公式替代每一个格点的导数。
逼近求解,换而言之,这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程。
从原则上说,这种方法仍然可以达到任意满意的计算精度。
因为方程的连续数值解可以通过减小独立变量离散取值的间格,或者通过离散点上的函数值进行插值计算来近似得到。
理论上,当网格步长趋近于零时,差分方程组的解应该收敛于精确解,但由于机器字节的限制,网格步长不可能也没有必要取得无限小,那么差分法的收敛性或者说算法的稳定性就显得至关重要。
因此,在运用有限差分法时,除了要保证精度外,还必须要保证其收敛性。
3.构造差分法的几种形式:主要草用的是泰勒级数展开的方法。
其基本差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等。
其中前两种形式为一阶计算精度,后一种为二阶计算精度。
4.有限差分法的应用:4.1抛物线形的差分法中的一维常系数抛物线型方程 考虑最简单的以为常系数抛物线型方程22()u uLu a f x t x∂∂=-=∂∂ (,)x t ∈Ω 其中Ω是(x.t )平面内的给定区域,可以是有节区域或无解区域;a>0是常数,L 是微分算子。
有限差分法在工程数学中的应用研究工程数学是一门研究工程问题的数学学科,它主要应用于解决工程实际问题中的数学模型。
而有限差分法是工程数学中的一种常用数值计算方法,通过将连续问题离散化为离散问题,从而求得问题的近似解。
本文将探讨有限差分法在工程数学中的应用研究。
一、有限差分法的基本原理有限差分法是一种基于差分逼近的数值计算方法,其基本原理是将连续问题离散化为离散问题,通过求解离散问题的近似解来获得原问题的近似解。
具体而言,有限差分法将求解区域划分为若干个小区域,然后在每个小区域内选取一些离散点,通过近似代替微分和积分算子,将原问题转化为一个线性代数方程组或一个差分方程组,进而求解得到近似解。
二、有限差分法在偏微分方程求解中的应用偏微分方程是工程数学中常见的数学模型,它描述了许多实际问题中的变化规律。
有限差分法在偏微分方程的求解中得到了广泛应用。
以二维热传导方程为例,假设一个矩形区域内的温度分布满足热传导方程,可以通过有限差分法将该方程离散化,然后求解离散化后的差分方程组,最终得到温度分布的近似解。
三、有限差分法在结构力学中的应用结构力学是研究结构物受力和变形规律的学科,它在工程领域中具有重要的应用价值。
有限差分法在结构力学中的应用主要体现在求解结构物的静力和动力问题上。
例如,在求解梁的挠度和应力分布时,可以通过有限差分法将梁的微分方程离散化,然后求解离散化后的差分方程组,从而得到梁的近似挠度和应力分布。
四、有限差分法在流体力学中的应用流体力学是研究流体运动规律的学科,它在工程领域中具有广泛的应用。
有限差分法在流体力学中的应用主要体现在求解流体流动的速度场和压力场上。
以二维不可压缩流体的流动为例,可以通过有限差分法将连续方程和动量方程离散化,然后求解离散化后的差分方程组,最终得到流体流动的速度场和压力场的近似解。
五、有限差分法的优缺点及发展趋势有限差分法作为一种常用的数值计算方法,具有一些优点和缺点。
有限差分法推导摘要:一、有限差分法简介1.有限差分法的概念2.有限差分法在数值计算中的应用二、有限差分法的推导1.差分法的定义2.有限差分法的推导过程3.有限差分法的性质三、有限差分法的应用1.微分方程的数值解法2.有限差分法在数值积分中的应用四、有限差分法的优缺点1.优点2.缺点正文:一、有限差分法简介有限差分法是一种数值计算方法,通过将连续函数离散化,用差分代替微分,从而实现对微分方程或积分方程的求解。
有限差分法广泛应用于科学、工程和金融领域,例如,在天气预报、海洋学、生物学、经济学等方面都有重要作用。
二、有限差分法的推导1.差分法的定义差分法是一种将函数在某一点上的值与该点附近点的值相减的方法,用于近似计算函数在该点处的导数或变化率。
给定一个函数f(x),在x=a 处求导,可以得到差分算子Df(a,h),其中h 为差分步长。
2.有限差分法的推导过程有限差分法是将差分法应用于离散点集,通过有限个差分算子来近似表示函数在某一点的值。
设函数f(x) 在区间[x0, x1] 上可导,离散点集为{x0,x0+h, x0+2h, ..., x1},有限差分法的表达式为:Df(x0+k h) ≈ (h/(k+1)) * [f(x0+k h) - f(x0+(k-1) h)] (k=1,2,3,...,n-1)3.有限差分法的性质有限差分法具有以下性质:(1) 线性性质:Df(x) + Dg(x) = D(f(x) + g(x))(2) 移位性质:Df(x+h) = Df(x) + h * df(x)/dx(3) 微分性质:Df(x) * (x - x0) = f"(x) * (x - x0) + O(h^2)三、有限差分法的应用1.微分方程的数值解法有限差分法可以用于求解微分方程,例如,对于一阶线性微分方程:df(x)/dx + p(x) * f(x) = q(x)可以用有限差分法将其离散化为一个线性代数方程组,从而求解离散解。
电磁学的数值计算方法电磁学是研究电场和磁场相互作用的学科,它在日常生活和科学研究中起着重要的作用。
随着计算机技术的快速发展,数值计算方法在电磁学中的应用也越来越广泛。
本文将介绍几种常用的电磁学数值计算方法,并探讨其原理和应用。
一、有限差分法(Finite Difference Method)有限差分法是一种基于离散化空间和时间的数值计算方法,常用于求解求解具有边值条件的偏微分方程。
在电磁学中,有限差分法可以用来求解电磁场的静电场、静磁场以及时变电磁场等问题。
该方法通过将空间和时间进行网格离散化,将偏微分方程转化为差分方程,并用迭代方法求解得到数值解。
二、有限元法(Finite Element Method)有限元法是一种广泛应用于各种物理问题求解的数值计算方法,电磁学也不例外。
该方法通过将求解区域划分为有限的小元素,并在局部内部逼近真实场量的变化。
在电磁学中,有限元法可以用来求解电场、磁场以及电磁波传播等问题。
通过选择合适的元素类型和插值函数,以及建立元素之间的边界条件,可以得到电磁场的数值解。
三、时域积分法(Time Domain Integral Method)时域积分法是一种基于格林函数的数值计算方法,通过积分形式表示电磁场的边界条件和过渡条件,进而求解电磁场。
时域积分法广泛应用于求解电磁波的辐射和散射问题,如天线辐射和散射、电磁波在介质中的传播等。
该方法通过离散化电磁场的源和观测点,并利用格林函数的性质进行数值积分,得到电磁场的数值解。
四、有限时域差分法(Finite-Difference Time-Domain Method)有限时域差分法是一种基于电磁场的离散化网格和时间的有限差分法,是求解各种电磁问题最常用的数值计算方法之一。
有限时域差分法通过离散化时空域,将麦克斯韦方程组转化为差分方程组,并通过时间步进的方式求解得到电磁场的数值解。
该方法适用于求解各种电磁波传播、辐射和散射等问题。
偏微分方程的有限差分法及地球物理应用有限差分法是一种常用的数值求解偏微分方程的方法。
它将连续的偏微分方程转化为离散的差分方程,通过近似求解差分方程,得到偏微分方程的数值解。
这种方法在地球物理学中有着广泛的应用,如地震波传播模拟、电磁场分布计算等领域。
首先,假设我们要研究地震波在地下介质中的传播,可以采用波动方程来描述地震波的传播过程。
波动方程可以写成:∂^2u/∂t^2 = c^2∇^2u其中,u是地震波场,c是地下介质中的波速。
为了用有限差分法求解波动方程,我们需要将连续的空间和时间离散化。
假设我们将空间离散化为网格点(i,j,k),其中i,j,k分别代表空间的x,y,z方向,将时间离散化为时间步长Δt。
对波动方程进行近似,我们可以得到:(u(i,j,k,t+Δt) - 2u(i,j,k,t) + u(i,j,k,t-Δt))/Δt^2 = c^2(u(i+1,j,k,t) + u(i-1,j,k,t) + u(i,j+1,k,t) + u(i,j-1,k,t) +u(i,j,k+1,t) + u(i,j,k-1,t) - 6u(i,j,k,t))/Δx^2将此差分方程应用于地震波传播模拟,我们可以得到地震波场在空间和时间上的离散解。
有限差分法在地球物理中有着广泛的应用。
例如,它可以用于模拟地震波在地下介质中的传播,帮助研究地震灾害的发生机制和地下构造的特征。
通过调整网格的大小和时间步长,可以模拟不同频率的地震波传播过程,从而了解地震波在不同介质中的传播规律。
此外,有限差分法还可以应用于电磁场的计算。
例如,在电磁勘探中,可以利用有限差分法求解麦克斯韦方程,计算电磁场在地下介质中的传播和散射过程。
通过模拟电磁场的分布情况,可以帮助研究地下矿产资源的寻找和勘探。
需要注意的是,有限差分法在应用过程中还需要考虑边界条件的处理。
通常情况下,边界条件是已知的,例如地震波在地表的边界条件可以假设为自由表面,电磁场计算中的边界条件可以假设为电场和磁场的边界条件等。
有限差分法的原理及应用1. 前言有限差分法(Finite Difference Method)是一种常见的数值计算方法,用于求解偏微分方程(Partial Differential Equations,简称PDE)。
它通过在求解域中采用离散点来逼近微分算子,将连续的微分方程转换为离散的代数方程,从而实现对PDE的数值求解。
有限差分法具有简单易懂、易于实现的优点,被广泛应用于科学计算、工程分析等领域。
2. 原理有限差分法的原理基于以下两个基本思想: - 寻找定义域上的离散点,并通过这些离散点来近似表示原方程中的未知函数。
- 使用差分格式来近似微分算子,从而将偏微分方程转化为代数方程组。
具体而言,有限差分法将定义域按照均匀的网格划分为一个个网格点,这些点被称为节点。
同时,有限差分法还使用网格点上的函数值来近似表示原方程中的未知函数。
通过将对原方程中的微商用差商来近似表示,然后将差商带入到原方程中,得到离散的代数方程。
3. 应用有限差分法广泛应用于各个科学领域和工程领域中的数值计算问题。
以下列举几个常见的应用领域:3.1 流体力学在流体力学中,有限差分法被用来模拟流体的运动。
通过将流体领域离散化,将流体的速度、压力等参数表示为离散点上的函数值,可以使用有限差分法求解Navier-Stokes方程,从而得到流体的流动行为。
3.2 热传导有限差分法可以用于求解热传导方程。
通过将传热领域离散化,并将温度表示为离散点上的函数值,可以使用有限差分法求解热传导方程,从而得到材料内的温度分布。
3.3 结构力学有限差分法也被广泛用于求解结构力学中的问题。
例如,在弹性力学中,可以通过将结构域离散化,并将结构的位移、应力等参数表示为离散点上的函数值,使用有限差分法求解相应的弹性方程,从而得到结构的应力分布和变形情况。
3.4 电磁场分析在电磁场分析中,有限差分法被用来求解麦克斯韦方程组。
通过将电磁场的定义域离散化,并将电场、磁场等参数表示为离散点上的函数值,可以使用有限差分法求解麦克斯韦方程组,从而得到电磁场的分布情况。
有限差分法在数值计算中的应用有限差分法是一种常用的数值计算方法,广泛应用于各个领域,包括物理学、工程学、金融学等。
本文将介绍有限差分法的基本原理,以及其在数值计算中的应用。
一、有限差分法的基本原理有限差分法是通过近似计算导数、积分等运算的一种方法,其基本思想是将函数在某一点处展开成一个泰勒级数,然后用有限个点处的函数值来逼近原函数。
有限差分法的核心是将连续的函数转化为离散的数据点,然后通过有限个离散点之间的差分来近似原函数的性质。
有限差分法的主要步骤包括以下几个:1. 网格划分:将计算区域划分为均匀的网格,即将连续的空间划分为一系列离散的点。
2. 逼近函数:将原函数在每个网格点处做泰勒级数展开,得到对应的近似函数。
3. 差分近似:根据泰勒级数展开的结果,利用有限个网格点之间的差分,来近似计算导数、积分等运算。
4. 求解方程:根据差分结果,可以得到离散的代数方程组,通过求解这个方程组得到数值解。
二、1. 偏微分方程求解:有限差分法可以用来求解各种类型的偏微分方程,包括抛物型、椭圆型和双曲型方程。
通过将偏微分方程离散化为代数方程组,再通过求解方程组得到数值解。
2. 数值积分:有限差分法可以用来近似计算函数的积分。
通过将积分区间划分为一系列小区间,并用离散点上的函数值来近似替代原函数,可以得到积分的数值结果。
3. 非线性方程求解:有限差分法也可以用来求解非线性方程。
通过将非线性方程转化为离散的代数方程组,并利用迭代方法求解方程组,可以得到非线性方程的数值解。
4. 边值问题求解:有限差分法可以应用于求解各类边值问题,如求解热传导方程的边值问题、求解电场分布的边值问题等。
通过将边值问题离散化为代数方程组,再通过求解方程组得到边值问题的数值解。
5. 优化问题求解:有限差分法可以用来求解各种类型的优化问题。
通过将优化问题转化为非线性方程组,并利用有限差分法求解方程组,可以得到优化问题的数值解。
总结:有限差分法作为一种常用的数值计算方法,在各个领域中有着广泛的应用。
有限差分法及其应用
1有限差分法简介
有限差分法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。
该方程将解域划分为差分网格,用有限个网络节点代替连续的求解域。
有限差分法通过泰勒级数展开等方法,把控制方程中的导数用网格节点上的函数值得差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
2有限差分法的数学基础
有限差分法的数学基础是用差分代替微分,用差商代替微商而用差商代替微商的意义是用函数在某区域内的平均变化率来代替函数的真是变化率。
而根据泰勒级数展开可以看出,用差商代替微商必然会带来阶段误差,相应的用差分方程代替微分方程也会带来误差,因此,在应用有限差分法进行计算的时候,必须注意差分方程的形式,建立方法及由此产生的误差。
3有限差分解题基本步骤
有限差分法的主要解题步骤如下:
1)建立微分方程
根据问题的性质选择计算区域,建立微分方程式,写出初始条件和边界条件。
2)构建差分格式
首先对求解域进行离散化,确定计算节点,选择网格布局,差分形式和步长;然后以有限差分代替无线微分,以差商代替微商,以差分方程代替微分方程及边界条件。
3)求解差分方程
差分方程通常是一组数量较多的线性代数方程,其求解方法主要包括两种:精确法和近似法。
其中精确法又称直接发,主要包括矩阵法,高斯消元法及主元素消元法等;近似法又称间接法,以迭代法为主,主要包括直接迭代法,间接迭代法以及超松弛迭代法。
4)精度分析和检验
对所得到的数值进行精度与收敛性分析和检验。
4商用有限差分软件简介
商用有限差分软件主要包括FLAC、UDEC/3DEC和PFC程序,其中,FLAC是一个基于显式有限差分法的连续介质程序,主要用来进行土质、岩石和其他材料的三维结构受力特性模拟和塑性流动分析;UDEC/3DEC是针对岩体不连续问题开发,用于模拟非连续介质在静,动态载荷作用下的反应;PFC是利用显式差分算法和离散元理论开发的微、细观力学程序,它是从介质的基本粒子结构的角度考虑介质的基本力学特性,并认为给定介质在不同应力条件下的基本特征主要取决于粒子之间接粗状态的变化,适用于研究粒状集合体的破裂和破裂发展问题,以及颗粒的流动(大位移)问题。
5有限差分法的应用
1)光子晶体光纤色散的有限差分法研究
光子晶体光纤的横截面结构和特征参数如图1所示。
对结构如此复杂的光子晶体光纤,采用有限差分法加以研究。
有限差分法是光波导分析中广为采用的数值方法之一,它利用台劳展开用有限差分式代替波动方程中的微分式,得到关于场分量的有限差分方程。
在纵向(z方向)均匀分布的光子晶体光纤中,电场的横向分量Ex和Ey 满足如下耦合方程组
其中,Ex和Ey为电场的两个偏振分量,εr为相对介电常数,k0为真空中的波数,β=koneff为传播常数,neff为模式有效折射率。
在光子晶体光纤分析中,通常求解基模(类比于传统光纤,也称
为HE11模)的两个偏振分量。
文献中已经证明了横向结构对称分布的光子晶体光纤基模的线偏振特性,因而(1),(2)式中的电场偏振分量Ex和Ey可以认为不再耦合,此时(1),(2)式将变为两个独立的分别关于两个偏振分量Ex和Ey的半矢量形式的波动方程
(3),(4)式中保留左边的第二项对于研究空气填充率较高的光子晶体光纤是必要的,因为这时标量波动方程给出的解会产生较大的误差。
在x,y方向采用相同的网格间距h,利用五点差分格式可以建立基于半矢量波动方程(3),(4)的差分方程。
其中E为Ex或Ey,p,q为网格点编号。
对于Ex,各项系数分别为
对于Ey,各项系数分别为
(5)式可化为特征值方程
其中A为系数矩阵,E为由各网格点构成的向量。
计算结果
为将基于半矢量波动方程的差分法得到的结果与文献中全矢量方法得到的结果比较,取与文献中相同的求解参数:光纤参数为空气孔间距Λ=2.3μm,空气孔直径d=1.0μm,整个求解区域取为6Λ×6Λ,石英的折射率为1.45,波长为1.5μm,采用电壁(electric wall)边界条件,x,y方向格点数均为240。
基于半矢量波动方程的差分法得到结果为neff=1.42806,与文献中给出的结果1.42868以及其他方法得到的结果较为接近。
考虑到全矢量方法计算结果表明在求解某一偏振态时另一偏振分量并不为零,但相差几个数量级,而半矢量波动方程完全忽略了另一偏振分量的影响,以及网格剖分的方法不同,这一结果是较为精确的。
2)抗滑桩全桩内力计算“m-k”法的有限差分法
如图1所示,设抗滑桩全长为H,其中滑动而以上,即受荷段桩长为h1;滑动而以下,即锚固段桩长为h2。
设滑坡推力ET按梯形分布,地而和滑动而处的分布荷载集度分别为q0和qa。
桩前滑体的剩余抗滑力为E'。
桩的截而尺寸:长x宽为axb,桩的抗弯刚度为EI。
为便于公式的推导,设桩前滑体地基系数为K1 (z )=mZ,滑而以下锚固体地基系数为K2 (z ) = K(常数)。
由于抗滑桩有刚性桩和弹性桩(或称柔性桩)之分,它们的受力特点是不同的,因而应注意区分,即应先按下式计算出桩的变形系数:
式中b0为桩的计算宽度的单位为m^-1 当h2<=1时,为刚性桩; h2>=1时,为弹性桩。
如图1所示,抗滑桩在滑坡推力和土体抗力作用下产生弹性挠曲变形,其受荷段.A)的挠曲微分方程为
采用等量分段h由桩顶往下至滑动面将桩长离散化,如图2(a)(b)所示。
则任意节点i处的控制差分方程为
即
由于桩顶自由,故在桩顶处弯矩和剪力均为0,即由
得
及由
得
令式(3)中i=0,并将式(4)和式(5)代入,可得
令
则式(4)和(6)可以表达为
仿此,可得
将(8)(9)代入(3)整理可得
上式中令
则式(10)仍可写成式(7)的形式。
滑动面以下,即锚固段可视为桩顶受水平荷载(集中荷载———剪力和弯矩)的侧荷桩,其挠曲微分方程为
从桩底由下往上至滑动面也采用等量分段h见图2(c),则可得该段任意节点j 处的控制差分方程为
即
桩底按自由端考虑,由弯矩和剪力均为0的边界条件,同理可得到式(4)和式(5)。
同样地式(12)中令j=0,可得
将式(4)和(5)代入上式,整理可得
令
则
防上2式可得
将(14)(15)代入(12)整理可得
由于抗滑桩在滑动面处应满足位移、转角、弯矩和剪力的连续条件。
参照图2(b),(c),并仿照式(7)~(9)和式(13)~(15),可得到下面的线性方程组
解此线性方程组,可求得滑动面处及附近节点(上、下段各5 个)处的位移。
从而采用迭代法可求得桩身各节点处的位移及其内力。
计算结果及准确度检验
计算结果几乎完全一致,说明本文方法是正确可靠的。
传统的查表手算法在桩底处出现了与边界条件不完全吻合的残余误差,这是由于传统方法所用表中系数截断误差所致,而差分法则没有出现这样的问题。
6结论
可以看出在实际问题中利用有限差分法对实际问题分析求解,再利用相应计算机程序如FLAC、UDEC/3DEC和PFC程序,有时也可以根据自己需要根据算法编制一个计算程序求出方程组的解从而得出的结果与查表结果和实际结果对比,精确对很高,从而证明了其正确性与实用性。
7参考文献:
1《计算机在材料科学与工程中的应用》-中南大学出版社、2 《抗滑桩全桩内力计算“m-k”法的有限差分法》戴自航,彭振斌3《光子晶体光纤色散的有限差分法研究》
栗岩锋,刘博文,王子涵,胡明列,王专,王清月(天津大学精密仪器与光电子工程学院光电信息技术科学教育部重点实验室,天津300072)。