法律状态公告日
法律状态信息
法律状态
2020-05-22
未缴年费专利权终止 IPC(主分类):G06F19/00 授权公告日:20171010 终止日期:20190529 申请日:20140529
专利权的终止
2017-10-10
授权
授权
2014-10-08
实质审查的生效 IPC(主分类):G06F19/00 申请日:20140529
实质审查的生效
2014-09-03
公开
公开
技术领域
本发明涉及柔性零件装配尺寸偏差分析领域,具体是一种柔性装配统计公差分析方法。
背景技术
在装配过程中,柔性零件受到各种复杂装配力的作用,零件受力下产生变形以及力释放后产生回弹,其装配尺寸公差的形成机理不同于刚性零件装配。刚性零件由于公差的存在,在装配过程中会产生装配间隙或干涉,而柔性零件可通过受力变形克服由于公差产生的间隙和干涉,但同时引入的装配应力将导致装配体发生变形。因此需要运用有效的方法分析柔性零件装配过程中的变形和偏差,预测最终的装配偏差。
在柔性装配系统中,除了存在零件制造、夹具定位、连接定位等几何误差源之外,还存在材料固有特性分散性和载荷随机性等不确定性因素。在基于线弹性变形假设的柔性件装配偏差分析中,弹性模量、泊松比等材料属性表征了零件的弹性变形能力,此类材料属性的随机误差直接引起装配闭合力的偏差,进而导致随机的柔性件装配尺寸。基于确定性的有限元分析方法假设材料属性在结构分析中保持不变,但从实际情况来看,不同批次(甚至同一批次)的材料间其弹性模量、泊松比等材料性能参数具有一定的随机变化特性。2002年,Mortensen研究了柔性件装配偏差统计分析方法,提出了继承刚体和柔性体几何协方差的装配尺寸公差的统一表示法,但是没有考虑材料属性随机误差对装配偏差的影响。2008年,于奎刚等采用影响系数法,将板料的厚度偏差这项材料属性纳入传统基于几何偏差分析的柔性装配偏差模型,给出了装配偏差和材料厚度偏差的定量关系,但是也没有考虑材料弹性模量和泊松比随机误差对装配偏差的影响。因此,建立一种考虑材料的弹性模量和泊松比随机误差的柔性装配统计公差分析方法,将有助于科学、准确地确定材料属性误差对柔性装配偏差的影响,提高工程应用的科学性。
发明内容
本发明为了解决现有柔性装配统计公差分析中不考虑材料的弹性模量和泊松比随机误差问题,提供了一种柔性装配统计公差分析方法,使柔性件装配的偏差预测更加准确,更加符合工程实际。
本发明包括如下步骤:
1)采集数据,得到柔性零件几何尺寸和偏差,运用统计公差分析方法(STA)得到零件闭合间隙的均值{μδ}j和方差
2)结合柔性件表面连续性模型,进行协方差(COV)计算,得到闭合间隙协方差矩阵[Σδ]j=[S]j[ΣS]j[S]Tj,式中,[S]j是根据表面连续性条件推导出的每个零件j的敏感度矩阵,[ΣS]j是由
3)根据零件几何尺寸和材料属性的均值,运用有限元分析方法(FEM)和超元刚度矩阵理论,得到柔性零件的缩减刚度矩阵[Kred]j;
4)建立考虑材料弹性模量和泊松比随机误差的零件刚度矩阵计算模型,在零件材料弹性模量和泊松比均值处,对[Kred]j进行一阶泰勒级数展开。
5)建立零件装配力的均值{μF}j和协方差矩阵[ΣF]j以及装配体的回弹力均值{μF}asm,S和协方差矩阵[ΣF]asm,S计算模型,根据弹性力学基本方程{F}=[K]{δ}可知弹性模量Ej与泊松比υj相互独立,由材料参数与几何参数的相互独立性,可知[Kred]j与{μδ}j相互独立,用符号E[]表示求数学期望,得到:
其中
同理
然而
根据>
>
得到>
>式中,“-”表示均值,
6)根据零件几何尺寸和材料属性的均值,运用有限元分析方法和超元刚度矩阵理论,得到装配体的缩减刚度矩阵[Kred]sam;
7)建立考虑材料弹性模量和泊松比随机误差的装配体刚度矩阵计算模型,在零件
8)建立装配回弹偏差均值{μδ}asm和协方差矩阵[Σδ]asm计算模型,根据{μF}asm,S=[Kred]asm{μδ}asm,得到{μδ}asm=[Kred]asm-1{μF}asm,S。
>
由于
其中,>
>
式中只有[Σδ]asm是未知量,是一个线性矩阵方程,可以利用矩阵的Kronecker积和拉直求得,令
[Q]vec([Σδ]asm)=vec([Y])
其中,>
>
符号
所以,vec([Σδ]asm)=[Q]-1vec([Y])
vec([Σδ]asm)的解即为[Σδ]asm中各元素的值,且[Q]非奇异,[Σδ]asm有唯一解。
所述的步骤4)中对[Kred]j进行一阶泰勒级数展开,得到:
>,其中,>式中,
所述步骤7)中对[Kred]sam进行一阶泰勒级数展开,得到:
>
,其中,>式中,
本发明有益效果在于:建立一种考虑材料的弹性模量和泊松比随机误差的柔性装配统计公差分析方法,有助于科学、准确地确定材料属性误差对柔性装配偏差的影响,提高工程应用的科学性。
附图说明
图1为考虑材料属性随机误差的柔性装配统计公差分析方法流程图;
图2为柔性薄板件装配定位阶段截面图;
图3为柔性薄板件装配夹紧阶段截面图;
图4为柔性薄板件装配连接阶段截面图;
图5为柔性薄板件装配释放阶段截面图。
具体实施方式
下面结合附图对本发明作进一步说明。
本发明的考虑材料属性随机误差的柔性装配统计公差分析方法,流程图如图1所示,具体实施按以下步骤实施:
步骤1:如图2所示柔性薄板件装配定位阶段,用三坐标测量机对薄板件采集数据,得到零件A、B的几何尺寸和偏差δA、δB,运用统计公差分析方法(STA)得到零件闭合间隙的均值{μδ}A、{μδ}B和方差
步骤2:结合柔性件表面连续性模型,进行协方差(COV)计算,得到闭合间隙协方差矩阵[Σδ]A、[Σδ]B;
[Σδ]A=[S]A[ΣS]A[S]TA (1)
[Σδ]B=[S]B[ΣS]B[S]TB (2)
式中,[S]A、[S]B是根据表面连续性条件推导出的零件A、B的敏感度矩阵;[ΣS]A、[ΣS]B是分别由
步骤3:根据零件几何尺寸和材料属性的均值,运用有限元分析方法(FEM)和超元刚度矩阵理论,得到柔性零件的缩减刚度矩阵[Kred]A、[Kred]B;
步骤4:建立考虑材料弹性模量和泊松比随机误差的零件刚度矩阵计算模型;
在零件材料弹性模量和泊松比均值处,对[Kred]A、[Kred]B进行一阶泰勒级数展开,即
>
>
其中,>
>
步骤5:如图3和图4所示,在薄板件夹紧阶段和连接阶段,建立零件装配力均值{μF}A、{μF}B和协方差矩阵[ΣF]A、[ΣF]B以及装配体回弹力均值{μF}asm,S和协方差矩阵[ΣF]asm,S计算模型;
其中,>
>
>
>
装配回弹力的均值和协方差分别为
{μF}asm,S=-{μF}A-{μF}B (9)
[ΣF]asm,S=[ΣF]A+[ΣF]B (10)
步骤6:运用有限元分析方法和超元刚度矩阵理论,得到装配体的缩减刚度矩阵[Kred]asm;
步骤7:建立考虑材料弹性模量和泊松比随机误差的装配体刚度矩阵计算模型;
>
其中,>
>
步骤8:如图5所示,在薄板件装配释放阶段,建立装配回弹偏差均值{μδ}asm和协方差矩阵[Σδ]asm计算模型。
{μδ}asm=[Kred]asm-1{μF}asm,S (12)
>
其中,>
>
>
>
令
[Q]vec([Σδ]asm)=vec([Y]) (14)
其中,
>
[Y]=[ΣF]asm-([D]A,asm+[H]A,asm)-([D]B,asm+[H]B,asm)
符号表示Kronecker积,vec()表示矩阵的列拉直(列展开);所以,
vec(Σδ]asm)=[Q]-1vec([Y]) (15)
装配回弹偏差均值{μδ}asm和协方差矩阵[Σδ]asm可以分别由式(12)和式(15)求得。
本发明具体应用途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进,这些改进也应视为本发明的保护范围。
机译: Insolation概率分布分析方法,Insolation概率分布分析系统,Insolation概率分布分析程序产品,Insolation归一化统计分析方法,Isolation归一化统计分析系统和缺失归一化统计分析计划产品
机译: Insolation概率分布分析方法,Insolation概率分布分析系统,Insolation概率分布分析程序产品,Insolation归一化统计分析方法,Isolation归一化统计分析系统和缺失归一化统计分析计划产品
机译: 统计时序分析的细胞库准备装置,统计时序分析器,统计时序分析的细胞库制备方法和统计时序分析方法