宿迁学院课程设计说明书设计题目:基于VB数字积分法一、二象限逆圆插补计算与仿真系(部)：机电工程系专业:自动化(数控技术）班级: ０9数控本(２)姓名：李伟学号: 20０9０7052指导老师（签名）:刘萍起止时间:20 1２年12月3日至2０1２年１2月7 日共１周20 １２年12 月6 日正文:数控原理与系统课程设计说明书一、课程设计的目的1）了解连续轨迹控制数控系统的组成原理。
2) 掌握逐点比较法插补的基本原理。
３）握逐点比较法插补的软件实现方法。
二．课程设计的任务逐点比较法插补是最简单的脉冲增量式插补算法之一,其过程清晰，速度平稳,但一般只用于一个平面内两个坐标轴的插补运算。
其基本原理是在刀具按要求轨迹运动加工零件轮廓的过程中，不断比较刀具与被加工零件轮廓之间的相对位置,并根据比较结果决定下一步的进给方向,使刀具向减小偏差的方向进给,且只有一个方向的进给。
也就是说,逐点比较法每一步均要比较加工点瞬时坐标与规定零件轮廓之间的距离,依此决定下一步的走向。
如果加工点走到轮廓外面去了，则下一步要朝着轮廓内部走;如果加工点处在轮廓的内部，则下一步要向轮廓外面走,以缩小偏差，这样周而复始，直至全部结束,从而获得一个非常接近于数控加工程序规定轮廓的轨迹。
逐点比较法插补过程中的每进给一步都要经过偏差判别、坐标进给、偏差计算和终点判别四个节拍的处理,其工作流程图如图所示。
三.逐点比较法基本原理逐点比较法I象限逆圆插补在加工圆弧过程中，人们很容易联想到使用动点到圆心的距离与该圆弧的名义半径进行比较来反映加工偏差。
假设被加工零件的轮廓为第Ⅰ象限逆走向圆弧SE,,圆心在O(０，０）,半径为Ｒ，起点为S(ＸS,YS ）,终点为Ｅ（Xe,Ye),圆弧上任意加工动点为N （X ｉ,Yi)。
当比较该加工动点到圆心的距离ON 与圆弧半径Ｒ的大小时，可获得刀具与圆弧轮廓之间的相对位置关系。
当动点Ｎ(Ｘｉ,Yi ）正好落在圆弧上时,则有下式成立22222R Y X Y X e e i i =+=+当动点N(Xi ，Y ｉ)落在圆弧外侧时,则有下式成立22222R Y X Y X e e i i =+>+当动点N （X ｉ,Y ｉ）落在圆弧内侧时,则有下式成立22222R Y X Y X e e i i =+<+由此可见，取逐点比较法圆弧插补的偏差函数表达式为222R Y X F i i -+=当动点落在圆外时，为了减小加工误差，应向圆内进给，即向(-X ）轴方向走一步;当动点落在圆内时,应向圆外进给，即向(+Y ）轴方向走一步。
当动点正好落在圆弧上且尚未到达终点时,为了使加工继续下去，理论上向(+Y)轴或(－X)轴方向进给均可以,但一般情况下约定向(-X ）轴方向进给。
综上所述,现将逐点比较法第Ⅰ象限逆圆插补规则概括如下：当F>0时，即222R Y X F i i -+=>０，动点落在圆外,则向(－X ）轴方向进给一步；当F=0时，即222R Y X F i i -+=＝0，动点正好落在圆上,约定向（-Ｘ)轴方向进给一步;当F ＜０时，即222R Y X F i i -+=＜0,动点落在圆内，则向(+Ｙ)轴方向进给一步。
由偏差函数表达式可知,计算偏差Ｆ值,就必须进行动点坐标、圆弧半径的平方运算。
显然，在用硬件或汇编语言实现时不太方便。
为了简化这些计算,按逐点比较法直线插补的思路，也可以推导出逐点比较法圆弧插补过程中偏差函数计算的递推公式。
假设第ｉ次插补后，动点坐标为N （Xi,Ｙi ），其对应偏差函数为222R Y X F i i i -+=当Ｆi ≥0，向(-X)轴方向进给一步，则新的动点坐标值为 X ｉ＋１=Xi-1, Yi ＋1=Yi 因此，新的偏差函数为()2222212111R Y X R Y X F i i i i i -+-=-+=+++∴ ﻩﻩ Fi+1＝Ｆi －2X ｉ+1 同理，当F ｉ<0,则向(＋Y)轴方向进给一步，则新的动点坐标值为 X ｉ+1＝Xi ， Yi ＋１=Ｙi ＋1 因此，可求得新的偏差函数为()2222212111R Y X R Y X F i i i i i -++=-+=+++∴ Fi+1＝Fi+2Yi ＋1将上式进行比较,可以看出两点不同:第一,递推形式的偏差计算公式中仅有加/减法以及乘２运算，而乘2可等效成该二进制数左移一位,这显然比平方运算来得简单。
第二，进给后新的偏差函数值与前一点的偏差值以及动点坐标N(Xi,Yi)均有关系。
由于动点坐标值随着插补过程的进行而不断变化,因此,每插补一次,动点坐标就必须修正一次,以便为下一步的偏差计算作好准备。
至此,将第Ⅰ象限逆圆弧插补的规则和计算公式汇总,见表（表2-1)第Ⅰ象限逆圆弧插补计算公式表2-1和直线插补一样，圆弧插补过程也有终点判别问题。
当圆弧轮廓仅在一个象限区域内，其终点判别仍可借用直线终点判别的三种方法进行,只是计算公式略不同。
Σ=｜Xe-Xｓ｜+|Ｙe－Ｙs|ﻩΣ=mａx｛｜Xｅ－Xs｜，|Ye-Yｓ|｝Σ1=|Ｘe－Xs| ，Σ2=｜Ye-Yｓ|式中ＸＳ、Ys ——被插补圆弧轮廓的起点坐标;Xe、Ye ——被插补圆弧轮廓的终点坐标。
ｂ、插补象限和圆弧走向前面所讨论的逐点比较法直线和圆弧插补，均是针对第一象限直线和逆圆插补这种特定情况进行的。
然而,任何数控机床都应具备处理不同象限、不同走向曲线的能力。
四个象限中圆弧插补圆弧插补情况比直线插补复杂，不仅有象限问题，而且还有圆弧走向问题。
现以第Ⅰ象限顺圆ＳＲ１插补为例，介绍圆弧插补的特性。
假设圆弧SE 起点为S （ＸS ，YS)，终点为Ｅ（Xe ，Ｙe)，圆心在坐标原点上。
与逆圆插补相似，当某一时刻动点Ｎ（Xi ，Yi ）在圆弧的外侧时，有F ≥0成立,应向－Y)轴方向进给一步,以减小误差;若动点N(Xi,Yi ）在圆弧内侧，则应向＋X ）轴方向进给一步。
由此可推导出第Ⅰ象限顺圆插补偏差函数的递推公式如下:当Fi ≥0时,向-Y ）轴方向进给一步,则新的动点坐标为Xi+1=Xi , Yi ＋1=Yi-1 新动点的偏差函数为()2222212111R Y X R Y X F i i i i i --+=-+=+++∴ Fi ＋1=Ｆi －2Yi ＋１ 当Fi ＜0时,向＋X)轴方向进给一步，则新的动点坐标为Xi ＋1＝Ｘi ＋1， Yi ＋1＝Ｙi新动点的偏差函数为()2222212111R Y X R Y X F i i i i i -++=-+=+++∴ Fi ＋1＝Ｆi+2X ｉ+１现将上式比较,可以看出它们有两点不同： ﻩ1）当Fi ≥0或Fi<０时，对应的进给方向不同;2）插补计算公式中动点坐标的修正也不同,以至于偏差计算公式也不相同。
进一步还可根据上述方法推导出其他象限不同走向圆弧的插补公式。
现将各种相应偏差计算见表(表2-2)四个象限圆弧插补偏差计算与进给方向表2-2从表可以看出,当按第Ⅰ象限逆圆NR1进行插补运算时，现若有意将X轴进给反向,则可以走出第Ⅱ象限顺圆ＳＲ2来;或者若将Y轴进给反向，则可以走出SR4来；或者将X轴和Ｙ轴的进给均反向,则可以走出NR3来;并且这四种线型（NR1、SＲ2、NR3、SＲ4)使用的偏差计算公式都相同，无须改变。
进一步还可以看出，当按第Ⅰ象限逆圆NＲ1线型插补时，现若将计算公式坐标X与Y对调，即把X当作Y，把Y当作X,那么就可得到SＲ1的走向。
类似地通过改变进给方向，利用SR1的公式就可获得其余三种线型(NＲ2、ＳR3、NR４）的走向。
下面,我们对圆弧逐点比较法作一个简要的介绍。
四.算法描述(逐点比较法在ＶB中的具体实现)根据上述基本原理，我们可以知道逐点比较法圆弧插补需要设置两个终点计器J∑X=|Ｘe – Xs|和J∑Y=｜Yｅ－Ys|,分别对X轴和Ｙ轴进行终点监控。
每当X轴或Y轴产生一个溢出脉冲,相应的终点计数器就作减1修正，直到为零,表明该坐标已到终点，并停止其坐标的累加运算。
只有当两个坐标轴均到达终点时，圆弧插补才结束。
如下图所示,圆弧起点Ｓ（4，-1）,终点（－1,４),且寄存器位数N＝3,当插补开始时，被积函数寄存器初值分别为ＪVＸ=Ｙs=-1和JＶY＝Xｓ＝4,终点判别寄存器J∑Ｘ＝｜Xｅ–Xs|=-１和J∑Y=|Yｅ－Ys｜=5.该圆弧插补运算过程如下表所示,插补轨迹如下图折线所示。
根据上述基本原理以及对预ＶB软件编程的掌握,查阅参考相关文献资料,编制了本课程设计的程序。
Prｉvaｔe Sub Ｃommand1_Cｌｉck()ＤimｒAs Ｓingｌｅ,ＸＡｓIｎtｅger, YＡs Integeｒ， c AｓＩnｔeｇeｒ, d AｓIｎｔeger, a AｓＳｉngle，b As Sinｇｌｅ, pi As Ｓingle，mＡs Iｎtegeｒ，n As ＩntegerX= Int(Teｘｔ1)Ｙ= Inｔ(Text2)ｃ= Ｉnｔ(Ｔｅxt3)d = Int(Texｔ4)ａ＝Atｎ(Aｂｓ(Y /X)）ｂ= Atn(Ａbs(d / c))pi = ３．14１５９2６5ｍ= Ｘ* X + Y* Yn = c * c +ｄ* dＰictｕre1．ＦoreCｏlｏr ＝vbBｌacｋPictｕre１.DｒawWiｄth = ２IｆＯｐtion1．Vａlue = Ｔrｕe ThｅnＩf X <= 0 OrＹ< ０Oｒ c <= 0 Ｏr ｄ< 0 Or m ＜＞n TheｎPrint MsgＢｏx("非法数据！＂，48, "提示信息")ElsｅPiｃture1.Linｅ（5０0，5００）－（５00, ６500）Picｔｕre1.Lｉne(５00, 6500)-(5５０0, 6５０0)Picture1.Ｌｉne (5４00, 6５30)-(5５００,６5０0)Piｃｔure1.Line（5400，64７0)-（５500, 6500)Ｐicturｅ1.Ｌine (5３0,60０）-(５00，500)Pｉｃｔｕre1.Liｎe (470, 60０)-(5０0，5００)Pictｕre１.FｏrｅＣolor= vbRｅdPiｃture1.DｒawWidth ＝5ｒ=Sqｒ（X *X ＋Y * Ｙ) * 400Ｐｉcturｅ1.Cirｃle （500, 650０), r, , ａ, bEnｄIfElsｅIf Option２．Vaｌuｅ＝True TheｎIｆX ＞= 0 Or Y < ０Or c >= ０Or d＜０Or m <＞n TheｎＰrinｔＭsｇＢoｘ("非法数据!＂, ４8, "提示信息")EｌsePicｔｕre1.Lｉｎｅ(550０, 50０）-(5500, ６500)Picture１.Line(５０0, ６5０0）－(５7０0,65０0)Ｐicｔurｅ1.Line (５６00，6530）-（５7０0, 6500)Ｐicｔure１．Line (5６０0, 6470)-(５700, 6500）Ｐｉcｔure1．Linｅ(553０，600)-(5５０0, ５00)Pｉｃｔure1．Lｉne（5470,600)-(5500, 500)Picｔurｅ1.FｏｒeColor = vbReｄPictuｒe1.ＤｒａwWｉdtｈ= 5r = Ｓqr(X ＊X+Y *Y) ＊4０0Ｐiｃture1.Cｉrcle (5500，650０)，r,,ｐi－a,pi－ｂEnd ＩｆEnd ＩｆＥnd SubPｒｉvatｅＳｕb Comｍanｄ2_Clｉck()Ｄim k，ｍ, j，l,n, F(3０), X（30）, Y(３0）As Integer，a As Ｉnteｇｅr，b As Integer, c As Inｔeger，ｄＡs Inｔegera =Int(Text1)b = Ｉnt(Tｅｘｔ2)c = Ｉnt（Texｔ3)d ＝Ｉnt（Ｔeｘt4)m ＝0l = 0ｋ= ０F（m)＝0Ｘ(m)= ａY(m) ＝ bPｉcture1.ForeColor = vbGｒeｅnPiｃture1．DrawWidth = 3ｊ= Abs(Abs（ａ) - Aｂs(c)）＋Ａbs(Abs(b)-Ａbs（ｄ)）Fｏrｍ1.ＣurreｎtX = 200Forｍ1．CｕrreｎtＹ= 200Pｒint "初始＂, "进给方向"，＂F(0)=０＂,＂X(0) ＝" ＆Ｉnt(Tｅxt1), "Y(0)=" & Iｎｔ(Texｔ２）, " Xe ＝"& Ｉnｔ(Ｔeｘｔ３), " Ｙe ＝" & Int(Texｔ4)，" ∑＝" & jIｆOｐtｉｏｎ1.Ｖalue ＝Ｔｒue ThenFor n = 1 TｏjIｆＦ(m) >= 0Aｎdｊ＞０Ｔhenｍ= m + 1l = l＋1F(m) ＝Ｆ(m- 1) － 2 * Abs(Ｘ（m - 1）) + 1X（m) ＝X(m - 1)- 1Y(m) = Y（m- 1)Picｔure1.Linｅ（500 + 400 *(a- ｌ), 6５00 - 4０0 * (ｂ+k))－(50０+ ４00 ＊(a - ｌ＋1), ６500 -400 *(b +ｋ))Form1．ＣuｒrｅntX = 2０0Form1．CurrｅntＹ= 20０+ m* 30０Ｐｒｉnt "第" &ｍ&"步＂,"-△X＂, "Ｆ(" ＆m & ＂)＝"& F(m),＂X（" & m & ")=" & X（m）, "Ｙ(" & m &")=" &Y（ｍ), " Ｘe = " &Inｔ(Text３），" Ye = " ＆Int(Text４)," ∑＝＂＆j - nElsek＝k＋1m= m +１F(ｍ) = Ｆ(m －1) + ２* Aｂs(Y(m - 1)）+１Y(m) = Y(m - 1) + 1X（ｍ) ＝X(ｍ－1)Ｐictｕrｅ1．Line （50０＋4０0 *（ａ- l）,650０- 4０0 *(b + k）)-(50０+ 40０* (a－ｌ),6500 - ４00*(b +k - 1))Forｍ1.CurｒenｔX = 200Forｍ1.CurｒentY = 200+m * 300Print "第" & m&"步＂, "＋△Y ", "F(" & ｍ& ")="& Ｆ(m）," X（＂& ｍ＆")＝" & Ｘ(m）," Ｙ（" & ｍ＆"）=" & Ｙ（m）," Xe ="& Iｎｔ(Text ３), ＂Ye =" &Iｎt(Teｘt4)，"∑= " & j -n;"＂EnｄIfNext nＥlsｅＩf Optiｏn2.Vａlue ＝ＴｒｕｅThenＦoｒn ＝1To jＩf F(m) >＝0And j ＞0 Tｈｅｎm =m＋１ｋ=k + 1F(ｍ）＝Ｆ(m- １) － 2 * Ａbs（Y（ｍ- 1)）+ 1Ｘ(ｍ) = X(m- 1)Y(m）= Ｙ(m - 1）- 1Pictｕre1.Line (５500 +40０*(a - l）,65００-４00 * (b - k ＋1))-(5５00 ＋４00 * (a -l),６500 - ４00* (b - ｋ))Forｍ１.CurrenｔX = 200Ｆorm1.CurreｎtY ＝2００+ m * 300Prｉnt ＂第＂& ｍ& "步", " -△Ｙ","F(" ＆m & ＂)=" &F(m), ＂Ｘ（" ＆ｍ＆")=" &X(m)，" Ｙ(" &m & "）=＂＆Y（m), ＂Xe =＂& Int(Teｘt3)，＂Ye=" & Int(Teｘｔ4), " ∑＝" & ｊ- nEｌsel =ｌ+ 1m= m ＋1F（m) = F(m －1) + 2 * Abｓ（X(m -1)) + １Ｙ(ｍ）=Y(m - 1)Ｘ(m) =X(m - １) －１Pｉctｕｒe1.Linｅ(55０0＋400 ＊(a －l）, 6500- 400 * （b - ｋ）)-(5500＋4０0 * (ａ- l ＋１), 6500 - ４００*（ｂ－k))Form1.CuｒｒeｎｔX = 200Form1.CuｒrｅntY = 2００＋ｍ* 300Print＂第" & ｍ＆"步", ＂-△X ", "Ｆ(" & m &")=" & F（m）, " X(＂＆ｍ& "）=" & X(m), " Ｙ(" & ｍ& ")="& Y(m)，＂Ｘｅ= " & Int(Ｔext3), "Ｙe = "＆Inｔ(Text4），" ∑＝＂＆j - nEnd IｆNext ｎEnd IfEnd SubPriｖａte Sub Command3_Ｃｌick()Tｅxt1.Texｔ= ""Texｔ２.Tｅxt =""Ｔexｔ３.Teｘｔ＝""Ｔｅｘt4.Text= ""Ｐiｃture1.ClsFoｒｍ1.ClｓEnｄSuｂ六、各象限选例展示一．参考文献【1】汪木兰主编.数控原理与系统.机械工业出版社，2０04【2】李金泽主编.基于ViｓuaｌＢａｓic编程软件的数控插补计算与设计方法.机械工业出版社,2０0７年。