CN111444644A - 一种基于有限元技术的工业机器人整机刚度优化方法 - Google Patents
一种基于有限元技术的工业机器人整机刚度优化方法 Download PDFInfo
- Publication number
- CN111444644A CN111444644A CN202010141357.7A CN202010141357A CN111444644A CN 111444644 A CN111444644 A CN 111444644A CN 202010141357 A CN202010141357 A CN 202010141357A CN 111444644 A CN111444644 A CN 111444644A
- Authority
- CN
- China
- Prior art keywords
- rigidity
- finite element
- substructure
- whole
- design
- 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.)
- Granted
Links
Images
Landscapes
- Numerical Control (AREA)
- Manipulator (AREA)
Abstract
本发明涉及机器人刚度优化设计领域,具体是一种基于有限元技术的工业机器人整机刚度优化方法,其具体步骤如下:S1:初始模型有限元建模;S2:子结构划分;S3:调整子结构刚度参数;S4:计算整机刚度变化量;S5:判断设计是否合理;S6:结构优化;利用有限元手段建立整机模型对整机刚度和零部件刚度的设计合理性进行分析,大大提升了分析效率,缩短了机器人设计开发周期;通过有限元技术和归一化方法提前暴露设计问题和缺陷,避免在样机阶段发生致命性问题,降低了开发成本,可与任何工业机器人设计开发,提高机器人的整机刚度水平的同时还能够保证整机的轻量化设计。
Description
技术领域
本发明涉及机器人刚度优化设计领域,具体是一种基于有限元技术的工业机器人整机刚度优化方法。
背景技术
机器人结构静刚度是机器人重要的性能标准之一,在国标GBT 12642中有明确的规定和测试方法介绍。目前国内的工业机器人设计多为逆仿,缺乏在设计过程中对机器人刚度的考虑,通常设计人员为了规避风险会在设计过程中相对保守,这必然导致了机器人的整机质量过重、各零部件的刚度过设计以及零部件之间的刚度匹配性不好。
如中国专利申请号为201811074857.2的基于有限元的六轴工业机器人动态性能优化方法中,涉及六轴工业机器人动态性能优化方法。本发明为了解决机器人在高速运行过程中定位精度低的问题。本发明包括:一:确定每个零件的低阶模态频率以及刚度分布;二:对每个零件刚度最小的部位进行优化;三:建立六轴工业机器人整机的有限元模型;四:确定六轴工业机器人在三种关节角配置时的前六阶模态频率;五:进行六轴工业机器人刚柔耦合动力学建模;六:得到六轴工业机器人每个零件在作业过程中的最大应力节点、出现最大应力的时刻、应力分布图以及末端弹性偏移量;七:与六轴工业机器人的设计要求中的对比,若符合设计要求则结束,否则重新执行步骤一至六,对每个部位单独进行优化方法,而非在整体的加载情况下进行刚度优化,同时在三维建模软件中按照关节2和关节3的不同配置将六轴工业机器人的零件进行装配,并生成x_t中间文件,步骤繁多,不利于刚度优化。
发明内容
为了解决上述问题,本发明提出一种基于有限元技术的工业机器人整机刚度优化方法。
一种基于有限元技术的工业机器人整机刚度优化方法,其具体步骤如下:
S1:初始模型有限元建模:建立机器人整机有限元模型,计算整机综合刚度Kw以及三方向刚度Kx、Ky、Kz;
S2:子结构划分:将机器人按照杆件和关节分割成i各子结构;
S3:调整子结构刚度参数:在有限元模型中调整各子结构刚度参数,计算调整后的整机综合刚度K’w以及三方向刚度K′x、K′y、K′z;
S4:计算整机刚度变化量:基于以上结果并计算各子结构对于整机刚度的贡献度系数以及三方向刚度的贡献度系数,得到DKw及DK;
S5:判断设计是否合理:根据已有评价标准,判断各子结构刚度是可否有轻量化空间或设计风险;
S6:结构优化:若有可轻量化空间可进行结构有限元优化,减小刚度过设计,针对设计风险子结构需要进一步的结构有限元校核。
所述的步骤S1中的建模是在选定任一工作位姿进行整机有限元建模。
所述的步骤S1中整机综合刚度以及三方向刚度均为机器人末端点在单位力的加载条件下的位移响应。
所述的步骤S2中子结构在有限元模型中需要为独立的零部件。
所述的步骤S3子结构刚度参数针对不同零部件有不同表征形式,即连杆子结构为连杆的弹性模量,关节子结构为关节的扭转刚度。
所述的步骤S5贡献度系数为利用归一化方法计算得出,即先计算出各子结构质量在整机中的质量占比系数,然后计算各子结构单位质量材料对整机的刚度贡献度系数。
所述的步骤S6结构有限元优化为拓扑优化技术。
本发明的有益效果是:利用有限元手段建立整机模型对整机刚度和零部件刚度的设计合理性进行分析,大大提升了分析效率,缩短了机器人设计开发周期;通过有限元技术和归一化方法提前暴露设计问题和缺陷,避免在样机阶段发生致命性问题,降低了开发成本,可与任何工业机器人设计开发,提高机器人的整机刚度水平的同时还能够保证整机的轻量化设计。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1为本发明的流程结构示意图;
图2为本发明的工业机器人结构示意图;
图3为本发明的刚度优化前转座结构示意图;
图4为本发明的刚度优化后转座结构示意图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面对本发明进一步阐述。
如图1至图4所示,一种基于有限元技术的工业机器人整机刚度优化方法,其具体步骤如下:
S1:初始模型有限元建模:建立机器人整机有限元模型,计算整机综合刚度Kw以及三方向刚度Kx、Ky、Kz;
S2:子结构划分:将机器人按照杆件和关节分割成i各子结构,子结构包括底座1、转座2、大臂3、连杆一4、连杆二5、连杆三6、小臂7、手腕体8;
S3:调整子结构刚度参数:通过有限元模型中调整各子结构刚度参数,在机器人末端加载三方向分力:Fx=200N、Fy=-500N以及Fz=-200N,计算调整后的整机综合刚度K’w以及三方向刚度K′x、K′y、K′z,如表1所示;
表1调整前后刚度对照表
S4:a:分别求取各子结构占整机的质量占比系数,得出如下质量占比系数矩阵θ:
θ(1~8)=[0.2109 0.2157 0.0784 0.0193 0.0237 0.0176 0.0421 0.0332]
b:计算整机刚度变化量:基于以上结果并计算各子结构对于整机刚度的贡献度系数以及三方向刚度的贡献度系数,得到DKw及DK;
DKw=[0.0804 0.3901 2.0078 0.0346 0.0988 0.2428 1.0439 0.0424]
S5:判断设计是否合理:根据已有评价标准,判断各子结构刚度是可否有轻量化空间或设计风险,对整机刚度贡献度系数DKw项量中各值求取平均值,可得为0.4926,该值小于1,即目前整机综合刚度过设计,有轻量化空间;
S6:结构优化:若有可轻量化空间可进行结构有限元优化,减小刚度过设计,针对设计风险子结构需要进一步的结构有限元校核:
a:进一步观察三方向刚度贡献度系数DK,在Z方向上刚度贡献度系数最大的转座,其贡献系数为0.6191,该值小于1,因此可考虑进行轻量化设计,同理可针对其余部件分别进行结构有限元优化和校核;
b:所有子结构优化完毕后重新计算新模型的,重新计算各子结构对整机综合刚度以及三方向刚度的贡献度系数D′Kw及D′K;
D′Kw=[0.1766 0.3891 2.8714 0.1256 0.0831 0.1383 0.7629 0.1117]
c:对整机刚度贡献度系数D′Kw项量中各值求取平均值,可得为0.5823,相对于原有的0.4926提升了约18.2%,即优化后整机刚度过设计现象有所改善;
d:进一步观察D′Kw,Z方向上刚度贡献度系数最大的转座,其贡献系数为0.6658,相对于优化前的0.6191有所增加,提升了7.5%,表明拓扑优化后转座的过设计有所改善。
综合以上步骤,即可表明,经过本发明方法的实施后,该型机器人的整机刚度设计更加合理。
所述的步骤S1中的建模是在选定任一工作位姿进行整机有限元建模。
通过有限元技术手段在机器人设计初期针对整机及零部件的刚度设计合理性进行校验和识别,从而有针对性的对相应零部件进行刚度优化或风险校核。本发明的实现能为机器人设计提供新的方法和思路,大大提升机器人产品开发速度和质量,基于此提升国产机器人的市场竞争力。
所述的步骤S1中整机综合刚度以及三方向刚度均为机器人末端点在单位力的加载条件下的位移响应。
所述的步骤S2中子结构在有限元模型中需要为独立的零部件。
所述的步骤S3子结构刚度参数针对不同零部件有不同表征形式,即连杆子结构为连杆的弹性模量,关节子结构为关节的扭转刚度。
所述的步骤S5贡献度系数为利用归一化方法计算得出,即先计算出各子结构质量在整机中的质量占比系数,然后计算各子结构单位质量材料对整机的刚度贡献度系数。
所述的步骤S6结构有限元优化为拓扑优化技术。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (7)
1.一种基于有限元技术的工业机器人整机刚度优化方法,其特征在于:其具体步骤如下:
S1:初始模型有限元建模:建立机器人整机有限元模型,计算整机综合刚度Kw以及三方向刚度Kx、Ky、Kz;
S2:子结构划分:将机器人按照杆件和关节分割成i各子结构;
S3:调整子结构刚度参数:在有限元模型中调整各子结构刚度参数,计算调整后的整机综合刚度K’w以及三方向刚度K′x、K′y、K′z;
S4:计算整机刚度变化量:基于以上结果并计算各子结构对于整机刚度的贡献度系数以及三方向刚度的贡献度系数,得到DKw及DK;
S5:判断设计是否合理:根据已有评价标准,判断各子结构刚度是可否有轻量化空间或设计风险;
S6:结构优化:若有可轻量化空间可进行结构有限元优化,减小刚度过设计,针对设计风险子结构需要进一步的结构有限元校核。
2.根据权利要求1所述的一种基于有限元技术的工业机器人整机刚度优化方法,其特征在于:所述的步骤S1中的建模是在选定任一工作位姿进行整机有限元建模。
3.根据权利要求1所述的一种基于有限元技术的工业机器人整机刚度优化方法,其特征在于:所述的步骤S1中整机综合刚度以及三方向刚度均为机器人末端点在单位力的加载条件下的位移响应。
4.根据权利要求1所述的一种基于有限元技术的工业机器人整机刚度优化方法,其特征在于:所述的步骤S2中子结构在有限元模型中需要为独立的零部件。
5.根据权利要求1所述的一种基于有限元技术的工业机器人整机刚度优化方法,其特征在于:所述的步骤S3子结构刚度参数针对不同零部件有不同表征形式,即连杆子结构为连杆的弹性模量,关节子结构为关节的扭转刚度。
6.根据权利要求1所述的一种基于有限元技术的工业机器人整机刚度优化方法,其特征在于:所述的步骤S5贡献度系数为利用归一化方法计算得出,即先计算出各子结构质量在整机中的质量占比系数,然后计算各子结构单位质量材料对整机的刚度贡献度系数。
7.根据权利要求1所述的一种基于有限元技术的工业机器人整机刚度优化方法,其特征在于:所述的步骤S6结构有限元优化为拓扑优化技术。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010141357.7A CN111444644B (zh) | 2020-03-04 | 2020-03-04 | 一种基于有限元技术的工业机器人整机刚度优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010141357.7A CN111444644B (zh) | 2020-03-04 | 2020-03-04 | 一种基于有限元技术的工业机器人整机刚度优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111444644A true CN111444644A (zh) | 2020-07-24 |
CN111444644B CN111444644B (zh) | 2023-05-23 |
Family
ID=71627264
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010141357.7A Active CN111444644B (zh) | 2020-03-04 | 2020-03-04 | 一种基于有限元技术的工业机器人整机刚度优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111444644B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112571452A (zh) * | 2020-11-13 | 2021-03-30 | 中国船舶重工集团公司第七一六研究所 | 一种基于刚度模型的机器人刚度设计方法与终端 |
CN112861287A (zh) * | 2021-03-05 | 2021-05-28 | 重庆大学 | 一种机器人轻量化效果评价方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005141645A (ja) * | 2003-11-10 | 2005-06-02 | Ohbayashi Corp | 非線形有限要素解析装置及び方法、コンピュータプログラム、記録媒体 |
CN102063548A (zh) * | 2011-01-07 | 2011-05-18 | 西安交通大学 | 一种机床整机动态性能优化设计方法 |
CN105844039A (zh) * | 2016-03-31 | 2016-08-10 | 吉林大学 | 一种车身结构各局部刚度贡献度的评价方法 |
US20160357893A1 (en) * | 2016-08-15 | 2016-12-08 | Xianwu Ling | Contact stiffness estimation based on structural frequency responses |
CN106815407A (zh) * | 2016-12-22 | 2017-06-09 | 四川大学 | 一种数控机床整机动态性能优化方法 |
CN108647405A (zh) * | 2018-04-24 | 2018-10-12 | 华中科技大学 | 多层级点阵结构拓扑优化设计的子结构插值模型建模方法 |
CN109255188A (zh) * | 2018-09-14 | 2019-01-22 | 哈尔滨工业大学 | 基于有限元的六轴工业机器人动态性能优化方法 |
-
2020
- 2020-03-04 CN CN202010141357.7A patent/CN111444644B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005141645A (ja) * | 2003-11-10 | 2005-06-02 | Ohbayashi Corp | 非線形有限要素解析装置及び方法、コンピュータプログラム、記録媒体 |
CN102063548A (zh) * | 2011-01-07 | 2011-05-18 | 西安交通大学 | 一种机床整机动态性能优化设计方法 |
CN105844039A (zh) * | 2016-03-31 | 2016-08-10 | 吉林大学 | 一种车身结构各局部刚度贡献度的评价方法 |
US20160357893A1 (en) * | 2016-08-15 | 2016-12-08 | Xianwu Ling | Contact stiffness estimation based on structural frequency responses |
CN106815407A (zh) * | 2016-12-22 | 2017-06-09 | 四川大学 | 一种数控机床整机动态性能优化方法 |
CN108647405A (zh) * | 2018-04-24 | 2018-10-12 | 华中科技大学 | 多层级点阵结构拓扑优化设计的子结构插值模型建模方法 |
CN109255188A (zh) * | 2018-09-14 | 2019-01-22 | 哈尔滨工业大学 | 基于有限元的六轴工业机器人动态性能优化方法 |
Non-Patent Citations (3)
Title |
---|
姜彪;谢小平;卿宏军;: "相对灵敏度与综合优化策略下的某商用车驾驶室轻量化设计" * |
赵永丰;王均涛;: "基于有限元的天平校准系统基座结构优化设计" * |
那景新;袁正;贺午阳;慕文龙;: "一种车身局部模块刚度贡献度的评价方法" * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112571452A (zh) * | 2020-11-13 | 2021-03-30 | 中国船舶重工集团公司第七一六研究所 | 一种基于刚度模型的机器人刚度设计方法与终端 |
CN112571452B (zh) * | 2020-11-13 | 2022-02-18 | 中国船舶重工集团公司第七一六研究所 | 一种基于刚度模型的机器人刚度设计方法与终端 |
CN112861287A (zh) * | 2021-03-05 | 2021-05-28 | 重庆大学 | 一种机器人轻量化效果评价方法 |
CN112861287B (zh) * | 2021-03-05 | 2022-09-16 | 重庆大学 | 一种机器人轻量化效果评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111444644B (zh) | 2023-05-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106919767A (zh) | 汽车白车身轻量化分析方法 | |
CN109255188B (zh) | 基于有限元的六轴工业机器人动态性能优化方法 | |
CN111444644A (zh) | 一种基于有限元技术的工业机器人整机刚度优化方法 | |
CN109726437B (zh) | 一种舱门气动载荷等效节点力处理方法 | |
CN103921954B (zh) | 基于三轴数控定位器的飞机壁板装配变形的数字化校正方法 | |
CN103950552B (zh) | 基于六轴数控定位器的飞机壁板装配变形的数字化校正方法 | |
CN110704953B (zh) | 一种大展弦比机翼静气弹性能设计敏度的分析方法 | |
CN111310380A (zh) | 一种电动车动力总成悬置橡胶衬套结构的设计开发方法 | |
CN112414668B (zh) | 一种风洞试验数据静气弹修正方法、装置、设备及介质 | |
CN107480357B (zh) | 一种多杆系结构有限元模型边界模拟方法 | |
Meheut et al. | Gradient-based single and multi-points aerodynamic optimizations with the elsA software | |
CN109255141B (zh) | 一种汽车车身正向概念设计截面形状优化方法 | |
Maheri et al. | Combined analytical/FEA-based coupled aero structure simulation of a wind turbine with bend–twist adaptive blades | |
CN109724767B (zh) | 一种基于仿生设计的风洞翼型动态试验模型 | |
CN110837682A (zh) | 一种基于正交试验的工业机器人大臂的结构优化方法 | |
CN108595769B (zh) | 一种基于优化算法的减振器刚度模拟方法 | |
CN108875188B (zh) | 汽车车身接头的优化方法与装置 | |
CN111797549A (zh) | 一种基于超单元和有限元混合模型的nvh分析方法 | |
CN107992709A (zh) | 基于中间函数的热结构模型修正方法 | |
CN114925563A (zh) | 车门开闭系统关键零部件极限载荷的仿真计算方法和系统 | |
CN115293069A (zh) | 一种用于飞行器外流场仿真控制参数智能优化的系统 | |
CN111177855B (zh) | 一种全局气动弹性优化中的气动结构求解方法及系统 | |
CN114912192A (zh) | 一种面向叶盘高保真动力学建模的界面减缩方法 | |
CN117077295B (zh) | 一种基于预连接数字孪生模型的预连接件数量优选方法 | |
CN115329505B (zh) | 具有三维特性的整机二维外部载荷及变形仿真方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |