基于准均匀三次B样条曲线的S试件建模仿真加工.docx

上传人:b****6 文档编号:5251613 上传时间:2022-12-14 格式:DOCX 页数:13 大小:790.54KB
下载 相关 举报
基于准均匀三次B样条曲线的S试件建模仿真加工.docx_第1页
第1页 / 共13页
基于准均匀三次B样条曲线的S试件建模仿真加工.docx_第2页
第2页 / 共13页
基于准均匀三次B样条曲线的S试件建模仿真加工.docx_第3页
第3页 / 共13页
基于准均匀三次B样条曲线的S试件建模仿真加工.docx_第4页
第4页 / 共13页
基于准均匀三次B样条曲线的S试件建模仿真加工.docx_第5页
第5页 / 共13页
点击查看更多>>
下载资源
资源描述

基于准均匀三次B样条曲线的S试件建模仿真加工.docx

《基于准均匀三次B样条曲线的S试件建模仿真加工.docx》由会员分享,可在线阅读,更多相关《基于准均匀三次B样条曲线的S试件建模仿真加工.docx(13页珍藏版)》请在冰豆网上搜索。

基于准均匀三次B样条曲线的S试件建模仿真加工.docx

基于准均匀三次B样条曲线的S试件建模仿真加工

基于准均匀三次B样条曲线的“S”试件建模仿真加工

1.背景与需求

国家在高端装备制造业“十二五”发展规划中明确提出了让智能制造装备产业基本满足国民经济重点领域和国防建设需求的发展目标。

近年来,数控机床,特别是高速五轴联动数控机床越来越多的应用在航空、航天、汽车、轨道交通、海洋工程、模具成型、化工原材料加工等领域。

相比普通机床,五轴数控机床具有非正交运动轴,可以连续改变刀具的空间姿态,能够加工自由曲面,具有加工精度高、零件加工一致性好、生产效率高等特点,适合加工复杂异形零件。

在一些关键、复杂的零部件加工中,常常对零件的尺寸精度提出较高要求,尤其在数控机床的高速联动加工中,不容易保证其加工精度。

由于现阶段高端制造装备的核心技术大多由发达国家所垄断,出于政治、军事、经济的战略角度考虑,其研发、制造以及保障技术一直被严格控制。

我国正处于加快转变经济发展方式、实现由制造业大国向强国转变的工业转型升级期,国产高档数控机床的性能严重制约着我国高端装备制造业发展。

数控机床在出厂安装、调试以及后续的生产、保养、维修中都需要着重关注其加工性能,尤其是其加工精度。

机床的各运动传动部件以及控制环节产生的误差会逐渐累积,最终反映在零件的加工误差上。

由于多轴数控机床在结构上的复杂性和多样性,现有的国内外机床检验通则和规范都还仅限于三轴数控机床的静态精度和低速加工性能,满足不了高速五轴数控机床的检验要求。

当前,在我国数控机床市场中,高档数控机床,特别是五坐标联动、大型数控机床的市场份额为国外品牌所占领。

机床工业的现状致使我国航空结构件数控加工装备长期以来以进口国外设备为主,尤其是五轴加工等高技术关键装备接近100%依赖进口。

中航工业成都飞机工业(集团)有限责任公司基于航空结构件的薄壁特性,提出了一种“S”形机床检验试件,如图1所示。

由沈阳机床(集团)有限责任公司、天津大学、成都飞机工业(集团)有限责任公司产学研用各方联合实施,研制出了国产首台用于复杂钛合金航空结构件加工的五轴联动加工中心。

同时,通过“S”试件切削工艺技术研究,完成了“S”试件试切实验。

该试件在实际生产中能够有效的检验五轴数控机床对航空结构件的加工性能。

通过三坐标测量机对“S”形检验试件的轮廓误差进行测量,评判试件的加工精度,从而判断机床的加工性能。

图1“S”形机床检验试件

 

2.专题作业

2.1建模过程

第一步:

通过翻阅施法中老师的《计算机辅助几何设计与非均匀有理B样条》,得知准均匀三次B样条曲线方程;通过matlab编程,将4组控制顶点拟合出4条B样条曲线,如图2所示。

由于每条B样条共有12个控制顶点,共能拟合9段曲线,通过给定参量步长,每段得到曲线上的101个点,从一百个点中的前、中、后分别选取一个点,又由于每段最后的一个点与下一段的第一个点相同,因此每条B样条共采集19个点。

相关程序见附录1。

图2

第二步:

将每条B样条的19个顶点导入到UG中,生成出4条B样条曲线,绘制出两个直纹面,并构建“S”形实体模型,如图3-4所示;

图3

图4

第三步:

根据毛坯尺寸,构建毛坯实体模型,如图5所示。

并将加工后的“S”形实体模型同毛坯模型组合在一起,如图6所示。

建模完成,等待进一步的仿真加工。

图5

图6

2.2工艺方案

2.2.1切削刀具和切削参数

切削刀具和对应的切削参数由机床厂家和用户协商确定并作出记录。

选用

下面为建议采用的刀具和切削参数。

切削刀具:

-刀具类型:

立铣刀

-20mm刀具直径:

20毫米

-刀具工作长度:

不小于50毫米

切削刀具参数设置如图7所示。

图7

切削参数:

推荐的切削深度为5-10毫米。

针对不同类型的机床的进给速度要求推荐如下:

-主轴转速大于18000转/分:

5000毫米/分

-主轴转速大于等于10000转/分并小于18000转/分:

3000毫米/分

-主轴转速大于等于6000转/分并小于10000转/分:

2000毫米/分

-主轴转速小于6000转/分:

1000毫米/分

主轴转速选取8000转/分,进给率为2000毫米/分,如图8所示。

图8

2.2.2加工过程

加工坐标系如图9所示

图9加工坐标系

步骤1:

直纹面的粗加工。

粗加工采用五轴侧铣加工,保证S形曲面每边各留1mm的加工余量。

步骤2:

直纹面的半精加工。

半精加工采用五轴侧铣加工,保证S形曲面每边各留0.8mm的加工余量,如图10所示。

步骤3:

直纹面的精加工。

精加工采用五轴侧铣加工来获得最终的S件曲面。

图10半精加工余量设计参数

2.3数控编程及仿真

仿真

粗加工轨迹生成如图11所示:

图11

半精加工轨迹生成图12所示:

图12

精加工轨迹生成图13所示:

图13

仿真界面如图14所示:

图14

仿真过程如图15-16所示:

图15

图16

数控编程

粗加工程序、半精加工程序及精加工程序,见附录2。

3课程小结

本实验课程主要有两部分组成,分别为理论部分和实践部分。

理论部分主要以老师上课讲授并操作示范为主。

可是由于出差在外,我错过了仅有的一节宝贵的理论课。

实践部分,则要求我们按照要求建立一个“S”试件模型,并在给定毛坯参数的情况下,建立毛坯模型,然后用软件对其仿真,生成刀轨和G代码,加工出“S”试件。

由于自己之前对CAD/CAM软件学习不深,用的不透,所以几乎是从零开始学习,遇到了很多的问题。

通过翻阅相关书籍、看视频以及请教别人,我一点点地建立起了模型,并学着去进行仿真加工,生成刀轨和G代码,终于完成了此次的作业。

虽然很艰难,但是收获颇多。

在这里,非常感谢帮助我的同学师姐和老师。

附录1

Main1:

clear;

clc;

loadctrlpz1.txt;%¶ÁÈ¡textÎļþ¡°¿ØÖƶ¥µã¡±Êý¾Ý

D=ctrlpz1;%½«textÎļþctrlptsµÄÊý¾Ý¸³¸øDÐÐÏòÁ¿

n=length(D);%Çó³öDÐÐÏòÁ¿ÔªËظöÊý

k=n/3;%DÐÐÏòÁ¿ÔªËØÈý¸öÊý¾Ý×é³ÉÒ»¸öÈýά¿ØÖƶ¥µã

p=zeros(3,k);%¶¨ÒåpΪ3ÐÐkÁÐÏòÁ¿£¬Òò´ËÒ»Áдú±íÒ»¸ö¿ØÖƶ¥µã

fori=1:

k

p(:

i)=D(1,3*i-2:

3*i);%¶ÔpÏòÁ¿½øÐи³Öµ

end

LyCAD(p);%µ÷Óú¯Êý

holdon;

Main2:

clear;

clc;

loadctrlpz2.txt;%¶ÁÈ¡textÎļþ¡°¿ØÖƶ¥µã¡±Êý¾Ý

D=ctrlpz2;%½«textÎļþctrlptsµÄÊý¾Ý¸³¸øDÐÐÏòÁ¿

n=length(D);%Çó³öDÐÐÏòÁ¿ÔªËظöÊý

k=n/3;%DÐÐÏòÁ¿ÔªËØÈý¸öÊý¾Ý×é³ÉÒ»¸öÈýά¿ØÖƶ¥µã

p=zeros(3,k);%¶¨ÒåpΪ3ÐÐkÁÐÏòÁ¿£¬Òò´ËÒ»Áдú±íÒ»¸ö¿ØÖƶ¥µã

fori=1:

k

p(:

i)=D(1,3*i-2:

3*i);%¶ÔpÏòÁ¿½øÐи³Öµ

end

LyCAD(p);%µ÷Óú¯Êý

holdon;

FUNCTION:

functionLyCAD(p)

t=0:

0.01:

1;%²ÎÊýtµÄ²½³¤Îª0.001£¬¾ö¶¨ÇúÏ߹⻬¶È

n=size(p,2);%½«ÏòÁ¿pµÄÐÐÊýºÍÁÐÊý·Ö±ð¸³¸ønºÍm

fori=1:

(n/2-3)

ifi==1

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½x²¿·Ö

x1=p(1,i)*(-t.^3+3*t.^2-3*t+1)+p(1,i+1)*(7/4*t.^3-9/2*t.^2+3*t)...

+p(1,i+2)*(-11/12*t.^3+3/2*t.^2)+p(1,i+3)*1/6*t.^3;

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½y²¿·Ö

y1=p(2,i)*(-t.^3+3*t.^2-3*t+1)+p(2,i+1)*(7/4*t.^3-9/2*t.^2+3*t)...

+p(2,i+2)*(-11/12*t.^3+3/2*t.^2)+p(2,i+3)*1/6*t.^3;

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½z²¿·Ö

z1=p(3,i)*(-t.^3+3*t.^2-3*t+1)+p(3,i+1)*(7/4*t.^3-9/2*t.^2+3*t)...

+p(3,i+2)*(-11/12*t.^3+3/2*t.^2)+p(3,i+3)*1/6*t.^3;

%»æÖÆ3´Î×¼¾ùÔÈBÑùÌõÇúÏß

end

ifi==2

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½x²¿·Ö

x1=p(1,i)*(-1/4*t.^3+3/4*t.^2-3/4*t+1/4)+p(1,i+1)*(7/12*t.^3-4/5*t.^2+1/4*t+7/12)...

+p(1,i+2)*(-1/2*t.^3+1/2*t.^2+1/2*t+1/6)+p(1,i+3)*1/6*t.^3;

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½y²¿·Ö

y1=p(2,i)*(-1/4*t.^3+3/4*t.^2-3/4*t+1/4)+p(2,i+1)*(7/12*t.^3-4/5*t.^2+1/4*t+7/12)...

+p(2,i+2)*(-1/2*t.^3+1/2*t.^2+1/2*t+1/6)+p(2,i+3)*1/6*t.^3;

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½z²¿·Ö

z1=p(3,i)*(-1/4*t.^3+3/4*t.^2-3/4*t+1/4)+p(3,i+1)*(7/12*t.^3-4/5*t.^2+1/4*t+7/12)...

+p(3,i+2)*(-1/2*t.^3+1/2*t.^2+1/2*t+1/6)+p(3,i+3)*1/6*t.^3;

end

if2

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½x²¿·Ö

x1=p(1,i)*(-1/6*t.^3+1/2*t.^2-1/2*t+1/6)+p(1,i+1)*(1/2*t.^3-1*t.^2+2/3)...

+p(1,i+2)*(-1/2*t.^3+1/2*t.^2+1/2*t+1/6)+p(1,i+3)*1/6*t.^3;

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½y²¿·Ö

y1=p(2,i)*(-1/6*t.^3+1/2*t.^2-1/2*t+1/6)+p(2,i+1)*(1/2*t.^3-1*t.^2+2/3)...

+p(2,i+2)*(-1/2*t.^3+1/2*t.^2+1/2*t+1/6)+p(2,i+3)*1/6*t.^3;

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½z²¿·Ö

z1=p(3,i)*(-1/6*t.^3+1/2*t.^2-1/2*t+1/6)+p(3,i+1)*(1/2*t.^3-1*t.^2+2/3)...

+p(3,i+2)*(-1/2*t.^3+1/2*t.^2+1/2*t+1/6)+p(3,i+3)*1/6*t.^3;

end

%»æÖÆ3´Î×¼¾ùÔÈBÑùÌõÇúÏß

plot3(x1,y1,z1,'r');

gridon;

axissquare;

holdon;

end

fori=(n/2-3)+4:

n-3

ifi==1

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½x²¿·Ö

x1=p(1,i)*(-t.^3+3*t.^2-3*t+1)+p(1,i+1)*(7/4*t.^3-9/2*t.^2+3*t)...

+p(1,i+2)*(-11/12*t.^3+3/2*t.^2)+p(1,i+3)*1/6*t.^3;

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½y²¿·Ö

y1=p(2,i)*(-t.^3+3*t.^2-3*t+1)+p(2,i+1)*(7/4*t.^3-9/2*t.^2+3*t)...

+p(2,i+2)*(-11/12*t.^3+3/2*t.^2)+p(2,i+3)*1/6*t.^3;

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½z²¿·Ö

z1=p(3,i)*(-t.^3+3*t.^2-3*t+1)+p(3,i+1)*(7/4*t.^3-9/2*t.^2+3*t)...

+p(3,i+2)*(-11/12*t.^3+3/2*t.^2)+p(3,i+3)*1/6*t.^3;

%»æÖÆ3´Î×¼¾ùÔÈBÑùÌõÇúÏß

end

ifi==2

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½x²¿·Ö

x1=p(1,i)*(-1/4*t.^3+3/4*t.^2-3/4*t+1/4)+p(1,i+1)*(7/12*t.^3-4/5*t.^2+1/4*t+7/12)...

+p(1,i+2)*(-1/2*t.^3+1/2*t.^2+1/2*t+1/6)+p(1,i+3)*1/6*t.^3;

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½y²¿·Ö

y1=p(2,i)*(-1/4*t.^3+3/4*t.^2-3/4*t+1/4)+p(2,i+1)*(7/12*t.^3-4/5*t.^2+1/4*t+7/12)...

+p(2,i+2)*(-1/2*t.^3+1/2*t.^2+1/2*t+1/6)+p(2,i+3)*1/6*t.^3;

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½z²¿·Ö

z1=p(3,i)*(-1/4*t.^3+3/4*t.^2-3/4*t+1/4)+p(3,i+1)*(7/12*t.^3-4/5*t.^2+1/4*t+7/12)...

+p(3,i+2)*(-1/2*t.^3+1/2*t.^2+1/2*t+1/6)+p(3,i+3)*1/6*t.^3;

end

if2

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½x²¿·Ö

x1=p(1,i)*(-1/6*t.^3+1/2*t.^2-1/2*t+1/6)+p(1,i+1)*(1/2*t.^3-1*t.^2+2/3)...

+p(1,i+2)*(-1/2*t.^3+1/2*t.^2+1/2*t+1/6)+p(1,i+3)*1/6*t.^3;

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½y²¿·Ö

y1=p(2,i)*(-1/6*t.^3+1/2*t.^2-1/2*t+1/6)+p(2,i+1)*(1/2*t.^3-1*t.^2+2/3)...

+p(2,i+2)*(-1/2*t.^3+1/2*t.^2+1/2*t+1/6)+p(2,i+3)*1/6*t.^3;

%3´Î×¼¾ùÔÈBÑùÌõ¹«Ê½z²¿·Ö

z1=p(3,i)*(-1/6*t.^3+1/2*t.^2-1/2*t+1/6)+p(3,i+1)*(1/2*t.^3-1*t.^2+2/3)...

+p(3,i+2)*(-1/2*t.^3+1/2*t.^2+1/2*t+1/6)+p(3,i+3)*1/6*t.^3;

end

%»æÖÆ3´Î×¼¾ùÔÈBÑùÌõÇúÏß

plot3(x1,y1,z1,'r');

gridon;

axissquare;

holdon;

end

forj=1:

n/2

plot3(p(1,j),p(2,j),p(3,j),'g+');%»æÖÆ¿ØÖƶ¥µã

ifj>1

plot3([p(1,j-1),p(1,j)],[p(2,j-1),p(2,j)],[p(3,j-1),p(3,j)]);%Á¬½Ó¿ØÖƶ¥µã

end

end

forj=n/2+1:

n

plot3(p(1,j),p(2,j),p(3,j),'g+');%»æÖÆ¿ØÖƶ¥µã

ifj>n/2+1

plot3([p(1,j-1),p(1,j)],[p(2,j-1),p(2,j)],[p(3,j-1),p(3,j)]);%Á¬½Ó¿ØÖƶ¥µã

end

end

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 工程科技 > 能源化工

copyright@ 2008-2022 冰豆网网站版权所有

经营许可证编号:鄂ICP备2022015515号-1