发明名称 一种VTL型数控蒙皮拉形机加载轨迹设计方法
摘要 本发明涉及一种VTL型数控蒙皮拉形机加载轨迹设计方法,该方法包括下列步骤:步骤一根据模具几何外形计算执行钳口运动轨迹;步骤二、驱动机构运动参数转化。依据机床使用限制,该算法可验证执行钳口加载轨迹数据以及转化数据是否超出机床运动范围。该方法可以通过理论计算预先给出合理的加载轨迹,进而减少实际试拉次数,提高零件成形质量。
申请公布号 CN101894181B 申请公布日期 2012.09.26
申请号 CN201010213784.8 申请日期 2010.06.30
申请人 上海飞机制造有限公司;北京航空航天大学 发明人 李小强;朱明华;李东升;吴志敏;李慧;桑阳;徐应强
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京市金杜律师事务所 11256 代理人 楼仙英;徐年康
主权项 1.一种VTL型数控蒙皮拉形机加载轨迹设计方法,该方法包括下列步骤:步骤一根据模具几何外形计算执行钳口运动轨迹在二维坐标系中描述模面特征线,原点建立在模具底面中心位置,Xc轴为拉形方向,Zc轴垂直于底面指向模面,此坐标系将作为执行钳口加载轨迹边界点的参考坐标系;找到模面最高点,将模面特征线左右各离散成n段线段单元[...,CD,...],以模面最高点右侧的模面特征线为例,用[x(i),z(i)]表示组成n段线段单元的n+1个离散点的坐标,i=0,1,2,...,n,初始板料用线段oo′表示,其长度为L<sub>r</sub>,预拉后包覆初始板料用线段OO′,其长度为L<sub>p</sub>;这样,拉形过程共分成n+2步,第1步属于预拉阶段,第2步至第n+1步为包覆阶段,最后一步是补拉;(一)、确定变形量的分配方式包覆阶段根据实际情况可以采用四种不同的变形量分配方式,方程如下所示:均匀分配:<maths num="0001"><![CDATA[<math><mrow><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><mn>2</mn><mi>SM</mi></mrow><msup><mi>XM</mi><mn>2</mn></msup></mfrac><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo></mrow></math>]]></maths>其中,0<x<sub>1</sub>≤XM    (1)抛物线函数分配:<maths num="0002"><![CDATA[<math><mrow><msub><mi>y</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mn>1.5</mn><mi>SM</mi></mrow><mi>XM</mi></mfrac><mo>-</mo><mfrac><mrow><mn>1.5</mn><mi>SM</mi></mrow><msup><mi>XM</mi><mn>3</mn></msup></mfrac><msup><msub><mi>x</mi><mn>2</mn></msub><mn>2</mn></msup><mo>,</mo></mrow></math>]]></maths>其中,0<x<sub>2</sub>≤XM    (2)余弦函数分配:<maths num="0003"><![CDATA[<math><mrow><msub><mi>y</mi><mn>3</mn></msub><mo>=</mo><mfrac><mrow><mn>4</mn><mi>SM</mi></mrow><mi>&pi;XM</mi></mfrac><mi>cos</mi><mfrac><msub><mi>&pi;x</mi><mn>3</mn></msub><mrow><mn>2</mn><mi>XM</mi></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>其中,0<x<sub>3</sub>≤XM    (3)负斜率直线分配:<maths num="0004"><![CDATA[<math><mrow><msub><mi>y</mi><mn>4</mn></msub><mo>=</mo><mfrac><mrow><mn>2</mn><mi>SM</mi></mrow><msup><mi>XM</mi><mn>2</mn></msup></mfrac><mrow><mo>(</mo><mi>XM</mi><mo>-</mo><msub><mi>x</mi><mn>4</mn></msub><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中,0<x<sub>4</sub>≤XM    (4)四个方程满足条件:<maths num="0005"><![CDATA[<math><mrow><msubsup><mo>&Integral;</mo><mn>0</mn><mi>XM</mi></msubsup><msub><mi>y</mi><mn>1</mn></msub><mi>dx</mi><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mi>XM</mi></msubsup><msub><mi>y</mi><mn>2</mn></msub><mi>dx</mi><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mi>XM</mi></msubsup><msub><mi>y</mi><mn>3</mn></msub><mi>dx</mi><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mi>XM</mi></msubsup><msub><mi>y</mi><mn>4</mn></msub><mi>dx</mi><mo>=</mo><mi>SM</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>公式(1-5)中:SM为包覆过程中的总变形量,XM为总离散步数;(二)、计算每条切线的长度通过离散点做模面特征线的切线,依据不同的拉形变形量计算每条切线的长度:L<sub>i</sub>=L<sub>i-1</sub>-L<sub>d</sub>+Δ<sub>j</sub>L<sub>i</sub>表示第i个离散点处切线的长度,L<sub>d</sub>表示线段单元的长度,Δ<sub>j</sub>表示包覆阶段中第j步板料应变量分配值;(三)、计算切线的端点即执行钳口加载轨迹边界点忽略摩擦和已包覆板料的变形,假设每步板料处于单向拉形变形模式下,j表示拉形的步数,j=1,2,...,n+2,当j=1时,为预拉阶段;当1<j<n+2时,为包覆阶段,第j步表示已经包覆了第1、2直到i=j-1段线段单元;当j=n+2时为补拉阶段;ε<sub>pre</sub>表示预拉阶段板料的伸长量,ε<sub>post</sub>表示补拉阶段板料的伸长量;用[X(j),Z(j)]表示n+2个执行钳口加载轨迹边界点,j=1,2,...,n+2;包覆阶段中,经过第j步,板料包覆第i段线段单元两端点分别为A、B,坐标分别为[x(i-1),z(i-1)],[x(i),z(i)],板料由线段组合[O..AA′]变为线段组合[O..ABB′],则B′的坐标和包覆角:当1<j<n+2时:<img file="FSB00000731186200021.GIF" wi="1467" he="239" />其中:j=i+1,1≤i≤n;当j=1时:<img file="FSB00000731186200022.GIF" wi="1476" he="241" />当j=n+2时:<img file="FSB00000731186200023.GIF" wi="1483" he="240" />这里,L<sub>i</sub>表示第i个离散点处切线的长度,L<sub>d</sub>表示线段单元的长度,L<sub>r</sub>为初始板料长度;步骤二、驱动机构运动参数转化(一)、检验模具初始高度是否能满足拉形条件模具固定在工作台上,模具安装初始位置必须满足一定的条件:模面特征线最高点必须位于拉形线下方,这样才能保证板料在无接触摩擦的情况下水平预拉形;判断条件:如果H<H′,H′为板料安装位置距工作台最低高度,则满足;(二)计算边界点外延已知执行钳口加载轨迹边界点的坐标和包覆角,欲求外延边界点相对于模具的坐标信息;实现方法有两种:一是定位拉形杆伸长量,控制托架位移;二是定位托架,控制拉形杆动作;方法的选择依据机床的行程能力;方法一的转化公式(9)如下:<maths num="0006"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>X</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>X</mi><mi>i</mi></msub><mo>+</mo><msub><mi>&delta;L</mi><mn>0</mn></msub><msub><mrow><mi>cos</mi><mi>&alpha;</mi></mrow><mi>i</mi></msub></mtd></mtr><mtr><mtd><msubsup><mi>Z</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>Z</mi><mi>i</mi></msub><mo>-</mo><msub><mi>L</mi><mn>0</mn></msub><msub><mrow><mi>sin</mi><mi>&alpha;</mi></mrow><mi>i</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>方法二的转化公式(10)如下:<maths num="0007"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>X</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>X</mi><mn>0</mn></msub><mo>+</mo><msub><mi>&delta;L</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><msubsup><mi>Z</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>Z</mi><mi>i</mi></msub><mo>-</mo><mo>[</mo><mo>|</mo><msub><mi>X</mi><mi>i</mi></msub><mo>-</mo><msubsup><mi>X</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>|</mo><mo>]</mo><mi>tan</mi><msub><mi>&alpha;</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>公式(9、10)中:当外延左边界点时,δ=-1;当外延右边界点时,δ=1;(X<sub>i</sub>,Z<sub>i</sub>)-第i-1个执行钳口加载轨迹边界点坐标;(X′<sub>i</sub>,Z′<sub>i</sub>)-第i-1个外延边界点坐标;α<sub>i</sub>-第i-1个执行钳口加载轨迹边界点处切线仰角;L<sub>0</sub>-初始外延量;(三)检验托架初置是否超出机床使用极限由于托架组件在机床中具有初始位置和行程范围,在给定L<sub>0</sub>数值后,需要检验外延后的外延边界点位置是否在机床使用极限范围内;应用方式一外延边界点时,在拉形过程中,每一步左右侧外延边界点X坐标绝对值须在托架组件行程范围内,判断条件:如果|X′<sub>min</sub>|>Xmin&amp;|X′<sub>max</sub>|<X<sub>max</sub>,X<sub>min</sub>为托架在X方向最小极限值;X<sub>max</sub>为托架在X方向最大极限值,则可以进行下一步,如果不满足,则需要重新预设L<sub>0</sub>数值;应用方式二外延边界点时,预拉形之前,首先要检查外延边界点X坐标绝对值是否在托架组件行程范围内;同时,还需要保证在包覆过程中,外延边界点始终位于执行钳口加载轨迹边界点外侧,判断条件如下:<maths num="0008"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>X</mi><mi>min</mi></msub><mo>&lt;</mo><mo>|</mo><msubsup><mi>X</mi><mn>0</mn><mo>&prime;</mo></msubsup><mo>|</mo><mo>&lt;</mo><msub><mi>X</mi><mi>max</mi></msub></mtd></mtr><mtr><mtd><mo>|</mo><msubsup><mi>X</mi><mn>0</mn><mo>&prime;</mo></msubsup><mo>|</mo><mo>></mo><msub><mrow><mo>|</mo><mi>X</mi><mo>|</mo></mrow><mi>max</mi></msub></mtd></mtr></mtable></mfenced></math>]]></maths>若不满足,则需要重新预设L<sub>0</sub>数值;(四)一次高度匹配第一次外延操作后的左右首末边界点存在高度差,需要以一边高度为基准,重新计算另一边的外延边界点,保证包覆终了状态左右外延边界点处于同一高度;一次高度匹配后,边界点的坐标由如下公式获得,其中公式(11)应用于外延方法一,公式(12)应用于外延方法二;<maths num="0009"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>X</mi><mi>i</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mo>=</mo><msubsup><mi>X</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>+</mo><msub><mrow><mi>&delta;</mi><mi>l</mi><mi>cos</mi><mi>&alpha;</mi></mrow><mi>i</mi></msub></mtd></mtr><mtr><mtd><msubsup><mi>Z</mi><mi>i</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mo>=</mo><msubsup><mi>Z</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>+</mo><msub><mrow><mi>l</mi><mi>sin</mi><mi>&alpha;</mi></mrow><mi>i</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0010"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>X</mi><mi>i</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mo>=</mo><msubsup><mi>X</mi><mi>i</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>Z</mi><mi>i</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mo>=</mo><msubsup><mi>Z</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>+</mo><mo>[</mo><mo>|</mo><msubsup><mi>X</mi><mi>i</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mo>-</mo><msubsup><mi>X</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>|</mo><msub><mrow><mi>tan</mi><mi>&alpha;</mi></mrow><mi>i</mi></msub><mo>]</mo></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow></math>]]></maths>公式11、12中:当外延左边界点时,δ=-1;当外延右边界点时,δ=1;(X′<sub>i</sub>,Z′<sub>i</sub>)-第i-1个外延边界点坐标;(X″<sub>i</sub>,Z″<sub>i</sub>)-第i-1个一次高度匹配后外延边界点坐标;α<sub>i</sub>-第i-1个执行钳口加载轨迹边界点处切线仰角;l-外延量修正值,由一次高度匹配后的拉形末了状态获得;X″<sub>t</sub>-由一次高度匹配后的包覆末了状态获得;(五)二次匹配前的检验在进行数据点的二次高度匹配之前,需要做以下检验:应用外延方法一时,需要进行二次匹配的钳口轨迹线必须位于二次匹配外延边界点轨迹线的内部;判断条件:如果|X′<sub>t</sub>|>|X<sub>t</sub>|,则可以进行二次高度匹配,如果不满足,则需要重新预设L<sub>0</sub>数值;应用外延方法二时,同样要保证需要进行二次匹配的钳口轨迹线必须位于二次匹配外延边界点轨迹线的内部;判断条件:如果|X′<sub>t</sub>|>|X|<sub>max</sub>,则可以进行二次高度匹配,如果不满足,则需要重新预设L<sub>0</sub>数值;(六)二次高度匹配从边界点数较少的一边出发,在相对边逐个匹配出等高边界点;执行钳口加载轨迹边界点从位置A运动到位置B时,切线与模面特征截面线的交点由位置m运动到位置k,当模面离散点足够多时,假设m与k重合;根据已知点A、B、C、D的坐标和切线俯仰角,可以求得点X、Z的坐标和切线俯仰角;(七)检验上顶量是否超出机床使用极限模具固定在工作台上,工作台上升到适用于板料预拉的位置,模具上顶量必须满足一定的条件:在保证模具安装初始位置的前提下,工作台上升的最高台面必须在行程范围之内;判断条件:如果h-H<Hmax,Hmax为工作台上顶最大极限,则满足。
地址 200436 上海市场中路3115号