国产精品乱码一区-性开放网站-少妇又紧又爽视频-西西大胆午夜人体视频-国产极品一区-欧美成人tv-四虎av在线-国产无遮挡无码视频免费软件-中文字幕亚洲乱码熟女一区二区-日产精品一区二区三区在线观看-亚洲国产亚综合在线区-五月婷婷综合色-亚洲日本视频在线观看-97精品人人妻人人-久久久久久一区二区三区四区别墅-www.免费av-波多野结衣绝顶大高潮-日本在线a一区视频高清视频-强美女免费网站在线视频-亚洲永久免费

機械社區(qū)

標題: 參數(shù)化計算原始不對稱型線 [打印本頁]

作者: shouce 時間: 2015-12-22 16:32
標題: 參數(shù)化計算原始不對稱型線
%CalpaMEF.m
%原始不對稱型線計算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齒數(shù)比
R1=(Z1/(Z1+Z2))*A; %陽轉(zhuǎn)子節(jié)圓半徑
R2=(Z2/(Z1+Z2))*A; %陰轉(zhuǎn)子節(jié)圓半徑
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數(shù)范圍在等腰三角形中求
%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
%plot(x1,y1)

%第二曲線方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個地方第一次弄錯了
%t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%陽轉(zhuǎn)子方程
%t1=linspace(0,t1,100);
%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉(zhuǎn)角參數(shù)
%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%轉(zhuǎn)角參數(shù)
%q=linspace(q1,q2,100);
k=i+1;
%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
%plot(x22,y22)

%t21=acos((2*R2^2-R^2)/(2*R2^2));
%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C點橫坐標
%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C點縱坐標
%cp=((x0-R2)^2+y0^2)^(1/2);%計算線段長度
%t22=acos((2*R2^2-cp^2)/(2*R2^2));
%P001=(A^2+R2^2-2*A*R2)^(1/2);
%P002=b1;
%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號輸錯
%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
%qm=linspace(qm01,qm02,100);
%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
%plot(x11,y11)

%第二曲線方程 EF  EF  EF
t21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21); %有點問題%  為什么是這個樣的？
p004=R2;
%PP=linspace(p003,p004,100);
qm03=1/i*(acos(k*p003/A)-t21);
qm04=1/i*(acos(k*p004/A)-t21);
qm1=linspace(qm03,qm04,100);
x12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
z12=0*qm1;
plot(x12,y12)

EF=[x12',y12',z12']
%save('EF.txt')
end

%CalpaMFG.m
%原始不對稱型線計算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
function [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齒數(shù)比
R1=(Z1/(Z1+Z2))*A; %陽轉(zhuǎn)子節(jié)圓半徑
R2=(Z2/(Z1+Z2))*A; %陰轉(zhuǎn)子節(jié)圓半徑
t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數(shù)范圍在等腰三角形中求
t=linspace(0,t,200);
x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
z1=0*t;
plot(x1,y1)
%
FG=[x1',y1',z1']
%save('FG.txt')
end

% [ x1, y1] = CalpMFG(110, 5, 6,30)

%CalpaMGH.m
%原始不對稱型線計算程序
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
i=Z1/Z2; %齒數(shù)比
R1=(Z1/(Z1+Z2))*A; %陽轉(zhuǎn)子節(jié)圓半徑
R2=(Z2/(Z1+Z2))*A; %陰轉(zhuǎn)子節(jié)圓半徑
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數(shù)范圍在等腰三角形中求
%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
%plot(x1,y1)

%第二曲線方程 GH  GH GH
b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個地方第一次弄錯了
t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%陽轉(zhuǎn)子方程
%t1=linspace(0,t1,100);
q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉(zhuǎn)角參數(shù)
%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%轉(zhuǎn)角參數(shù)
%q=linspace(q1,q2,100);
k=i+1;
%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
%plot(x22,y22)

%第三段曲線
%t21=acos((2*R2^2-R^2)/(2*R2^2));
x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C點橫坐標
y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C點縱坐標
cp=((x0-R2)^2+y0^2)^(1/2);%計算線段長度
t22=acos((2*R2^2-cp^2)/(2*R2^2));
P001=(A^2+R2^2-2*A*R2)^(1/2);
P002=b1;
qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號輸錯
qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
qm=linspace(qm01,qm02,100);
x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
z11=0*qm;
plot(x11,y11)
%
GH=[x11',y11',z11']
%save('GH.txt')
end

歡迎光臨機械社區(qū) (http://www.whclglass.com.cn/)