发明名称 一种二维变曲率型材零件工艺模型的建模方法
摘要 本发明属于飞机制造技术领域,涉及一种二维变曲率型材零件工艺模型的建模方法,通过对轮廓线的离散,将轮廓线离散为多个直线和圆弧段。采用有限元模拟,计算每个离散段的回弹量,通过多次迭代补偿得补偿后的零件轮廓曲线,用于零件的工艺模型定义。有益效果:本发明所采用的方法能够快速、精确设计型材零件的工艺模型,提供给工艺人员用于拉弯模具的设计。
申请公布号 CN102968524B 申请公布日期 2016.02.24
申请号 CN201210437515.9 申请日期 2012.11.06
申请人 西北工业大学 发明人 雷湘衡;王俊彪;刘闯;杨忆湄;谭浩
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种二维变曲率型材零件工艺模型的建模方法,其特征在于步骤如下:步骤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="FDA0000663585340000011.GIF" wi="251" he="124" />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="FDA0000663585340000012.GIF" wi="83" he="85" />步骤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="FDA0000663585340000013.GIF" wi="309" he="93" />离散点<img file="FDA0000663585340000014.GIF" wi="64" he="78" />的曲率为<img file="FDA0000663585340000015.GIF" wi="103" he="78" />切线方向为<img file="FDA0000663585340000016.GIF" wi="105" he="94" />步骤1c根据离散点<img file="FDA00006635853400000121.GIF" wi="67" he="76" />的曲率<img file="FDA00006635853400000122.GIF" wi="69" he="70" />将曲线C<sub>1</sub>划分为直线段和圆弧段:对i进行循环(i=2,…,N<sup>Z</sup>),若<img file="FDA0000663585340000017.GIF" wi="280" he="83" />成立,则确定一个直线段L,起点为<img file="FDA0000663585340000018.GIF" wi="97" he="85" />终点为<img file="FDA0000663585340000019.GIF" wi="91" he="78" />直线长度为<img file="FDA00006635853400000110.GIF" wi="186" he="110" />否则以一个圆弧来逼近<img file="FDA00006635853400000111.GIF" wi="286" he="84" />不成立的相邻离散点<img file="FDA00006635853400000112.GIF" wi="300" he="91" />使得逼近精度<img file="FDA00006635853400000113.GIF" wi="311" he="88" />而以圆弧逼近相邻离散点<img file="FDA00006635853400000114.GIF" wi="327" he="91" />的精度<maths num="0001" id="cmaths0001"><math><![CDATA[<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><img file="FDA00006635853400000115.GIF" wi="330" he="88" /></maths>圆弧逼近精度<img file="FDA00006635853400000116.GIF" wi="130" he="88" />的计算步骤如下:步骤a计算离散点<img file="FDA00006635853400000117.GIF" wi="288" he="92" />到逼近圆弧段距离的算术平均值<img file="FDA00006635853400000118.GIF" wi="144" he="88" />计算公式为<img file="FDA00006635853400000119.GIF" wi="458" he="212" />其中圆心为O,半径计算公式为<img file="FDA00006635853400000120.GIF" wi="429" he="198" />圆心角计算公式<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&theta;</mi><mo>=</mo><mi>arccos</mi><mrow><mo>(</mo><mfrac><mrow><mover><msubsup><mi>V</mi><mi>i</mi><mi>Z</mi></msubsup><mo>&RightArrow;</mo></mover><mo>&CenterDot;</mo><mover><msubsup><mi>V</mi><mrow><mi>i</mi><mo>+</mo><mi>j</mi></mrow><mi>Z</mi></msubsup><mo>&RightArrow;</mo></mover></mrow><mrow><mo>|</mo><mover><msubsup><mi>V</mi><mi>i</mi><mi>Z</mi></msubsup><mo>&RightArrow;</mo></mover><mo>|</mo><mo>|</mo><mover><msubsup><mi>V</mi><mrow><mi>i</mi><mo>+</mo><mi>j</mi></mrow><mi>Z</mi></msubsup><mo>&RightArrow;</mo></mover><mo>|</mo></mrow></mfrac><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000663585340000021.GIF" wi="458" he="205" /></maths>保留<img file="FDA0000663585340000022.GIF" wi="284" he="91" />中的两个端点<img file="FDA0000663585340000023.GIF" wi="60" he="82" />和<img file="FDA0000663585340000024.GIF" wi="114" he="91" />删除<img file="FDA0000663585340000025.GIF" wi="60" he="82" />和<img file="FDA0000663585340000026.GIF" wi="78" he="91" />之间的点;得到分为a<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="FDA0000663585340000027.GIF" wi="226" he="147" />圆弧角度为θ<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="FDA0000663585340000028.GIF" wi="329" he="93" />;对曲线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="FDA0000663585340000029.GIF" wi="247" he="129" />P<sub>N</sub>为离散曲线C<sub>2</sub>的另一端点;当Seg<sub>2v</sub>为圆弧段时,圆弧曲率为K<sub>2v</sub>,圆弧半径大小为<img file="FDA00006635853400000210.GIF" wi="242" he="140" />圆弧角度为θ<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="FDA00006635853400000211.GIF" wi="330" he="94" />;将分段结果Seg<sub>1u</sub>,Seg<sub>2v</sub>作为工艺模型原始模型;步骤2:采用数值模拟确定轮廓曲线C<sub>1</sub>和C<sub>2</sub>回弹后曲线<img file="FDA00006635853400000213.GIF" wi="82" he="81" />和<img file="FDA00006635853400000212.GIF" wi="108" he="81" />根据弧长相等确定工艺模型回弹数据模型<img file="FDA00006635853400000214.GIF" wi="296" he="78" />步骤如下:步骤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="FDA0000663585340000031.GIF" wi="81" he="86" />和<img file="FDA0000663585340000032.GIF" wi="108" he="88" />根据弧长相等,确定Seg<sub>1u</sub>,Seg<sub>2v</sub>在曲线<img file="FDA0000663585340000033.GIF" wi="82" he="89" />和<img file="FDA0000663585340000034.GIF" wi="80" he="93" />中相对应分段<img file="FDA0000663585340000035.GIF" wi="300" he="77" />如果Seg<sub>1u</sub>,Seg<sub>2v</sub>是一个圆弧段,其圆弧角度为<img file="FDA0000663585340000036.GIF" wi="206" he="81" />同时根据<img file="FDA0000663585340000037.GIF" wi="328" he="92" /><img file="FDA0000663585340000038.GIF" wi="410" he="87" />计算回弹后的近似圆弧半径<img file="FDA0000663585340000039.GIF" wi="220" he="78" />得到工艺模型回弹数据模型;<img file="FDA00006635853400000310.GIF" wi="294" he="84" />的圆弧半径为<img file="FDA00006635853400000311.GIF" wi="218" he="80" />圆弧角度为<img file="FDA00006635853400000312.GIF" wi="320" he="94" />和<img file="FDA00006635853400000313.GIF" wi="81" he="86" />对应的离散点分别为<img file="FDA00006635853400000314.GIF" wi="330" he="93" />和<img file="FDA00006635853400000315.GIF" wi="358" he="98" />步骤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="FDA00006635853400000316.GIF" wi="535" he="77" /><maths num="0003" id="cmaths0003"><math><![CDATA[<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><img file="FDA00006635853400000317.GIF" wi="552" he="83" /></maths>λ为补偿系数,取值为1;根据<maths num="0004" id="cmaths0004"><math><![CDATA[<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><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>v</mi></mrow><mi>C</mi></msubsup><msubsup><mi>&theta;</mi><mrow><mn>2</mn><mi>v</mi></mrow><mi>C</mi></msubsup><mo>,</mo></mrow>]]></math><img file="FDA00006635853400000318.GIF" wi="772" he="91" /></maths>计算补偿后圆弧角度<img file="FDA00006635853400000319.GIF" wi="217" he="85" />得到Seg<sub>1u</sub>,Seg<sub>2v</sub>补偿后的离散段<img file="FDA00006635853400000320.GIF" wi="316" he="83" />其圆弧半径为<img file="FDA00006635853400000321.GIF" wi="230" he="80" />圆弧角度为<img file="FDA00006635853400000322.GIF" wi="216" he="78" />补偿后对应曲线<img file="FDA00006635853400000323.GIF" wi="86" he="88" />和<img file="FDA00006635853400000324.GIF" wi="82" he="90" />的离散点分别为<img file="FDA00006635853400000325.GIF" wi="333" he="98" />和<img file="FDA00006635853400000326.GIF" wi="355" he="87" />步骤3b各离散段的补偿:轮廓曲线C<sub>1</sub>和C<sub>2</sub>的补偿以起点P<sub>1</sub>为基准,使补偿后曲线<img file="FDA00006635853400000327.GIF" wi="98" he="86" />和<img file="FDA00006635853400000328.GIF" wi="88" he="86" />的起点<img file="FDA00006635853400000329.GIF" wi="58" he="84" />和<img file="FDA00006635853400000330.GIF" wi="64" he="74" />与P<sub>1</sub>是同一个点,且在P<sub>1</sub>点处的切线方向<img file="FDA00006635853400000331.GIF" wi="42" he="84" />相同;以任意连续的两个离散段<img file="FDA00006635853400000332.GIF" wi="328" he="90" />在连接点<img file="FDA00006635853400000333.GIF" wi="119" he="87" />处一阶连续,且曲率方向一致,得到补偿曲线为工艺模型建模参考模型;步骤3c工艺模型设计:用型材的截面形状在CAD软件中,以工艺模型建模参考模型为轮廓线扫略得到零件的工艺模型;步骤4工艺模型优化:以<img file="FDA00006635853400000334.GIF" wi="86" he="85" />和<img file="FDA00006635853400000335.GIF" wi="80" he="90" />取代C<sub>1</sub>和C<sub>2</sub>,重复步骤2中的步骤2a~步骤2b,dis是轮廓曲线C上的点到<img file="FDA0000663585340000041.GIF" wi="84" he="85" />和<img file="FDA0000663585340000042.GIF" wi="91" he="85" />的最大距离,如果dis≤0.5mm则结束,如果dis>0.5mm则重复步骤3使得dis≤0.5mm;步骤5:如果重复步骤3后仍为dis>0.5mm,则重复步骤4直到dis≤0.5mm结束。
地址 710072 陕西省西安市友谊西路127号