1.二次不等幅非圆齿轮传动行星轮系取苗机构的设计方法，其特征在于：该方法具体步骤如下：步骤一、根据取苗要求，确定取苗轨迹的五个关键点的坐标以及秧针在各关键点的姿态角：关键点A处秧针开始进入钵苗盘，设定秧针在A点的姿态角为γ1，秧针的姿态角定义为秧针与水平面的夹角；关键点B处秧针夹紧钵苗，设定秧针在B点的姿态角为γ2；关键点C为取苗结束点，设定秧针在C点的姿态角为γ3；关键点D为持苗阶段点，设定秧针在D点的姿态角为γ4；关键点E落苗点，设定秧针在E点的姿态角为γ5；

步骤二、建立五位姿二杆开链机构的数学模型，求解五位姿二杆开链机构的两杆杆长；

所述的五位姿二杆开链机构由杆Ⅰ和杆Ⅱ组成，其中FG表示杆Ⅰ，GP表示杆Ⅱ，F为固定铰链点，G为动铰链点，P为杆Ⅱ自由端端点；

2.1 建立杆Ⅱ从关键点A运动到另外的第i个关键点时的坐标及位姿转换关系矩阵，i＝2,3,4,5：其中，关键点A对应杆Ⅱ的位姿为q1，位姿包括坐标和姿态角；中间变量D11i＝cosγ1i,D12i＝-sinγ1i,D13i＝xi-x1cosγ1i+y1sinγ1i,D21i＝sinγ1i,D22i＝cosγ1i,D23i＝yi-x1sinγ1i-y1cosγ1i,D31i＝0,D32i＝0,D33i＝1；

γ1i＝γi-γ1，γ1i表示杆Ⅱ的位姿qi相对位姿q1的角位移，第i个关键点对应杆Ⅱ的位姿为qi，γi表示杆Ⅱ在第i个关键点处的姿态角，γ1表示杆Ⅱ在关键点A处的姿态角，xi表示第i个关键点的x坐标，yi表示第i个关键点的y坐标，x1表示关键点A的x坐标，y1表示关键点A的y坐标；

2.2 建立矩阵坐标方程：

其中矢量G1＝(xg1,yg1)T中xg1和yg1为位姿q1对应的G点坐标，Gi＝(xgi,ygi)T中xgi和ygi为位姿qi对应的G点坐标，矢量G1＝(xg1,yg1)T为未知定量，Gi＝(xgi,ygi)T为未知变量；

2.3 根据FG杆长不变，建立约束方程：

[Gi-F]T[Gi-F]＝[G1-F]T[G1-F]        (3)其中矢量F＝(xf,yf)T为未知定量，xf和yf为固定铰链点F坐标；

2.4 将矩阵坐标方程(2)代入约束方程(3)，整理得：Hixf+Jiyf+Ki＝0(i＝2,3,4,5)         (4)其中

Hi＝Gi1xg1-Gi2yg1+Gi3,Ji＝Gi2xg1+Gi1yg1+Gi4,Ki＝Gi5xg1+Gi6yg1+Gi7,Gi1＝1-D11i,Gi2＝D12i,Gi3＝-D13i,Gi4＝-D23i,Gi5＝D11iD13i+D21iD23i,

2.5 在四个方程Hixf+Jiyf+Ki＝0(i＝2,3,4,5)中任选三个构成方程组，共构成四个方程组，对每个方程组进行消元整理，共得到四个关于xg1和yg1的布尔梅斯特曲线方程：其中，Lkj，k＝1,2,···,8分别为第j个布尔梅斯特曲线方程中从左至右依次排列的各个系数；

2.6 联立四个布尔梅斯特曲线方程(5)，消去xg1化简后得到：其中，a、b、c、d和e均是系数；

一元四次方程(6)的实数解为0个、2个或4个，当出现0个解时，调整秧针在五个关键点的姿态角值，每次调整只针对一个关键点处进行，每次调整时姿态角值增大或减小10％，然后回到步骤2.1，直到一元四次方程(6)存在实数解时，执行步骤2.7；

2.7 将求得的yg1代入布尔梅斯特曲线方程(5)得到xg1，然后将G1的各个实数解与预设的取苗轨迹一同绘制在一副图中，选取该图中x坐标最小的G1作为最终的G1位置，从而得到xg1和yg1的解；

2.8 将xg1和yg1的值代入任意两个方程Hixf+Jiyf+Ki＝0(i＝2,3,4,5)中，求得：xf1代表杆Ⅱ处于位姿q1时固定铰链点F的x坐标，yf1代表杆Ⅱ处于位姿q1时固定铰链点F的y坐标；

步骤三、计算两杆间角位移函数关系；

3.1 取苗轨迹上给定n-5个杆Ⅱ自由端端点P，n≥360，五个关键点也均选作杆Ⅱ自由端端点P，则共有n个杆Ⅱ自由端端点P，保证相邻两个点P的x坐标差值小于0.1；然后，计算杆Ⅰ和杆Ⅱ杆长：LFG+LGP＝max(LFP)

-LFG+LGP＝min(LFP)

其中，LFG、LGP分别是杆Ⅰ和杆Ⅱ的杆长；max(LFP)是杆Ⅰ和杆Ⅱ共线但不重叠时，两杆长度之和；min(LFP)是杆Ⅰ和杆Ⅱ重叠时两杆长度之差；

3.2 建立两自由度角位移关系模型：

①根据每个杆Ⅱ自由端端点P的坐标值计算α和β：

杆Ⅱ自由端端点P由距离固定铰链点F最远位置运动至最近位置过程中，α和β计算式如下：α＝α1+α2

β＝π-β1

其中α1为PF连线与水平方向的夹角；α2为杆Ⅰ与PF连线的夹角；β1为杆Ⅰ与杆Ⅱ的夹角；α表示杆Ⅰ的角位移，即杆Ⅰ与水平方向的夹角；β表示杆Ⅱ相对杆Ⅰ的角位移，即杆Ⅰ与杆Ⅱ夹角的互补角； 点P坐标(xp,yp)，点P到固定铰链点F的长度

杆Ⅱ自由端端点P由距离固定铰链点F最近位置运动至最远位置过程中，α和β计算式如下：α＝α1-α2+2π

β＝π-β1

②求解出n个杆Ⅱ自由端端点P对应的α和β后，便得到杆Ⅰ和杆Ⅱ的两自由度角位移函数关系式：β＝F(α)                           (8)步骤四、总传动比的计算及两级传动比的分配；

4.1 将两自由度机构转化为单自由度的两级齿轮传动行星轮系机构；

计算总传动比i总：

其中ω1是杆Ⅰ相对固定铰链点F的角速度，ω2是杆Ⅱ相对杆Ⅰ的角速度；

4.2 两级传动比的分配：

建立太阳轮和第一中间轮的节曲线方程：

其中太阳轮和第一中间轮的节曲线均采用椭圆方程拟合曲线；r1为太阳轮的曲率半径；

为杆Ⅰ的转角，也就是行星架的转角；A为椭圆的长半轴，k1为椭圆的偏心率，A、k1均为给定值；r2为第一中间轮的曲率半径； 为第一中间轮相对太阳轮的转角； 和 的变化周期都是[0,2π)；i12为太阳轮与第一中间轮的传动比；

4.3 添加初始相位角 使传动比i12曲线上波峰位置对应的杆Ⅰ角位移与总传动比曲线上第二个波峰位置对应的杆Ⅰ角位移调到一致，并计算第二中间轮与行星轮的传动比以及第二中间轮和行星轮的节曲线方程：i34＝i总/i′12                          (13)其中 为添加初始相位角 后行星架的转角， 为添加初始相位角后第一中间轮相对太阳轮的转角，i′12为添加初始相位角 后太阳轮与第一中间轮的传动比，i34为第二中间轮与行星轮的传动比， 为第二中间轮的转角， 为行星轮相对第二中间轮的转角；

给定第二中间轮与行星轮的中心距a34，求得：

第二中间轮的曲率半径

行星轮的曲率半径r4＝a34-r3             (17)步骤五、求得太阳轮、第一中间轮、第二中间轮以及行星轮的齿廓数据；

步骤六、取苗机构的装配：太阳轮、第一中间轮、第二中间轮、行星轮分别与行星架铰接，第一中间轮与第二中间轮固连，行星轮与移栽臂的壳体固连，移栽臂的移栽凸轮与行星架固连；太阳轮与第一中间轮啮合，第二中间轮与行星轮啮合。

2.根据权利要求1所述的二次不等幅非圆齿轮传动行星轮系取苗机构的设计方法，其特征在于：秧针在A点的姿态角γ1设为171度，秧针在B点的姿态角γ2设为143度，秧针在C点的姿态角γ3设为124.5度，秧针在D点的姿态角γ4设为130.50度，秧针在E点的姿态角γ5设为115度。

3.根据权利要求1所述的二次不等幅非圆齿轮传动行星轮系取苗机构的设计方法，其特征在于：步骤四与步骤五之间还有如下步骤：以传动比i34的两个波峰之比最小为优化函数，通过增加或减小k1的数值并通过公式(10)-(13)进行迭代，得到第二级传动比的两个波峰之比最小时的最优解k1，其中，k1增加或减小的步长为0.0001；然后将k1的最优解代入公式(10)-(17)进行计算。

4.根据权利要求1所述的二次不等幅非圆齿轮传动行星轮系取苗机构的设计方法，其特征在于：太阳轮、第一中间轮、第二中间轮以及行星轮的齿廓采用齿形法线法求得。

5.根据权利要求1所述的二次不等幅非圆齿轮传动行星轮系取苗机构的设计方法，其特征在于：若传动比i34曲线具有两个不可导的尖点，则将第二中间轮的齿廓设计为固定在一起的中间不完全齿轮和凹齿轮，将行星轮的齿廓设计为为固定在一起的行星不完全齿轮和凸齿轮；第二中间轮的中间不完全齿轮与行星轮的行星不完全齿轮啮合时，第二中间轮的凹齿轮与行星轮的凸齿轮分离，第二中间轮的凹齿轮与行星轮的凸齿轮啮合时，第二中间轮的中间不完全齿轮与行星轮的行星不完全齿轮分离。