发明名称 一种体现制造几何缺陷的三维有限元模型修调方法
摘要 本发明属于航空航天结构主承力构件设计技术领域,涉及一种体现制造几何缺陷的三维有限元模型修调方法,包括以下步骤:1)采用光学测量方法测量得到三维点云信息;2)对三维点云信息进行滤波处理,消除噪声;3)建立不考虑三维点云信息的完美结构三维有限元模型;4)在三维点云信息的基础上,或者对得到的三维点云信息进行放大、缩小处理后,采用插值方法对完美结构三维有限元模型进行修调,体现制造几何缺陷。本发明能够增强有限元预测模型的真实性,进一步提高有限元预测的准确性,同时能够保证对已有有限元模型的充分利用,实现自适应缺陷幅值缩放和位置选择,工作量极小,有望成为航空航天领域中有限元模型修正的关键技术之一。
申请公布号 CN106372282A 申请公布日期 2017.02.01
申请号 CN201610712399.5 申请日期 2016.08.24
申请人 大连理工大学 发明人 王博;郝鹏;毕祥军;杜凯繁;任明法;黄诚;蒋亮亮
分类号 G06F17/50(2006.01)I;G06T19/20(2011.01)I 主分类号 G06F17/50(2006.01)I
代理机构 大连理工大学专利中心 21200 代理人 梅洪玉;潘迅
主权项 一种体现制造几何缺陷的三维有限元模型修调方法,其特征在于以下步骤:第一步,采用光学测量方法,对网格加筋筒壳结构试验件的三维形貌进行测量,得到网格加筋筒壳结构试验件的三维点云信息;第二步,将第一步三维点云信息采用滤波方法进行滤波处理,消除噪声;所述的滤波方法包括标准高斯平均、中值滤波算法或超限邻域平均滤波算法;所述的超限邻域平均滤波算法为:设置门限,当三维形貌点云数据与其邻域平均值的差值超过门限时,令其等于该邻域的平均值,反之,数据值不变;所述的超限邻域平均算法的数学表示为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>p</mi><mi>i</mi></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mover><msub><mi>q</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover></mtd><mtd><mrow><mo>|</mo><mover><msub><mi>f</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>-</mo><msub><mi>f</mi><mi>i</mi></msub><mo>|</mo><mo>&gt;</mo><mi>M</mi></mrow></mtd></mtr><mtr><mtd><msub><mi>q</mi><mi>i</mi></msub></mtd><mtd><mrow><mo>|</mo><mover><msub><mi>f</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>-</mo><msub><mi>f</mi><mi>i</mi></msub><mo>|</mo><mo>&le;</mo><mi>M</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001089401130000011.GIF" wi="1301" he="166" /></maths><maths num="0002"><math><![CDATA[<mrow><mover><msub><mi>q</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msub><mi>q</mi><mi>j</mi></msub><msub><mi>v</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>/</mo><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msub><mi>v</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001089401130000012.GIF" wi="1261" he="93" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>v</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>|</mo><mo>|</mo><msub><mi>A</mi><mi>i</mi></msub><mo>-</mo><msub><mi>A</mi><mi>j</mi></msub><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub><mo>&gt;</mo><mi>L</mi></mrow></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mo>|</mo><mo>|</mo><msub><mi>A</mi><mi>i</mi></msub><mo>-</mo><msub><mi>A</mi><mi>j</mi></msub><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub><mo>&lt;</mo><mi>L</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001089401130000013.GIF" wi="1318" he="174" /></maths>其中,p<sub>i</sub>为测点i数据滤波后的值;<img file="FDA0001089401130000014.GIF" wi="43" he="77" />为测点i的邻域平均值;q<sub>i</sub>为测点i的原始半径;M为门限;N为测点数目;v<sub>ij</sub>为测点j对测点i邻域平均值权系数;L为测点i邻域半径;A<sub>i</sub>为测点i的三维坐标向量;第三步,采用有限元分析软件建立不考虑三维点云信息的完美结构三维有限元模型;第四步,在第二步滤波处理后得到的三维点云信息的基础上,或者对得到的三维点云信息进行放大、缩小处理后,采用插值方法对第三步得到的完美结构三维有限元模型整体或局部的节点坐标进行修调,体现制造几何缺陷;4.1)根据插值公式(4),基于三维点云信息和完美结构三维有限元模型,计算得到实测的三维点云信息的权系数;<maths num="0004"><math><![CDATA[<mrow><msub><mi>v</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mfrac><mn>1</mn><msup><mrow><mo>(</mo><mo>|</mo><mo>|</mo><msub><mi>Q</mi><mi>i</mi></msub><mo>-</mo><msub><mi>P</mi><mi>j</mi></msub><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub><mo>)</mo></mrow><mn>4</mn></msup></mfrac></mtd><mtd><mrow><mo>|</mo><mo>|</mo><msub><mi>Q</mi><mi>i</mi></msub><mo>-</mo><msub><mi>P</mi><mi>j</mi></msub><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub><mo>&le;</mo><mi>F</mi></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>|</mo><mo>|</mo><msub><mi>Q</mi><mi>i</mi></msub><mo>-</mo><msub><mi>P</mi><mi>j</mi></msub><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub><mo>&gt;</mo><mi>F</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001089401130000021.GIF" wi="1510" he="270" /></maths>其中,v<sub>ij</sub>是实测的三维点云信息第j个节点对完美结构三维有限元模型中第i个节点修正的权系数;P<sub>j</sub>为测点j的三维坐标向量;Q<sub>i</sub>为完美结构三维有限元模型第i个节点的三维坐标向量;F是完美结构三维有限元模型单元节点的修调距离,根据实测的三维点云信息中的样本点分布密度取值;4.2)根据插值公式(5),基于4.1)得到的三维点云信息的权系数和第二步滤波处理后得到的三维点云信息中的测点半径,计算得到完美结构三维有限元模型的插值半径;<maths num="0005"><math><![CDATA[<mrow><msub><mi>r</mi><mrow><msub><mi>im</mi><mi>i</mi></msub></mrow></msub><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msub><mi>r</mi><msub><mi>p</mi><mi>j</mi></msub></msub><msub><mi>v</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>/</mo><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msub><mi>v</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001089401130000022.GIF" wi="1405" he="126" /></maths>其中,<img file="FDA0001089401130000023.GIF" wi="54" he="54" />是完美结构三维有限元模型中第i个节点的插值半径;<img file="FDA0001089401130000024.GIF" wi="51" he="61" />是实测的三维点云信息中第j个测点的半径;N是实测的三维点云信息中测点的个数;4.3)根据插值公式(6),基于4.2)得到的完美结构三维有限元模型的插值半径和第三步得到的完美结构三维有限元模型,计算得到完美结构三维有限元模型的修调偏移量;基于修调偏移量,移动完美结构三维有限元模型节点坐标,实现完美结构三维有限元模型修调,将真实缺陷引入完美结构三维有限元模型中,获得基于真实缺陷的结构有限元模型,体现制造几何缺陷;<maths num="0006"><math><![CDATA[<mrow><msub><mi>r</mi><mi>i</mi></msub><mo>=</mo><msub><mi>r</mi><mrow><msub><mi>im</mi><mi>i</mi></msub></mrow></msub><mo>-</mo><msub><mi>r</mi><mrow><msub><mi>per</mi><mi>i</mi></msub></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001089401130000025.GIF" wi="1211" he="78" /></maths>其中,r<sub>i</sub>是完美的三维有限元模型中第i个节点的修调偏移量;<img file="FDA0001089401130000026.GIF" wi="56" he="54" />是完美的三维有限元模型中第i个节点的插值半径;<img file="FDA0001089401130000027.GIF" wi="67" he="55" />是完美的三维有限元模型中第i个节点修调前的半径。
地址 116024 辽宁省大连市甘井子区凌工路2号
您可能感兴趣的专利