叶轮强度计算中二次算法的改进

叶轮强度计算中二次算法的改进

武汉化工学院 胡家顺

摘 要 本文根据离心机械中叶轮轮盘应力计算的基本方程式，导出了用迭代法对叶轮轮盘进行应

力计算的迭代公式和求解方法。该方法克服了二次计算法的不足，特别适合编程用计算机计算求解。 关键词：叶轮 强度 计算 迭代

1 引言

工程上对离心式机械中具有复杂形状叶轮轮盘的强度计算，通常都是将其分为若干段，使每一段成为简单的等厚段或锥形段。于是只要已知轮盘内孔处应力时，即可从第一段开始采用求解等厚或锥形轮盘的方法逐段计算出各个截面的应力，从而得到整个轮盘的应力分布。

轮盘内孔处的径向应力可根据设计要求由轮盘的过盈配合选取，故可认为是已知的，但是切向应力实际上无法事先知道。这样也就无法用求解等厚或锥形轮盘的公式进行应力计算。为解决这个问题，工程中通常采用二次计算法求解如下形式的方程方程组：

T

T

c to t ro r t c to t ro r r βσβσβσασασασ++=++= （1）

在求解过程中，所要用到的系数,,,,r c t r βααα…均已绘成大量计算曲线或制成表格，即“尺寸系数”。这

种二次计算法，为求解复杂形状轮盘的应力提供了卓有成效效的手段。但在实际使用过程中需多次作出假设反复对各截面应力进行计算，过程十分繁复。为此本文给出了一种计算叶轮二向应力的迭代法，克服了上述二次计算法的不足。只要给出一定的边界条件，就很容易编程在计算机上计算出叶轮轮盘各截面上的应力。

2 递推公式及求解方程

若将轮盘沿半径方向分割为k-1个截面，则可从轮盘的内孔处向外周方向递推计算各截面处的应力。对于旋转的轮盘，外径处的径向应力ra σ为定值（对于汽轮机轮盘，外径处装有叶片，由叶片离心力确定；对于离心式压缩机或离心泵轮盘，0=ra σ），而内孔处的径向应力0r σ为已知值（可由设计要求根据装配和驱动的情况确定），切向应力0t σ可假设为Q 。

于是可得初始条件：Q t r r ===0t101 σσσσ

当i=1时，

1

111111211

11111121V Q Z F E D Y Q X C B A t r t ti t r r ri +=++==+=++==++σσσσσσσσ

式中， 11111111111111,,,,,T F E D T C B A c t r t r βββααα====== 111111111111 , , ,F D V E Z C A Y B X r r +==+==σσ

当i=2时，

222121212122

2222312

22121212122

222231)(= )(= V Q Z F V E Y D Q Z E X D F E D Y Q X C V B Y A Q Z B X A C B A t r t i t t r r ri +=++++++==+=++++++==++σσσσσσσσ

当i=k-1时，

1

11212121211

11111

11212121211

1111)( )( --------------------------------+=++++=++=+=++++=++=k k k k k k k k k k k k tk k rk k k t k k k k k k k k k k k k rk k tk k rk V Q Z F V E Y D Q Z E X D F E D Y Q X C Y B V A Q X B Z A C B A σσσσσσ

从以上递推过程可得：

i i i i i i i i i i i i i i i i i i i i i i F V E Y D V Z E X D Z C V B Y

A Y Z

B X A X ++=+=++=+=--------11111111

,

, （i=2,3,…，k-1） （2）

因此， σσri i i t i i i X Q Y Z Q V ++=+=+11

（i=2,3,…k-1） （3）

其中Q 由外径处（即i=k-1处）的边界条件确定。

: , = , 1

1

11-----=

+=k k ra ra k k ra rk X Y Q Y Q X σσσσ得即于是由 （4）

对于汽轮机轮盘，外径处装有叶片，ra α由叶片离心力确定；对于离心式压缩机或离心泵轮盘，0=ra σ 将所得Q 代入式（3）便可求得不同截面上的应力值ti ri σσ,。

3 计算方法与步骤

3.1根据轮盘的实际形状将轮盘沿半径方向取k-1个截面，将其分为等厚、锥形等若干段，得各段i i y D ,，

据此计算（或查图）得各截面的强度计算系数)(,,),(,,'

'ci ci ti ri ci ci ti ri ββββαααα；

3.2 依式（2）并由 A B F i i i ，，…，计算得X Y Z V i i i i ，，，； 3.3 根据式（4）计算Q ；

3.4 将Q 代入式（3）求得不同半径截面上的应力σσri ti ，。

对于截面突变的情况，如图1所示可将截面突变处视为两个截面，两截面

的应力ti ri σσ、有如下关系[1]

：

)

(111

1ri ri ti ti ri i i

ri y y σσμσσσσ-+==

++++

为便于迭代，将上式改写为： 图1

i

i i ti i ri i ti ri i ri ri i ti ri ri ti ti i

i i ti i ri i ri i ri V Q Z F E D m m Y Q X C B A m +=++=+-=-+=-+=+=++==+++σσσσμσσμσσσμσσσσσσ)1( )()(111

式中，0,1),1(,0,0,,1

==-=====

+i i i i i i i i i i

i F E m D C B m A y y m μ，i i i t V Z Y X ,,,同式（2）。

4计算示例

为与二次计算法进行比较，现以文[1]叶轮强度二次计算法实例中所提供的参数进行计算。已知某离心式压缩机叶轮的基本参数为：单圆弧叶型叶片曲率R K =959.09mm,叶片圆弧所在圆半径R C =718.15mm ，叶片进口角β130A =?出口角β247A =?，叶片厚度δ=45.mm ，叶片进口宽度b 176=mm ，出口宽度b a =39mm ，

叶片（Z 型）折边宽度221=?mm ，5.272=?mm ，叶片数Z=22，叶轮转速n=6000rpm,轮盘材料为35CrNi3Mo 。

轮盘和轴的配合过盈量所引起的轮盘内孔处的应力16.10-=ro σMpa 。试计算轮盘的应力。 解：将轮盘截面划分为6段，除第Ⅰ段为等厚段外其余各段均为锥形段，图2所示。各段截面的数据

列于表1。第Ⅰ段和第Ⅱ～Ⅵ段分别按等厚和锥形段计算（或查图）得强度计算系数，另外按文[1]的处理方法考虑Ⅴ、Ⅵ段叶片离心力的影响，分别计算：

i

轮盘各截面处的叶片安装角βi βi K C i K i

R R R R R =-+-cos 1

222

2

轮盘各截面处的叶片宽度b i b b R R tg tg b R R i a i i a =+-=

-+()()

1142θ

θ

轮盘的截面处叶片截面积f i f b i i =++()??12δ 轮盘各截面附加厚度y i ' y f z

D i i i i

'.

sin =05πβ

表1

图 2

轮盘各截面处相对附加厚度S i

S y y i i i

='

轮盘各截面处系数p q i i ,

p S R S R R R i i i i i

i i =

--+++111

q S S R R i i i

i i

=--++11

折合重度 γγxi i i i p q R =++()12

由此确定该两段的强度计算系数'

',ci

ci βα，然后按上述计算方法与步骤可求得轮盘各截面处的应力，如表2所示，其迭代计算结果与二次计算法的结果基本一致，个别地方的些许出入，仅仅是由于两种计算方法在确定强度计算系数时所产生的误差而引起的。

5结 语

本文所给出的算法求解叶轮强度，只需根据已知的边界条件，通过简单的有限次迭代计算即可求得叶轮的应力．克服了二次算法过于繁复的不足，具有运算过程简单明了，能很好适应不同截面形状叶轮的应力计算等特点，特别适合用计算机编程计算．

参 考 文 献

[1] 潘永密 李斯特主编，化工机器（下册），化学工业出版社，北京，1981年

[2] 西安交通大学透平压缩机教研室编著，离心式压缩机强度，机械工业出版社，北京1980年

a a

5

54

43

322

11

04

15.4129.3544

72.2169169

ⅥⅤ

Ⅳ

Ⅲ

ⅡⅠ

111.5

计算电场强度的基本方法

计算电场强度的基本方法 电场强度是静电学中最基本最重要的概念之一，是历年高考考查的热点。高考中将静电学与力学、磁学等问题放在一起作为综合题考查在每年是必不可少的。这些题目中往往涉及有电场力、电势和电势能等参数，这些参数与静电场最基本的物理性质参数——电场强度是紧密相关的。因此，要解决好这些问题，我们首先必须熟练掌握计算电场强度的方法。 在这里，我们首先介绍一下计算电场强度的基本方法。结合所分析的静电场的特点，很多求解电场强度的问题都可以用它来解决。对于一些比较特殊的电场，我们将在下一节介绍一些特殊的方法，那些特殊的方法也是由这些基本方法衍生而来的，因此，我们需要掌握好这些基本方法。下面来看一看这些基本方法。 方法特点 电场强度的定义是检验电荷在电场中某点受到的电场力F 与电荷q 的比值，用E 表示。因此，我们可以利用这一定义去求电场中某点的电场强度。想办法求出电荷q 在某点所受的电场力，使用公式F q E =，即可求出电场强度。在这里需要注意两点：（1）这里q 代表 电量，如果带正电则值为正，此时E 的方向与F 相同；如果带负电则值为负，此时E 的方向与F 相反。（2）由于E 有方向，是矢量，因此我们可以使用矢量的运算法则（正交分解法、平行四边形法则、矢量三角形法则等）求几个不同的电场在某一点所产生的合场强。 根据这一定义，点电荷Q 在周围某点所产生的场强为22 Qq F r q k Q E k q r ===。根据这一定义以及匀强电场中电场力做功与电势能的关系有W F d qE d q U === ，因此匀强电场的场强为U d E =。 从定义引出来的方法是最基本的方法，下面我们来看一看具体该怎么用。 经典体验（1） 如图所示，带正电小球质量为m=1×10-2kg ，带电量为q=1.6×10-6 C 。置于光滑绝缘水平面上的A 点，当空间存在着斜向上的匀强电场时，该小 球从静止开始始终沿水平面做匀加速直线 运动，当运动到B 点时，测得速度v B =1.5m/s ， 此时小球的位移为s=0.15m ，求此匀强电场 的场强E 的取值范围（g=10m/s 2 ）。 体验思路： 要求E 的取值范围，我们已知电量q ，根据上面的定义，即是要求电场力的

§10-怎样计算电场强度

§10 怎样计算电场强度？ 静电场的电场强度计算，一般有三种方法： 1、 从点电荷场强公式出发进行叠加； 2、 用高斯定理求解； 3、 从电场强度和电势的微分关系求解。 这三种方法各有优点： 从点电荷的场强公式出发，通过叠加原理来计算，在原则上，是没有不可应用的。但是，叠加是矢量的叠加，因此计算往往十分麻烦。 用高斯定理求电场强度，方法简单，演算方便，它有较大的局限性，只适宜于某些电荷对称分布的场强的计算，或者场强不是对称的，但为几种能用高斯定理求解折场的合成。 用场电势的微分关系求场强也有普遍性，而且叠加是代数叠加。这一种方法也简便，不过还比不上高斯定理。 所以求场强时，一般首先考虑是琐能用高斯定理，其次考虑是否能用场强与电势的微分关系去求。下面分别加以讨论。 一、从点电荷的场强公式出发通过叠加原理进行计算 点电荷的场强公式： 301 (1)4i i i q E r r πε= ∑r r 当电荷连续分布时： ()() 303 0301(2) 4134144r E dl r r E ds r r E d r λπεσπερτπε===???r r r r r r 式中 λ－电荷的线密度； σ－电荷的面密度； ρ－电荷的体密度。 式（2）、（3）、（4）中，积分应普遍一切有电荷分布的地方。计算时，还必须注意这是矢量和。 1、 善于积分变量的统一问题

如果积分上包含有几个相关的变量，只有将它们用同一变量来表示，积分才能积得结果。 这在应用点电荷的场强公式求带电体的场强时，或者应用毕－沙－拉定律求B r 时，常常遇到。 因此，要积分必须先解决积分变量的统一问题。 积分上包含有几个变量，相互之间存在一定的关系。因此，任一变量都可选作自变量，而将其他变量用该变量来统一表示。必须指出，不但可以将积分号中包含的变量选作自变量，而且也可选择不包含在积分号中但与积分号中的变量都有关的量作为自变量，要根据具体情况而定。 现以图2－10－1所示均匀带电直线的场强计算为例来讨论积分变量的统一问题。 由图可知： 2 0cos 4x dl dE r λθπε= 2 0sin 4y dl dE r λθπε= 202 0cos (5) 4sin (6) 4x x y y dl E dE r dl E dE r λθπελθπε∴====?? ?? 上述三个变量中，共有三个相关变量：θ、l 、r 。为了把积分计算出来，必须把三个变量统一用某一个变量，可以θ、l 、r 中的任一个，或者用它的相关变量来表示。究竟选哪 一个好呢？ 如果选择θ为自变量，则应把l 、r 都化作θ的函数来表示。由图示几何关系可得： 2222cot l a dl acse d r a cse θθθθ =-== 于是得： ()()2 12 1 21002100cos sin sin 44sin cos cos 44x y E a a E a a θθθθλλ θθθπεπελλ θθθπεπε==-==-? ? x 图2－10－1

MAAB计算方法迭代法牛顿法二分法实验报告

姓名 实验报告成绩 评语： 指导教师（签名） 年 月 日 说明：指导教师评分后，实验报告交院（系）办公室保存。 实验一 方程求根 一、 实验目的 用各种方法求任意实函数方程0)(=x f 在自变量区间[a ，b]上，或某一点附近的实根。并比较方法的优劣。 二、 实验原理 (1)、二分法 对方程0)(=x f 在[a ，b]内求根。将所给区间二分，在分点 2a b x -=判断是否0)(=x f ；若是，则有根2a b x -=。否则，继续判断是否0)()(

+)(0x f 0))(('0=-x x x f 设0)('0≠x f ,则=x -0x )(') (00x f x f 。取x 作为原方程新的近似根1x ，然后将1x 作为0x 代入上式。迭代公式为：=+1 k x -0x )(')(k k x f x f 。 三、 实验设备：MATLAB 软件 四、 结果预测 （1）11x = （2）5x = （3）2x =0,09052 五、 实验内容 （1）、在区间[0,1]上用二分法求方程0210=-+x e x 的近似根，要求误差不超 过3105.0-?。 （2）、取初值00=x ，用迭代公式=+1 k x -0x )(') (k k x f x f ，求方程0210=-+x e x 的近似根。要求误差不超过3105.0-?。 （3）、取初值00=x ，用牛顿迭代法求方程0210=-+x e x 的近似根。要求误差 不超过3105.0-?。 六、 实验步骤与实验程序 （1） 二分法 第一步：在MATLAB 软件，建立一个实现二分法的MATLAB 函数文件如下： function x=agui_bisect(fname,a,b,e) %fname 为函数名，a,b 为区间端点，e 为精度 fa=feval(fname,a); %把a 端点代入函数，求fa fb=feval(fname,b); %把b 端点代入函数，求fb if fa*fb>0 error('两端函数值为同号'); end

电场强度地计算

电场力的性质之考点一（电场强度的理解及计算） 班级：：编写：熠 学习目标：1、理解电场强度的矢量性；2、掌握电场强度的计算方法。 自主学习：一、三个公式的比较 二、 (1)电场叠加：多个电荷在空间某处产生的电场的电场强度为各电荷在该处所产生的电场场强的矢量和． (2)计算法则：平行四边形定则． 题型一、点电荷产生的电场 正点电荷电场方向背离电荷负点电荷电场方向指向电荷中心 1、如图所示，真空中有两个点电荷Q1 =+3.0×10-8C和Q2 =-3.0×10-8C，它们相距0.1m，A点与两个点电荷的距离r相等，r=0.1m 。求：电场中A点的场强。 2、如图，A、B两点放有均带电量为+2×10-8C两个点电荷，相距60cm，试求：

（1）AB 连线中点O 的场强； （2）AB 连线的垂直平分线上离开O 点距离为30cm 处的P 点的场强。 合作学习： 【拓展训练】：3、(2013·重点中学联考)如图所示，一个均匀的带电圆环， 带电荷量为＋Q ，半径为R ，放在绝缘水平桌面上．圆心为O 点，过O 点作一竖直线，在此线上取一点A ，使A 到O 点的距离为d 。求A 点处的电场强度。 方法归纳： 【变式训练】：4、在某平面上有一个半径为r 的绝缘带电圆环： (1)若在圆周上等间距地分布n (n ≥2)个相同的点电荷，则圆心处的合场强为多少？ (2)若有一半径同样为r ，单位长度带电荷量为q (q ＞0)的均匀带电圆环上有一个很小的缺口Δl (且Δl r )，如图所示，则圆心处的场强又为多少？ 方法归纳：补偿法。 解题关键：把带有缺口的带电圆环―――→转化为 点电荷 解析： (1)当n 分别取2、3、4时圆心处的场强均为零，结合点电荷电场的对称性可知，n 个相同的点电荷在圆心处的合场强为零． (2)可以把均匀带电圆环视为由很多点电荷组成，若将缺口补上，再根据电荷分布的对称性可得，圆心O 处的合场强为零，由于有缺口的存在，圆心O 处的电场即为缺口相对圆心O 的对称点产生的电场，其电场强度为该处电荷(可视为点电荷)在O 点的电场强度(包括 大小和方向)．其电场强度的大小为E ＝k q Δl r 2，方向由圆心O 指向缺口． 答案： (1)合场强为零 (2) k q Δl r 2，方向由圆心O 指向缺口 分析电场叠加问题的一般步骤 电场强度是矢量，叠加时应遵从平行四边形定则，分析电场的叠加问题的一般步骤是： (1)确定分析计算的空间位置； (2)分析该处有几个分电场，先计算出各个分电场在该点的电场强度的大小和方向； (3)依次利用平行四边形定则求出矢量和． 题型二特殊带电体产生的电场

直齿圆柱齿轮强度计算

4.5 直齿圆柱齿轮强度计算 一、轮齿的失效 齿轮传动就装置形式来说，有开式、半开式及闭式之分；就使用情况来说有低速、高速及轻载、重载之别；就齿轮材料的性能及热处理工艺的不同，轮齿有较脆（如经整体淬火、齿面硬度较高的钢齿轮或铸铁齿轮）或较韧（如经调质、常化的优质钢材及合金钢齿轮），齿面有较硬（轮齿工作面的硬度大于350HBS或38HRC，并称为硬齿面齿轮）或较软（轮齿工作面的硬度小于或等于350HBS或38HRC，并称为软齿面齿轮）的差别等。由于上述条件的不同，齿轮传动也就出现了不同的失效形式。一般地说，齿轮传动的失效主要是轮齿的失效，而轮齿的失效形式又是多种多样的，这里只就较为常见的轮齿折断和工作面磨损、点蚀,胶合及塑性变形等略作介绍，其余的轮齿失效形式请参看有关标准。至于齿轮的其它部分（如齿圈、轮辐、轮毂等），除了对齿轮的质量大小需加严格限制外，通常只需按经验设计，所定的尺寸对强度及刚度均较富裕，实践中也极少失效。 轮齿折断

轮齿折断有多种形式，在正常情况下，主要是齿根弯曲疲劳折断，因为在轮齿受载时，齿根处产生的弯曲应力最大，再加上齿根过渡部分的截面突变及加工刀痕等引起的应力集中作用，当轮齿重复受载后，齿根处就会产生疲劳裂纹，并逐步扩展，致使轮齿疲劳折断（见图1 图2 图3）。此外，在轮齿受到突然过载时，也可能出现过载折断或剪断；在轮齿受到严重磨损后齿厚过分减薄时，也会在正常载荷作用下发生折断。在斜齿圆柱齿轮（简称斜齿轮）传动中，轮齿工作面上的接触线为一斜线（参看），轮齿受载后，如有载荷集中时，就会发生局部折断。 若制造或安装不良或轴的弯曲变形过大，轮齿局部受载过大时，即使是直齿圆柱齿轮（简称直齿轮），也会发生局部折断。 为了提高齿轮的抗折断能力，可采取下列措施：1）用增加齿根过渡圆角半径及消除加工刀痕的方法来减小齿根应力集中；2）增大轴及支承的刚性，使轮齿接触线上受载较为均匀；3）采用合适的热处理方法使齿芯材料具有足够的韧性；4）采用喷丸、滚压等工艺措施对齿根表层进行强化处理。 齿面磨损 在齿轮传动中，齿面随着工作条件的不同会出现不同的磨损形式。例如当啮合齿面间落入磨料性物质（如砂粒、铁屑等）时，齿面即被逐渐磨损而至报废。这种磨损称为磨粒磨损（见图4、图5、图6）。它

电场强度的几种计算方法

电场强度的几种求法 一．公式法 1．q F E =是电场强度的定义式：适用于任何电场，电场中某点的场强是确定值，其大小和方向与试探电荷无关，试探电荷q 充当“测量工具”的作用。 2．2 r k Q E =是真空中点电荷电场强度的决定式，E 由场源电荷Q 和某点到场源电荷的距离r 决定。 3．d U E =是场强与电势差的关系式，只适用于匀强电场，注意式中的d 为两点间的距离在场强方向的投影。 二．对称叠加法 当空间的电场由几个点电荷共同激发的时候，空间某点的电场强度等于每个点电荷单独存在时所激发的电场在该点的场强的矢量和，其合成遵守矢量合成的平行四边形定则。 例：如图，带电量为+q 的点电荷与均匀带电。 例：如图，带电量为+q 的点电荷与均匀带

电薄板相距为2d ，点电荷到带电薄板的垂线通过板的几何中心，如图中a 点处的场强为零，求图中b 点处的场强多大 例：一均匀带负电的半球壳，球心为O 点，AB 为其对称轴，平面L 垂直AB 把半球壳一分为二，L 与AB 相交于M 点，对称轴AB 上的N 点和M 点关于O 点对称。已知一均匀带电球壳内部任一点的电场强度为零，点电荷q 在距离其为r 处的电势为r q k =?。假设左侧部分在M 点的电场强度为 E 1，电势为1?；右侧部分在M 点的电场强 度为E 2，电势为2?；整个半球壳在M 点的电场强度为E 3，在N 点的电场强度为E 4，下列说法中正确的是（ ） A ．若左右两部分的表面积相等，有E 1＞E 2，1?＞2 ?

B ．若左右两部分的表面积相等，有E 1＜E 2，1?＜2 ? C ．只有左右两部分的表面积相等，才有E 1＞E 2，E 3=E 4 D ．不论左右两部分的表面积是否相等，总有 E 1＞E 2，E 3=E 4 答案：D 例：ab 是长为L 的均匀带电细杆，P1、P2是位于ab 所在直线上的两点，位置如图所示．ab 上电荷产生的静电场在P1处的场强大小为E 1，在P2处的场强大小为E2。则以下说法正确的是( ) A ．两处的电场方向相同， E1＞E2 B ．两处的电场方向相反， E1＞E2 C ．两处的电场方向相同，E1＜E2 D ．两处的电场方向相反，E1＜E2 A B M O N L

电场强度的几种计算方法

电场强度的几种求法 一． 公式法 1．q F E = 是电场强度的定义式：适用于任何电场，电场中某点的场强是确定值，其大小和方向与试探电荷无关，试探电荷q 充当“测量工具”的作用。 2．2r k Q E =是真空中点电荷电场强度的决定式，E 由场源电荷Q 和某点到场源电荷的距离r 决定。 3．d U E = 是场强与电势差的关系式，只适用于匀强电场，注意式中的d 为两点间的距离在场强方向的投影。 二．对称叠加法 当空间的电场由几个点电荷共同激发的时候，空间某点的电场强度等于每个点电荷单独存在时所激发的电场在该点的场强的矢量和，其合成遵守矢量合成的平行四边形定则。 例：如图，带电量为+q 的点电荷与均匀带电。 例：如图，带电量为+q 的点电荷与均匀带电薄板相距为2d ，点电荷到带电薄板的垂线通过板的几何中心，如图中a 点处的场强为零，求图中b 点处的场强多大？ 例：一均匀带负电的半球壳，球心为O 点，AB 为其对称轴，平面L 垂直AB 把半球壳一分为二，L 与AB 相交于M 点，对称轴AB 上的N 点和M 点关于O 点对称。已知一均匀带电球壳内部任一点的电场强度为零，点电荷q 在距离其为r 处的电势为r q k =?。假设左侧部分在M 点的电场强度为E 1，电势为1?；右侧部分在M 点的电场强度为E 2，电势为2?；整个半球壳在M 点的电场强度为E 3，在N 点的电场强度为E 4，下列说法中正确的是（ ） A ．若左右两部分的表面积相等，有E 1＞E 2，1?＞2? B ．若左右两部分的表面积相等，有E 1＜E 2，1?＜2?

C ．只有左右两部分的表面积相等，才有E 1＞E 2，E 3=E 4 D ．不论左右两部分的表面积是否相等，总有 E 1＞E 2，E 3=E 4 答案：D 例：ab 是长为L 的均匀带电细杆，P1、P2是位于ab 所在直线上的两点，位置如图所示．ab 上电荷产生的静电场在P1处的场强大小为E 1，在P2处的场强大小为E2。则以下说法正确的是( ) A ．两处的电场方向相同，E1＞E2 B ．两处的电场方向相反，E1＞E2 C ．两处的电场方向相同，E1＜E2 D ．两处的电场方向相反，E1＜E2 三．等效替代法 例：均匀带电的球壳在球外空间产生的电场等效于电荷集中于球心处产生的电场，如图，在半球面A 、B 上均匀分布正电荷，总电荷量为q ，球面半径为R ，CD 为通过半球顶点与球心O 的轴线，在轴线上有M 、N 两点，OM=ON=2R ，已知M 点的场强大小为E ，则N 点场强大小为（ ） A ．E R -22kq B ．24kq R C ．E R -24kq D ．E R +2 4kq 答案：A 例：【2013安徽20】如图所示，xOy 平面是无穷大导体的表面，该导体充满0z <的空间， 0z >的空间为真空。将电荷为q 的点电荷置于z 轴上z=h 处，则在xOy 平面上会产生感应 电荷。空间任意一点处的电场皆是由点电荷q 和导体表面上的感应电荷共同激发的。已知静电平衡时导体内部场强处处为零，则在z 轴上2 h z = 处的场强大小为（k 为静电力常量） A ．24q k h B ．249q k h C ．2329q k h D ．2 409q k h 【答案】D C D A B

KISSsoft关于齿轮强度的计算中文版

3. 强度计算 输入你自己的材料数据 在Kisssoft的数据库中已经包含了一些塑料的数据，如果你想在kisssoft中储存你的一些关于塑料齿轮的数据，你可以使用以下方法： 这里我们用已经做好的POM表 首先点击“Extras”->“Data base tool”，选择相应的数据然后进行计算，如图3-1。或者输入自己的数据，点击“material basic base”并在对话框的底部点击“+”，就会出现一个对话框，在这个对话框中就可以输入数据。如图3-2 （图3-1）

（图3-2） 结合有效的齿型计算强度 在KISSsoft系统中如何激活“graphical method(图解法)”。当你输入强度时，在对话框的右下方点击“Details”按钮，然后在“Form factor Yf and Ys”的下拉菜单中选择“using graphical method”如图所示

现在，计算时首先计算出的是齿轮的齿形系数Yf和它的应力修整系数Ys. 你也可以在KISSsoft系统中显示齿根应变系数，点击“Path of contact”输入你所需的设置参数，并进行运算。如下图： “Path of contact”的设置版面 然后你点击“Graphics”->“Path of contact”, 选择你所需要的图表，例如选择应力强度曲线（stress curve）的2D形式。

Tooth root stresses and Hertzian pressure

Tooth root stresses, progression in the tooth root

电场强度的几种计算方法

微专题训练16 电场强度的几种计算方法 1．(公式法)(单选)如图1所示，真空中O 点有一点电荷，在它产生的电场中有a 、 b 两点，a 点的场强大小为E a ，方向与ab 连线成60°角，b 点的场强大小为E b ，方向与ab 连线成30°角．关于a 、b 两点场强大小E a 、E b 的关系，以下结论正确的是 ( )． 图1 A ．E a ＝33E b B ．E a ＝13E b C ．E a ＝3E b D ． E a ＝3E b 解析 由题图可知，r b ＝3r a ，再由E ＝kQ r 2可知，E a E b ＝r 2b r 2a ＝31，故D 正确． 答案 D 2．(图象斜率法)(多选)如图2甲所示，在x 轴上有一个点电荷Q (图中未画出)，Q 、 A 、 B 为轴上三点，放在A 、B 两点的试探电荷受到的电场力跟试探电荷所带电荷量的关系如图乙所示，则 ( )． 图2 A ．A 点的电场强度大小为2×103 N/C B ．B 点的电场强度大小为2×103 N/C C ．点电荷Q 在A 、B 之间 D ．点电荷Q 在A 、O 之间 解析 对于电场中任意一点而言，放在该处的试探电荷的电荷量q 不同，其受

到的电场力F的大小也不同，但比值F q是相同的，即该处的电场强度．所以F-q 图象是一条过原点的直线，斜率越大则场强越大．由题图可知A点的电场强度 E A＝2×103 N/C，B点的电场强度的大小为E B＝0.6×103 N/C，A正确，B错误．A、 B两点放正、负不同的电荷，受力方向总为正，说明A、B的场强方向相反，点电荷Q只能在A、B之间，C正确． 答案AC 3．(叠加法)(多选)如图3所示，在x轴坐标为＋1的点上固定一个电荷量为4Q的正点电荷，坐标原点O处固定一个电荷量为Q的负点电荷，那么在x坐标轴上，电场强度方向沿x轴负方向的点所在区域应是()． 图3 A．(0,1)B．(－1,0) C．(－∞，－1)D．(1，＋∞) 解析在区域(0,1)中4Q和－Q的电场的电场强度方向都向左，合场强仍向左， A对；在－Q左侧距－Q为x处场强为零，由k Q x2＝k 4Q (1＋x)2 得x＝1，所以区域(－∞，－1)内合场强向左，C对． 答案AC 4．(叠加法)(单选)如图4所示，中子内有一个电荷量为＋2e 3的上夸克和两个电荷量 为－e 3的下夸克，3个夸克都分布在半径为r的同一圆周上，则3个夸克在其圆 心处产生的电场强度大小为() 图4

电场强度的叠加原理及电场强度的计算

第二讲：电场强度的叠加原理及电场强度的计算 容：§9－3 电场强度的求法 要求： 1．理解场强叠加原理； 2．掌握用积分的方法计算电场强度。 重点与难点： 1．电场强度及其计算。 作业： 习题：P37：9，11 预习：电场强度的叠加原理

四、电场强度叠加原理 1．点电荷的场强：电荷Q ，空间r 处 2 04r r Q q F E ＝ 2．点电荷系： 在点电荷系Q 1，Q 2，…，Q n 的电场中，在P 点放一试验电荷q 0，根据库仑力的叠加原理，可知试验电荷受到的作用力为 i F F ，因而P 点的电场强度为 i i i E q F q F q F E ＝ 即 3 04r r Q E E i i ＝＝ 点电荷系电场中某点的场强等于各个点电荷单独存在时在该点的场强的矢量和。这就是电场强度的叠加原理。 3．连续分布电荷激发的场强 将带电区域分成许多电荷元d q ，则 2 04r r dq E d E ＝ 其中，对于电荷体分布，d q =ρd v ， v r r dv E 0 204 ＝ 对于电荷面分布，d q =σds ，02 04r r ds E s ＝ 对于电荷线分布，d q =λd l ， l r r dl E 0 204 ＝ 其中体密度 dV dQ V Q V lim 0＝ 单位C/m 3 ； 面密度 dS dQ S Q S lim 0＝ 单位C/m 2；

线密度 dl dQ l Q l lim 0＝ 单位C/m 。 五、 电场强度的计算： 1．离散型的： 3 04r r Q E E i i ＝＝ 2．连续型的： 2 04r r dq E d E ＝ 空间各点的电场强度完全取决于电荷在空间的分布情况。如果给定电荷的分布，原则上就可以计算出任意点的电场强度。计算的方法是利用点电荷在其周围激发场强的表达式与场强叠加原理。计算的步骤大致如下： ● 任取电荷元d q ，写出d q 在待求点的场强的表达式； ● 选取适当的坐标系，将场强的表达式分解为标量表示式； ● 进行积分计算； ● 写出总的电场强度的矢量表达式，或求出电场强度的大小和方向； ● 在计算过程中，要根据对称性来简化计算过程。 例1． 电偶极子（Electric Dipole ）的场强。 1. 几个概念： （1）两个电量相等、符合相反、相距为l 的点电荷+q 和-q ，若场点到这两个电荷的距离比l 大得多时，这两个点电荷系称为电偶极子。 （2）从-q 指向+q 的矢量l 称为电偶极子的轴。 （3）l q p 称为电偶极子的电偶极矩 2. 电偶极子的电场强度 （1）电偶极子轴线延长线上一点的电场强度 如图所示，取电偶极子轴线的中点为坐标原点O ，沿极轴的延长线为O x 轴，轴上任意点A 距原点的距离为x ，则正负电荷在点A 产生的场强为 i l x q E 2 02/41 i l x q E 2 02/41 由叠加原理可知点A 的总场强为 i l x xl q i l x q l x q E E E 2220 2204/242/2/41 ＝＋－＋＝－ 当x >>l 时，2 224/x l x

标准齿轮模数齿数计算公式

齿轮的直径计算方法： 齿顶圆直径=（齿数+2）*模数 分度圆直径=齿数*模数 齿根圆直径=齿顶圆直径-(4.5×模数) 比如：M4 32齿34*3.5 齿顶圆直径=（32+2）*4=136 分度圆直径=32*4=128 齿根圆直径=136-4.5*4=118 7M 12齿 中心距（分度圆直径1+分度圆直径2）/2 就是（12+2）*7=98 这种计算方法针对所有的模数齿轮（不包括变位齿轮）。 模数表示齿轮牙的大小。 齿轮模数=分度圆直径÷齿数 =齿轮外径÷（齿数-2） 齿轮模数是有国家标准的（1357-78） 模数标准系列（优先选用）1、1.25、1.5、2、2.5、3、4、5、6、8、10、12、14、16、20、25、32、40、50 模数标准系列（可以选用）1.75，2.25，2.75，3.5，4.5，5.5，7，9，14，18，22，28，36，45 模数标准系列（尽可能不用）3.25，3.75，6.5，11，30

上面数值以外为非标准齿轮，不要采用！ 塑胶齿轮注塑后要不要入水除应力 精确测定斜齿轮螺旋角的新方法 （）周节 齿轮分度圆直径d的大小可以用模数(m)、径节()或周节()与齿数(z)表示 径节P()是指按齿轮分度圆直径(以英寸计算)每英寸上所占有的齿数而言

径节与模数有这样的关系: 25.4 1/8模=25.48=3.175 3.175/3.1416(π)=1.0106模 1) 什么是「模数」？ 模数表示轮齿的大小。 R模数是分度圆齿距与圆周率（π）之比,单位为毫米()。 除模数外,表示轮齿大小的还有ＣＰ（周节：）与ＤＰ（径节：）。【参考】齿距是相邻两齿上相当点间的分度圆弧长。 2) 什么是「分度圆直径」？ 分度圆直径是齿轮的基准直径。 决定齿轮大小的两大要素是模数和齿数、 分度圆直径等于齿数与模数(端面)的乘积。 过去,分度圆直径被称为基准节径。最近,按标准,统一称为分度圆直径。 3) 什么是「压力角」？ 齿形与分度圆交点的径向线与该点的齿形切线所夹的锐角被称为分度圆压力角。一般所说的压力角,都是指分度圆压力角。 最为普遍地使用的压力角为２０°,但是,也有使用１４．５°、１５°、１７．５°、２２．５°压力角的齿轮。 4) 单头与双头蜗杆的不同是什么？ 蜗杆的螺旋齿数被称为「头数」,相当于齿轮的轮齿数。

齿轮强度计算公式

第7节 标准斜齿圆柱齿轮的强度计算 一． 齿面接触疲劳强度计算 1. 斜齿轮接触方式 2. 计算公式 校核式： 设计式： 3. 参数取值说明 1) Z E ---弹性系数 2) Z H ---节点区域系数 3) εα---斜齿轮端面重合度 4) β---螺旋角。斜齿轮：β=80～250；人字齿轮β=200～350 5) 许用应力：[σH ]=([σH1]+[σH2])/2≤1.23[σH2] 6) 分度圆直径的初步计算 在设计式中，K 等与齿轮尺寸参数有关，故需初步估算： a) 初取K=K t b) 计算d t c) 修正d t 二． 齿根弯曲疲劳强度计算 1. 轮齿断裂 2. 计算公式校核式： 设计式： 3. 参数取值说明 1) Y F a 、Y Sa ---齿形系数和应力修正系数。Z v =Z/cos 3β→Y Fa 、Y Fa 2) Y β---螺旋角系数。 3) 初步设计计算 在设计式中，K 等与齿轮尺寸参数有关，故需初步估算： d) 初取K=K t e) 计算m nt f) 修正m n [] H t H E H u u bd KF Z Z σεσα≤±=1 1[] 3 2 1112??? ? ??±≥H H E d Z Z u u KT d σεψα[]3 2 1112 ??? ? ??±≥H H E d t t Z Z u u T K d σψ311t t K K d d ≥[]F n sa Fa t F bm Y Y Y KF σεσα β ≤=[]3 2121cos 2F sa Fa d n Y Y z Y KT m σεψβα β≥3t t n n K K m m ≥[] 3 2121cos 2F sa Fa d t nt Y Y z Y T K m σεψβα β≥

MATLAB计算方法迭代法牛顿法二分法实验报告要点

姓名实验报告成绩 评语： 指导教师（签名） 年月日

说明：指导教师评分后，实验报告交院（系）办公室保存。 实验一 方程求根 一、 实验目的 用各种方法求任意实函数方程0)(=x f 在自变量区间[a ，b]上，或某一点附近的实根。并比较方法的优劣。 二、 实验原理 (1)、二分法 对方程0)(=x f 在[a ，b]内求根。将所给区间二分，在分点2a b x -= 判 断是否0)(=x f ；若是，则有根 2a b x -= 。否则，继续判断是否0)()(

（1）11x =0.09033 （2）5x =0.09052 （3）2x =0,09052 五、 实验内容 （1）、在区间[0,1]上用二分法求方程0210=-+x e x 的近似根，要求误差不 超过 3 105.0-?。 （2）、取初值00=x ，用迭代公式=+1k x -0x )(') (k k x f x f ，求方程0210=-+x e x 的 近似根。要求误差不超过 3 105.0-?。 （3）、取初值00=x ，用牛顿迭代法求方程0210=-+x e x 的近似根。要求误 差不超过 3 105.0-?。 六、 实验步骤与实验程序 （1） 二分法 第一步：在MATLAB 7.0软件，建立一个实现二分法的MATLAB 函数文件agui_bisect.m 如下： function x=agui_bisect(fname,a,b,e) %fname 为函数名，a,b 为区间端点，e 为精度 fa=feval(fname,a); %把a 端点代入函数，求fa fb=feval(fname,b); %把b 端点代入函数，求fb if fa*fb>0 error('两端函数值为同号'); end %如果fa*fb>0，则输出两端函数值为同号 k=0 x=(a+b)/2 while(b-a)>(2*e) %循环条件的限制

齿轮强度计算公式

第7节 标准斜齿圆柱齿轮的强度计算 一． 令狐采学 二． 齿面接触疲劳强度计算 1. 斜齿轮接触方式 2. 计算公式 校核式： 设计式： 3. 参数取值说明 1) Z E---弹性系数 2) Z H---节点区域系数 3) ---斜齿轮端面重合度 4) ---螺旋角。斜齿轮：=80～250；人字齿轮=200～350 5) 许用应力：[H]=([H1]+[H2])/2 1.23[H2] 6) 分度圆直径的初步计算 在设计式中，K 等与齿轮尺寸参数有关，故需初步估算： a) 初取K=Kt b) 计算dt c) 修正dt 三． 齿根弯曲疲劳强度计算 1. 轮齿断裂 2. 计算公式校核式： 设计式： 3. 参数取值说明 1) Y Fa 、YSa---齿形系数和应力修正系数。Zv=Z/cos3YFa 、YFa 2) Y ---螺旋角系数。 3) 初步设计计算 在设计式中，K 等与齿轮尺寸参数有关，故需初步估算： d) 初取K=Kt e) 计算mnt [] H t H E H u u bd KF Z Z σεσα≤±=1 1[]32 1112 ??? ? ??±≥H H E d t t Z Z u u T K d σψ[]3 2121cos 2F sa Fa d n Y Y z Y KT m σεψβα β≥[] 32 121cos 2F sa Fa d t nt Y Y z Y T K m σεψβα β≥

f) 修正mn 第8节 标准圆锥齿轮传动的强度计算 一． 作用：用于传递相交轴之间的运动和动力。 二． 几何计算 1. 锥齿轮设计计算简化 2. 锥距 3. 齿数比： u=Z2/Z1=d2/d1=tan 2=cot 1 4. 齿宽中点分度圆直径 dm/d=(R-0.5b)/R=1-0.5b/R 记R=b/R---齿宽系数R=0.25～0.3 dm=(1-0.5R)d 5. 齿宽中点模数 mn=m(1-0.5R) 三． 受力分析 大小： Ft1=2T1/dm1(=Ft2) Fr1=Ft1tan cos Fa2) Fa1=Ft1tan sin 1(=Fr2) 方向： 四． 强度计算 1. 齿面接触疲劳强度计算 1)计算公式： 按齿宽中点当量直齿圆柱齿轮计算，并取齿宽为0.85b ，则： 以齿轮大端参数代替齿宽中点当量直齿圆柱齿轮参数，代入 n 1 n 2 相交轴 n 2 两轴夹角900 n 1 2 2 2122212 21Z Z m d d R +=+= d 1 d m b R d m2 d 2 δ1 δ2 O C 2 C 1 A 2 A 1 q Fr α δ Fa Fn Ft Fa1 Fr 2 2 1 n 1 Fa2 Fr 1 Ft 1 Ft 2 []H v v v v H E H u u bd KT Z Z σσ≤+=1 85.023 1 1

16 电场强度的几种计算方法

微专题训练16　电场强度的几种计算方法1．(公式法)(单选)如图1所示，真空中O点有一点电荷，在它产生的电场中有a、b两点，a点的场强大小为E a，方向与ab连线成60°角，b点的场强大小为E b，方向与ab连线成30°角．关于a、b两点场强大小E a、 E b的关系，以下结论正确的 是 ( )． 图1 A．E a＝E b B．E a＝E b C．E a＝E b D．E a＝3E b 解析　由题图可知，r b＝r a，再由E＝可知，＝＝，故D正确． 答案　D 2．(图象斜率法)(多选)如图2甲所示，在x轴上有一个点电荷Q(图中未画出)，Q、A、B为轴上三点，放在A、B两点的试探电荷受到的电场力跟试探电荷所带电荷量的关系如图乙所示， 则 ( )．

图2 A．A点的电场强度大小为2×103 N/C B．B点的电场强度大小为2×103 N/C C．点电荷Q在A、B之间 D．点电荷Q在A、O之间 解析　对于电场中任意一点而言，放在该处的试探电荷的电荷量q不同，其受到的电场力F的大小也不同，但比值是相同的，即该处的电场强度．所以Fq图象是一条过原点的直线，斜率越大则场强越大．由题图可知A点的电场强度E A＝2×103N/C，B点的电场强度的大小为E B ＝0.6×103 N/C，A正确，B错误．A、B两点放正、负不同的电荷，受力方向总为正，说明A、B的场强方向相反，点电荷Q只能在A、B之间，C正确． 答案　AC 3．(叠加法)(多选)如图3所示，在x轴坐标为＋1的点上固定一个电荷量为4Q的正点电荷，坐标原点O处固定一个电荷量为Q的负点电荷，那么在x坐标轴上，电场强度方向沿x轴负方向的点所在区域应 是 ( )．

MATLAB计算方法迭代法牛顿法二分法实验报告

姓名 _______________ 实验报告成绩 __________________________ 评语： 指导教师（签名）

说明：指导教师评分后，实验报告交院(系)办公室保存 实验一方程求根 一、实验目的 用各种方法求任意实函数方程f(x)0在自变量区间［a，b］上，或某一点附 近的实根。并比较方法的优劣。 二、实验原理 (1)、二分法 b a x --------- 对方程f(x)0在［a，b］内求根。将所给区间二分，在分点2判 b a x --------- 断是否f(x) 0；若是，贝y有根2。否则，继续判断是否f(a)?f(x) 0, 若是，则令b x,否则令a x。否则令a x。重复此过程直至求出方程f(x) 0在［a,b］中的 近似根为止。 (2)、迭代法 将方程f(x) 0等价变换为x = ?( x )形式，并建立相应的迭代公式xk 1 J x)。 (3)、牛顿法 若已知方程的一个近似根X。，则函数在点X。附近可用一阶泰勒多项 式pi(x) f(X0) f'(X0)(X X0)来近似，因此方程f(x) 0可近似表示为f(X0) f(X0) f'(X0)(X X)0设f'(X0) 0,则x X0 f'(X0)。取x 作为原方程新的近 f (X k) 似根X1，然后将X1作为X0代入上式。迭代公式为：Xk 1 X0 f'(X k)。 三、实验设备：MATLAB 7.0软件 四、结果预测

(1) 心=0.09033 (2) x 5=o.o9O52 (3) x 2=0,09052 五、实验内容 (1)、在区间［0,1］上用二分法求方程e x 10x 2 0的近似根，要求误差不 六、实验步骤与实验程序 (1)二分法 第一步：在MATLAB 7.0软件，建立一个实现二分法的 MATLAB 函数 文 件 agui_bisect.m 女口下： fun cti on x=agui_bisect(fname,a,b,e) %fname 为函数名，a,b 为区间端点，e 为精度 fa=feval(fname,a); % 把a 端点代入函数，求fa fb=feval(fname,b); % 把b 端点代入函数，求fb if fa*fb>0 error(' 两端函数值为同号'); end %如果fa*fb>0，则输出两端函数值为同号 k=0 x=(a+b)/2 while(b-a)>(2*e) % 循环条件的限制 fx 二feval(fname,x);% 把x 代入代入函数，求fx 超过 0.5 10 3 (2)、取初值 xo 0 ，用迭代公式 Xk1 f (X k ) f'(X k )，求方程 e x 10x 2 0 的 近似根。要求误差不超过 0.5 10 (3)、取初值x 。 ，用牛顿迭代法求方程e x 10x 2 0的近似根。要求误 差不超过 0.5 10

标准直齿圆柱齿轮传动强度计算

§8-5 标准直齿圆柱齿轮传动的强度计算 一．齿轮传动承载能力计算依据 轮辐、轮缘、轮毂等设计时，由经验公式确定尺寸。若设计新齿，可参《工程手册》20、22篇，用有限元法进行设计。 轮齿的强度计算： 1．齿根弯曲强度计算：应用材料力学弯曲强度公式W M b = σ进行计算。数学模型：将轮齿看成悬臂梁，对齿根进行计算，针对齿根折断失效。

险截面上，γcos ca p --产生剪应力τ，γsin ca p 产生压应力σc ，γcos .h p M ca =产生弯曲应力σF 。分析表明，σF 起主要作用，若只用σF 计算齿根弯曲疲劳强度，误差很小（<5%），在工程计算允许范围内，所以危险剖面上只考虑σF 。 单位齿宽（b=1）时齿根危险截面的理论弯曲应力为 2 20cos .66 *1cos .S h p S h p W M ca ca F γγσ=== 令α cos ,,b KF L KF p m K S m K h t n ca S h = ===，代入上式，得 ()αγαγσcos cos 6.cos cos ..622 0S h t S h t F K K bm KF m K b m K KF == 令 αγc o s c o s 62 S h Fa K K Y = Fa Y --齿形系数，表示齿轮齿形对σF 的影响。Fa Y 的大小只与轮齿形状有关（z 、h *a 、c *、

α）而与模数无关，其值查表10-5。 齿根危险截面理论弯曲应力为 bm Y KF Fa t F = 0σ 实际计算时，应计入载荷系数及齿根危险剖面处的齿根过渡曲线引起的应力集中的影响。 bm Y Y KF Sa Fa t F = σ 式中：Sa Y --考虑齿根过渡曲线引起的应力集中系数，其影响因素同Fa Y ，其值可查表10-5。 2．齿根弯曲疲劳强度计算 校核公式 []F Fa Sa Sa Fa t F Y Y bmd KT bm Y Y KF σσ≤== 1 1 2 MPa 令1 d b d = φ，d φ--齿宽系数。 将111,mz d d b d ==φ代入上式 设计公式 [])(.23 211mm Y Y z KT m F Sa Fa d σφ≥

数值计算方法-简单迭代

《数值计算方法》实验2报告 班级： 学号： 姓名： 成绩： 1. 实验名称 实验2 非线性方程的迭代解法(之简单迭代法) 2. 实验题目 用简单迭代法求方程010423=-+x x 在区间[1,2]内的一个实根，取绝对误差限为410-． 3. 实验目的 掌握非线性方程的简单迭代法． 4. 基础理论 简单迭代法：将方程0)(=x f 改写成等价形式)(x x ?=，从初值0x 开始，使用迭代公式)(1k k x x ?=+可以得到一个数列，若该数列收敛，则其极限即为原方程的解．取数列中适当的项可作为近似解． 5. 实验环境 操作系统：Windows xp ； 程序设计语言：Matlab 6. 实验过程 (1)、输入初值x0，将方程改写等价形式； (2)、构造迭代公式 ，k=0,1,2.......； (3)、采用for 循环实现反复迭代； (4)、以误差410-终止循环； (5)、输出结果。 7. 结果与分析 x0=0 del=0.0001 N=100 k x(k) 0 0.000000 1 1.581139 2 1.229548

3 1.426638 4 1.331951 5 1.381756 6 1.356640 7 1.369593 8 1.362987 9 1.366376 10 1.364643 11 1.365531 12 1.365076 13 1.365309 14 1.365190 近似解=1.365251 100次迭代后未达到精度要求. 分析：计算结果与理论一致。 8. 附录：程序清单 phi=inline('0.5*sqrt((10-x^3))'); x0=input('x0='); del=input('del='); N=input('N='); n=1; fprintf('\n k x(k) '); fprintf('\n %4d %f',0,x0); while n