保角变换和曲线坐标
- 格式:doc
- 大小:244.00 KB
- 文档页数:10
一、基础知识 1 定义在自变量域我们对同一个点从两个方向趋近,这两个趋近方向的夹角与在因变量上趋近的方向夹角一致,称为保角变换 2泊松方程与拉普拉斯方程对于泊松方程:20ρϕε∇=(在静电场中,可以表示电势与电荷的分布关系) 同时在没有电荷分布的地方满足拉普拉斯方程:20ϕ∇=3将在原来复杂的区域上的表达式通过一个变换,折射到宁一个区域上,使得某一分布函数得到简化变换的条件是泊松方程与拉普拉斯方程仍然成立22222x y∂∂∇=+∂∂,同时,我们定义x 、y 为ξ、η的函数:(,)x ξη、(,)y ξη 则x x x ξηξη∂∂∂∂∂=+∂∂∂∂∂2222222()x x x x x x x x x x ξηξξηηξηξξηη∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂=+=+++ ∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂ 其中:222x x x x x ξηξηξξξηξξηξ ∂∂∂∂∂∂∂∂∂∂∂∂=⋅+=+⋅∂∂∂∂∂∂∂∂∂∂∂∂∂ 同理:222x x x x x ηξηξηηηξηηξη∂∂∂∂∂∂∂∂∂∂∂∂=⋅+=+⋅ ∂∂∂∂∂∂∂∂∂∂∂∂∂ 所以:222222222222x x x x x x x ξηξηξηξηξηξη∂∂∂∂∂∂∂∂∂∂∂∂ =++++∂∂∂∂∂∂∂∂∂∂∂∂∂ 同理:222222222222y y y y y y y ξηξηξηξηξηξη∂∂∂∂∂∂∂∂∂∂∂∂=++++ ∂∂∂∂∂∂∂∂∂∂∂∂∂ 所以拉普拉斯方程变换为:22222222222222222222222x y x y x y xy xy x y y x ξξηηξξηηξηξηξηξηξη ∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂ ∇=+=+++++++++∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂要满足保角变换,其实部与虚部都需要满足拉普拉斯方程:20ξ∇=、20η∇= 将实部与虚部要满足的拉普拉斯方程代入上式:2222222222222x y x y x y ξξηηξη∂∂∂∂∂∂∂∂ ∇=+=+++∂∂∂∂∂∂∂∂ ()'f z ix xξη∂∂=+∂∂(对于趋近方向为:0,0x y ∆→∆=) 222222"()f z x x x y y x ξηξξηη ∂∂∂∂∂∂=+=+=+ ∂∂∂∂∂∂将其代入:22222222'()'()'f z f z ξη ∂∂∇=+=∇∂∂也就是说,原坐标下的拉普拉斯方程与泊松方程变换为:220'0ϕϕ∇=⇒∇=222001''()f z ρρϕϕεε∇=⇒∇= 那么对于一个线段,在原坐标系下长度为1,其在新的坐标系下长度为'()f z 二、常用的保角变换1. 线性变换f az b =+,显然'f a =,其几何效果如下:线性变换一般不单独使用:仅对原来的二位分布做了位似2.幂和根式n xn f z = i n in z Ae f A e ϕϕ=⇒=用来处理过原点的射线,原来的射线的长度ρ的取值范围为(0,+∞),求幂或根还是(0,+∞)将原来的自变量求幂次积,几何效果如下:假设有变换3f z =,其效果为:将原来的60°夹角变为180°,并且其中的点的分布也随之扩大角度,假设原来的函数为电势分布函数,求p 点的电势,则通过变换之后,在新的复平面得到了一个平行分布的电势图,设新的电势分布图中,边界上的电势为V 0,则空间中的电势分布为0u V C η=+⋅,其中,C 为常数,C 与介质表面的面密度σ相关,其正负与σ的正负相反我们在新的复平面中求出电势的表达式之后,再求逆变换得到在原来的复平面上的电势表达式:0u V C η=+⋅中,由原来的变换:()()()32332322333(3)(3)i x iy x x iy x iy iy x xy i x y y ξη+=+=+⋅++=−+−由实部对实部,虚部对虚部,得:233x y y η=− 将η代入电势表达式中:()2303u V C x y y =+⋅−得到电势关于x 、y 的表达式同理可以得到将原来的复平面上的表达式开根得到将原来的夹角缩小相应的倍数的变换方法3. 指对数变换(一)、对于指数函数:()z x iy x iy f e e e e +===⋅此处需要注意,这里使用了复变函数的幅角表示法,即:i z Ae ϕ=,所以此处的x e 为幅值,iy e 为幅角其几何空间意义如下: (1),复平面中平行于实轴的直线,其变换后的图像为过原点的射线对于原空间有一条平行于实轴的直线((,)y const x ∈−∞+∞,),原来的x 的值为幅角,y 的值为幅值。
§ 8.7 保角变换和曲线坐标学习思路:弹性力学问题的求解有赖于边界条件的简化。
对于复杂的边界形状,如果利用空间的变换,将是简化问题求解的最好途径。
保角变换就是充分发挥复变函数的特长,将孔口问题映射到平面的单位圆。
这一节将介绍保角变换和曲线坐标的概念。
由于应用保角变换,矢量-位移,张量-应力公式以及K-M 函数等均必须做出曲线坐标描述。
保角变换使得问题的公式复杂,但是边界条件的简化,以及柯西积分的应用将简化问题的分析。
在本节学习之前,请你先学习附录2,(有关保角变换的知识)学习要点:1.保角变换和曲线坐标;2.矢量的保角变换;3.位移分量的曲线坐标表达式;4.应力分量的曲线坐标表达式。
为了便于根据边界条件确定 K-M 函数,采取保角变换z = ()将物体在z 平面上所占的区域变为在平面所占的区域。
一般的说,通过保角变换可以将非圆边界映射为圆边界,使得问题得以简化。
假设将z 平面上的有限区域或者无限区域S 映射为平面的单位圆内的区域,并且将z 平面上的区域S的边界l 映射为单位圆,对应的关系如下表:由于平面上的任一点可以表示为,。
和是点的极坐标。
而根据保角变换公式z = ( ),则z平面任意一点也可以通过和表示。
因此,和又称为曲线坐标。
对于某些问题的描述中,采用曲线坐标形式表示位移和应力有利于问题的分析。
曲线坐标的概念:平面的一个圆周=const 和一条径向直线=const分别对应于z 平面的两条曲线,这两条曲线就记作=const和=const。
于是和可以看作z 平面上一点的曲线坐标。
由于变换的保角性,这个曲线坐标总是正交的,而且坐标轴和的相对位置和坐标轴Ox和Oy的相对位置相同,如图所示。
首先讨论矢量的保角变换。
设曲线坐标,即=const 与x 轴夹角,如果A 为z 平面上的任一矢量,设A 与曲线坐标夹角。
设A x, A y 分别表示矢量A 在x,y 轴的投影;A , A 表示在=const 和=const 上的投影,则上式的几何意义为,将矢量A 绕z 点顺时针方向转动角后,其在Oxy 坐标系的位置,相当于A 在曲线坐标系( , )中的位置,如图所示。
1 应用原理及特点在矿场水力压裂中,如何针对有效渗透率和厚 度不等的特定储层,设计出缝长和导流能力的优化 方案, 是应考虑的首要问题之一。
另外需要一种计算裂缝井产能的简易方法。
应用保角变换方法研究压裂井产能,其原理及特点是:①能将 z 平面上特别复杂的渗流问题转化为平面上一相对简单和易于求解的渗流问题;② 可准确地描述井筒附近较为复杂的流动型态( 裂缝 内流动和非裂缝区域拟径向流动) 对压裂后产能的贡献,而且能对不同导流能力造成的复杂流线型态 统一转化,因而具有广泛的适应性;③经过保角变换后假设的缝端封闭边界条件更符合实际,因保角变换后, 裂缝端部位于主流线上。
以此为基础,应用质量守衡定律和达西运动方程,推导出了裂缝内原油 流动所满足的压力二阶微分方程, 并进行了产量的 求解,与现有的典型曲线对比,一致性程度较好。
2 数学模型2、1模拟的假设条件 模拟的假设条件是: ①垂直裂缝 , 且对称分布于油井的两边; ②假设裂缝剖面为矩形, 高度恒定, 并等于油层厚度 ; ③裂缝宽度相对油藏的供给半径来 说非常小,即在进行保角变换时可忽略不记; ④裂缝 内导流能力可以是有限导流, 也可以是无限导流; ⑤油藏及裂缝内为单相流动,且符合达西线性定律; ⑥稳态渗流,且不考虑地层的垂向流动; ⑦不考虑地层和裂缝内的污染。
2、2模型 的建立在 z 平面上建立 一 Y 坐标系,保角变换转化为平面 r — s 坐标系( 图1 )图一 保角变换示意图取保角变换为:chw L z f =2ww e e chw -+=式中:z 为Z 平面上的复变函数,i y x z +=,f L 为裂缝半长,m;w 为变换后的W 平面,''i y x w +=。
裂缝井的渗流问题从而演变为带状地层向中心 线A 的单向渗流问题。
由于对称性 , 只研究 平 面中图示阴影部分的单向渗流问题。
其中'O 为''B A 的中点 , 即2''π=A O 。
保角变换能计算力摘要:1.保角变换的定义和作用2.保角变换在计算力中的应用3.保角变换在实际问题的应用案例4.保角变换的局限性和发展前景正文:保角变换是一种数学变换,它在数学、物理等领域具有广泛的应用。
保角变换能够保持角度不变,仅改变长度和面积的比值。
在计算力方面,保角变换能够简化复杂的计算问题,提高计算效率。
保角变换在计算力中的应用主要体现在以下几个方面:1.坐标变换:在平面上,保角变换可以将一个复杂的图形变换到一个简单的图形,从而降低问题的复杂度。
例如,将极坐标变换为直角坐标,可以简化计算过程。
2.微积分:在求解微分方程、积分等问题时,保角变换可以将复杂的问题转化为简单的三角函数问题。
例如,在研究波动方程时,利用保角变换可以将空间坐标变换为复数坐标,从而简化问题的求解。
3.数值计算:在数值计算中,保角变换可以提高计算精度和稳定性。
例如,在求解非线性方程时,采用保角变换可以将方程变为易于求解的形式。
4.信号处理:在信号处理领域,保角变换被广泛应用于信号分析、滤波和信号重建。
例如,傅里叶变换和拉普拉斯变换就是两种常见的保角变换方法,它们能够将时域信号转换为频域信号,从而方便对信号进行分析和处理。
在实际问题中,保角变换的应用案例众多。
例如,在地震勘探、无线通信、图像处理等领域,保角变换技术都发挥着重要作用。
然而,保角变换也存在一定的局限性,如在处理奇异值问题时,保角变换可能失效。
因此,在未来发展中,我们需要不断探索新的变换方法,以应对更为复杂的问题。
总之,保角变换在计算力领域具有重要的应用价值。
通过简化复杂问题、提高计算效率,保角变换为科学研究和实际工程带来极大的便利。
[平衡微分方程的适用范围]平衡微分方程平衡微分方程的适用范围弹性力学、塑性力学、弹塑性力学。
张量:怎样判断?商判则:和任意矢量点积为K-1阶张量的量一定为K 阶张量。
能否满足分量转换规律是判断某个数的集合是否表示一个张量的基本准则。
3、n 维张量的举例标量零阶张量,矢量为一阶张量,应力、应变为二阶张量,应力、应变之间的弹性关系可用四阶张量表示。
4、▽的意义?▽为一个梯度,▽2为调和算子,▽4为重调和算子。
5、柯西应变张量与格林应变张量的区别?柯西应变张量适用于线弹性小变形,格林应变张量适用于任何情况。
6、任意斜面上的应力的本质是?平衡微分方程和转轴公式。
7、如何描述正应变,剪应变,体积应变,应力的球张量,应力的偏张量?对于各向同性材料,正应力引起正应变,引起线元长度变化;剪应力引起剪应变,引起角度的变化;应力的球张量,只引起体积变化,不会引起形状的变化;应力的偏张量,只引起形状变化,不会引起体积的变化。
动力学的平衡微分方程如何表示?根据达朗贝尔原理,把惯性力当作体力来满足力平衡和力矩平衡条件。
9、转轴公式的理论依据:柯西公式。
10、等效应力、等效应变物理意义、公式:等效应力将6个应力分量的对变形体的作用,等效于一个单向拉伸力的作用;等效应变将6个应变分量等效于一个单向拉伸力所产生的应变。
利用实验,就可以直接建立等效应变与等效应力的数值关系11、体积不可压:从体积弹性模量来看,当时,K 趋向于无穷大,也就是说体积变化无限小,即表示体积不可压缩。
12、为什么等值拉压是纯剪切等值拉压时,线元只有角度发生变化,长度有发生变化,故等值拉压是纯剪切。
13、里茨和伽辽金法的物理思想均是利用利用最小势能原理,寻找满足约束边界条件的试验函数。
14、弹性力学为什么可用逆解法、半逆解法:解的唯一性定理表明,无论用什么方法求得的解,只要能满足全部基本方程和边界条件,就一定是问题的真解。
15、叠加原理建立在什么条件下:基本方程和边界条件满足线弹性条件,举例:在线弹性条件下,复杂问题可通过简单叠加处理。
§3.3 保角变换通过保角变换,把物理平面上的复杂几何域映射成像平面上的单位圆、半无限平面等简单规则域;同时把物理平面上的基本关系也用像平面上的复变量表示。
先在像平面的规则域上寻找满足这些基本关系的解,然后把结果返回物理平面就得到实际问题的解。
这种保角变换技术在下面介绍的级数展开法,柯西积分法以及解析延拓法中均能采用。
3.3.1 保角变换与曲线坐标 采用保角变换()ζω=z ,把弹性体在z 平面上所占的区域变换为ζ平面上的区域。
数学家已经进行了大量的研究,各种相应区域的保角变换解析函数)(ζω可从保角变换手册中查到。
在ζ平面上令θρθθρζi e i =+=)sin (cos , (3-17)式中ρ和θ是ζ点的极坐标(不是z 点的极坐标)。
ζ平面上的一个圆周const.ρ=和一根径向线const.θ=分别对应于z 平面上的一根曲线。
这两根曲线也就可以用const.ρ=和const.θ=来表示,如图3-3所示。
于是,ρ和θ是z 平面上一点的曲线坐标。
由于变换的保角性,这两组曲线总是正交的,相应的切线ρ和θ叫曲线坐标轴,它们的相对方向与坐标轴x 和y 相同。
设z 平面上有一个矢量F ,它的起点在()()i z e θωζωρ==。
F x 及F y 为这矢量在x 及y 轴上的投影,ρF 及θF 是它在ρ及θ轴上的投影。
设ρ轴与x 轴成角λ,则由几何关系有cos sin ,sin cos x y F F F F F F ρθρθλλλλ=-=+.于是可得()i x y F iF F iF e λρθ+=+即()i x y F iF F iF e λρθ-+=+ (1)为了求得λi e -,设想沿ρ轴方向给z 点以位移d z ,因而对应点ζ得径向位移d ζ,且d d , d d i i ze z e λθζζ==。
故()()()d ()d ()()d d i i ze e z λθωζζωζζωζρωζζωζωζ'''===='''⋅. (2)上式两边取共轭,得i e λ-,于是(1)式变为()()()x y F iF F iF ρθζωζρωζ'+=+' (3)3.3.2 保角变换后的位移与应力公式首先把其中z 的函数变换为ζ的函数。
§1 复变函数的定义由两个实数x,y确定的数z=x+i y称为复数。
x,y分别称为复数z的实部和虚部,记作x=Re z 和y =Im z。
Re和Im分别为表示复数实部和虚部的符号。
其中称为虚数单位。
显然z可以用直角坐标系(x,y)表示,x称为实轴,y称为虚轴。
坐标平面称为复平面,或者z平面。
因此,z平面上的任一点可记作称为复数z的模,称为z的幅角,其在[0,2 ]之间的值称为主幅角。
显然,复数可以写作极坐标表达形式。
设有一个复数z=x+i y的集合g。
对于集合g中的每一个复数z都有对应的复数值,w=u+i v,则称w是z的复变函数,记作w = f (z)。
给定一个复变函数就是在点(x,y)与(u,v)之间给出了一一对应关系。
因此,u,v均随x,y而确定,这就是说给定了一个复变函数和给定两个实变函数u=u(x,y),v=v(x,y)是等价的。
而且w=u(x,y)+i v(x,y)复变函数和实变函数同样有单值函数和多值函数,应该注意到实变函数的性质对于复变函数可能是不成立的。
例如复变函数中的对数函数w=ln z是多值的。
为了便于理解,以对数函数为例。
设。
上式对于z的所有不等于零的复数值定义了函数ln z。
在公式中包含一个任意的整数k,这就是说ln z是一个多值函数。
对于k的任一整数值,就有函数ln z的一个分支。
通常取k=0的那一支叫做的主值,即如果z的一个值对应着w的一个值,那么函数f(z)是单值函数;如果z的一个值对应着两个或两个以上的w值,则f(z)是多值函数。
集合g称为f(z)的定义集合。
§2 解析函数--复变函数的可导性复变函数的导数与实变函数的导数定义是相同的。
因此,关于实变函数的一系列微分公式与法则,可以完全照搬到复变函数上。
不过应该注意的是,复变函数的变量是复变量,不是实变量。
值得指出的是,实变函数的可导性要求当x=x0+∆x 由左右两方趋近x0时,∆y/∆x的极限都存在而且相等。
§8.7 保角变换和曲线坐标
学习思路:
弹性力学问题的求解有赖于边界条件的简化。
对于复杂的边界形状,如果利用空间的变换,将是简化问题求解的最好途径。
保角变换就是充分发挥复变函数的特长,将孔口问题映射到ξ 平面的单位圆。
这一节将介绍保角变换和曲线坐标的概念。
由于应用保角变换,矢量-位移,张量-应力公式以及K-M函数等均必须做出曲线坐标描述。
保角变换使得问题的公式复杂,但是边界条件的简化,以及柯西积分的应用将简化问题的分析。
在本节学习之前,请你先学习附录2,(有关保角变换的知识)
学习要点:
1. 保角变换和曲线坐标;
2. 矢量的保角变换;
3. 位移分量的曲线坐标表达式;
4. 应力分量的曲线坐标表达式。
为了便于根据边界条件确定K-M函数,采取保角变换
z = ω (ξ)
将物体在z平面上所占的区域变为在ξ平面所占的区域。
一般的说,通过保角变换可以将非圆边界映射为圆边界,使得问题得以简化。
假设将z平面上的有限区域或者无限区域S映射为ξ平面的单位圆内的区域∑,并且将z平面上的区域S的边界l 映射为单位圆γ,对应的关系如下表:
由于ξ 平面上的任一点可以表示为,。
ρ和ϕ是点ξ 的极坐标。
而根据保角变换公式z = ω (ξ),则z平面任意一点也可以通过ρ和ϕ表示。
因此,ρ 和ϕ 又称为曲线坐标。
对于某些问题的描述中,采用曲线坐标形式表示位移和应力有利于问题的分析。
曲线坐标的概念:ξ平面的一个圆周ρ =const和一条径向直线ϕ =const分别对应于z平面的两条曲线,这两条曲线就记作ρ =const和ϕ =const。
于是ρ和ϕ可以看作z平面上一点的曲线坐标。
由于变换的保角性,这个曲线坐标总是正交的,而且坐标轴ρ 和ϕ 的相对位置和坐标轴Ox和Oy的相对位置相同,如图所示。
首先讨论矢量的保角变换。
设曲线坐标ρ,即ϕ =const与x轴夹α角,如果A 为z平面上的任一矢量,设A与曲线坐标ρ 夹β角。
设A x, A y分别表示矢量A 在x,y轴的投影;Aρ ,Aϕ 表示在ρ=const和ϕ =const上的投影,则
上式的几何意义为,将矢量A绕z点顺时针方向转动α角后,其在Oxy坐标系的位置,相当于A在曲线坐标系(ρ,ϕ)中的位置,如图所示。
如果用uρ , uϕ 分别表示曲线坐标下的位移矢量分量,则
根据保角变换,有
所以
沿曲线(ρ)取微分线段d z,则在ξ平面对应的有dξ,由于
所以,取其共轭可得。
将上式回代到公式,可得
下面通过保角变换对弹性力学的公式作对应的转换。
首先,设K-M函数和ψ (z)分别使用和ψ 1(z)代替,同时令
根据位移表达式,有
在z 平面上,将位移矢量向曲线坐标ρ和ϕ投影。
由公式
可得
上式两边同时乘以2G,可得
上式是ξ平面上的曲线坐标系表达的位移表达式。
下面建立曲线坐标中应力分量的复变函数表达式。
如果用σρ, σϕ , τρϕ表示物体在曲线坐标中的应力分量。
则
因为和,而由公式
所以
上式为经过保角变换后,z平面上的曲线坐标系中的应力分量的复变函数表达式。
§8.8 无限大薄板的孔口问题
学习思路:
本节的主要任务是将保角变换用于无限大薄板的孔口问题,确定K-M 函数的基本求解公式。
推导中首先确定无限大板孔口问题的保角变换公式,将K-M 函数转换为曲线坐标形式。
采用的方法仍然是将K-M 函数分解为以级数表达的解析函数和对数表达的多值函数两部份。
对于K-M 函数的级数形式,通过孔口面力边界条件可以确定级数函数的求解方程。
这个求解过程,利用保角变换后孔口边界的特殊性质,使用柯西积分使得计算简化。
学习要点:
1. 保角变换公式与K-M 函数;
2. 利用孔口边界条件确定K-M 函数求解公式;
3. 柯西积分确定K-M 函数的级数形式。
保角变换的目标是:将z平面上的孔口边界l映射为ξ 平面上的单位圆γ,将l 以外的无限区域S 映射为ξ 平面上的单位圆内的有限区域∑,将z平面上的无穷远点映射为ξ平面的坐标原点,如图所示。
保角变换公
式:
是将l
以外的无限区域映射为单位圆γ 内(|ξ|<1)的普遍变换式,公式中R为实数,C k为复数,而且<1。
保角变换公式确定以后,可以确定K-M函数和ψ(ξ),即将K-M函数和ψ1(z)变换到曲线坐标中去。
其
中
因
为
由于<1 ,将上式展开,有
所以,ln z = ln ξ +单位圆内部ξ的解析函数。
另
外。
根据上述分析,的各项都转变为单位圆内ξ 的单值解析函数。
因此
其
中,。
讨论边界条件确定K-M 函数和ψ 0(ξ)。
根据面力边界条件
,经过保角变换后,可得
在单位圆的圆周上,。
所以上述面力边界条件可以表示为
根据公式
,则在边界即单位圆周上
将上述K-M函数的边界值回代面力边界条件,并且将已知函数与需要确定的未知函数分开,可得
其中已知函数为
因为和ψ0(ξ)是单位圆内的泰勒级数,它们是从z平面上l R之外无穷区域的罗伦级数转化而来的。
因此对于公式
幂级数求解时,由于方程两边都含有σk=e i kϕ的各个项(k由-∞到∞),比较各个同类项的系数,即可求得a k,b k的
值。
不过这样作太麻烦了,由于和ψ0(ξ )在单位圆内是解析的,而且在
圆内和圆周上是连续的,因此可以直接采用柯西积分计算。
将边界条件的第一式两边乘以,积分可得
由于在单位圆内是解析的,因此公式的第一个积分即等于
,它是级数之和。
对于公式第三项的积分函数,由于
在单位圆外是解析的,在圆外和圆周上是连续的,所以。
因此,边界条件的第一式就成为
同理,边界条件的第二式成为
上述公式就是边界条件通过柯西积分所推导出的计算和ψ0(ξ)表达式。
其中是边界的已知函数。