CN102495964A - 一种二维多工位装配成功率计算方法 - Google Patents
一种二维多工位装配成功率计算方法 Download PDFInfo
- Publication number
- CN102495964A CN102495964A CN2011104066413A CN201110406641A CN102495964A CN 102495964 A CN102495964 A CN 102495964A CN 2011104066413 A CN2011104066413 A CN 2011104066413A CN 201110406641 A CN201110406641 A CN 201110406641A CN 102495964 A CN102495964 A CN 102495964A
- Authority
- CN
- China
- Prior art keywords
- tolerance
- deviation
- deviations
- model
- alignment pin
- 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
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种二维多工位装配成功率计算方法,包括以下步骤:1)建立由夹具定位销公差、零件孔公差所引起的定位偏差模型;2)采用数论网格方法对定位销公差、零件孔公差进行均匀采样,得到公差合适样本空间,并代入定位偏差模型,得到定位偏差样本空间;3)根据偏差流状态空间模型,计算得到偏差传递矩阵,并将定位偏差作为状态空间模型的输入偏差,并计算得到测点输出偏差;4)将计算结果和测点允许偏差进行比较,统计出合格样本数,再将其除以总样本数,则计算得到二维多工位装配成功率。本发明提出了一种面向二维多工位装配,基于偏差流状态空间模型的二维装配成功率计算方法,为二维多工位装配成功率预测提供了一种新的方法。
Description
技术领域
本发明涉及机械制造技术领域,更具体的说是一种二维多工位装配成功率计算方法。
背景技术
装配是产品制造过程中的最后一个环节,是产品生命周期的重要组成部分,也是实现产品功能的主要环节。相关研究表明,产品生产过程中大约三分之一的人力在从事与产品装配有关的活动,装配成本占到整个产品制造成本的40%~50%。由于装配过程占用的劳动量大、费用高,因此,提高装配成功率带来的经济效益将会非常显著。近年来,随着计算机技术的快速发展,利用计算机进行基于装配成功率的机械产品尺寸公差设计,受到国内外学术界越来越广泛的重视。在进行尺寸公差设计时,装配成功率一般作为约束条件。国内外学者就明确提出将装配成功率作为公差优化设计的约束条件。由于当变量维数较高时,直接计算装配成功率非常复杂,所以国内外学者一般采用简化方法或计算机模拟方法来求得装配成功率,并就此开展了大量的研究工作。经过查阅有关技术文献,发现最常用的方法有如下两种:一是基于可靠度β法的装配成功率计算方法;二是基于蒙特卡洛模拟法的装配成功率计算方法。另外亦有学者提出基于全因子试验法、统计法和拉丁方抽样法等装配成功率计算方法。
但是以上研究重点集中在一维装配所形成的一维尺寸链公差优化设计中,在公差优化设计时以成本最低为目标,装配成功率为约束条件,对组成环尺寸公差进行优化。但随着科学技术的发展和顾客对产品质量要求的日益强烈,产品的功能和结构日趋复杂,二维多工位装配已经成为制造企业的重要特征。二维多工位装配过程尺寸偏差传递和累积相当复杂,查找二维多工位装配过程的尺寸链非常困难,若仍根据尺寸链计算装配成功率将不可取。
发明内容
为了解决上述技术问题,本发明提供了一种面向二维多工位装配,基于偏差流状态空间模型的二维装配成功率计算方法。
本发明解决上述的技术问题的技术方案包括以下步骤:
1)建立由夹具定位销公差、零件孔公差所引起的定位偏差模型;
2) 采用数论网格方法对定位销公差、零件孔公差进行均匀采样,得到公差合适样本空间,并代入定位偏差模型,得到定位偏差样本空间;
3) 根据偏差流状态空间模型,计算得到偏差传递矩阵,并将定位偏差作为状态空间模型的输入偏差,并计算得到测点输出偏差;
4) 将计算结果和测点允许偏差进行比较,统计出合格样本数,再将其除以总样本数,则计算得到二维多工位装配成功率。
进一步,所述的步骤1)中由夹具定位销公差、零件孔公差所引起的定位偏差模型的具体方法为,根据夹具定位布局原理及公差“入体原则”,分析四向定位销与零件孔、二向定位销与零件槽的位置关系,建立定位偏差模型,其中模型参数包括,定位销公差T P、零件孔公差T H、四向定位销与零件孔接触定位角α、二向定位销与零件槽接触定位角β,定位销X方向、Z方向偏差。
进一步,所述的步骤2)中采用数论网格方法对定位销公差、零件孔公差进行均匀采样,得到公差合适样本空间,并代入定位偏差模型,得到定位偏差样本空间的具体方法为,先将定位销公差T P和零件孔公差值T H换算成对称的上下偏差 、,在Matlab软件中编程,利用格点点集表中合适的生成矢量或者分圆域法产生好格子点集合,利用对公差进行均匀采样,并将公差样本值代入定位偏差模型,得到定位偏差样本空间。
进一步,所述的步骤3)具体步骤为:根据状态空间模型计算得到偏差传递矩阵M,并将状态方程转化为Y=MU,将定位偏差作为输入偏差U,根据Y=MU计算得到测点输出偏差Y。
进一步,所述的步骤4)的具体步骤为:零件上测点包括X方向和Z方向两个方向的偏差,仿真得到输出偏差,t为测量工位上测点个数,k=1,…,n,n为数论网格法采样样本总数,规定测点允许偏差为[ ],判断是否在范围[ ]中,统计出合格样本数,记为h,则二维多工位装配成功率为P=h/n。
由于采用上述技术方案,本发明的技术效果是:本发明提出了一种面向二维多工位装配,基于偏差流状态空间模型的二维装配成功率计算方法,为二维多工位装配成功率预测提供了一种新的方法。
附图说明
图1为本发明的流程图。
图2为本发明中四向定位销与零件孔位置及坐标方向偏差示意图,该图中(a)为四向定位销与零件孔位置图;(b)为坐标方向偏差图。
图3为本发明中二向定位销与零件槽位置及坐标方向偏差示意图,该图中(a)二向定位销与零件槽位置图;(b)坐标方向偏差图。
图4为本发明车身地板装配定位销平面布局图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示,一种二维多工位装配成功率计算方法,包括以下步骤:
1)建立由夹具定位销公差、零件孔公差所引起的定位偏差模型;
2) 采用数论网格方法对定位销公差、零件孔公差进行均匀采样,得到公差合适样本空间,并代入定位偏差模型,得到定位偏差样本空间;
3) 根据偏差流状态空间模型,计算得到偏差传递矩阵,并将定位偏差作为状态空间模型的输入偏差,并计算得到测点输出偏差;
4) 将计算结果和测点允许偏差进行比较,统计出合格样本数,再将其除以总样本数,则计算得到二维多工位装配成功率。
具体地,所述步骤(1)包括:
根据附图2所示,四向定位销与零件孔位置情况及坐标方向偏差,则定位偏差为:
根据附图3所示,二向定位销与零件槽位置情况及坐标方向偏差,则定位偏差为:
(3)
具体地,所述步骤(2)包括:
式中为好格子点集合中一个好点。
k=1,…,n (7)
将其代入(1)—(4)式得到定位偏差样本空间。
具体地,所述步骤(3)包括:
根据田兆青《多工位薄板装配偏差流传递的状态空间模型》文献可得偏差流状态空间模型表示为:
X(m)=A(m-1)X(m-1)+B(m)U(m)+V(m) (10)
Y(m)=C(m)X(m)+W(m), m=1,2,…N (11)
式中m为工位号,N为总工位数,X(m)和X(m-1)分别表示工位m和m-1上所有偏差组合矢量,U(m)为工位m上夹具定位偏差组合矢量,Y(m)为工位m上与所有测量点相关的观测矢量,V(m)为装配过程中随机噪声因素,W(m)为测量过程中随机噪声因素,A(m-1)为装配系统动态矩阵,B(m)为装配输入矩阵,C(m)为观测矩阵。
则模型可以转化为:
式中U为输入偏差,M为偏差传递矩阵,Y为零件上测点输出偏差。
(13)
t为测量工位上测点个数,k=1,…,n,n为数论网格法采样样本总数。
具体地,所述步骤(4)包括:
规定测点允许上下偏差如下式所示
P=h/n (14)
实施例一
以某汽车企业车身地板总成装配为例,根据附图4所示,零件1为后地板,零件2为中地板,零件3为前地板。,,为四向定位销。,,为两向定位销。,,,为测点。整个装配过程分为三个工位,各工位的定位布局分布如下:
各定位销和测点的名义坐标值如表1、2所示。
(1) 建立由夹具定位销公差、零件孔(槽)公差所引起的定位偏差模型
工位Ⅰ时:
工位Ⅱ时:
工位Ⅲ时:
(2) 采用数论网格方法对定位销公差、零件孔(槽)公差进行均匀采样,得到公差合适样本空间,并代入定位偏差模型,得到定位偏差样本空间
已知定位销公差和零件孔(槽)公差值为0.08mm,换算成对称上下偏差为±0.04mm,定位销和零件孔(槽)分别采样,此时变量个数为6,所以应采用生成矢量产生好格子点集合。根据格点点集表,选取生成矢量(2129;1,41,1681,793,578,279)。应用Matlab软件,编程产生相应的样本矩阵,。将公差样本矩阵代入定位偏差模型中得到定位偏差样本矩阵。
(3) 将定位偏差作为偏差流状态空间模型的输入偏差,并计算得到测点输出偏差。
(4)将计算结果和测点允许的偏差进行比较,统计出合格样本数,再将其除以总样本数,则计算得到二维多工位装配成功率。
规定零件上观测点允许偏差为:
Claims (5)
1.一种二维多工位装配成功率计算方法,包括以下步骤:
1)建立由夹具定位销公差、零件孔公差所引起的定位偏差模型;
2) 采用数论网格方法对定位销公差、零件孔公差进行均匀采样,得到公差合适样本空间,并代入定位偏差模型,得到定位偏差样本空间;
3) 根据偏差流状态空间模型,计算得到偏差传递矩阵,并将定位偏差作为状态空间模型的输入偏差,并计算得到测点输出偏差;
4) 将计算结果和测点允许偏差进行比较,统计出合格样本数,再将其除以总样本数,则计算得到二维多工位装配成功率。
2.根据权利要求1所述的二维多工位装配成功率计算方法,所述的步骤1)中由夹具定位销公差、零件孔公差所引起的定位偏差模型的具体步骤为,根据夹具定位布局原理及公差“入体原则”,分析四向定位销与零件孔、二向定位销与零件槽的位置关系,建立定位偏差模型,其中模型参数包括,定位销公差T P、零件孔公差T H、四向定位销与零件孔接触定位角α、二向定位销与零件槽接触定位角β,定位销X方向、Z方向偏差。
4.根据权利要求1所述的二维多工位装配成功率计算方法,其特征在于:所述的步骤3)的具体步骤为:根据状态空间模型计算得到偏差传递矩阵M,并将状态方程转化为Y=MU,将定位偏差作为输入偏差U,根据Y=MU计算得到测点输出偏差Y。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110406641.3A CN102495964B (zh) | 2011-12-09 | 2011-12-09 | 一种二维多工位装配成功率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110406641.3A CN102495964B (zh) | 2011-12-09 | 2011-12-09 | 一种二维多工位装配成功率计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102495964A true CN102495964A (zh) | 2012-06-13 |
CN102495964B CN102495964B (zh) | 2014-09-17 |
Family
ID=46187789
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110406641.3A Expired - Fee Related CN102495964B (zh) | 2011-12-09 | 2011-12-09 | 一种二维多工位装配成功率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102495964B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413024A (zh) * | 2013-07-18 | 2013-11-27 | 合肥工业大学 | 一种再制造零部件装配尺寸分级选配方法 |
CN105608270A (zh) * | 2015-12-20 | 2016-05-25 | 湖南科技大学 | 一种二维装配工艺系统可靠度计算方法 |
CN111891260A (zh) * | 2020-08-11 | 2020-11-06 | 上汽大众汽车有限公司 | 一种汽车底盘整体合装方法 |
CN112307568A (zh) * | 2020-11-27 | 2021-02-02 | 中国航发沈阳发动机研究所 | 一种基于调整件组尺寸关系的装配可靠性分析方法 |
CN112697041A (zh) * | 2020-12-10 | 2021-04-23 | 北京新风航天装备有限公司 | 一种基于蒙特卡洛法的装配位姿测量精度预评估方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6349237B1 (en) * | 1997-12-23 | 2002-02-19 | The Regents Of The University Of Michigan | Reconfigurable manufacturing system having a production capacity method for designing same and method for changing its production capacity |
CN102128589A (zh) * | 2010-01-20 | 2011-07-20 | 中国科学院自动化研究所 | 一种在轴孔装配中零件内孔方位误差的校正方法 |
-
2011
- 2011-12-09 CN CN201110406641.3A patent/CN102495964B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6349237B1 (en) * | 1997-12-23 | 2002-02-19 | The Regents Of The University Of Michigan | Reconfigurable manufacturing system having a production capacity method for designing same and method for changing its production capacity |
CN102128589A (zh) * | 2010-01-20 | 2011-07-20 | 中国科学院自动化研究所 | 一种在轴孔装配中零件内孔方位误差的校正方法 |
Non-Patent Citations (3)
Title |
---|
WENZHEN HUANG等: "Tolerance Analysis for Design of Multistage manufacturing processes using number-theoretical net method (NT-net)", 《THE INTERNATIONAL JOURNAL OF FLEXIBLE MANUFACTURING SYSTEM》 * |
文泽军等: "基于偏差流分析法(SOVA)的二维多工位装配顺序规划的研究", 《中国机械工程》 * |
文泽军等: "多工位装配过程夹具系统公差和维护综合优化设计", 《机械工程学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413024A (zh) * | 2013-07-18 | 2013-11-27 | 合肥工业大学 | 一种再制造零部件装配尺寸分级选配方法 |
CN105608270A (zh) * | 2015-12-20 | 2016-05-25 | 湖南科技大学 | 一种二维装配工艺系统可靠度计算方法 |
CN105608270B (zh) * | 2015-12-20 | 2018-09-21 | 湖南科技大学 | 一种二维装配工艺系统可靠度计算方法 |
CN111891260A (zh) * | 2020-08-11 | 2020-11-06 | 上汽大众汽车有限公司 | 一种汽车底盘整体合装方法 |
CN111891260B (zh) * | 2020-08-11 | 2021-08-20 | 上汽大众汽车有限公司 | 一种汽车底盘整体合装方法 |
CN112307568A (zh) * | 2020-11-27 | 2021-02-02 | 中国航发沈阳发动机研究所 | 一种基于调整件组尺寸关系的装配可靠性分析方法 |
CN112307568B (zh) * | 2020-11-27 | 2024-06-11 | 中国航发沈阳发动机研究所 | 一种基于调整件组尺寸关系的装配可靠性分析方法 |
CN112697041A (zh) * | 2020-12-10 | 2021-04-23 | 北京新风航天装备有限公司 | 一种基于蒙特卡洛法的装配位姿测量精度预评估方法 |
CN112697041B (zh) * | 2020-12-10 | 2022-03-08 | 北京新风航天装备有限公司 | 一种基于蒙特卡洛法的装配位姿测量精度预评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102495964B (zh) | 2014-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102495964A (zh) | 一种二维多工位装配成功率计算方法 | |
Liu et al. | The spatial distribution of forest carbon sinks and sources in China | |
CN109445305A (zh) | 一种基于数字孪生的装配精度仿真分析方法与系统 | |
CN109145471B (zh) | 基于cad与实测数据共融模型的虚拟装配系统及方法 | |
Huang et al. | Estimates of methane emissions from Chinese rice paddies by linking a model to GIS database | |
CN106354920B (zh) | 风电机组中主轴承轴套强度分析方法及过盈量分析方法 | |
CN104298814B (zh) | 一种基于参数误差累积的齿轮系统性能可靠度计算方法 | |
Luo et al. | A method for axis straightness error evaluation based on improved artificial bee colony algorithm | |
CN104461688B (zh) | 一种详细燃烧化学反应机理骨架简化的方法 | |
CN114996892A (zh) | 一种基于着色Petri网的电力碳足迹计算方法及装置 | |
CN103870658A (zh) | 一种基于动态规划与遗传算法的装配序列规划方法及装置 | |
CN104794309A (zh) | 一种基于仿真技术的农机设计方法 | |
CN112287484B (zh) | 一种基于矢量代理模型的复杂工程系统可靠性设计方法 | |
CN103605849B (zh) | 一种产品设计分析开发环境中联动计算的实现方法 | |
CN117933095A (zh) | 一种基于机器学习的地表发射率实时反演及同化方法 | |
Zhi-hui et al. | CLDAS drive land surface model to simulate latent heat flux in China | |
Fritzsche et al. | Interactive production planning and ergonomic assessment with Digital Human Models–Introducing the Editor for Manual Work Activities (ema) | |
CN103390251A (zh) | 一种电力系统状态估计的量测量权重设置方法 | |
CN109684649A (zh) | 一种基于地形的风速修订方法及系统 | |
CN109753687A (zh) | 一种复合式气动弹性动力学建模方法 | |
CN105740585B (zh) | 一种面向三维夹具设计的公差迭代式优化方法 | |
CN117852353B (zh) | 一种变截面结构的结构应力计算方法 | |
CN103455689B (zh) | 基于泰勒级数展开的区间模型修正方法 | |
CN116663334B (zh) | 一种监测站点大气污染物本地排放贡献的量化识别方法 | |
Mahowald | Collaborative Project: Building improved optimized parameter estimation algorithms to improve methane and nitrogen fluxes in a climate model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140917 Termination date: 20211209 |