发明名称 一种二维变曲率型材零件工艺模型的建模方法
摘要 本发明属于飞机制造技术领域,涉及一种二维变曲率型材零件工艺模型的建模方法,通过对轮廓线的离散,将轮廓线离散为多个直线和圆弧段。采用有限元模拟,计算每个离散段的回弹量,通过多次迭代补偿得补偿后的零件轮廓曲线,用于零件的工艺模型定义。有益效果:本发明所采用的方法能够快速、精确设计型材零件的工艺模型,提供给工艺人员用于拉弯模具的设计。
申请公布号 CN102968524A 申请公布日期 2013.03.13
申请号 CN201210437515.9 申请日期 2012.11.06
申请人 西北工业大学 发明人 雷湘衡;王俊彪;刘闯;杨忆湄;谭浩
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种二维变曲率型材零件工艺模型的建模方法,其特征在于步骤如下:步骤1将型材零件轮廓线分割后离散为直线段和圆弧段:取T截面形状、L截面形状、U截面形状和Z截面形状的型材零件成形时与模具外缘面相接触的一条轮廓曲线C,C的起点和终点分别为P<sup>Start</sup>和P<sup>End</sup>,C总长度为l;连接点P<sup>Start</sup>和P<sup>End</sup>得到直线段L;以轮廓曲线C上到直线L距离最大的点P<sup>Top</sup>将轮廓曲线C分割为两个曲线C<sub>1</sub>和C<sub>2</sub>,曲线C<sub>1</sub>起点为P<sup>Top</sup>,终点为P<sup>Start</sup>,长度为l<sub>1</sub>,曲线C<sub>2</sub>起点为P<sup>Top</sup>,终点为P<sup>End</sup>,长度为l<sub>2</sub>;采用下述步骤将曲线C<sub>1</sub>和C<sub>2</sub>离散为直线和圆弧段:步骤1a:以间距Δl等长度离散曲线C<sub>1</sub>,得到C<sub>1</sub>曲线上的N个离散点P<sub>i</sub>(i=1,2,…,N),其中Δl=(5~10)mm,<img file="FDA00002360165300011.GIF" wi="227" he="106" />P<sub>1</sub>为P<sup>Top</sup>,P<sub>n</sub>为离散曲线C<sub>1</sub>的另一端点,离散点P<sub>i</sub>的曲率为K<sub>i</sub>,切线方向为<img file="FDA00002360165300012.GIF" wi="60" he="61" />步骤1b去除P<sub>1</sub>,P<sub>2</sub>,…,P<sub>N</sub>中曲率相同的中间点:对i进行循环(i=2,…,N-1),若K<sub>i-1</sub>=K<sub>i</sub>=K<sub>i+1</sub>成立,则去掉P<sub>i</sub>点;设去除曲率相同的中间点后剩下N<sup>Z</sup>个点,对剩下的N<sup>Z</sup>个点重新编号得<img file="FDA00002360165300013.GIF" wi="293" he="68" />离散点<img file="FDA00002360165300014.GIF" wi="52" he="55" />的曲率为<img file="FDA00002360165300015.GIF" wi="86" he="55" />切线方向为<img file="FDA00002360165300016.GIF" wi="81" he="67" />步骤1c根据离散点<img file="FDA00002360165300017.GIF" wi="52" he="55" />的曲率<img file="FDA00002360165300018.GIF" wi="60" he="55" />将曲线C<sub>1</sub>划分为直线段和圆弧段:对i进行循环(i=2,…,N<sup>Z</sup>),若<img file="FDA00002360165300019.GIF" wi="269" he="54" />成立,则确定一个直线段L,起点为<img file="FDA000023601653000110.GIF" wi="82" he="54" />终点为<img file="FDA000023601653000111.GIF" wi="78" he="54" />直线长度为<img file="FDA000023601653000112.GIF" wi="167" he="92" />否则以一个圆弧来逼近<img file="FDA000023601653000113.GIF" wi="269" he="55" />不成立的相邻离散点<img file="FDA000023601653000114.GIF" wi="285" he="66" />使得逼近精度<img file="FDA000023601653000115.GIF" wi="275" he="68" />而以圆弧逼近相邻离散点<img file="FDA000023601653000116.GIF" wi="304" he="67" />的精度<maths num="0001"><![CDATA[<math><mrow><mover><msub><mi>p</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>i</mi><mo>+</mo><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&OverBar;</mo></mover><mo>></mo><mn>0.1</mn><mo>;</mo></mrow></math>]]></maths>所述圆弧逼近精度<img file="FDA000023601653000118.GIF" wi="111" he="67" />的计算步骤如下:步骤a计算离散点<img file="FDA000023601653000119.GIF" wi="276" he="68" />到逼近圆弧段距离的算术平均值<img file="FDA000023601653000120.GIF" wi="120" he="67" />计算公式为<maths num="0002"><![CDATA[<math><mrow><mover><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>i</mi><mo>+</mo><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mi>i</mi></mrow><mrow><mi>i</mi><mo>+</mo><mi>j</mi></mrow></munderover><mo>|</mo><mi>R</mi><mo>-</mo><mo>|</mo><mover><mrow><msubsup><mi>P</mi><mi>k</mi><mi>Z</mi></msubsup><mi>O</mi></mrow><mo>&RightArrow;</mo></mover><mo>|</mo><mo>|</mo></mrow><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>其中圆心为O,半径计算公式为<maths num="0003"><![CDATA[<math><mrow><mi>R</mi><mo>=</mo><mfrac><mrow><mo>|</mo><mover><mrow><msubsup><mi>P</mi><mi>i</mi><mi>Z</mi></msubsup><msubsup><mi>P</mi><mrow><mi>i</mi><mo>+</mo><mi>j</mi></mrow><mi>Z</mi></msubsup></mrow><mo>&RightArrow;</mo></mover><mo>|</mo></mrow><msqrt><mn>2</mn><mo>&times;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>cos</mi><mo>)</mo></mrow></msqrt></mfrac><mo>,</mo></mrow></math>]]></maths>圆心角计算公式<img file="FDA00002360165300021.GIF" wi="438" he="179" />保留<img file="FDA00002360165300022.GIF" wi="276" he="67" />中的两个端点<img file="FDA00002360165300023.GIF" wi="53" he="55" />和<img file="FDA00002360165300024.GIF" wi="91" he="60" />删除<img file="FDA00002360165300025.GIF" wi="52" he="55" />和<img file="FDA00002360165300026.GIF" wi="65" he="60" />之间的点;得到分为α<sub>1</sub>个离散段Seg<sub>1u</sub>(u=1,2,…,a<sub>1</sub>)的C<sub>1</sub>曲线,当Seg<sub>1u</sub>为圆弧段时,圆弧曲率为K<sub>1u</sub>,圆弧半径大小为<img file="FDA00002360165300027.GIF" wi="211" he="117" />圆弧角度为θ<sub>1u</sub>,圆弧长度l<sub>1u</sub>;当Seg<sub>1u</sub>为直线段时,曲率为K<sub>1u</sub>=0,定R<sub>1u</sub>=0角度为θ<sub>1u</sub>=0,直线长度l<sub>1u</sub>;a<sub>1</sub>+1个离散点<img file="FDA00002360165300028.GIF" wi="312" he="69" />;对曲线C<sub>2</sub>重复步骤1a~步骤1c,得到为a<sub>2</sub>个离散段Seg<sub>2v</sub>(v=1,2,…,a<sub>2</sub>)的C<sub>2</sub>曲线,其中:<img file="FDA00002360165300029.GIF" wi="227" he="106" />P<sub>n</sub>为离散曲线C<sub>2</sub>的另一端点;当Seg<sub>2v</sub>为圆弧段时,圆弧曲率为K<sub>2v</sub>,圆弧半径大小为<img file="FDA000023601653000210.GIF" wi="220" he="116" />圆弧角度为θ<sub>2v</sub>,圆弧长度l<sub>2v</sub>;当Seg<sub>2v</sub>为直线段时,曲率为K<sub>2v</sub>=0,定R<sub>2v</sub>=0角度为θ<sub>2v</sub>=0,直线长度l<sub>2v</sub>;a<sub>2</sub>+1个离散点<img file="FDA000023601653000211.GIF" wi="346" he="68" />将分段结果Seg<sub>1u</sub>,Seg<sub>2v</sub>作为工艺模型原始模型;步骤2:采用数值模拟确定轮廓曲线C<sub>1</sub>和C<sub>2</sub>回弹后曲线<img file="FDA000023601653000212.GIF" wi="66" he="60" />和<img file="FDA000023601653000213.GIF" wi="95" he="60" />根据弧长相等确定工艺模型回弹数据模型<img file="FDA000023601653000214.GIF" wi="286" he="54" />步骤如下:步骤2a型材拉弯数值模拟:根据零件工艺模型,采用有限元分析软件建立有限元分析模型,其中网格尺寸小于等于5mm,网格单元类型采用壳单元,设置C<sub>1</sub>,C<sub>2</sub>的工艺参数在预拉量分别为(0.3%~0.5%)×l<sub>1</sub>、(0.3%~0.5%)×l<sub>2</sub>,补拉量分别为(1%~4%)×l<sub>1</sub>、(1%~4%)×l<sub>2</sub>,弯曲速度为0.5deg/s~2.5deg/s;然后采用显式方法对拉弯过程进行数值模拟,采用隐式方法对回弹过程进行数值模拟;步骤2b工艺模型回弹数据模型确定:提取回弹后轮廓曲线C<sub>1</sub>和C<sub>2</sub>对应的网格节点集set<sub>1</sub>和set<sub>2</sub>,将这些节点在CAD软件中拟合形成样条线spline<sub>1</sub>和spline<sub>2</sub>,spline<sub>1</sub>和spline<sub>2</sub>表示C<sub>1</sub>和C<sub>2</sub>回弹后的曲线<img file="FDA00002360165300031.GIF" wi="66" he="61" />和<img file="FDA00002360165300032.GIF" wi="95" he="61" />根据弧长相等,确定Seg<sub>1u</sub>,Seg<sub>2v</sub>在曲线<img file="FDA00002360165300033.GIF" wi="67" he="60" />和<img file="FDA00002360165300034.GIF" wi="67" he="60" />中相对应分段<img file="FDA00002360165300035.GIF" wi="283" he="55" />如果Seg<sub>1u</sub>,Seg<sub>2v</sub>是一个圆弧段,其圆弧角度为<img file="FDA00002360165300036.GIF" wi="191" he="54" />同时根据,<img file="FDA00002360165300037.GIF" wi="314" he="64" /><img file="FDA00002360165300038.GIF" wi="332" he="64" />计算回弹后的近似圆弧半径<img file="FDA00002360165300039.GIF" wi="213" he="48" />得到工艺模型回弹数据模型;<img file="FDA000023601653000310.GIF" wi="269" he="54" />的圆弧半径为<img file="FDA000023601653000311.GIF" wi="207" he="54" />圆弧角度为<img file="FDA000023601653000312.GIF" wi="307" he="60" />和<img file="FDA000023601653000313.GIF" wi="66" he="60" />对应的离散点分别为<img file="FDA000023601653000314.GIF" wi="308" he="69" />和<img file="FDA000023601653000315.GIF" wi="342" he="69" />步骤3工艺模型计算:步骤3a对离散段Seg<sub>1u</sub>和Seg<sub>2v</sub>进行回弹补偿:当离散段Seg<sub>1u</sub>,Seg<sub>2v</sub>为直线段时不需要补偿;当Seg<sub>1u</sub>,Seg<sub>2v</sub>为圆弧段时,补偿后的圆弧半径为<img file="FDA000023601653000316.GIF" wi="519" he="54" /><maths num="0004"><![CDATA[<math><mrow><msubsup><mi>R</mi><mrow><mn>2</mn><mi>v</mi></mrow><mi>C</mi></msubsup><mo>=</mo><msub><mi>R</mi><mrow><mn>2</mn><mi>v</mi></mrow></msub><mo>-</mo><mi>&lambda;</mi><mo>&times;</mo><mrow><mo>(</mo><msubsup><mi>R</mi><mrow><mn>2</mn><mi>v</mi></mrow><mi>S</mi></msubsup><mo>-</mo><msub><mi>R</mi><mrow><mn>2</mn><mi>v</mi></mrow></msub><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>λ为补偿系数,取值为1;根据<maths num="0005"><![CDATA[<math><mrow><msub><mi>R</mi><mrow><mn>1</mn><mi>u</mi></mrow></msub><msub><mi>&theta;</mi><mrow><mn>1</mn><mi>u</mi></mrow></msub><mo>=</mo><msubsup><mi>R</mi><mrow><mn>1</mn><mi>u</mi></mrow><mi>C</mi></msubsup><msubsup><mi>&theta;</mi><mrow><mn>1</mn><mi>u</mi></mrow><mi>C</mi></msubsup><mo>,</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mi>R</mi><mrow><mn>2</mn><mi>v</mi></mrow></msub><msub><mi>&theta;</mi><mrow><mn>2</mn><mi>v</mi></mrow></msub><mo>=</mo><msubsup><mi>R</mi><mrow><mn>2</mn><mi>u</mi></mrow><mi>C</mi></msubsup><msubsup><mi>&theta;</mi><mrow><mn>2</mn><mi>u</mi></mrow><mi>C</mi></msubsup><mo>,</mo></mrow></math>]]></maths>计算补偿后圆弧角度<img file="FDA000023601653000320.GIF" wi="198" he="54" />得到Seg<sub>1u</sub>,Seg<sub>2v</sub>补偿后的离散段<img file="FDA000023601653000321.GIF" wi="298" he="54" />其圆弧半径为<img file="FDA000023601653000322.GIF" wi="212" he="54" />圆弧角度为<img file="FDA000023601653000323.GIF" wi="227" he="55" />补偿后对应曲线<img file="FDA000023601653000324.GIF" wi="70" he="60" />和<img file="FDA000023601653000325.GIF" wi="70" he="60" />的离散点分别为<img file="FDA000023601653000326.GIF" wi="298" he="68" />和<img file="FDA000023601653000327.GIF" wi="346" he="68" />步骤3b各离散段的补偿:轮廓曲线C<sub>1</sub>和C<sub>2</sub>的补偿以起点P<sub>1</sub>为基准,使补偿后曲线<img file="FDA000023601653000328.GIF" wi="70" he="60" />和<img file="FDA000023601653000329.GIF" wi="70" he="60" />的起点<img file="FDA000023601653000330.GIF" wi="54" he="54" />和<img file="FDA000023601653000331.GIF" wi="54" he="54" />与P<sub>1</sub>是同一个点,且在P<sub>1</sub>点处的切线方向<img file="FDA000023601653000332.GIF" wi="35" he="62" />相同;以任意连续的两个离散段<img file="FDA000023601653000333.GIF" wi="323" he="60" />在连接点<img file="FDA000023601653000334.GIF" wi="102" he="61" />处一阶连续,且曲率方向一致,得到补偿曲线为工艺模型建模参考模型;步骤3c工艺模型设计:用型材的截面形状在CAD软件中,以工艺模型建模参考模型为轮廓线扫略得零件的工艺模型;步骤4工艺模型优化:以<img file="FDA000023601653000335.GIF" wi="70" he="60" />和<img file="FDA000023601653000336.GIF" wi="70" he="60" />取代C<sub>1</sub>和C<sub>2</sub>,重复步骤2中的步骤2a~步骤2b,dis是轮廓曲线C上的点到<img file="FDA000023601653000337.GIF" wi="67" he="60" />和<img file="FDA000023601653000338.GIF" wi="66" he="60" />的最大距离,如果dis≤0.5mm则结束,如果dis>0.5mm则重复步骤3使得dis≤0.5mm;步骤5:如果重复步骤3后仍为dis>0.5mm,则重复步骤4直到dis≤0.5mm结束。
地址 710072 陕西省西安市友谊西路127号