主权项 |
跌落高度拟合方法,其特征在于,包括以下步骤:S1、有效测试点列表维护:跌落测试机进行冲击测试,得到冲击测试点的跌落数据并生成测试点列表,将新的冲击测试点的跌落数据加入已经生成的测试点列表时进行冲突分析,若不冲突,生成新的有效测试点列表后进入步骤S2,若冲突,则按照冲突判定规则删除冲突点,生成新的有效测试点列表后进入步骤S2;S2、确定是否可以三次B样条拟合:根据S1中新的有效测试点列表进行判断是否包含四个以上的有效测试点,若是,则进入步骤S3,若否,则进行跌落测试引导后进入步骤S1;S3、三次B样条拟合:利用三次B样条基函数,根据四个相邻的有效测试点的跌落数据计算出每一段拟合曲线的若干个型值点,然后将密集的型值点连接起来形成拟合的三次B样条曲线;三次B样条基函数为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mo>[</mo><msub><mi>N</mi><mn>0,4</mn></msub><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>,</mo><mo>,</mo><msub><mi>N</mi><mn>1,4</mn></msub><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>N</mi><mn>2,4</mn></msub><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>N</mi><mn>3,4</mn></msub><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>]</mo><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>u</mi></mtd><mtd><msup><mi>u</mi><mn>2</mn></msup></mtd><mtd><msup><mi>u</mi><mn>3</mn></msup></mtd></mtr></mtable></mfenced><mfrac><mn>1</mn><mrow><mn>3</mn><mo>!</mo></mrow></mfrac><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>4</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>-</mo><mn>3</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>3</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>3</mn></mtd><mtd><mo>-</mo><mn>6</mn></mtd><mtd><mn>3</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>-</mo><mn>1</mn></mtd><mtd><mn>3</mn></mtd><mtd><mo>-</mo><mn>3</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo><mrow><mo>(</mo><mn>0</mn><mo>≤</mo><mi>u</mi><mo>≤</mo><mn>1</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000734311930000011.GIF" wi="1688" he="313" /></maths>三次B样条曲线为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>r</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mrow><mo>[</mo><msub><mi>N</mi><mn>0,4</mn></msub><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>N</mi><mn>1,4</mn></msub><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>N</mi><mn>2,4</mn></msub><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>N</mi><mn>3,4</mn></msub><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>]</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>V</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>V</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>V</mi><mrow><mi>i</mi><mo>+</mo><mn>2</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>V</mi><mrow><mi>i</mi><mo>+</mo><mn>3</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mrow><mo>(</mo><mn>0</mn><mo>≤</mo><mi>u</mi><mo>≤</mo><mn>1</mn><mo>)</mo></mrow></mrow></mrow>]]></math><img file="FDA0000734311930000012.GIF" wi="1398" he="319" /></maths>其中,r<sub>i</sub>(u)表示三次B样条曲线第i段的参数方程,N<sub>0,4</sub>(u)表示三次B样条曲线的第0个基函数,N<sub>1,4</sub>(u)表示三次B样条曲线的第1个基函数,N<sub>2,4</sub>(u)表示三次B样条曲线的第2个基函数,N<sub>3,4</sub>(u)表示三次B样条曲线的第3个基函数,V<sub>i</sub>表示三次B样条曲线特征多边形第i个顶点矢量,V<sub>i+1</sub>表示三次B样条曲线特征多边形第i+1个顶点矢量,V<sub>i+2</sub>表示三次B样条曲线特征多边形第i+2个顶点矢量,V<sub>i+3</sub>表示三次B样条曲线特征多边形第i+3个顶点矢量;S4、临界跌落高度求解:根据步骤S3得到的三次B样条曲线,求解三次B样条曲线上HIC=1000对应点的跌落高度;S5、确定是否取得有效临界跌落高度:根据步骤S4得到的跌落高度,在步骤S1中有效测试点列表中寻找,在区间【跌落高度‑500mm,跌落高度】内是否存在二个有效测试点,以及在区间【跌落高度,跌落高度+500mm】内是否存在二个有效测试点,若存在,则进入步骤S6,若不存在,则进行跌落测试引导后进入步骤S1;S6、输出结果:输出得到的三次B样条曲线以及临界跌落高度在得到的三次B样条曲线上HIC=1000对应点的跌落高度。 |