径向辛体系一种新的差分格式.pdf

径向辛体系一种新的差分格式1 1 2 1 任谊滨 ，朱炳麒 ，卓家寿 1 河海大学土木工程学院，江苏南京（210098 ） 2 河海大学机电工程学院，江苏常州（213022 ） E-mail：renyibin@ 摘 要：通过对Hellinger-Reissner 变分原理进行坐标变换，将径向模拟为时间，导向辛体系， 得到Hamilton 对偶方程组。将微分形式的有限差分法引入弹性力学极坐标系下径向辛体系， 把对偶方程组中的微分方程直接改用差分方程代替，推导出极坐标系下问题的辛差分方程， 从而得到一种全新的径向辛体系差分格式。求解方程组，可直接得到位移和应力。编程计算 曲梁等算例，结果表明该辛差分格式是有效的，丰富了弹性力学辛体系差分法的内容。 关键词：弹性力学，辛体系，有限差分法，新差分格式 中图分类号：CN11-5484/N 1. 引 言 极坐标系下典型的区域是扇形域，通过分别将径向及环向模拟为时间坐标，引入对偶变 量，由弹性力学的变分原理可将极坐标问题导向两种不同形式的辛体系，从而给出了圆形及 环扇形域平面弹性问题的一个解析求解方法[1,2]，在环向辛体系下，求得了极坐标弹性力学 [3] 问题的一个新解，利用这个新解可求解一类有实际意义的弹性力学问题 ，在直角坐标辛体 系下，建立了辛差分格式，可以求解具有应力边界的直角坐标平面问题[4]，还有将积分形式 的差分法和哈密顿体系应用到极坐标下，求解了一类问题[5] 。本文是将哈密顿体系与微分形 式的有限差分法应用到极坐标下，把径向模拟为时间，建立了径向辛体系一种新的差分格式， 并编程对算例进行计算，结果令人满意，丰富了弹性力学辛体系差分法的内容。 2. 极坐标哈密顿对偶体系 如图1 所示的顶角为2 α 的弹性扇型域。在极坐标（ρ ，ϕ ）中，记uρ ，uϕ 分别为径 向和环向位移分量，而σ ，σ ，τ 为相应的应力分量。 ρ ϕ ρϕ 下面扼要介绍弹性平面极坐标下的哈密顿体系，有关详细内容及推导过程参见文献 [1,2]。从Hellinger-Reissner 变分原理出发，引入变换ξ ln ρ ，即ρ exp(ξ) ，ξ ln R ， 1 1 ξ ln R ，且引入变量S ρσ ，S ρτ ，S ρσ ，使扇型域变成如图2 所示的 2 2 ρ ρ ρϕ ρϕ ϕ ϕ 矩形区域，并将径向ξ模拟为时间，记(i) ∂/ ∂ξ ，可将原问题导入哈密顿体系。此时的 变分方程为： α ξ2 ⎧ ∂u ∂u 1 ⎫ ϕ ρ 2 2 2 （1） δ + + + − − − + − ν + +ν ξ ϕ ⎨S u S (u ) S (u u ) [S S 2 S S 2(1 )]S ⎬d d 0 ∫∫ ρ ρ ϕ ρ ρϕ ϕ ϕ ρ ϕ ρ ϕ ρϕ −α ξ1 ⎩ ∂ϕ ∂ϕ 2E ⎭ 1本课题得到国家自然科学基金（项目名称：Hamilton 体系下新的数值方法，项目编号：2002504343 ）的资 助。 - 1 -