CN103605893B - 基于单叶双曲面各向同性并联机构全局优化设计方法 - Google Patents
基于单叶双曲面各向同性并联机构全局优化设计方法 Download PDFInfo
- Publication number
- CN103605893B CN103605893B CN201310608105.0A CN201310608105A CN103605893B CN 103605893 B CN103605893 B CN 103605893B CN 201310608105 A CN201310608105 A CN 201310608105A CN 103605893 B CN103605893 B CN 103605893B
- Authority
- CN
- China
- Prior art keywords
- parallel institution
- parameter
- design
- hyperboloid
- battle array
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于单叶双曲面各向同性并联机构全局优化设计方法,采用复合单叶双曲面型式的并联机构,消除了负载惯性参数的约束,实现了完全各向同性。在此基础上提出用于衡量全域工作空间模态变化程度的全局模态灵敏度指标,采用该指标作为优化目标,通过优化结构参数实现了并联机构的具体参数设计。采用该发明设计的并联机构,不仅实现了控制中心的完全解耦及动态各向同性,而且保证了全域工作空间内的最优性能。本方法消除了传统六自由度并联机构无法实现完全各向同性的缺陷,实现了并联机构的完全解耦及优越的动态特性,使得其在工业应用中无需为了提升控制性能而研究复杂的控制策略,降低了工业应用成本。
Description
技术领域
本发明涉及结构设计及优化领域,具体是一种基于单叶双曲面各向同性并联机构全局优化设计方法。
背景技术
并联机构以其刚度大、结构稳定、承载能力强的优点在工业领域获得了广泛应用。然而并联机构各自由度间存在的强耦合特性会导致控制性能的严重降低。完全动态各向同性意味着并联机构不仅实现了解耦并且其各阶模态完全相等,由于受到实际负载惯性参数的约束,标准Stewart并联机构无法实现完全动态各向同性。为了获得高性能的控制品质,国内外学者提出了能够满足完全动态各向同性的复合单叶双曲面型式的并联机构,符合该型式的并联机构消除了实际负载惯性参数的约束条件,可实现柔顺控制中心的局部解耦且具有良好的控制性能。然而满足复合单叶双曲面形式的并联机构的解空间并不唯一,考虑全局特性并确定并联机构的具体设计参数在工程设计中具有重要意义。
发明内容
基于以上不足之处,本发明的目的在于提供一种基于单叶双曲面各向同性并联机构全局优化设计方法。本方法基于复合单叶双曲面的满足完全动态各向同性且以全局模态灵敏度为优化目标,根据负载特性设计满足完全各向同性的具体并联机构结构参数。
本发明采用以下技术方案予以实现:
步骤1:确定系统构型
根据负载特性Mt,计算满足完全各向同性的中位雅可比矩阵Jlx0:
式(1)中:
p1(α)=[-ka1sinαka1cosαkc1-a1zka1cosα+r1kc1sinα-α1zka1sinα-r1kc1cosαr1ka1]T
p2(α)=[-ka2sinα-ka2cosαkc2-a2zka2cosα-r2kc2sinαa2zka2sinα-r2kc2cosα-r2ka2]T
复合单叶双曲面的特征参数包括:喉部半径r1及r2,双曲面中心距a1z及a2z,系数c1及c2。
负载特性Mt=[mxmymzIxxIyyIzz]应满足:
mx=my=mz=m,Ixx=Iyy。
各参数计算过程如下:
定义喉部半径比
(1)当n=1时:
(2)当n≠1时:
根据设计要求选取喉部半径比n及角度α,运用公式(1)计算中位雅可比矩阵Jlx0=Jlx(α),其中0≤n≤1,
步骤2:选择优化参数
工程设计中一般固定支腿长度l并通过调节支腿布置方式来获得具体的并联机构结构设计参数。
支腿长度l:lscale为特征尺度。
可优化参数包括:平台中心距d、上平台高度h。
各参数选取原则为:
平台中心距d:位于两组双曲面中的上平台中心与解耦中心距离,当n=1时,可选d=0。
平移高度h:上下平台沿z向平移的高度。
本发明中采用单参数优化,每次优化取上述各参数其中之一。
步骤3:结构参数计算
对并联机构的结构参数设计实质为根据步骤1中获得的中位雅可比矩阵Jlx0求取上下平台铰点空间阵A和B。
上平台铰点空间阵:A=[a1a2a3a4a5a6]
下平台铰点空间阵:B=[b1b2b3b4b5b6]
ai为上平台各铰点空间矢量,bi为下平台各铰点空间矢量i=1,2…6。
具体设计过程如下:
(1)采用步骤1求得的中位雅可比矩阵Jlx0,代入公式(2)中,取出中位各支腿单位矢量Ini0和矢量矩vi0。
(2)计算上下平台铰点空间阵A和B:
当Ini0(3)=0时,k2i=0。
ai=pi+(k1i+k2i)Ini0
B=A-L0
步骤4:建立模态灵敏度函数
模态灵敏度矩阵:
其中:
ε为一微小摄动量,一般情况下可取ε=10-5
雅可比矩阵Jlx(sx):
其中:
Ini为各支腿空间单位矢量i=1,2…6。
c=[xyz]T
c表示cos,s表示sin。
步骤5:计算最小灵敏度
运用公知的黄金分割法寻找最小值:ρmin=minρ(d,h)。
步骤6:生成优化曲线
根据优化变量的不同,绘制优化曲线
确定优化目标阈值,一般选取0≤f≤3db,选取符合目的的优化参数。
步骤7:校核
检验设计的结构参数是否存在干涉,若存在,返回步骤2重新修改设计参数进行优化。
步骤8:结束。
本发明的优点是:
本发明采用新型的基于复合单叶双曲面的并联机构结构形式,给出了一种满足局部完全动态各向同性且在全局工作空间内保证动态特性最优的结构参数优化方法,从设计角度消除了传统六自由度并联机构无法实现完全各向同性的缺陷,实现了并联机构的完全解耦及优越的动态特性,使得其在工业应用中无需为了提升控制性能而研究复杂的控制策略,降低了工业应用成本。
附图说明
图1为典型的复合单叶双曲面型式的并联机构的立体图;
图2为图1的俯视图;
图3为典型的复合单叶双曲面型式的并联机构的数学描述立体图;
图4为图3的主视图;
图5为并联机构构造方法示意图;
图6为优化算法流程图;
图7为实施例1优化曲线图;
图8为实施例1优化结构图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明:
图1-2为典型的复合单叶双曲面型式的并联机构,其各支腿的上下铰点分别位于上下平台的两组圆ra1与ra2、rb1与rb2内,六条支腿按奇偶序号分为两组:a1b1、a3b3、a5b5为一组,a2b2、a4b4、a6b6为一组,两组支腿分别位于双曲面S1和S2中,如图3-4所示,并通过线列Γ1和Γ2分别逆时针依次旋转120°生成。
其中S1的数学描述公式为:
S2的数学描述公式为:
图5为并联机构具体结构参数设计方法,选取平台中心距d及平移高度h,则可确定上平台平面1及上平台平面2,根据支腿长度可确定整体平台结构参数。
图6为优化算法流程图,下面结合具体实施例对其进行说明。
实施例1
步骤1:确定系统构型
根据负载特性Mt,计算满足完全动态各向同性的中位雅可比矩阵Jlx0。
负载特性Mt=[430043004300410041006700]。
各参数的计算过程如下:
选取喉部半径比n=1。
r1=r2=0.8826,ka1=ka2=0.8165,kc1=kc2=0.5774,a1z=a2z=0.7510
选取角度
步骤2:选择优化参数
支腿长度l:取特征尺度lscale=3,
平台中心距d:取d=0。
平移高度h:作为优化变量。
步骤3:结构参数计算
对并联机构的结构参数设计实质为根据步骤1中获得的中位雅可比矩阵Jlx0求取上下平台铰点空间阵A和B。
上平台铰点空间阵:A=[a1a2a3a4a5a6]
下平台铰点空间阵:B=[b1b2b3b4b5b6]
ai为上平台各铰点空间矢量,bi为下平台各铰点空间矢量i=1,2…6。
具体设计过程如下:
(1)采用步骤1求得的中位雅可比矩阵Jlx0,代入公式(2)中,取出中位各支腿单位矢量Ini0和矢量矩vi0。
(2)计算上下平台铰点空间阵A和B:
ai=pi+(k1i+k2i)Ini0
B=A-L0
步骤4:建立模态灵敏度函数
步骤5:计算最小灵敏度
运用公知的黄金分割法寻找最小值:ρmin=minρ(h)=51.5951。
步骤6:生成优化曲线
根据优化变量的不同,绘制优化曲线如附图7所示。
确定优化目标阈值,选取0≤f≤3db,选取优化参数h=0.4m。
步骤7:校核
设计完成的并联机构结构如附图8所示,经检验不存在干涉,优化完成。
步骤8:结束。
本发明的设计方法,采用复合单叶双曲面型式的并联机构,消除了负载惯性参数的约束,实现了完全各向同性。在此基础上提出用于衡量全域工作空间模态变化程度的全局模态灵敏度指标,采用该指标作为优化目标,通过优化结构参数实现了并联机构的具体参数设计。采用该发明设计的并联机构,不仅实现了控制中心的完全解耦及动态各向同性,而且保证了全域工作空间内的最优性能。
Claims (1)
1.基于单叶双曲面各向同性并联机构全局优化方法,其特征在于,方法如下:
步骤1:确定系统构型
根据负载特性Mt,计算满足完全各向同性的中位雅可比矩阵Jlx0:
式(1)中:
p1(α)=[-ka1sinαka1cosαkc1-a1zka1cosα+r1kc1sinα-a1zka1sinα-r1kc1cosαr1ka1]T
p2(α)=[-ka2sinα-ka2cosαkc2-a2zka2cosα-r2kc2sinαa2zka2sinα-r2kc2cosα-r2ka2]T
复合单叶双曲面的特征参数包括:喉部半径r1及r2,双曲面中心距a1z及a2z,系数c1及c2;
负载特性Mt=[mxmymzIxxIyyIzz]应满足:
mx=my=mz=m,Ixx=Iyy;
各参数计算过程如下:
定义喉部半径比
(1)当n=1时:
(2)当n≠1时:
根据设计要求选取喉部半径比n及角度α,运用公式(1)计算中位雅可比矩阵Jlx0=Jlx(α),其中0≤n≤1,
步骤2:选择优化参数
工程设计中固定支腿长度l,并通过调节支腿布置方式来获得具体的并联机构结构设计参数;
支腿长度l:lscale为特征尺度;
可优化参数包括:平台中心距d、上平台高度h;
各参数选取原则为:
平台中心距d:位于两组双曲面中的上平台中心与解耦中心距离,当n=1时,d=0;
平移高度h:上下平台沿z向平移的高度;
本发明中采用单参数优化,每次优化取上述各参数其中之一;
步骤3:结构参数计算
对并联机构的结构参数设计实质为根据步骤1中获得的中位雅可比矩阵Jlx0求取上下平台铰点空间阵A和B;
上平台铰点空间阵:A=[a1a2a3a4a5a6]
下平台铰点空间阵:B=[b1b2b3b4b5b6]
ai为上平台各铰点空间矢量,bi为下平台各铰点空间矢量i=1,2…6;
具体设计过程如下:
(1)采用步骤1求得的中位雅可比矩阵Jlx0,代入公式(2)中,取出中位各支腿单位矢量lni0和矢量矩Vi0;
(2)计算上下平台铰点空间阵A和B:
当lni0(3)=0时,k1i=0;
当lni0(3)=0时,k2i=0;
ai=pi+(k1i+k2i)lni0
L0=l·[1n10ln20ln30ln40ln50ln60]
B=A-L0
步骤4:建立模态灵敏度函数
模态灵敏度矩阵:
其中:
ε为一微小摄动量,取ε=10-5;
GT(sx)为广义频率阵,U(sx)i为模态解耦阵U(sx)的第i列,U(sx)为广义频率阵GT(sx)的模态解耦阵,
雅可比矩阵Jlx(sx):
其中:
lni为各支腿空间单位矢量i=1,2…6;
c=[xyz]T
步骤5:计算最小灵敏度
运用公知的黄金分割法寻找最小值:ρmin=minρ(d,h);
步骤6:生成优化曲线
根据优化变量的不同,绘制优化曲线
确定优化目标阈值,0≤f≤3db,选取符合目的的优化参数;
步骤7:校核
’检验设计的结构参数是否存在干涉,若存在,返回步骤2重新修改设计参数进行优化;
步骤8:结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310608105.0A CN103605893B (zh) | 2013-11-19 | 2013-11-19 | 基于单叶双曲面各向同性并联机构全局优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310608105.0A CN103605893B (zh) | 2013-11-19 | 2013-11-19 | 基于单叶双曲面各向同性并联机构全局优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103605893A CN103605893A (zh) | 2014-02-26 |
CN103605893B true CN103605893B (zh) | 2016-03-23 |
Family
ID=50124114
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310608105.0A Expired - Fee Related CN103605893B (zh) | 2013-11-19 | 2013-11-19 | 基于单叶双曲面各向同性并联机构全局优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103605893B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107127738A (zh) * | 2017-06-19 | 2017-09-05 | 哈尔滨工业大学 | 刚柔混合线解耦变刚度的可重构冗余驱动广义并联机构 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109555812B (zh) * | 2018-07-18 | 2023-04-28 | 中国航空工业集团公司洛阳电光设备研究所 | 一种压电驱动的各向同性多自由度隔振平台 |
CN110398332A (zh) * | 2019-07-25 | 2019-11-01 | 北京航宇振控科技有限责任公司 | 一种异面构型的六自由度振动模拟装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103009375A (zh) * | 2012-11-30 | 2013-04-03 | 哈尔滨工业大学 | 动态各向同性广义Stewart容错并联机构及其设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7880478B2 (en) * | 2008-08-04 | 2011-02-01 | Chung-Yuan Christian University | Sensing device for measuring a position of nanoscale motion apparatus |
-
2013
- 2013-11-19 CN CN201310608105.0A patent/CN103605893B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103009375A (zh) * | 2012-11-30 | 2013-04-03 | 哈尔滨工业大学 | 动态各向同性广义Stewart容错并联机构及其设计方法 |
Non-Patent Citations (4)
Title |
---|
Dynamic isotropic design of a class of Gough-Stewart parallel manipulators lying on a circular hyperboloid of one sheet;Jiang Hong-zhou等;《Mechanism and Machine Theory》;20101119;第46卷(第3期);第358-374页 * |
Influence of Passive Joint Damping on Modal Space Decoupling for a Class of Symmetric Spatial Parallel Mechanisms;TX Tian等;《Applied Mechanics and Material》;20130213;第1152-1157页 * |
基于加速度各向同性的空间光学仪器主动隔振Stewart平台设计;佟志忠 等;《振动与冲击》;20120128;第31卷(第2期);第109-114页 * |
基于复合单叶双曲面的广义Gough-Stewart结构六维力传感器各向同性优化设计;佟志忠 等;《航空学报》;20130319;第2607-2615 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107127738A (zh) * | 2017-06-19 | 2017-09-05 | 哈尔滨工业大学 | 刚柔混合线解耦变刚度的可重构冗余驱动广义并联机构 |
Also Published As
Publication number | Publication date |
---|---|
CN103605893A (zh) | 2014-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103605893B (zh) | 基于单叶双曲面各向同性并联机构全局优化设计方法 | |
CN102521435B (zh) | 一种对称层合板的铺层制作方法 | |
CN106570823B (zh) | 基于平面特征匹配的点云粗拼接方法 | |
CN104866666A (zh) | 一种基于有限元模型修正的索网反射面形面精度的鲁棒调整方法 | |
CN102493569B (zh) | 一种建筑结构基于抗震性能的优化方法和系统 | |
Vescovini et al. | Efficient post-buckling analysis of variable-stiffness plates using a perturbation approach | |
CN102607330A (zh) | 惯导系统传递对准中基准信息的优化处理方法 | |
CN109657289A (zh) | 基于灵敏度分析的参数分类-分步模型修正方法 | |
CN103615054A (zh) | 一种基于区格剪切变形的屈曲约束支撑布置方法 | |
Banichuk et al. | Structural optimization with uncertainties | |
Broujerdian et al. | Effect of curvature and aspect ratio on shear resistance of unstiffened plates | |
CN107656227B (zh) | 基于Levenberg-Marquardt算法的磁力计校准方法 | |
Zhou et al. | An approximation method for computing the dynamic responses and equivalent static wind loads of large-span roof structures | |
CN105806572B (zh) | 一种单层网壳结构节点半刚性测定方法 | |
CN102729249B (zh) | 适用于模态空间控制的六自由度并联机构参数优化方法 | |
Yan et al. | Calculation theory and experimental study of the K6 single-layer reticulated shell | |
CN102354391B (zh) | 基于全向搜索机制变异函数的最优遥感尺度选择方法 | |
CN116306178B (zh) | 基于自适应形函数和等效中性层的结构应变反演方法 | |
Song et al. | Nonlinear investigation of Gol’denveizer’s problem of a circular and elliptic elastic torus | |
CN108229055B (zh) | 一种质刚重合的三向等刚度气浮式隔振系统优化设计方法 | |
CN103761367B (zh) | 一种基于位移控制目标的弦支梁结构内力的解析计算方法 | |
CN102982206B (zh) | 评价汽轮发电机隔振基础静位移和变形的设备和方法 | |
Zhou et al. | Wind loads and wind-induced responses of Guangzhou New TV Tower | |
Koo et al. | A study on the stability of the single-layer latticed dome during erection using the step-up method | |
Blaise et al. | Principal static wind loads on a large roof structure |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160323 Termination date: 20161119 |
|
CF01 | Termination of patent right due to non-payment of annual fee |