关于共形映射法求解挠度.docVIP

应用复分析课程设计 ——基于共形映射法求解挠度 基于共形映射法求解挠度 摘要 为了克服传统的分析求解挠度的方法-----几何方法（面矩法，共轭梁法）和能量法（虚功法，卡氏定理），因此我们采用了基于保角映射法求解挠度的方法。此种方法更加精确。首先用解析方法（包括傅里叶展开式、留数定理）求解出保角映射函数，再用此映射函数求得渐开线直齿轮轮齿挠度的解，对影响挠度的因素进行了简要分析。 关键词：保角映射函数； 渐开线直齿轮； 轮齿挠度。 引言 应用平面弹性理论的复变函数保角映射法求解齿轮轮齿挠度，由日本学者首先提出，由试算法得到，计算复杂，映射精度难以保证。后国内学者程乃士等用计算机求得挠度，但求解方法仍属于数值回归法，且求解受到初值影响。本文沿袭盛行的解析法，计算简单映射精度更高。 相关词条的物理解释 应力 当材料在外力作用下不能产生位移时，它的几何形状和尺寸将发生变化，这种形变称为应变（Strain）。材料发生形变时内部产生了大小相等但方向相反的反作用力抵抗外力，结构构件的轴线或中面由于弯曲引起垂直于轴线或中面方向的线位移。 定义为理想材料有小形变时应力与相应的应变之比。泊松比泊松比是材料横向应变与纵向应变的比值，也叫横向变形系数，它是反映材料横向变形的弹性常数。 () 为将Z平面齿轮边界围成的域D映射为平面的下半平面域 的映射函数，因为用复变函数方法解平面弹性力学问题的时候，Airy应力函数的表示为 U=Re[ψ1(z)+θ1(z)] 式中ψ1(z)和θ1(z)是z(z=x+iy, =x-iy)的两个解析函数 因此轮齿位移的复变函数表示为 u+iv=ψ1(z)-[z + ] （1） 式中ψ1(z)= ；E为弹性模量；为泊松比，对于平面应变问题，式中E用E/（1- ）代替，用/（1-）代替。 所以，此时问题集中于求应力函数的两个解析函数和使其满足边界条件。 下面根据集中力作用点分情况进行讨论 当集中力作用于半平面域的边界上，且作用点为原点时 =-lnz，=lnz X=Pcos，Y=Psin 当集中力作用于边界上任意点Z0时，由坐标平移的 =-ln（z-z0）， =ln（z-z0）+* 保角映射 把持性边界映射C映射为虚轴，采用下面形式的映射函数 Z==- 其中适当的选择m0,mj,aj,就能得到精确度很高的齿形。 与的映射关系 ==[()],==[()] 因此得 =-ln[()-(r)] =ln[()-(r)]+ * 因为工程上是没有理论上的集中力的，所以我们需要把和分成非解析部分和及解析部分和两部分 由文献[1]，通过求极限我们有 设A和B到的距离为rA和rB，则 [()+(-)+(-)] + [-]=0 下面的解满足上面方程 (2) 利用在复分析中的Cauchy积分 当函数f(z)在扩充复平面上仅有有限个孤立奇点： 用此分析上式的每一项，故 即 再通过对(2)式取共轭 对上式取Cauchy积分，经过化简得 将上面的式子代入 (1) 得到位移表达式为 (3) 式中u为x方向的位移，v为y方向的位移 （2）为求解保角公式，我们先作如下推导： 图中左半平面相交两线为两啮合齿轮的外围轮廓线。 如图1a所示，渐开线齿轮齿廓有渐开线和过度曲线两部分组成。 设渐开线齿廓上任一点坐标为（Xk ,Yk）为 = (－)－[＋ －－＋] = ±[＋－ －＋] 式中 ——齿顶圆半径 ——基圆半径 ——齿形角 ——齿数 ——变位系数 ——含义如图la所示 ——齿廓上任意点处压力角。在渐开线部分。的取值范围为： ——渐开线与过渡曲线的界限点的压力角 =[－] ——刀具齿顶高系数 ——刀具齿顶圆角半径 齿根过度曲线部分任一点坐标设为（Xk ,Yk） 式中 r－－分度圆半径 －－相应于齿根过渡曲线上计算点的道具齿顶圆角计算点的展开角。 的取值范围为。这里取。 （）、（）为满足边界条件的两个解析函数。由文献[2]中的分析可知：在力的作用点附近， 和 应分解为非解析部分与以及解析部分：与： （