并联式太阳跟踪机构设计与尺寸优化.pdf

制造使用改进 ：收稿日期：；修回日期：基金项目：国家自然科学基金项目（）；浙江省 人才培养计划项目和浙江省高校中青年学科带头人培养项目；衢州市科技计划项目（）；衢江区科技计划项目（）。第一作者简介：赵世豪（），男，河南周口人，硕士研究生，主要研究方向为并联机构。通信作者：杜小强（），男，福建福清人，博士，教授，博士 硕士研究生导师，主要从事机构设计及理论研究。：并联式太阳跟踪机构设计与尺寸优化赵世豪，杜小强（浙江理工大学 机械工程学院，浙江 杭州 ）摘要：为实现聚光光伏系统的高效跟踪作业，课题组提出了一种可用于太阳跟踪，且具有刚度高、解耦性好及结构简单特点的 并联机构。首先，运用螺旋理论对并联机构的自由度进行了验证；其次，根据并联机构的几何约束条件构建了运动学逆解模型，并通过对速度雅克比矩阵进行分析可知该机构为解耦机构；然后，基于太阳方向角与机构欧拉角的映射关系，采用数值离散搜索法得到了并联机构的工作空间，并提出了一种工作空间评价指标；最后，基于聚光光伏系统对并联跟踪机构的性能要求，以机构的灵巧度和工作空间指标为尺寸优化的目标，采用遗传算法对主动杆、从动杆和动平台的尺寸进行多目标优化。结果表明优化后并联机构的灵巧度和工作空间都有所提升。虽然优化后并联机构的工作空间在高度角为 的区域并未完全覆盖太阳轨迹区域，但未覆盖的面积集中在一天中辐照度较低的区域。关键词：聚光光伏；太阳跟踪机构；工作空间；数值离散搜索法；灵巧度中图分类号：文献标志码：文章编号：（），（，）：，：；第 卷 第 期 年 月轻工机械 为提高光伏（，）电池的发电效率，通常采用太阳能跟踪机构对太阳进行跟踪 。常见的跟踪机构按驱动方式可以分为单轴式和双轴式（其本质是串联机构）。相较于双轴式机构，并联式机构具有刚度大、承载能力强、结构稳定性好和精度高等优点 ，可 用 于 跟 踪 精 度 要 求 较 高 的 聚 光 光 伏（，）系统，但也有工作空间较小的缺点。近些年来，国内外为并联机构应用于太阳能跟踪做了很多研究，但大多数机构都存在结构复杂、动力源较多、运动不解耦和伴随有额外运动的问题。等 提出了一种 的并联机构，并将其应用于 跟踪，该机构的跟踪误差可以限制在 以内。等 提出了一种冗余的 并联式太阳能跟踪机构，该机构具有工作空间大、能耗低的特点。等 提出了一种 并联太阳能跟踪机构，被用来进行塔式发电。技术是通过透镜或反射镜等光学元件将太阳光聚集到一个极小的 电池上使光能转化为电能 。相较于传统的 电池发电，技术具有更高的单位面积输出功率和更低的发电成本 。但由于其结构特点，该发电方式对机构的跟踪精度要求较高 。因此，为提高 系统的发电效率，首先要对用于 发电的并联机构的综合性能进行优化。尺寸优化是进行并联机构结构设计的重要环节，涉及到运动学分析、工作空间分析和动力学分析等 ，但具体优化指标的选择要依据该机构的实际工作要求。课题组针对现有并联式太阳跟踪机构存在的问题，提出了一种 旋转并联机构，该并联机构具有刚度高、解耦性好、结构简单和工作空间大的特点；并针对 系统对并联机构跟踪精度、工作空间的性能要求，采用灵巧度和工作空间指标作为优化目标，进行尺寸优化研究。太阳轨迹区域分析 太阳的高度角与方位角为确定优化结构是否满足对太阳的跟踪区域要求，在进行优化分析前要对太阳的轨迹进行分析。太阳在天空中的具体位置可以通过计算太阳的高度角和方位角求得，太阳的高度角和方位角的计算公式如下：。（）；。（）式中：为太阳时角，为高度角，为方位角，为赤纬角度，为地区纬度，为真太阳时。杜春旭等 对各种赤纬角的求解方法进行了对比分析，得出 提出的算法精度最高，的赤纬角算法如下：（）（）（）（）；（）其中，（）。式中：为日子数，为由世界时零时开始算起。计算与结果修正基于真太阳时与太阳时角的关系可求得太阳时角的变化范围。以杭州地区为例，其经度为 ，纬度为 ，结合式（）、（）和（）即可得全年中杭州地区太阳的轨迹区域。由式（）的第 个公式中太阳方位角与真太阳时的关系可知，在时间段 ：和 ：内太阳方位角分别存在一个拐点，该拐点可通过对该公式进行时间求导得到。存在拐点说明方位角存在突变，与真实的太阳运动不符，因此结合式（）的第 个公式对其进行修正，太阳方位角的表达式为：（），（）；（），（）；（），（）。（）其中，（）（）。（）根据太阳高度角、方位角与真太阳时的关系，可求得太阳高度角与方位角的关系。由太阳的运动规律可知：只有当高度角的范围为 ，太阳方位角的范 制造使用改进赵世豪，等：并联式太阳跟踪机构设计与尺寸优化围为 时，对太阳的跟踪才是有效的 。此跟踪范围对应一天中的白天时间段，此时太阳的运动轨迹区域如图 所示。图 杭州地区太阳轨迹区域 并联机构简介与自由度分析 并联机构的结构简图和运动螺旋分布如图 所示。该机构由定平台（）、动平台（）和 条支链组成；条支链分别为 支链（）、支链（）和 支链（）；动平台和定平台都为直角三角形，条支链通过运动副分别与动平台、定平台的 个角相连，固定支链一端固定于定平台的直角点，另一端通过万向副与动平台相连。在固定支链中，万向副（）的静旋转轴线与定平台的直角边 平行，动旋转轴线与定平台的直角边位于同一平面；在 支链中，万向副（）的静旋转轴线与静平台的直角边 平行，且 个转动副（）的旋转轴线与万向副（）的静旋转轴线相互平行；在 支链中，万向副（）的静旋转轴线（固定在连杆 上）与定平台的直角边 平行，且转动副（）的旋转轴线与万向副（）的静旋转轴线相互平行。在定平台直角点 和动平台直角点 分别建立系统坐标系 和动坐标系 。在系统坐标系 下，固定支链 的运动螺旋系可表示为：（；）；（；）。（）图 并联机构 式中：为运动副中点 在系统坐标系下 轴方向的坐标值，为第 条支链第 个运动副的轴线在坐标轴 上的方向余弦，为第 条支链第 个运动副的轴线在坐标轴 上的方向余弦。对式（）求互易积，可得反螺旋系为：（；）；（；）；（；）；（；）。（）式中：为过 点并沿 轴方向的约束力，为过点并沿 轴方向的约束力，为过系统坐标系原点 并沿 轴方向的约束力，为垂直于 和 轴的约束力偶。轻工机械 年第期因此，固定支链绕 方向的转动和沿 ，和 方向的移动被限制。支链的运动螺旋系可以表示为：（；）；（；）；（；）；（；）。（）同样的，可得其反螺旋系为：（；）；（；）。（）式中：为垂直于 和 轴的约束力偶；为过点并与 轴平行的约束力。因此，支链绕 轴的转动和沿 的移动被限制。支链的运动螺旋系可以表示为：（；）；（；）；（；）；（；）；（；）；（；）。（）式中，为第 条支链第 个运动副的轴线在坐标轴 上的方向余弦。通过对式（）求互易积发现，支链的运动螺旋系组成的矩阵为满秩矩阵 （），所以，支链在空间中具有 个独立自由度。结合式（）和（）可建立 并联机构的约束螺旋系：（；）；（；）；（；）；（；）；（；）；（；）。（）由固定支链与 支链的几何关系和相对运动可以看出：在系统坐标系 下，多边形 始终在 平面内，因此，动平台的关节点 也始终在 平面内运动，多边形 各关节点的位置关系如图 所示。图 固定支链和 支链的几何位置关系 基于固定支链和 支链的几何关系可得：（）。（）式中，为动坐标系 轴与系统坐标系 轴的夹角。同时，系统在运动过程中，螺旋 和螺旋 始终共轴线，即有 螺旋线沿 轴和 轴的余弦分量为：；。（）根据固定支链和 支链的几何位置关系，并考虑运动副的相对位姿，结合式（）和式（），对式（）进行分析，可以看到约束螺旋 和 线性相关，约束螺旋 和 线性相关。求式（）的反螺旋，可得 并联机构的运动螺旋系为：（；）；（；）。（）从式（）中可知：运动螺旋 为 个绕平行于 轴并过 点的旋转运动，即并联机构具有 个绕系统坐标系 轴转动自由度；运动螺旋 为 个绕平行于 轴并过 点的旋转运动，即并联机构具有 个绕动坐标系 轴的转动自由度。由 并联机构的约束螺旋系可知 个约束力偶具有相同的方向，因此为并联机构添加了 个公共约束，即 ，可得并联机构的阶数为 。个约束力偶相当于 个约束力偶，并限制了并联机构在 方向上的转动。基于自由度修正 制造使用改进赵世豪，等：并联式太阳跟踪机构设计与尺寸优化公式：（）。（）式中：为并联机构的阶数，为并联机构的冗余约束个数，表示包含机架的并联机构构件数，为并联机构的运动副数目，为第 个运动副的自由度数目。对于非对称并联机构，冗余约束可由以下公式计算求得：。（）式中：为第 个支链的独立反螺旋数，为支链个数，为除公共约束外剩余的独立约束数。基于上述螺旋分析，并结合式（）和（），可得 并联机构的自由度为：（）。（）运动学分析 逆运动学分析此并联机构的逆运动学问题本质上是在已知动平台旋转角度 和 时，通过位姿逆解公式求得驱动杆转角 和 的位移参数。通过对 并联机构的分析可知该机构在空间中拥有 个旋转自由度，即绕 轴和 轴的旋转。因此，动平台相对于静平台的欧拉角只有 个，即绕 轴旋转的欧拉角为 及绕 轴旋转的欧拉角为 ，于是可得旋转变换矩阵为：（，）（，）。（）在系统坐标系 下，由并联机构的几何关系可知：，。（）将旋转矩阵和各关节的位置矢量代入式（），可得关于主动关节转角的三角超越方程：，。（）可得 并联机构的运动学逆解公式：槡，（）式中：，和 为已知的并联机构结构参数。速度雅克比矩阵并联机构的速度雅克比矩阵反映了主动关节转角与输出位姿之间的映射关系，通过对运动学逆解公式（）进行全微分可得并联机构的速度映射关系：。（）式中：为驱动关节的雅克比矩阵；为反映了当驱动关节锁死时，被动关节的雅克比矩阵；为驱动关节位移导数（速度）的矢量；为动平台的角度导数（角速度）的矢量。对式（）进行变换，可得如下公式：。（）由于 ，机构的速度雅克比矩阵可表示为：。（）式中：和 为驱动关节的雅克比矩阵的主对角线元素，为被动关节的雅克比矩阵的元素。由式（）可知该机构的速度雅克比矩阵为上三角矩阵，因此，该机构为解耦机构。灵巧度指标为提高并联机构的跟踪精度和运动 力传递特性，在此引入灵巧度评价指标，灵巧度被定义为机构速度雅克比矩阵条件数的倒数 。矩阵的条件数越小，说明机构的灵巧度越大，机构的跟踪精度和运动 力传递特性越好。计算公式可表达为：（）。（）其中，（）。式中：（）为矩阵的条件数，为 雅克比矩阵在某一机构位姿下的数值。被称为局部灵巧度指标。为反映并联机构在整个工作空间内的运动 力传递性能，定义全局灵巧度指标，计算公式可表达为：。（）轻工机械 年第期式中 为并联机构的工作空间。工作空间分析 并联机构欧拉角与太阳高度角、方位角的关系为研究并联机构的工作空间对于太阳轨迹区域的覆盖情况，在进行工作空间求解前，要对并联机构欧拉角与太阳的方位角、高度角之间的转化关系进行分析，将机构的欧拉角下的工作空间转化到太阳轨迹区域衡量下的工作空间。并联机构的方位角 可被定义为系统坐标系 轴负方向矢量 与机构动平台法向量在系统坐标系 平面上投影矢量 之间的夹角；机构的高度角 可被定义为系统坐标系 轴正方向矢量 与动平台法向量 的夹角。可得如下表达式：；。（）把参数代入式（）得：（）槡 槡 （）（）槡 （）如图 所示，由动平台与太阳的几何关系，可知并联机构与太阳的高度角及方位角有如下转化公式：；。（）图 太阳位置角与动平台相应转角之间关系 工作空间求解机构工作空间的本质为机构动平台所有可达位姿点的集合，用以衡量在特定结构参数和约束条件下机构活动空间的大小。课题组采用数值离散搜索法 进行求解。求解流程如图 所示，此机构的约束条件有 个，即几何约束、球副转角限制、万向副转角限制和机构奇异性。在首次循环中给定欧拉角 ，的最小值，根据几何约束条件，判断运动学逆解是否非实数解，并判断在此欧拉角下是否满足关节转角约束和奇异条件，若存在 个约束条件不满足，则对驱动杆 进行迭代更新，若全部满足，输出对应的欧拉角 和 即为工作空间边界位姿，并将其转化为 和。对欧拉角 进行迭代，直到 。重复上述操作，对欧拉角 进行迭代，即可得机构的工作空间。图 并联机构工作空间求解流程 工