CN107194051B - 一种后桥系统非线性强度分析方法 - Google Patents
一种后桥系统非线性强度分析方法 Download PDFInfo
- Publication number
- CN107194051B CN107194051B CN201710330088.7A CN201710330088A CN107194051B CN 107194051 B CN107194051 B CN 107194051B CN 201710330088 A CN201710330088 A CN 201710330088A CN 107194051 B CN107194051 B CN 107194051B
- Authority
- CN
- China
- Prior art keywords
- point
- error
- data
- rear axle
- straight line
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Body Structure For Vehicles (AREA)
- Complex Calculations (AREA)
Abstract
一种后桥系统非线性强度分析方法,包括1、对后桥材料曲线进行处理;2、设置误差判定值eb;3、生成初始点计数器i=1,结束点计数器j=3;4、连接初始点和结束点,生成判定直线;5、计算初始点与结束点之间的原始数据与判定直线的误差e;6、对计算误差进行判定:如果都小于eb,则结束点计数器j=j+1,存储误差的数组清零,返回第4步;如果至少有一个大于eb,则初始点计数器i=j‑1,结束点计数器j=j+1,存储误差的数组清零,返回第4步;7、将每次生成判定直线的初始点及第1步处理后的材料曲线数据的最后一点作为保留点,其他点删除,得到优化的材料数据。本发明对材料应力应变曲线、衬套刚度曲线进行合理缩减,提高计算效率。
Description
技术领域
本发明涉及汽车CAE分析领域,具体涉及一种后桥系统非线性强度分析优化工具。
背景技术
汽车后桥系统有连接副多、结构承载高的特点。近年来,汽车后桥强度分析逐渐由单零部件分析向系统级分析发展;由线性分析向非线性分析发展。通过对整个后桥系统的分析,能够得到各个零部件间的连接关系及各个零部件的应力大小,合理地把控后桥系统性能。
在后桥CAE求解过程中,大量使用非线性迭代算法,非线性迭代算法主要是通过牛顿迭代、拟牛顿迭代等迭代法,将非线性数据分段线性化,再进行后桥系统线性方程或方程组求解。由于每一次的迭代计算就等于进行了一次后桥系统线性方程组的求解,所以非线性计算往往是线性计算时间的数倍甚至数百倍。在后桥系统分析过程中,材料应力应变曲线、衬套刚度曲线是主要的两类非线性数据。而这两类非线性数据量的大小,直接影响了后桥系统分析的计算效率。在进行求解前,需进行合理的数据压缩,尽可能的减少数据量,从而减少非线性迭代的计算时间,同时还需保证压缩数据相对原始数据的精度。
现有的材料应力应变曲线、衬套刚度曲线压缩方法一般采用多项式拟合缩减法,即将原始数据用选定了阶次的多项式拟合出一条曲线,再通过曲线上等间距取点,得到压缩后的非线性数据。
多项式拟合缩减法的使用过程中,有两个步骤引入与原始数据的误差:第一步在多项式拟合时,多项式的阶次一般可以取到5阶或6阶,依然会产生拟合曲线与原始数据的误差,且误差大小和位置不易计算,往往只能靠人眼观察识别;第二步在等间距取点时,等间距取点在拟合曲线曲率较高的位置误差会较大,且误差大小和位置不易计算,也只能靠人眼观察识别。
发明内容
本发明的目的是提供一种后桥系统非线性强度分析方法,通过对非线性数据进行优化,在保证误差控制的基础上,对材料应力应变曲线、衬套刚度曲线进行合理缩减,提高后桥系统非线性强度分析的计算效率。
本发明所述的后桥系统非线性强度分析方法,其步骤如下:
第一步,对测得的后桥材料曲线进行处理,删除纵轴数值下降的点,如果出现多个同一横坐标点,则仅保留纵轴数值最高的点。
第二步,根据材料曲线数据的精度需要,设置误差判定值eb。材料曲线误差一般设置eb为1-5Mpa。
第三步,生成初始点计数器i=1,结束点计数器j=3。
第四步,连接初始点结束点,生成判定直线。
第五步,计算初始点与结束点之间的原始数据与判定直线的误差e。如果初始点与结束点间有n个原始数据点,则计算每个原始数据点的误差e1,e2…en。误差ei的计算方法:过原始数据点作横坐标的垂线,与判定直线相交于一点,该点的纵轴坐标与原始数据点的纵轴坐标的差值即是误差e。
第六步,对第五步所述的计算误差e1,e2…en进行判定。
如果e1,e2…en都小于eb,则结束点计数器j=j+1,存储误差e1,e2…en的数组清零,并且返回到第四步。
如果e1,e2…en至少有一个大于eb,则初始点计数器i=j-1,结束点计数器j=j+1,存储误差e1,e2…en的数组清零,并返回到第四步。
第七步,所有数据均进行判定后,每次生成判定直线的初始点及数据的最后一点,作为保留点,原始数据其他点删除,得到优化后的材料数据。
本发明的有益效果:可以根据材料数据的精度需要,自行设置材料应力应变曲线、衬套刚度数据的误差eb大小,保证压缩后数据任意横轴坐标对应的纵轴坐标的误差在可控范围内。由于该方法是以误差判断进行迭代计算的,所以在曲率较大的位置会自动多保留一些原始数据,而在接近直线的线性区域便会保留较少的原始数据,实现了保证误差精度的同时,最大限度的压缩了非线性材料数据量,提高数值计算中,非线性求解的效率。本方法也不会出现误差积累现象,由于每经过一段时间迭代,由于e的逐渐增大,会导致上述第六步中判定生成新的判定直线,而新的判定直线初始点的取值来自于原始数据,即经过判定迭代过程后总会回归到初始数据上,在该点与原始数据的误差即为0,误差得到了修正。
附图说明
图1为本发明的分析逻辑流程图。
图2、图3、图4、图5为本发明的迭代例子示意图。
图6为优化后与优化前的材料应力应变数据对比。
具体实施方式
以下结合附图进一步说明本发明的技术实现:
参见图1,后桥系统非线性强度分析方法的逻辑步骤如下
第一步,对测得的后桥材料曲线进行处理,删除纵轴数值下降的点,如果出现多个同一横坐标点,则仅保留纵轴数值最高的点。
第二步,根据材料曲线数据的精度需要,设置误差判定值eb,按照经验,选择1-5Mpa。
第三步,生成初始点计数器i=1,结束点计数器j=3。
第四步,连接初始点结束点,生成判定直线。
第五步,计算初始点与结束点之间的原始数据与判定直线的误差e。如果初始点与结束点间有n个原始数据点,则计算每个原始数据点的误差e1,e2…en。误差e的计算方法:过原始数据点作横坐标的垂线,与判定直线相交于一点,该点的纵轴坐标与原始数据点的纵轴坐标的差值即是误差e。
第六步,对第五步所述的计算误差e1,e2…en进行判定。
如果e1,e2…en都小于eb,则结束点计数器j=j+1,存储误差e1,e2…en的数组清零,并且返回到第四步。
如果e1,e2…en至少有一个大于eb,则初始点计数器i=j-1,结束点计数器j=j+1,存储误差e1,e2…en的数组清零,并返回到第四步。
第七步,所有数据均进行判定后,每次生成判定直线的初始点及数据的最后一点,作为保留点,原始数据其他点删除,得到优化后的材料数据。
图2-图5是本发明方法的一个材料曲线缩减迭代例子示意图:
图2中判定直线的初始点与结束点分别为材料曲线的第一点与第三点,第二点与判定直线产生的误差为e1。对e1进行判定,假设e1小于判定误差值eb,那么认为第二点是可以删除的点。
图3中由于图2已判定第二点可以删除,所以判定直线的初始点与结束点更改为第一点和第四点,存储误差e1,e2…en的数组清零。第二点、第三点与新的判定直线产生新的误差为e1,e2。对e1,e2进行判定,假设此次判定结果为e1大于误差判定值eb,那么认为该判定直线误差过大,所以需保留第三点。
图4中由于图3判定结果为误差过大,所以判定直线的初始点需从新生成,即从第三点作为起始点,重新开始迭代生成判定直线及判定的过程。
图5为经过三轮迭代判定后,第四点、第五点、第六点、第七点与每次生成的判定直线的误差均小于误差判定值eb,那么现在判定直线的初始点依然为第三点,结束点为第八点。
图6为经过数轮迭代后,从材料曲线中选出了几个特征点保留,这几个特征点所构成的曲线与原始数据的误差在误差判定值eb以内。可以看到保留的曲线与原始数据相比,在曲率较大的地方相应保留较多的点,以保证数据精度,在接近线性段部分则保留较少的点,同样满足数据精度。
采用本发明优化方法处理后桥系统的材料及衬套刚度,计算后桥强度分析算例,计算结果效率对比如表1:
表1
算例分析精度对比如表2:
表2
Claims (3)
1.一种后桥系统非线性强度分析方法,其步骤如下:
第一步,对测得的后桥材料曲线进行处理,删除纵轴数值下降的点,如果出现多个同一横坐标点,则仅保留纵轴数值最高的点;
第二步,根据材料曲线数据的精度需要,设置误差判定值eb;
第三步,生成初始点计数器i=1,结束点计数器j=3;
第四步,连接初始点和结束点,生成判定直线;
第五步,计算初始点与结束点之间的原始数据与判定直线的误差e,如果初始点与结束点间有n个原始数据点,则计算每个原始数据点的误差e1,e2…en;
第六步,对第五步所述的计算误差e1,e2…en进行判定:
如果e1,e2…en都小于eb,则结束点计数器j=j+1,存储误差e1,e2…en的数组清零,并且返回到第四步;
如果e1,e2…en至少有一个大于eb,则初始点计数器i=j-1,结束点计数器j=j+1,存储误差e1,e2…en的数组清零,并返回到第四步;
第七步,所有数据均进行判定后,将每次生成判定直线的初始点及上述第一步处理后的材料曲线数据的最后一点,作为保留点,原始数据其他点删除,得到优化后的材料数据。
2.根据权利要求1所述的后桥系统非线性强度分析方法,其特征在于,误差e的计算方法为:过原始数据点作横坐标的垂线,与判定直线相交于一点,该点的纵轴坐标与原始数据点的纵轴坐标的差值即是误差e。
3.根据权利要求1或2所述的后桥系统非线性强度分析方法,其特征在于,误差判定值eb设置为1-5Mpa。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710330088.7A CN107194051B (zh) | 2017-05-11 | 2017-05-11 | 一种后桥系统非线性强度分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710330088.7A CN107194051B (zh) | 2017-05-11 | 2017-05-11 | 一种后桥系统非线性强度分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107194051A CN107194051A (zh) | 2017-09-22 |
CN107194051B true CN107194051B (zh) | 2020-10-16 |
Family
ID=59872473
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710330088.7A Active CN107194051B (zh) | 2017-05-11 | 2017-05-11 | 一种后桥系统非线性强度分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107194051B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112232002B (zh) * | 2020-12-09 | 2021-03-16 | 北京智芯仿真科技有限公司 | 一种基于误差估计的集成电路的电磁响应确定方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101303280A (zh) * | 2008-04-11 | 2008-11-12 | 北京联合大学 | 采用曲线平移法处理材料载荷变形曲线的方法 |
CN102214256A (zh) * | 2011-05-20 | 2011-10-12 | 中国汽车技术研究中心 | 一种汽车碰撞波形特征参数提取与梯形波构建方法 |
CN102607502A (zh) * | 2012-01-18 | 2012-07-25 | 中北大学 | 汽车后桥总成尺寸自动检测装置及检测方法 |
CN103603262A (zh) * | 2013-12-09 | 2014-02-26 | 中铁二十三局集团有限公司 | 斜高墩大跨径连续梁的线形控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7908112B2 (en) * | 2008-11-24 | 2011-03-15 | GM Global Technology Operations LLC | Dynamic observer for the estimation of vehicle lateral velocity |
-
2017
- 2017-05-11 CN CN201710330088.7A patent/CN107194051B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101303280A (zh) * | 2008-04-11 | 2008-11-12 | 北京联合大学 | 采用曲线平移法处理材料载荷变形曲线的方法 |
CN102214256A (zh) * | 2011-05-20 | 2011-10-12 | 中国汽车技术研究中心 | 一种汽车碰撞波形特征参数提取与梯形波构建方法 |
CN102607502A (zh) * | 2012-01-18 | 2012-07-25 | 中北大学 | 汽车后桥总成尺寸自动检测装置及检测方法 |
CN103603262A (zh) * | 2013-12-09 | 2014-02-26 | 中铁二十三局集团有限公司 | 斜高墩大跨径连续梁的线形控制方法 |
Non-Patent Citations (1)
Title |
---|
"矿用自卸车后桥强度有限元分析";王翠凤;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20051015(第6期);第C035-122页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107194051A (zh) | 2017-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10867083B2 (en) | Technique for generating approximate design solutions | |
US11468292B2 (en) | Techniques for analyzing vehicle design deviations using deep learning with neural networks | |
CN109326002B (zh) | 基于钻孔数据的矿体建模方法、装置、系统及存储介质 | |
CN107369140B (zh) | 非结构化环境下的高精密靶球中心提取方法 | |
US11030360B2 (en) | Aircraft designing method, aircraft designing program and aircraft designing apparatus | |
CN110781558B (zh) | 基于疲劳和侧倾性能的汽车稳定杆多学科优化设计方法 | |
CN112597610B (zh) | 机械臂结构轻量化设计的优化方法、装置及设备 | |
CN107194051B (zh) | 一种后桥系统非线性强度分析方法 | |
US20240046086A1 (en) | Quantization method and quantization apparatus for weight of neural network, and storage medium | |
CN105654460A (zh) | 点云拼接系统及方法 | |
CN104460429A (zh) | Can信号的解析方法和装置 | |
CN113987684A (zh) | 基于sn曲线的随机应力载荷谱分级方法、系统和存储介质 | |
CN110398942B (zh) | 一种用于工业生产过程控制的参数辨识方法 | |
CN113761684A (zh) | 一种结构件加工变形分段校正载荷优化方法及系统 | |
US8738678B2 (en) | Methods and systems for determining an enhanced rank order value of a data set | |
Khan et al. | Geometric correction of deformed chromosomes for automatic Karyotyping | |
CN112733471A (zh) | 用于分离两体非定常气动特性的方法 | |
CN109902417B (zh) | 基于可视化界面的火箭发动机法兰快速设计系统 | |
CN102928830A (zh) | 激光高度计无效数据剔除方法 | |
CN106001933A (zh) | 激光切割修边线的优化方法 | |
CN113032894B (zh) | 一种基于冯卡门形线的双锥整流罩形线优化设计方法 | |
US20180335645A1 (en) | Ophthalmic lens designing device and ophthalmic lens designing method | |
CN109214056A (zh) | 一种基于流动物理的气动优化设计变量选取方法 | |
CN105138807A (zh) | 均匀接触应力液压联轴节及其优化方法 | |
KR101910736B1 (ko) | 직접모사법 기반의 정보보존 방법에서의 이완함수를 이용한 입자 및 격자의 정보갱신단계의 보정방법 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |