%子函數(shù)：大輪輪齒凸面方程和法線方程                                                %程序功能：由大輪加工刀具的內(nèi)刀忍切削錐面產(chǎn)成大輪的凸面                                                %編程時間：2008年9月1日                                                %編 程 者：樊紅衛(wèi) 谷霽紅%u_g,theta_g:大輪切削錐面的曲面坐標%R_c2：刀尖半徑%alphi_c2：刀具齒形角%phi_cr2:搖臺轉角%r_c2_jxdw:徑向刀位%q_2：角向刀位%smallgama_2:大輪的機床安裝角（根錐角）%X_B2:大輪垂直輪位%phi_2:被加工大輪轉角phi_2=phi_cr2/m_cg%m_cg：加工大輪的切削滾比clear alltic;%………………開始計時語句，與最后的toc配合syms u_g  theta_g  phi_cr2R_c2=0.5*D_c2-0.5*W_2;%內(nèi)刀刀尖半徑r_c2=[(R_c2-u_g*sin(alphi_c2))*cos(theta_g);(R_c2-u_g*sin(alphi_c2))*sin(theta_g);-u_g*cos(alphi_c2)];%坐標系S_c2中大輪刀具內(nèi)刀刃切削錐面方程N_c2=cross(diff(r_c2,'theta_g'),diff(r_c2,'u_g'));%切削錐面的法矢量n_c2=N_c2/abs(N_c2);%大輪切削錐面的單位法矢量phi_2=phi_cr2/m_cg;%大輪切削滾比恒定時M_cr2c2=[1 0 0 r_c2_jxdw*cos(q_2);0 1 0 r_c2_jxdw*sin(q_2);0 0 1 0;0 0 0 1];%對徑矢從S_c2到S_cr2M_m2cr2=[cos(phi_cr2) sin(phi_cr2) 0 0;-sin(phi_cr2) cos(phi_cr2) 0 0;0 0 1 0;0 0 0 1];%對徑矢從S_cr2到S_m2M_am2=[cos(smallgama_2) 0 sin(smallgama_2) -X_B2*sin(smallgama_2);0 1 0 0;-sin(smallgama_2) 0 cos(smallgama_2) -X_B2*cos(smallgama_2);0 0 0 1];%對徑矢從S_m2到S_aM_2a=[1 0 0 0;0 cos(phi_2) -sin(phi_2) 0;0 sin(phi_2) cos(phi_2) 0;0 0 0 1];%對徑矢從S_a到S_2r_c2_yipie=[(R_c2-u_g*sin(alphi_c2))*cos(theta_g);(R_c2-u_g*sin(alphi_c2))*sin(theta_g);-u_g*cos(alphi_c2);1];%為了利用矩陣乘法做變換對r_c2的變形r_2(u_g,theta_g,phi_cr2)=M_2a*M_am2*M_m2cr2*M_cr2c2*r_c2_yipie;%坐標系S_2中大輪凸面方程M_cr2c2_yipie=[1 0 0;0 1 0;0 0 1];%對自由矢量從S_c2到S_cr2M_m2cr2_yipie=[cos(phi_cr2) sin(phi_cr2) 0;-sin(phi_cr2) cos(phi_cr2) 0;0 0 1];%對自由矢量從S_cr2到S_m2n_m2=M_m2cr2_yipie*M_cr2c2_yipie*n_c2;%固定坐標系S_m2中大輪切削錐面的單位法矢量w_2_m2=[-cos(smallgama_2);0;-sin(smallgama_2)];%固定坐標系S_m2中大輪的旋轉角速度w_cr2_m2=[0;0;-m_cg];%固定坐標系S_m2中搖臺的旋轉角速度r_m2=M_m2cr2*M_cr2c2*r_c2_yipie;%固定坐標系S_m2中大輪切削錐面上某一點的位置R_m2=[0;0;X_B2];%S_m2和S_2坐標原點之間的距離v_cg_m2=cross((w_cr2_m2-w_2_m2),r_m2)-cross(R_m2,w_2_m2);%固定坐標系S_m2中被加工大輪和產(chǎn)形輪之間的相對角速度