圆外弹性问题应力的自然边界元方法解法.docxVIP

圆外弹性问题应力的自然边界元方法解法 挖一条圆形道路后，道路边缘会发生新的位移。由于圆形混合，材料的位错是如图1所示的。这些问题可以追溯到圆外位移的边值问题。假设在内边界Γ上各点存在位移[ux(x,y),uy(x,y)]T|Γ,边界Γ变为Γ′,Ω上应力场必然也随之改变。研究因圆孔边界发生改变,而形成的应力场,在工程应用中有着重要的意义。 在圆外问题中,常规解析方法只能解得边界位移为常数情况下的应力场,对于含有任意边界位移的情况,可借助有限元和边界元等数值方法,但数值方法通常计算量较大,边界位移函数也只能用逼近的方法得到,精度难以保证。本文将自然边界元方法,作为解决该问题的主要方法,自然边界元法应用自然边界规化的思想将微分方程边值问题化为边界上的超奇异积分方程,特别适合求解无界区域,由这种方法可直接得到应力场解析或半解析解,相比于数值方法其具有更高的精度,计算量较小。本文根据文献给出的圆外弹性问题的自然边界元解答,来分析各类位移边界条件下的孔边应力场,在边界位移为常数时其结果与文献得到的解析解一致。 1 圆外调和方程解析函数 以位移表示的弹性平衡方程组为 在平面区域Ω上的任意一组位移解均可表示为复变函数的实部或虚部形式 其中φ(z)和ψ(z)为Ω的两个解析函数。因此对于以半径为R的圆外问题,其边界Γ上的位移为 代入圆外调和方程边值问题Poisson公式, 及自然积分方程 其中,为边界Γ上的外法向导数。 将弹性方程(4)代入到(5)、(6)中,即可得到当Ω为半径为R的圆外区域时,平面弹性问题的Poisson积分公式为 *为卷积符号,其中 自然积分方程为 其中 δ(θ)为Dirac函数。 2 边界上的应力 算例1:半径为R的圆孔,其边界Γ上各点都向圆心方向作径向位移s,求圆外区域Ω的应力分布。 边界Γ上的各点的位移为 将此代入到Poisson积分公式(7)式中,可以得圆外的位移场为 进一步求得应力场分布 利用自然边界积分方程(8)式可以求得边界上的应力为 文献给出了此算例的复变函数解法,两种方法的结果一致。 如果边界Γ上各点都有切向方向位移s,用同样的方法可以得到应力场为 可以看到,在极坐标系下,边界各点有相同法向位移时,切向力为0,两个正应力大小相等;若边界各点有相同切向位移,切向应力与r2成反比,法向应力与r3成反比,法向应力随半径的变化值要大于切向应力。其在x和y方向应力等值线图,分别如图2和图3所示。 对于较为复杂的位移边界条件,一般的解析方法已很难求得解析解,而利用自然边界元方法仍然可以直接快速求得应力场。 算例2:假设边界位移函数为 将此代入到Poisson积分公式(7)式中,可以得圆外的位移场为 由上述位移场可得到相应应力场。 3 自然边界元法 利用自然边界元方法求得任意位移边界条件下的圆外应力场。给出的算例表明,利用自然边界元法可以方便、准确地求得此类问题的解答,这种方法不受复杂位移边界条件的限制,可广泛用于含有复杂位移边界条件的圆外问题求解。