共轭曲面的数字化方法实现与计算-第一章节.docVIP

两个相互运动且保持相切接触的曲面互为共轭曲面。共轭曲面解问题是在已知一个由解析方程 式表达或一系列离散数字坐标值表达的曲面和已知两个曲面的共轭运动的前提下，求解另一个曲面的问题，在机械设计和机械加工中有着广泛的用途，诸如工具曲 面，轧钢辊面以及齿轮齿面设计、凸轮轮廓设计、机构的运动分析与综合、加工仿真等众多领域。基于MATLAB的数学化共轭曲面求解及仿真研究，旨在充分利 用MATLAB强大硬捷的数学运算和图形显示功能[84-92]，将共轭曲面求解中复杂的推导和计算过程由MATLAB编程实现，并增强计算机对于各种已知曲面（包括解析曲面和数字化曲面）和各种曲面共轭运动的适应能力，以减轻工程人员的负担，并增强计算的精度。 3.2??基于MATLAB的数字化共轭曲面求解实现 3.2.1??求解数字模型 用MATLAB解决数字共轭曲面求解问题，首先要处理已知曲面的数据，把已知曲面变换为适合MATLAB计算的矩阵形式。已知曲面如果由解析表达式表示，则就把解析式离散化为曲面坐标值矩阵；如果由不规则离散坐标值表达，则应把不规则离散坐标值重构为规则坐标值矩阵。 用一维数组与单个点的运动轨迹对应，用二维数组与曲面对应，用三维数组与曲面族对应，使得共轭曲面求解中所用到的各个要素，都由适合MATLAB计算的形式来表达。 共轭曲面求解的关键是求曲面上一点的法矢N和点的运动轨迹的切矢v12，基于MATLAB的数字化共轭曲面求解也是围绕这两个关键量展开。 根据以上论述的曲面共轭条件和MATLAB的运算特点，构造如下的共轭曲面求解数学模型： 式（3-1）中第一式表示曲面Σ1由坐标系S1变换到坐标系S2所形成的曲面族，Xmatrix，Ymatrix，Zmatrix分别表示曲面族中各点在x，y，z坐标方向上的坐标值；第二式中的，NXt，NYt，NZt表示曲面Σ1在点(u,v)对应于t的相对法矢分量，DXt，DYt，DZt表示点(u,v)运动轨迹对应于t的相对切矢分量。 3.2.2 求解模型算法 （1）对规则化后的坐标值进行坐标变换，求出已知曲面Σ1在坐标系S2中对应于序列t的位 置坐标和曲面各点的法矢，分别得到对应于曲面族坐标值的三个三维数组Xmatrix，Ymatrix，Zmatrix和对应于曲面族中曲面法向量的三个三 维数组NXmatrix，NYmatrix，NZmatrix。三维数组的各维分别对应不同的u，v，t； （2）给定uj值； （3）给定vk值； （4）分别将三维数组Xmatrix，Ymatrix，Zmatrix沿点（uj,vk）运动方向进行一维抽取，得到三个一维数组Xt，Yt，Zt，它们表示点（uj,vk）对应序列t的坐标值，也即点（uj,vk）的离散运动轨迹。将Xt，Yt，Zt分别对应于序列t进行三次样条曲线插值，然后分别求导，得到三个一维数组DXt，DYt，DZt，它们表示点（uj,vk）的运动轨迹对应于序列t的相对切矢分量； （5）分别将三维数组NXmatrix，NYmatrix，NZmatrix沿点（uj,vk）运动方向进行一维抽取，得到三个一维数组：NXt，NYt，NZt，它们表示曲面Σ1在点（uj,vk）对应于序列t的相对法矢分量； （6）取ψ（uj,vk,t）=NXt·DXt+NYt·DYt+NZt·DZt，那么ψ（uj,vk,t）是一维数组，它表示点（uj,vk）对应序列t的ψ（u,v,t）的值，然后利用tt=spline(ψ,t,0)实现优化的功能，得到ψ（uj,vk,t）=0时t的值tt； （7）利用计算分段多项式的方法求出点(uj,vk)对应于tt的坐标值Xtt，Ytt，Ztt，然后将它们分别赋给数组元素Xcon（uj,vk），Ycon（uj,vk），Zcon（uj,vk）； （8）vk←vk+△v，转至（3），直到v方向计算完毕； （9）uj←uj+△u，转至（2），直到u方向计算完毕。 如果点（uj,vk）在整个t的运动过程中都没有使ψ（uj,vk,t）=0，则可判断此点没有参与共轭，在共轭曲面的求解中剔除该点；如果有超过一次的使ψ（uj,vk,t）=0，则可判断此点在整个运动过程中多次参与了共轭，这种已知点和所求共轭点“一对多”的映射会导致程序判断的紊乱，对于这种情况，可利用使t自动分段和递归调用共轭求解函数的方法来处理，直至在每一个t的分段中（uj,vk）参与共轭的次数不超过一次。 通过上述计算，最后得到三个二维数组Xcon(u,v)，Ycon(u,v)，Zcon(u,v)，这就是共轭曲面求解的结果，它们表示与曲面Σ1相共轭的数字曲面Σ2的坐标值。在此基础上，调用MATLAB的绘图命令可以显示出共轭曲面的图形。 3.2.3 算法实现 根据上述求解模型与算法，编制出基于MATLAB共轭曲面求解流程图如图3-1所示。