飞剪计算程序代码
飞剪机构尺寸设计源程序代码
/**飞剪机构尺寸设计源程序代码**/
#include
#include
void main ()
{
/**参数定义**/
double
a0,b0,c0,d0=1.0,a,b,c,d,blc,blcx,e,f,h=0.25,ax,bx,cx,dx,ex,fx,deltah=0.005,deltaarf,deltafaiE,m,mx, k,k1,k1x,L=1.0,
fai,faiBD,fai2,fai3,fai2x,fai3x,faiE,gama1,gama2,thita,thita0,A,B,N,arf3,arf4,
LBD,LCE,LCF,LCEF,FZ2,FZ3,arf2x,arf3x,arf4x,jBEC,jDEC,jDAB,
jqj,jqjx,w1,w2,w3,VEt,VFt,Vd,Vt=2.0,lgxs0,lgxs,pai=3.141593
,Mb,Fc=98000.0,VEx,VFx;
/**参数输入**/
printf("请输入行程速比系数k:");
scanf("%lf",&k);
printf("请输入远极位夹角gama2:");
scanf("%lf",&gama2);
printf("请输入摇杆摆角fai:");
scanf("%lf",&fai);
printf("请输入机架角:");
scanf("%lf",&arf4);
printf("请输入拉钢系数lgxs:");
scanf("%lf",&lgxs);
printf("请输入k1:");
scanf("%lf",&k1);
/**相对尺寸设计**/
m=pai/180.0;
thita=pai*(k-1)/(k+1);
gama2*=m;
fai*=m;
arf4*=m;
gama1=fai+gama2-thita;
thita0=atan(sin(gama2)*sin(thita)/(sin(gama1)-sin(gama2)*cos(thita)));
A=cos(thita+thita0)*sin(gama2+thita0);
B=sin(gama2)+sin(thita0)*cos(gama1+thita+thita0);
N=2.0*sin(gama2)*cos(thita+thita0);
a0=(A-B)/N;b0=(A+B)/N;c0=sin(thita0)/sin(gama2);
/**绝对尺寸**/
a=k1*lgxs*L/(2.0*pai);
blc=a/a0;
b=blc*b0;c=blc*c0;d=blc;
/**刀刃位置确定**/
f=d*cos(arf4)-h;
e=sqrt((f-deltah)*(f-deltah)+(d-a)*(d-a)-2.0*(f-deltah)*(d-a)*cos(arf4)); deltaarf=acos((c*c+(d-a)*(d-a)-b*b)/(2.0*c*(d-a)));
arf3=arf4+deltaarf;
LCE=sqrt(c*c+(f-deltah)*(f-deltah)-2.0*c*(f-deltah)*cos(arf3));
LCF=sqrt(c*c+f*f-2*c*f*cos(arf3));
LCEF=(LCE+LCF)/2.0;
arf2x=acos((b*b+e*e-LCEF*LCEF)/(2.0*b*e));
arf3x=acos((c*c+f*f-LCEF*LCEF)/(2.0*c*f));
/**剪切角的确定**/
jBEC=acos((e*e+LCEF*LCEF-b*b)/(2.0*e*LCEF));
jDEC=acos((f*f+LCEF*LCEF-c*c)/(2.0*f*LCEF));
FZ2=jBEC+jDEC;
LBD=sqrt(e*e+f*f-2.0*e*f*cos(FZ2));
jDAB=acos((a*a+d*d-LBD*LBD)/(2.0*a*d));
jqj=arf4-jDAB;
/**调整上下刀刃水平速度误差**/
faiBD=atan((d*sin(arf4)-a*sin(jqj))/(d*cos(arf4)-a*cos(jqj)));
FZ3=acos((b*b-LBD*LBD-c*c)/(2.0*LBD*c));
FZ3=2.0*pai-FZ3;
fai3=faiBD+FZ3;
/**验证正负**/
fai2=atan((LBD*sin(faiBD)+c*sin(fai3))/(LBD*cos(faiBD)+c*cos(fai3)));
fai2+=2.0*pai;
faiE=atan((e*sin(fai2+arf2x)-b*sin(fai2))/(e*cos(fai2+arf2x)-b*cos(fai2))); /**验证正负**/
faiE+=pai;
deltafaiE=faiE-pai/2.0;
arf4x=arf4-deltafaiE;
jqjx=jqj-deltafaiE;
fai2x=fai2-deltafaiE;
fai3x=fai3-deltafaiE;
/**调整刀刃与钢材运动速度同步**/
w1=2.0*pai*Vt/L;
w2=(-a*w1*sin(jqjx-fai3x))/(b*sin(fai2x-fai3x));
w3=(a*w1*sin(jqjx-fai2x))/(c*sin(fai3x-fai2x));
VEt=a*w1*cos(jqjx)+e*w2*cos(fai2x+arf2x);
VFt=w3*f*cos(jDEC-pai/2.0);
VFt*=-1.0;/**改的没有理由，可能是计算过程中没有规定的正负**/
Vd=(VEt+VFt)/2.0;
lgxs0=Vd/Vt;
k1x=a*w1/Vd;
ax=k1x*lgxs*Vt/w1;
blcx=ax/a;
bx=blcx*b;cx=blcx*c;dx=blcx*d;ex=blcx*e;fx=blcx*f;
/**求曲柄的平衡力矩Mb**/
VEx=-ax*w1*sin(jqjx)-ex*w2*sin(fai2x+arf2x);
VFx=w3*fx*sin(jDEC-pai/2.0);
Mb=Fc*(VEx-VFx)/w1;
/**弧度角转度数**/
mx=180.0/pai;
jqjx*=mx;
arf2x*=mx;arf3x*=mx;arf4x*=mx;
faiE*=mx;deltafaiE*=mx;
thita0*=mx;deltaarf*=mx;
jBEC*=mx;jDEC*=mx;jDAB*=mx;jqj*=mx;
faiBD*=mx;fai3*=mx;fai2*=mx;fai2x*=mx;fai3x*=mx;
printf("机构的相对尺寸是：\na0=%lf b0=%lf c0=%lf d0=%lf\n",a0,b0,c0,d0);
printf("机构的绝对尺寸是：\na=%lf b=%lf c=%lf d=%lf e=%lf f=%lf\n",a,b,c,d,e,f); printf("机构的修正绝对尺寸是：\nax=%lf bx=%lf cx=%lf dx=%lf ex=%lf fx=%lf\n",ax,bx,cx,dx,ex,fx);
printf("其他尺寸分别是：\nthita=%lf gama1=%lf thita0=%lf k1x=%lf\n",thita,gama1,thita0,k1x);
printf("\nlgxs=%lf blc=%lf blcx=%lf deltaarf=%lf LCE=%lf\n",lgxs,blc,blcx,deltaarf,LCE); printf("\nLCF=%lf LCEF=%lf arf2x=%lf arf3x=%lf\n",LCF,LCEF,arf2x,arf3x);
printf("\njBEC=%lf jDEC=%lf LBD=%lf jDAB=%lf\n",jBEC,jDEC,LBD,jDAB);
printf("\njqj=%lf fai2=%lf faiE=%lf deltafaiE=%lf\n",jqj,fai2,faiE,deltafaiE);
printf("\narf4x=%lf jqjx=%lf fai2x=%lf fai3x=%lf\n",arf4x,jqjx,fai2x,fai3x);
printf("\nw2=%lf w3=%lf VEt=%lf VFt=%lf\n",w2,w3,VEt,VFt);
printf("\nVd=%lf lgxs0=%lf VEx=%lf VFx=%lf Mb=%lf\n",Vd,lgxs0,VEx,VFx,Mb);
}
剪切力的计算方法
第3章 剪切和挤压的实用计算 3.1 剪切的概念 在工程实际中，经常遇到剪切问题。剪切变形的主要受力特点是构件受到与其轴线相垂直的大小相等、方向相反、作用线相距很近的一对外力的作用(图3-1a)，构件的变形主要表现为沿着与外力作用线平行的剪切面(n m -面)发生相对错动(图3-1b)。 图3-1 工程中的一些联接件，如键、销钉、螺栓及铆钉等，都是主要承受剪切作用的构件。构件剪切面上的力可用截面法求得。将构件沿剪切面n m -假想地截开，保留一部分考虑其平衡。例如，由左部分的平衡，可知剪切面上必有与外力平行且与横截面相切的力Q F (图3-1c)的作用。Q F 称为剪力，根据平衡方程∑=0Y ，可求得F F Q =。 剪切破坏时，构件将沿剪切面(如图3-la 所示的n m -面)被剪断。只有一个剪切面的情况，称为单剪切。图3-1a 所示情况即为单剪切。 受剪构件除了承受剪切外，往往同时伴随着挤压、弯曲和拉伸等作用。在图3-1中没有完全给出构件所受的外力和剪切面上的全部力，而只是给出了主要的受力和力。实际受力和变形比较复杂，因而对这类构件的工作应力进行理论上的精确分析是困难的。工程中对这类构件的强度计算，一般采用在试验和经验基础上建立起来的比较简便的计算方法，称为剪切的实用计算或工程计算。 3.2 剪切和挤压的强度计算 3.2.1 剪切强度计算 剪切试验试件的受力情况应模拟零件的实际工作情况进行。图3-2a 为一种剪切试验装置的简图，试件的受力情况如图3-2b 所示，这是模拟某种销钉联接的工作情形。当载荷F 增大至破坏载荷b F 时，试件在剪切面m m -及n n -处被剪断。这种具有两个剪切面的情况，称为双剪切。由图3-2c 可求得剪切面上的剪力为 2 F F Q =
数值计算方法比较
有限差分方法(FDM:Finite Difference Method)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。有限差分法主要集中在依赖于时间的问题（双曲型和抛物型方程）。有限差分法方面的经典文献有Richtmeyer & Morton的《Difference Methods for Initial-Value Problems》；R. LeVeque《Finite Difference Method for Differential Equations》；《Numerical Methods for C onservation Laws》。 注：差分格式： （1）从格式的精度来划分，有一阶格式、二阶格式和高阶格式。 （2）从差分的空间形式来考虑，可分为中心格式和逆风格式。 （3）考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。 目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 构造差分的方法： 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 有限差分法的不足：由于采用的是直交网格，因此较难适应区域形状的任意性，而且区分不出场函数在区域中的轻重缓急之差异，缺乏统一有效的处理自然边值条件和内边值条件的方法，难以构造高精度(指收敛阶)差分格式，除非允许差分方程联系更多的节点(这又进一步增加处理边值条件韵困难)。另外它还有编制不出通用程序的困难。 有限差分法的优点：该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念 直观，表达简单，精度可选而且在一个时间步内，对于一个给定点来说其相关的空间点只是 与该相邻的几点，而不是全部的空间点。是发展较早且比较成熟的数值方法 广义差分法（有限体积法）（GDM:Generalized Difference Method）：1953年，Mac—Neal 利用积分插值法(也称积分均衡法)建立了三角网格上的差分格 式，这就是以后通称的不规划网格上的差分法．这种方法的几何误差小，特别是给出了处理自然边值条件(及内边值条件)的有效方法，堪称差分法的一大进步。1978年，李荣华利用有限元空间和对偶单元上特征函数的推广——局部Taylor展式的公项，将积分插值法改写成广义Galerkin法形式，从而将不规则网格差分法推广为广义差分法.其基本思路是，将计算区域划分为一系列不重复的控制体积，并使每个网格点周围有
高程布置参考—给水处理厂课程设计计算手册
给水处理厂课程设计计算书 12．高程布置 为了配合平面布置，我们首先应根据下表估计各构筑物之间连接管渠的大小及长度大致水头损失。然后在平面布置确定后，按水力学公式逐步计算各构筑物之间的水 构筑物 沉淀池~滤池0.3~0.5 快滤池内 2.0~3.0 虹吸、无阀滤池 1.5~2.0 滤池到清水池0.3~0.5 1.3.4高程布置设计计算
1.3.4.1水处理构筑物的高程布置设计计算 1.水头损失计算 在处理工艺流程中，各构筑物之间水流应为重力流。两构筑物之间水面高差即为流程中的水头损失，包括构筑物本身、连接管道、计量设备等水头损失在内。水头损失应通过计算确定，并留有 余地. （1）处理构筑物水头损失 处理构筑物中的水头损失与构筑物的型式和构造有关，具体根据设计手册第3册表15-13 g ——重力加速度，2/m s 。 ① 配水井至絮凝池连接管线水头损失 a ）沿程水头损失 配水井至絮凝池连接管采用800DN 钢管，管长15l m =。 考虑浑水的因素0.015n =，按0.013n =查设计手册第1册水力计算表得 1.8i =‰，换算成相当 于0.015n =时的i ： 浑水管长15m 算得沿程损失为：
b）局部水头损失 管路中，进口1个，局部阻力系数 10.50 ξ=；急转弯管1个， 20.90 ξ=；闸阀1个， 30.06 ξ=； 90o弯头1个， 41.05 ξ= ；出口1个，局部阻力系数 5 0.04 ξ=，则局部阻力系数总计为： 管内流速 1.11/ v m s =，则管路局部水头损失为： c）总水头损失 ②絮凝池至沉淀池 絮凝池与沉淀池合建，其损失取0.1m。 ③沉淀池至V a）沿程水头损失 沉淀池至V型滤池连接管采用900 DN钢管，管长l= 21.052 2.1 ξ=?＝； 闸阀2 43.0 ξ=；出口1个，V，按0.013 n=查设计手册第1册水力计算表得 2.4 i=‰，则V型滤池至清水池连接管沿程损失为： b）局部水头损失 管路中，进口1个，局部阻力系数 10.50 ξ=；90?弯头3个，局部阻力系数 21.053 3.15 ξ=?＝； 闸阀1个， 30.06 ξ=；出口1个，局部阻力系数 41.00 ξ=，则局部阻力系数总计为：管内流速 1.0/ v m s =，则管路局部水头损失为： c）总水头损失
数值分析(计算方法)总结
第一章绪论 误差来源：模型误差、观测误差、截断误差（方法误差）、舍入误差 是的绝对误差，是的误差，为的绝对误差限（或误差限） 为的相对误差，当较小时，令 相对误差绝对值得上限称为相对误差限记为：即： 绝对误差有量纲，而相对误差无量纲 若近似值的绝对误差限为某一位上的半个单位，且该位直到的第一位非零数字共有n位，则称近似值有n位有效数字，或说精确到该位。 例：设x==…那么,则有效数字为1位，即个位上的3，或说精确到个位。 科学计数法：记有n位有效数字，精确到。 由有效数字求相对误差限：设近似值有n位有效数字，则其相对误差限为 由相对误差限求有效数字：设近似值的相对误差限为为则它有n位有效数字 令 1.x+y近似值为和的误差（限）等于误差（限）的 和 2.x-y近似值为 3.xy近似值为 4. 1．避免两相近数相减 2．避免用绝对值很小的数作除数 3．避免大数吃小数 4．尽量减少计算工作量 第二章非线性方程求根 1.逐步搜索法 设f (a) <0, f (b)> 0，有根区间为(a, b)，从x0=a出发，按某个预定步长(例如h=(b-a)/N)
一步一步向右跨，每跨一步进行一次根的搜索，即判别f(x k)=f(a+kh)的符号，若f(x k)>0(而 f(x k-1)<0),则有根区间缩小为[x k-1,x k] (若f(x k)=0，x k即为所求根), 然后从x k-1出发，把搜索步长再缩小，重复上面步骤，直到满足精度：|x k-x k-1|
易飞ERP月结和成本操作手册
y 易飞ERP系统月结方案 1、月结第一项：存货管理系统-录入库存交易单-把那些没有货号的产品出库（999）-记入研发部门-进入高级 查询 月结流程描述： 1、相关子系统单据完整性检查（销售、采购、库存、生产） 销售部门：应该检查系统中的销货单和销退单，查看这些单据的完整性：未审核、单据数量、价格有无错误以及单据中的其他信息，如果是未审核又不要的单据，应作“作废”处理。另外，对于不再执行发货的客户订单，应作“指定结束”处理。 采购部门：应该检查系统中的进货单和退货单，查看这些单据的完整性：未审核、单据数量、价格有无错误以及单据中的其他信息，如果是未审核又不要的单据，应作“作废”处理。另外，对于不再执行进货的请购单和采购单，应作“指定结束”处理。 重点注意是否有进货单上进货价格为0的,为防止成本计算不准确,必须由采购员提供实际价格。 仓管部门：应该检查系统中的库存交易单、调拨单和成本调整单，查看这些单据的完整性：未审核、单据数量、价格有无错误以及单据中的其他信息，如果是未审核又不要的单据，应作“作废”处理。如果是要抛转到总帐系统的单据,如制费领料单,维修领料单等还要检查单头“部门”字段是否已填入正确信息。 生产部门：应该检查系统中的领料单,退料单,生产入库单,委外进货单,委外退货单, 查看这些单据的完整性：未审核、单据数量、价格有无错误以及单据中的其他信息，如果是未审核又不要的单据，应作“作废”处理。另外,对于不再执行生产的工单,应作”指定结束”处理，对于不需要发放的LRP生产计划，要及时用“清除批次计划”来清除。 财务部门：视需要抽查上述部门的单据。 以销货单为例，设定查询条件【高级设置：审核码等于未审核、销货日期小于或等于20110131、销货日期大于等于20110101】
数值计算方法大作业
目录 第一章非线性方程求根 (3) 1.1迭代法 (3) 1.2牛顿法 (4) 1.3弦截法 (5) 1.4二分法 (6) 第二章插值 (7) 2.1线性插值 (7) 2.2二次插值 (8) 2.3拉格朗日插值 (9) 2.4分段线性插值 (10) 2.5分段二次插值 (11) 第三章数值积分 (13) 3.1复化矩形积分法 (13) 3.2复化梯形积分法 (14) 3.3辛普森积分法 (15) 3.4变步长梯形积分法 (16) 第四章线性方程组数值法 (17) 4.1约当消去法 (17) 4.2高斯消去法 (18) 4.3三角分解法 (20)
4.4雅可比迭代法 (21) 4.5高斯—赛德尔迭代法 (23) 第五章常积分方程数值法 (25) 5.1显示欧拉公式法 (25) 5.2欧拉公式预测校正法 (26) 5.3改进欧拉公式法 (27) 5.4四阶龙格—库塔法 (28)
数值计算方法 第一章非线性方程求根 1.1迭代法 程序代码： Private Sub Command1_Click() x0 = Val(InputBox("请输入初始值x0")) ep = Val(InputBox(请输入误差限ep)) f = 0 While f = 0 X1 = (Exp(2 * x0) - x0) / 5 If Abs(X1 - x0) < ep Then Print X1 f = 1 Else x0 = X1 End If Wend End Sub 例：求f(x)=e2x-6x=0在x=0.5附近的根（ep=10-10）
1.2牛顿法 程序代码： Private Sub Command1_Click() b = Val(InputBox("请输入被开方数x0")) ep = Val(InputBox(请输入误差限ep)) f = 0 While f = 0 X1 = x0 - (x0 ^ 2 - b) / (2 * b) If Abs(X1 - x0) < ep Then Print X1 f = 1 Else x0 = X1 End If Wend End Sub 例：求56的值。（ep=10-10）
塔的水力学计算手册
塔的水力学计算手册
1.目的与适用范围 (1) 2.塔设备特性 (1) 3.名词术语和定义 (1) 4.浮阀/筛孔板式塔盘的设计 (1) 5.填料塔的设计 (1)
1.目的与适用范围 为提高工艺工程师的设计质量,推广计算机应用而编写本手册。 本手册是针对气液传质塔设备中的普遍性问题而编写。对于某些具体塔设备的数据(比如:某生产流程中针对某塔设备的板效率而采用的计算关联式,或者对于某吸收填料塔的传质单元高度或等板高度而采用的具体计算公式)则未予收入。本设计手册以应用为主，主要是指导性的计算方法和步骤，并配合相应的计算程序，具体公式及理论推阐可参考有关文献。 2.塔设备特性 作为气(汽)、液两相传质用的塔设备，首先必须能使气(汽)、液两相得到充分的接触，以得到较高的传质分离效率。 此外，塔设备还应具有以下一些特点： (1)当气(汽)、液处理量过大(超过设计值)时，仍不致于发生大量的雾 沫挟带或液泛等影响正常操作的现象。 (2)当操作波动(设计值的50％～120％)较大时，仍能维持在较高的传 质效率下稳定操作，并具有长期连续操作所必须具备的可靠性。 (3)塔压力降尽量小。 (4)结构简单、耗材少、制造和安装容易。 (5)耐腐蚀、不易堵塞。 (6)塔内的滞留液量要小。 3.名词术语和定义 3.1 塔径(tower diameter)，D T 塔筒体内壁直径，见图3.1-(a)。 3.2 板间距(tray spacing)，H T 塔内相邻两层塔盘间的距离，见图3.1-(a)。 3.3 降液管(downcomer)，DC 各层塔盘之间专供液相流体通过的组件，单溢流型塔盘为侧降液管，双溢流型塔盘有侧降液管和中央降液管，三或多溢流型塔盘有侧降液管、偏侧降液管、偏中央降液管及中央降液管。 3.4 降液管顶部宽度(DC top width)，Wd 弓形降液管面积的弦高。掠堰另有算法，见图3.1-(a),-(b)。 3.5 降液管底间隙(DC clearance)，ho 降液管底部边缘至塔盘(或受液盘)之间的距离，见图3.1-(a)。 3.6 溢流堰高度(weir height)，hw 降液管顶部边缘高出塔板的距离，见图3.1-(a)。 3.7 总的塔盘横截面积(total tower cross-section area)，A T
材料力学-切应力计算
第四章弹性杆横截面上的切应力分析 § 4-3梁横力弯曲时横截面上的切应力 梁受横弯曲时，虽然横截面上既有正应力，又有切应力。但一般情况下，切应力 对梁的强度和变形的影响属于次要因素，因此对由剪力引起的切应力，不再用变形、物理和静力关系进行推导，而是在承认正应力公式（6-2）仍然适用的基础上，假定剪应力在横截面 上的分布规律，然后根据平衡条件导出剪应力的计算公式。 1.矩形截面梁 对于图4-15所示的矩形截面梁，横截面上作用剪力F Q。现分析距中性轴z为y的横线aa1 上的剪应力分布情况。根据剪应力成对定理，横线aa1两端的剪应力必与截面两侧边相切， 即与剪力F Q的方向一致。由于对称的关系，横线aa i中点处的剪应力也必与F Q的方向相同。 根据这三点剪应力的方向，可以设想aa i线上各点切应力的方向皆平行于剪力F Q。又因截面高度h大于宽度b,切应力的数值沿横线aa i不可能有太大变化，可以认为是均匀分布的。基于上述分析，可作如下假设： 1）横截面上任一点处的切应力方向均平行于剪hj力F Q。 2）切应力沿截面宽度均匀分布。 图4-15 图4-16 基于上述假定得到的解，与精确解相比有足够的精确度。从图4-16a的横弯梁中截出dx 微段，其左右截面上的内力如图4-16b所示。梁的横截面尺寸如图4-16c所示，现欲求距中性 轴z为y的横线aa1处的切应力。过aa1用平行于中性层的纵截面aa2C1自dx微段中截出 一微块（图4-16d）。根据切应力成对定理，微块的纵截面上存在均匀分布的剪应力。微块左右侧面上正应力的合力分别为N1和N2，其中
y 1dA 。 A * 由微块沿x 方向的平衡条件 这样，式（4-32）可写成 N 1 I dA A * My 1 dA Ms ； z A * I z (4-29) N 2 II dA (M dM)y 1dA A * A * I z (M dM)。 * ^n^Sz (4-30) 式中，A 为微块的侧面面积， （ii ）为面积 A 中距中性轴为 y i 处的正应力， 将式 N 1 N 2 （4-29）和式（4-30）代入式 dM * nr S z bdx 0 4-31)，得 bdx 0 dM S ； dx bI z (4-31) 因 F Q , dx ，故求得横截面上距中性轴为 y 处横线上各点的剪应力 * F Q S Z bn (4-32) 式（4-32）也适用于其它截面形式的梁。式中， F Q 为截面上的剪力; I z 为整个截面 对中性轴z 的惯性矩；b 为横截面在所求应力点处的宽度; S y 为面积A *对中性轴的静矩。 对于矩形截面梁（图4-17），可取dA bdy i ，于是 * S z y i dA A 2(h y 2) 电( h! y 2) 上式表明，沿截面高度剪应力 4-17 ）。 按抛物线规律变化（图 在截面上、下边缘处，y= ± h , =0；在中性轴上，y=0， 2 切应力值最大，其值为 ■ 1 1 r 尸蛰 T *17 A" y 图 4-17 * S z 0,得
数值计算方法》试题集及答案
《计算方法》期中复习试题 一、填空题： 1、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案：2.367，0.25 2、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2 x 的系数为 ，拉 格朗日插值多项式为 。 答案：-1， )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 3、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字； 4、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( )； 答案 )(1)(1n n n n n x f x f x x x '--- =+ 5、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 )； 6、计算方法主要研究( 截断 )误差和( 舍入 )误差； 7、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时，二分n 次后的误差限为 ( 1 2+-n a b )； 8、已知f (1)＝2，f (2)＝3，f (4)＝5.9，则二次Newton 插值多项式中x 2系数为( 0.15 )； 11、 两点式高斯型求积公式?1 d )(x x f ≈( ?++-≈1 )] 321 3()3213([21d )(f f x x f )，代数精度 为( 5 )； 12、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少，应将该表达 式改写为 11 ,))64(3(10-= -++=x t t t t y ，为了减少舍入误差，应将表达式1999 2001-
工艺专业塔器水力学计算设计导则
1 塔器设计概述 1.1 石油化工装置中塔器占有很大的比重。几乎每种工艺流程都存在蒸馏或吸收等分离单元过程，因此塔器设计至关重要。往往塔器设计的优劣，决定着装置的先进性和经济性，必须给予重视。 1.2 塔器设计与工艺流程设计有着非常密切的关系，亦即塔器的选型和水力学计算与工艺流程的设计计算是结合在一起的。有时塔器设计影响着分离流程和操作条件的选择。例如减小蒸馏塔的回流比，能降低能耗，但塔板数增加，对塔器讲就是减小塔径和增加塔高，其中必有一个最经济条件的选择。又如真空塔或对釜温有要求的蒸馏塔均对压降要求较严，需要选择压降低的板式塔或填料塔，在塔器水力学计算后，压降数据要返回工艺作釜温核算。 1.3 一般工艺流程基本确定后，进行塔器的选型、设计等工作。塔器设计涉及到工艺、化学工程、设备、仪表、配管等专业。化学工程专业的任务及与各专业间关系另有说明。见化学工程专业工作手册H-P0101-96、H-P0301-96。 1.4 随着石油化工和科技的迅猛发展，蒸馏塔从一般的一股进料、二股产品的常规塔发展为多股进料、多侧线，有中间换热的复杂塔。要求塔的生产能力大、效率高、塔板数多，即大塔径、多程数、高效、低压降等，对塔器设计提出了更高的要求，并推动了塔器设计工作的发展。 1.5 近年来电子计算机的普及和发展，为工艺与塔器设计提供了有力的工具。我们可应用PROCESS或PRO/Ⅱ等工艺流程模拟软件进行计算，得到塔的最大和最小汽液负荷、密度等数据，以便进行分段的塔的水力学计算，使工艺和塔的水力学计算能同步进行，并作多方案比较，求得最佳设计。 1.6 设计中主要考虑的问题 1.6.1 确定工艺流程（尤其是分离流程） 通过工艺流程模拟电算，选定最佳切割方案，其中包括多股进料、侧线采出、进料状态和位置等方面的选择。 1.6.2 塔压的设定
易飞ERP高级成本计算
成本系统 (一)成本参数设定 1：成本分类方式 使用时机：打印[直接原料明细表]时，会按照此处定义的分类显示 2、半成品投入归成上阶材料成本 成本计算时，若领用的下阶材料有半成品，其半成品的材料成本、人工和制费直接视为领用的材料成本；打勾分部步结转法，不然就是平行结转法。 3、材料在制约量 ⑴按实际投入材料成本计算 计算步骤 1、计算工单生产一单位所需耗用的材料 每单位材料成本=(工单单身需领用量*月档单位成本)/工单的预计产量 2、计算出工单的约单产量 约当产量=工单实际领用材料金额/每单位材料成本 3、计算出工单的在制约量 在制约量=期初在制约量+约当产量-(生产入库量+委外进货量+报废量+破坏量) 注意事项： a：领料超领太多时，造成在制约量高估计 b：领料少领时，造成在制约量低估，若工单有入库则可能会呈现负数。 c：成会人员须定期检视工但用料的合理性或调整当期材料的在制约量 d：此种方式计算出来的材料在制约量会比较接近实际 ⑵按投入套数计算（领料套数) 适用对象：整批领料、制造前置时间短、领料日与投料日比较相近 计算逻辑：期初在制约量+已领套数– (生产入库量+委外进货量-委外退货量+报废数量+破坏数量) ⑶无在制约当产量 当期投入材料成本全数分摊至当期入库数量 适用于：当期投入当期完工者，或者在制成本非常低可以忽略不计 4、人工制费在制约量 ⑴按实际投入材料成本计算 适用成本要素中材料成本占极大部分，或标准成本无法收集者 ⑵按投入套数计算 ⑶无在制约量 当期投入的人工及制费成本全数分摊至当期入库数量 适用于：当期投入当期完工者 ⑷按实际投入工时计算 有设定标准工时才可预估生产一单位产品的工时 若无工时则需维护[录入品号信息]的[单位标准人工成本]与[单位标准制造费用]，才可预估生产一单位产品的人工制费成本。 若公司有标准公司为依据，还需用考虑[工资率与制费分摊] 标准：取[录入工作中心]的标准工资率及标准制费分摊率 实际：取[录入工作中心成本]的单位人工及单位制费 计算逻辑：半成品归上阶 If 工艺管理子系统->录入工单工艺->标准工时 <> 0 then Begin
剪切应力计算
拉伸、压缩与剪切 1 基本概念及知识要点 1.1 基本概念 轴力、拉（压）应力、力学性能、强度失效、拉压变形、胡克定律、应变、变形能、静不定问题、剪切、挤压。 以上概念是进行轴向拉压及剪切变形分析的基础，应准确掌握和理解这些基本概念。 1.2 轴向拉压的内力、应力及变形 1．横截面上的内力：由截面法求得横截面上内力的合力沿杆的轴线方向，故定义为轴力 F N ，符号规定：拉力为正，压力为负。工程上常以轴力图表示杆件轴 力沿杆长的变化。 2．轴力在横截面上均匀分布，引起了正应力，其值为 F A σ= N 正应力的符号规定：拉应力为正，压应力为负。常用的单位为MPa 、Pa 。 3．强度条件 强度计算是材料力学研究的主要问题之一。轴向拉压时，构件的强度条件是 []F A σσ= ≤N 可解决三个方面的工程问题，即强度校核、设计截面尺寸及确定许用载荷。 4．胡克定律 线弹性范围内，杆的变形量与杆截面上的轴力F N 、杆的长度l 成正比，与截面尺寸A 成反比；或描述为线弹性范围内，应力应变成正比，即 F l l E E A σε?= =N 式中的E 称为材料的弹性模量，EA 称为抗拉压刚度。胡克定律揭示在比例极限内，应力和应变成正比，是材料力学最基本的定律之一，一定要熟练掌握。 1.3 材料在拉压时的力学性能 材料的力学性能的研究是解决强度和刚度问题的一个重要方面。材料力学性能的研究一般是通过实验方法实现的，其中拉压试验是最主要、最基本的一种试验，由它所测定的材料性能指标有： E —材料抵抗弹性变形能力的指标；b s σσ,—材料的强度指标； ψδ, —材料的塑性指标。低碳钢的拉伸试验是一个典型的试验。
理正岩土使用手册-水力学
第一章 功能概述 理正工程水力学计算软件包含有五个计算内容：倒虹吸水力学计算、渠道水力学计算、水闸水力学计算、隧洞水力学计算和消能工水力学计算。 倒虹吸水力学计算模块可计算倒虹吸的过水能力、设计倒虹吸管径； 渠道水力学计算模块含有清水渠道均匀流的水力计算、清水渠道非均匀流的水力计算和挟沙水流渠道的水力计算； 水闸水力学计算模块适用于无坎宽顶堰、有坎宽顶堰、WES实用堰上的平板和弧形闸门，可计算水闸的泄流能力、设计闸孔宽度和确定闸门的开启度； 水工隧洞水力学计算模块适用于矩形、圆形、拱形断面隧洞的水力设计，对无压隧洞可计算洞的过流能力和设计断面尺寸，半有压隧洞可校核隧洞的过流能力，对于有压隧洞可计算隧洞在不同水位、不同闸门开度下的泄流量，并可在已知过流量条件下校核上游水位，还可绘制出总水头线和压坡线，形象的显示洞身各点有无负压； 消能工水力学计算模块适用于底流式消能工和挑流式消能工的水力设计。底流式消能工中包括下挖式消力池、突槛式消力池（消力墙）和综合式消力池三种基本型式，可进行消力池尺寸设计计算和校核消能能力。挑流式消能工可进行连续式挑流鼻坎的水力计算。 五个计算模块最后都给出计算的图形结果、文字结果及图文并茂的计算书。 第二章 快速操作指南 2.1 操作流程 理正工程水力学计算软件的操作流程如图2.1-1，每一步骤都有相对应的菜单操作。 图2.1-1 操作流程 2.2 快速操作指南
2.2.1 选择工作路径 设置工作路径，既可以调入已有的工作目录，也可在输入框中键入新的工作目录，后面操作中生成的所有文件（包括工程数据及计算书等）均保存在设置的工作目录下。 图2.2-1 指定工作路径 注意：此处指定的工作路径是所有岩土模块的工作路径。进入某单个计算模块后，还可以通过按钮【选工程】重新指定此模块的工作路径。 2.2.2 增加计算项目 工程水力学计算软件包含有五个计算内容：倒虹吸水力学计算、渠道水力学计算、水闸水力学计算、隧洞水力学计算和消能工水力学计算。用户可根据需要选择。 图2.2-2 当选好一个计算项目后，点击【工程操作】菜单中的“增加项目”或“增”按钮来新增一个计算项目（以水闸水力学计算为例）。
易飞成本计算步骤详解
成本会计操作手册 岗位说明：成本会计 岗位职责：负责成本计算数据的收集整理，检核并保证原始数据的正确性，计算每月的生产成本，做好成本分析的工作，向领导层及相关部门提供成本数据及改进建议 主要工作： 一、日常工作：成本计算的正确性，完全依赖各相关部门平时工作的正确与细致。因此，成本会计应负起督促与检核前端各部门的日常作业正确性的重大职责，主要包括以下内容： 1、日常存货异常之检视。 保证方法：避免存货异常的不二法门是保证日常交易数据的准确及时地录入，而这又完全仰仗所有操作人员按照流程的准确一致的行动。成本会计在平时应随时注意检查系统数据的准确与及时性。例如，成本会计应经常查询存货管理子系统中的“库存明细表”，查询是否有存货异常的情况（如库存数量为负数）。另外，成本会计也应不定期对重要的料件，进行帐实的核对。 2、所有生产性的领退料与生产入库必须在工单系统进行录入，并且必须设定核对工单。 保证方法：（1）在工单委外管理系统的设置单据性质中，将所有领退料单、生产入库单全设成核对工单。 （2）在库存管理中录入的库存交易单据，必须要求录入人员录入部门信息及在备注上写明用途，由成本会计随时抽查。 3、工单指定完工管理。 保证方法：已确定不再做的工单，应立即做指定完工的处理。可由生管人员在每月末月结前，通过工单委外管理子系统中的“工单明细表”，查询出未完工的工单，再逐一核对实际情况，进行指定完工的作业。 4、落实余料退库作业。 保证方法：指定完工的的工单，余料必须做退库的处理，同样可由生管通过工单委外管理子系统中的“工单明细表”，对指定完工的工单的用料情况进行核对，核对完后。由生管打印出当月未完工的“工单明细表”，在上面签字，传给成本会计，表示指定完工工单已处理完毕。 5、工单工时须正确。 保证方法：如果采用实际工时计算成本，车间管理人员必须保证每个工单工艺的工时录入正确，另外应注意：工时应是该工艺的有效工时，无效工时不应记录进去。 每月底，由生管或车间管理人员在工艺管理子系统中打印出“工作中心效率分析表”，将实际工时与标准工时进行核对，对于差异率较大的记录，应由车间管理人员仔细查核原因，并给出合理解释。 打印出的“工作中心效率分析表”由负责人员签字，转交给成本会计。 6、非工单用料之管理需明确。 保证方法：非生产性的用料，必须全部在库存管理子系统中用库存交易单录入，并注明领料部门及原因。该类的领退料单，每月末，由库管转交给成本会计，成本会计应检查单据是否符合规定。
易飞ERP月结和成本
目录： 存货子系统月结处理............................................................................... ２应收子系统月结 .................................................................................. １２应付子系统月结 .................................................................................. １３总账子系统月结 .................................................................................. １４自动分录 ........................................................................................... １５实际成本计算 ..................................................................................... １８
存货子系统月结处理月结流程图：
月结流程描述： 1、相关子系统单据完整性检查（销售、采购、库存、生产） 销售部门：应该检查系统中的销货单和销退单，查看这些单据的完整性：未审核、单据数量、价格有无错误以及单据中的其他信息，如果是未审核又不要的单据，应作“作废”处理。另外，对于不再执行发货的客户订单，应作“指定结束”处理。 采购部门：应该检查系统中的进货单和退货单，查看这些单据的完整性：未审核、单据数量、价格有无错误以及单据中的其他信息，如果是未审核又不要的单据，应作“作废”处理。另外，对于不再执行进货的请购单和采购单，应作“指定结束”处理。 重点注意是否有进货单上进货价格为0的,为防止成本计算不准确,必须由采购员提供实际价格。 仓管部门：应该检查系统中的库存交易单、调拨单和成本调整单，查看这些单据的完整性：未审核、单据数量、价格有无错误以及单据中的其他信息，如果是未审核又不要的单据，应作“作废”处理。如果是要抛转到总帐系统的单据,如制费领料单,维修领料单等还要检查单头“部门”字段是否已填入正确信息。 生产部门：应该检查系统中的领料单,退料单,生产入库单,委外进货单,委外退货单, 查看这些单据的完整性：未审核、单据数量、价格有无错误以及单据中的其他信息，如果是未审核又不要的单据，应作“作废”处理。另外,对于不再执行生产的工单,应作”指定结束”处理，对于不需要发放的LRP生产计划，要及时用“清除批次计划”来清除。 财务部门：视需要抽查上述部门的单据。 备注：ERP项目组或企管部应指定每月上述工作最终的完成确认时间。 2、重计现有库存 在各业务部门完成单据完整性检查后，执行“重计现有库存”作业，防止出现库存数量信息异常。 “重计年月”可以选择从结帐当月开始,也可以选择从”关帐年月”后的任意月份开始。 本作业的逻辑是从”维护品号每月统计信息”档中取得品号的月初数量和月初成本,从”维护库存交易明细”档取得品号当月的库存发生数量和成本,计算得到本月末和下月初的数量和成本,反复进行,一直到库存现行年月。 本作业要求在单一用户状态下进行,建议由系统管理员或财务部专岗（如材料会计）执行,不能随意随时执行。
剪切力的计算方法
第3章剪切和挤压的实用计算 3.1 剪切的概念 在工程实际中，经常遇到剪切问题。剪切变形的主要受力特点是构件受到与其轴线相垂直的大小相等、方向相反、作用线相距很近的一对外力的作用(图3-1a)，构件 m-面)发生相对错动(图3-1b)。的变形主要表现为沿着与外力作用线平行的剪切面(n 图3-1 工程中的一些联接件，如键、销钉、螺栓及铆钉等，都是主要承受剪切作用的构件。构件剪切面上的内力可用截面法求得。将构件沿剪切面n m-假想地截开，保留一部分考虑其平衡。例如，由左部分的平衡，可知剪切面上必有与外力平行且与横截面相切的内力Q F(图3-1c)的作用。Q F称为剪力，根据平衡方程∑=0 F Q=。 Y，可求得F 剪切破坏时，构件将沿剪切面(如图3-la所示的n m-面)被剪断。只有一个剪切面的情况，称为单剪切。图3-1a所示情况即为单剪切。 受剪构件除了承受剪切外，往往同时伴随着挤压、弯曲和拉伸等作用。在图3-1中没有完全给出构件所受的外力和剪切面上的全部内力，而只是给出了主要的受力和内力。实际受力和变形比较复杂，因而对这类构件的工作应力进行理论上的精确分析是困难的。工程中对这类构件的强度计算，一般采用在试验和经验基础上建立起来的比较简便的计算方法，称为剪切的实用计算或工程计算。 3.2 剪切和挤压的强度计算 3.2.1 剪切强度计算
剪切试验试件的受力情况应模拟零件的实际工作情况进行。图3-2a 为一种剪切试验装置的简图，试件的受力情况如图3-2b 所示，这是模拟某种销钉联接的工作情形。当载荷F 增大至破坏载荷b F 时，试件在剪切面m m -及n n -处被剪断。这种具有两个剪切面的情况，称为双剪切。由图3-2c 可求得剪切面上的剪力为 2 F F Q = 图3-2 由于受剪构件的变形及受力比较复杂，剪切面上的应力分布规律很难用理论方法确定，因而工程上一般采用实用计算方法来计算受剪构件的应力。在这种计算方法中，假设应力在剪切面内是均匀分布的。若以A 表示销钉横截面面积，则应力为 A F Q =τ (3-1) τ与剪切面相切故为切应力。以上计算是以假设“切应力在剪切面上均匀分布”为基础的，实际上它只是剪切面内的一个“平均切应力”，所以也称为名义切应力。 当F 达到b F 时的切应力称剪切极限应力，记为b τ。对于上述剪切试验，剪切极限应力为 A F b b 2= τ
数值分析计算方法试题集及答案
数值分析复习试题 第一章 绪论 一. 填空题 1.* x 为精确值 x 的近似值；() **x f y =为一元函数 ()x f y =1的近似值； ()**,*y x f y =为二元函数()y x f y ,2=的近似值，请写出下面的公式：**e x x =-： *** r x x e x -= ()()()*'1**y f x x εε≈? ()() () ()'***1**r r x f x y x f x εε≈ ? ()()()() ()* *,**,*2**f x y f x y y x y x y εεε??≈?+??? ()()()()() ** * *,***,**222r f x y e x f x y e y y x y y y ε??≈ ?+??? 2、 计算方法实际计算时，对数据只能取有限位表示，这时所产生的误差叫 舍入误 差 。 3、 分别用2.718281，2.718282作数e 的近似值，则其有效数字分别有 6 位和 7 位；又取 1.73≈-21 1.73 10 2 ≤?。 4、 设121.216, 3.654x x ==均具有3位有效数字，则12x x 的相对误差限为 0.0055 。 5、 设121.216, 3.654x x ==均具有3位有效数字，则12x x +的误差限为 0.01 。 6、 已知近似值 2.4560A x =是由真值T x 经四舍五入得 到,则相对误差限为 0.0000204 . 7、 递推公式,??? ? ?0n n-1y =y =10y -1,n =1,2, 如果取0 1.41y ≈作计算,则计算到10y 时,误 差为 81 10 2 ?;这个计算公式数值稳定不稳定 不稳定 . 8、 精确值 14159265.3* =π，则近似值141.3*1=π和1415.3*2=π分别有 3
剪切计算公式
2.剪切强度计算 (1) 剪切强度条件 剪切强度条件就是使构件的实际剪应力不超过材料的许用剪应力。 []s F A ττ= ≤ (5-6) 这里[τ]为许用剪应力，单价为Pa 或MPa 。 由于剪应力并非均匀分布，式(5-2)、(5-6)算出的只是剪切面上的平均剪应力，所以在使用实验的方式建立强度条件时，应使试件受力尽可能地接近实际联接件的情况，以确定试样失效时的极限载荷τ0，再除以安全系数n ，得许用剪应力[τ]。 []n ττ= (5-7) 各种材料的剪切许用应力应尽量从相关规范中查取。 一般来说，材料的剪切许用应力[τ]与材料的许用拉应力[σ]之间，存在如下关系： 对塑性材料： []0.60.8[]τσ= 对脆性材料： []0.8 1.0[]τσ= (2) 剪切实用计算 剪切计算相应地也可分为强度校核、截面设计、确定许可载荷等三类问题，这里就不展开论述了。但在剪切计算中要正确判断剪切面积，在铆钉联接中还要正确判断单剪切和双剪切。下面通过几个简单的例题来说明。 例5-1 图5-12(a)所示电瓶车挂钩中的销钉材料为20号钢，[τ]＝30MPa ，直径d=20mm 。挂钩及被连接板件的厚度分别为t ＝8mm 和t 1＝12mm 。牵引力F=15kN 。试校核销钉的剪切强度。 图5-12 电瓶车挂钩及其销钉受力分析示意图 解：销钉受力如图5-12(b)所示。根据受力情况，销钉中段相对于上、下两段沿m-m 和 n-n 两个面向左错动。所以有两个剪切面，是一个双剪切问题。由平衡方程容易求出： 2s F F = 销钉横截面上的剪应力为： 332151023.9MPa<[] 2(2010)4s F A ττπ-?===?? 故销钉满足剪切强度要求。 例5-2 如图5-13所示冲床，F max =400KN ，冲头[σ]=400MPa ，冲剪钢板的极限剪应力τb =360 MPa 。试设计冲头的最小直径及钢板最大厚度。
数值计算方法试题及答案
数值计算方法试题一 一、 填空题（每空1分，共17分） 1、如果用二分法求方程043=-+x x 在区间]2,1[内的根精确到三位小数，需对分（ ）次。 2、迭代格式 ) 2(2 1-+=+k k k x x x α局部收敛的充分条件是α取值在 （ ）。 3、已知?????≤≤+-+-+-≤≤=31)1()1()1(2110)(2 33x c x b x a x x x x S 是三次样条函数， 则 a =( )， b =（ ）， c =（ ）。 4、)(,),(),(10x l x l x l n Λ是以整数点n x x x ,,,10Λ为节点的Lagrange 插值基函数，则 ∑== n k k x l 0)(( )， ∑== n k k j k x l x 0 )(( )，当2≥n 时 = ++∑=)()3(20 4 x l x x k k n k k ( )。 5、设 1326)(247+++=x x x x f 和节点,,2,1,0,2/Λ==k k x k 则=],,,[10n x x x f Λ 和=?07f 。 6、5个节点的牛顿-柯特斯求积公式的代数精度为 ，5个节点的求积公式最高代数精度为 。 7、{}∞=0)(k k x ?是区间]1,0[上权函数x x =)(ρ的最高项系数为1的正交多项式族，其中1)(0=x ?，则?= 1 4)(dx x x ? 。
8、给定方程组?? ?=+-=-221121b x ax b ax x ，a 为实数，当a 满足 ，且 20<<ω时，SOR 迭代法收敛。 9、解初值问题00 (,)()y f x y y x y '=?? =?的改进欧拉法 ??? ??++=+=++++)],(),([2),(] 0[111] 0[1n n n n n n n n n n y x f y x f h y y y x hf y y 是 阶方法。 10、设 ?? ??? ?????=11001a a a a A ，当∈a （ ）时，必有分解式T LL A =，其中L 为下三角阵，当其对角线元素)3,2,1(=i l ii 满足（ ）条件时，这种分解是唯一的。 二、 二、选择题（每题2分） 1、解方程组b Ax =的简单迭代格式g Bx x k k +=+)() 1(收敛的充要条件是 （ ）。 （1）1)(A ρ, (4) 1)(>B ρ 2、在牛顿-柯特斯求积公式： ?∑=-≈b a n i i n i x f C a b dx x f 0 )() ()()(中，当系数 ) (n i C 是负值时，公式的稳定性不能保证，所以实际应用中，当（ ）时的牛顿-柯特斯求积公式不使用。 （1）8≥n ， （2）7≥n ， （3）10≥n ， （4）6≥n ， 3、有下列数表