机构设计凸轮曲线程序.doc

ClsDim t(7) as doubleDim p(7) as doubleDim f(7) as doubleDim c(7) as doubleDim b(7) as doubleDim jsd(7) as doubleDim yd(7) as doubleDim sd(7) as doubleDim wy(7) as doubleDim pi as doubleDim a1 as doubleDim a2 as doubleDim m as doubleDim i as doubleDim j as doubleDim mc as doubleDim ss as doubleJ=0Dim tt as doubleDim aa as doubleStar:Pi=3.1415926ss=”请选择通用间歇体形曲线的名称:”& chr(10)+chr(13)&”1、修正等速“& chr(10)+chr(13)&” 2、修正梯形”& chr(10)+chr(13) &” 3、修正正弦”mc=inputbox(ss)if mc=1 thent(0)=0t(1)=1/16t(2)=1/16t(3)=1/4t(4)=3/4t(5)=15/16t(6)=15/16t(7)=1if mc=2 thent(0)=0t(1)=1/8t(2)=3/8t(3)=1/2t(4)=1/2t(5)=5/8t(6)=7/8t(7)=1elseif mc=3 thent(0)=0t(1)=1/8t(2)=1/8t(3)=1/2t(4)=1/2t(5)=7/8t(6)=7/8t(7)=1elsemsgbox(“请输入1-3的数”)goto startend iffor i= 1 to7f(i)=2*(t(i)-f(i-1)/pinext im=(f(1)+f(2)*pi/2+f(3)/(f(5)+f(6)*pi/2+f(7)a1=(f(3)2+0.5*(t(2)-t(1)2)-f(1)2+f(3)*(1-t(3)+(t(2)-t(1)*(1-t(2)+f(1)-m*(f(7)2+0.5*(t(6)-t(5)2)- f(5)2+f(5)*(1-t(4)+(t(6)-t(5)*(1-t(6)(-1)a2=m*a1c(1)=f(1)*a1c(2)=-t(1)*a1+c(1)c(3)= t(2)*a1+c(2)c(4)= f(3)*a1+c(3)c(5)=-f(5)*a2+c(4)c(6)=t(5)*a2+c(5)c(7)=- t(6)*a2+c(6)b(1)=0b(2)=-a1*(f(1)2+0.5*t(1)2)+t(1)*(c(1)-c(2)+b(1)b(3)= a1*(f(3)2+0.5*t(2)2)+t(2)* (c(2)-c(3)+b(2)b(4)=t(3)*(c(3)-c(4)+b(3)b(5)=t(4)*(c(4)-c(5)+b(4)b(6)=a(2)*(f(5)2+0.5*t(5)2)+t(5)*(c(5)-c(6)+b(5)b(7)=-a(2)*(f(7)2+0.5*t(6)2)+t(6)*(c(6)-c(7)+b(6)for tt=0 to 1 step 0.01for i=1 to 7select case icase 1 and (tt=t(0) and ttt(1)and ttt(2)and ttt(3)and ttt(4)and ttt(5)and ttt(6)and tt=t(7)p(i)=(tt-t(i-1)/f(i)+(i-1)*pi/4j=j+1jsd(j)=a2*sin(p(i)yd(j)=a2*cos(p(i)/f(i)sd(j)=-a2*f(i)*cos(p(i)+c(i)wy(j)= -a2*(f(i)2)*sin(p(i)+c(i)*tt+b(i)open “d:加速度.txt”for append as #1print #1,jsd(j)closeopen “d:跃度.txt”for append as #1print #1,yd(j)closeopen “d:速度.txt”for append as #1print #1,sd(j)closeopen “d:位移.txt”for append as #1print #1,wy(j)closeend selectnext inext iimsgbox (“计算完毕”&chr(10)+chr(13)&”位移数据在 d:位移.txt” &chr(10)+chr(13)& ”速度数据在 d:速度.txt” &chr(10)+chr(13)&”加速度数据在 d:加速度.txt” &chr(10)+chr(13)& ” 跃度数据在 d:位移.txt”)msbox(“输出图形“)for i=1 to 1200 step 2form1.currentx=i*5+100form1.currenty=1500print”.”Nextfor i=1 to 101currentx=i*20+1500currenty=-wy(i)*600+1500print”.”Nextcurrentx=1500currenty=2500print”位移曲线”for i=1500 to 3000 step 2form1.currentx=i*5+100form1.currenty=1500print”.”Nextfor i=1 to 101currentx=i*20+8500currenty=-sd(i)*200+1500print”.”Nextcurrentx=9500currenty=2500print”速度曲线”for i=1 to 1200 step 2form1.currentx=i*5+100form1.currenty=5500print”.”Nextfor i=1 to 1000 step 2form1.currentx=i*2+5000form1.currenty=1500print”.”Nextfor i=1 to 101currentx=i*20+1500currenty=-jsd(i)*50+5500print”.”Nextcurrentx=2000currenty=7000print”加速度曲线”for i=1500 to 3000 step 2form1.currentx=i*5+100form1.currenty=5500print”.”Nextfor i=1 to 1000 step 2form1.currentx=i*2+5000form1.currenty=8500print”.”Nextfor i=1 to 101 step 0.8currentx=i*20+8500currenty=-yd(i)*6+5500print”.”Nextcurrentx=9000currenty=7000print”跃度曲线”printprintselect case mccase 1currentx=6000currenty=12000print “修正等速”;”的四个曲线”case 2curren