CN106096136B - 一种车辆安全性设计的聚类-多椭球不确定性度量方法 - Google Patents

一种车辆安全性设计的聚类-多椭球不确定性度量方法 Download PDF

Info

Publication number
CN106096136B
CN106096136B CN201610412294.8A CN201610412294A CN106096136B CN 106096136 B CN106096136 B CN 106096136B CN 201610412294 A CN201610412294 A CN 201610412294A CN 106096136 B CN106096136 B CN 106096136B
Authority
CN
China
Prior art keywords
cluster
ellipsoid
sample
vehicle safety
response
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
Application number
CN201610412294.8A
Other languages
English (en)
Other versions
CN106096136A (zh
Inventor
刘杰
刘浩
姜潮
韩旭
刘启明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hunan University
Original Assignee
Hunan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hunan University filed Critical Hunan University
Priority to CN201610412294.8A priority Critical patent/CN106096136B/zh
Publication of CN106096136A publication Critical patent/CN106096136A/zh
Application granted granted Critical
Publication of CN106096136B publication Critical patent/CN106096136B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Body Structure For Vehicles (AREA)

Abstract

本发明涉及一种车辆安全性设计的聚类‑多椭球不确定性度量方法，包括以下步骤：测量车身部分结构的尺寸参数样本；以混合高斯模型为基础，对车身部分结构的尺寸参数样本进行聚类分析；对每个聚集进行椭球不确定性建模；对不确定性参数进行空间转换；对车辆安全性吸能结构的响应面方程在每个聚类单位球空间中获取响应极值点；将聚类单位球空间中的极值点转换到聚类多椭球空间计算得到车辆安全性吸能结构响应的准确边界。所述方法仅在小样本下就能建立合理紧凑的车辆结构参数不确定模型，且该聚类‑多椭球模型度量了参数不确定性的多分布及相关性等特征，为不确定性环境下的车辆安全性设计提供了重要基础模型和依据。

Description

一种车辆安全性设计的聚类-多椭球不确定性度量方法
技术领域
本发明涉及结构不确定性分析技术领域，尤其涉及一种车辆安全性设计的聚类-多椭球不确定性度量方法。
背景技术
车身结构的安全性设计是整车设计的关键部分之一，但由于结构的几何特征、材料属性、加工测量误差等常伴有不确定性，这必然会使得车身结构实验碰撞响应结果不会是一个确定的值，而是一个范围，这给身结构安全性设计带来很大的困难及浪费。因此正确有效地估计和量化结构不确定性参数对结构响应的影响，对于工程结构的分析、设计等具有重要指导意义。然而，由于车身结构实验昂贵只能获取不确定性参数的小样本信息，且样本具有多分布形貌。虽然单椭球凸模型可度量参数的不确定性，但却忽略了样本分布可能存在的多聚类属性，致使可能会扩大参数的不确定性域，最终导致分析结果过保守。
因此，发展一种合理的不确定性度量方法对小样本多聚类情况下的参数不确定性进行度量显得非常必要。
发明内容
鉴于上述的分析，本发明旨在提供一种车辆安全性设计的聚类-多椭球不确定性度量方法，用以解决现有技术针对车身结构设计过程中，结构参数伴有不确定性的问题。
本发明的目的主要是通过以下技术方案实现的：
一种车辆安全性设计的聚类-多椭球不确定性度量方法，其特征在于，包括以下步骤：
步骤1：测量车身部分结构的尺寸参数样本；进行车辆碰撞实验；测量相应部位碰撞后的形变量参数样本；
步骤2：建立所述车辆碰撞系统的有限元模型，对车辆碰撞前后获得的参数样本进行仿真，获得车辆碰撞安全性吸能结构的响应面方程；
步骤3：以混合高斯模型为基础，对步骤1获得的车身部分结构的尺寸参数样本进行聚类分析，将具有相似属性的样本点进行聚集；
步骤4：对每个聚集进行椭球不确定性建模；
步骤5：对不确定性参数进行空间转换，将原先的聚类多椭球空间转换成聚类单位球空间；
步骤6：对步骤2所构建的车辆安全性吸能结构的响应面方程在每个聚类单位球空间中获取响应极值点；
步骤7：将聚类单位球空间中的极值点转换到聚类多椭球空间计算得到车辆安全性吸能结构响应的准确边界，对车身结构尺寸进行再设计。
步骤1中所述尺寸参数样本进一步包括：
观测m个样本X(m)＝(x1,x2,…,xp)满足混合高斯分布，该分布由K个(K≥1)成分构成，每个成分代表一个不同的类别，K个成分即K个类别、也即K个独立的高斯分布，其中每个X(m)是n维的(n≥1)。
所述步骤3进一步包括：
构造样本期望，利用期望最大化算法，根据小样本集合X和聚类数目K，通过迭代求解每一聚类集合的参数，即确定每一个聚类集合的概率，均值和协方差矩阵，实现对不确定性参数的类别划分。
所述步骤4进一步包括：
利用高斯模型的等高椭球特性对每个聚集进行椭球不确定性建模，获得不确定性参数的聚类多椭球模型。
所述步骤5进一步包括：
当进行空间转换后，原聚类-多椭球空间的每个聚集被转换成对应的单位球空间
所述步骤6进一步包括：
构造如下优化问题获取响应极值，进行迭代更新搜索方向直到迭代收敛，获取响应极值点。
所述步骤3进一步包括：
利用不确定性参数样本的联合概率函数构造样本期望。
本发明有益效果如下：
(1)在少样本量的情况下，对车身结构不确定性参数进行聚类-多椭球建模，即可获得合理紧凑的参数不确定性域。
(2)聚类-多椭球模型描述了不确定性参数的多分布形貌，能够获得合理的不确定性参数域，使得传播分析更为准确。
(3)根据车辆安全性结构吸能响应的准确边界值，在不确定性参数椭球域内选取结构参数值以对车身结构进行再设计，提高了车身结构的安全性。
本发明的其他特征和优点将在随后的说明书中阐述，并且，部分的从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
附图说明
附图仅用于示出具体实施例的目的，而并不认为是对本发明的限制，在整个附图中，相同的参考符号表示相同的部件。
图1为基于聚类-多椭球模型的车驾驶室不确定性分析流程图；
图2为驾驶室结构有限元模型示意图；
图3为驾驶室不确定性参数的聚类-多椭球模型示意图；
图4为不确定性参数空间转换图。
具体实施方式
下面结合附图来具体描述本发明的优选实施例，其中，附图构成本申请一部分，并与本发明的实施例一起用于阐释本发明的原理。
根据本发明的一个具体实施例，公开了一种车辆安全性设计的聚类-多椭球不确定性度量方法，具体包括以下步骤：
步骤1：按照乘用车正面碰撞的实验方法和程序进行车辆碰撞实验前，测量车身部分结构的尺寸参数样本，包括防撞梁厚度、加强板厚度、纵梁支柱板厚度、前支柱外板厚度；进行车辆碰撞实验，测量相应部位碰撞后的形变量参数样本；其中，碰撞实验前后所获得的参数样本为以获得具有不确定性的参数样本，且样本具有多分布形貌。
步骤2：按照车辆碰撞实验工况，建立所述车辆碰撞系统的有限元模型，利用车辆碰撞的有限元模型对车辆碰撞前后获得的参数样本进行仿真以获得车辆碰撞安全性吸能结构的响应面方程；
步骤3：以混合高斯模型为基础，采用基于EM(期望最大化)算法对步骤1获得的车身部分结构的尺寸参数样本进行聚类分析，将具有相似属性的样本点进行聚集。
具体地，观测m个样本X(m)＝(x1,x2,…,xp)满足混合高斯分布，该分布由K个(K≥1)成分构成，每个成分代表一个不同的类别，K个成分即K个类别、也即K个独立的高斯分布，其中每个X(m)是n维的(n≥1)。
以车辆驾驶室结构支柱板厚度为例，车辆驾驶室结构不确定性参数样本X＝(x1,x2)，如表1所示，其中x1为前支柱外板厚度，x2为纵梁支柱板厚度。
表1车驾驶室结构支柱板的厚度样本(mm)
X中样本数据的混合概率密度函数为
上式表示样本X是第k类的密度函数，其中，αk表示观察值属于第k类的概率，即观察值属于某一类别的可能性，这样0＜αk＜1、Φ(X|μk,∑k)表示X是第k类的高斯分布概率密度函数，μk、∑k表示第k类高斯分布的参数，μk是第k类分布的期望，∑k是第k类分布的协方差矩阵；K代表聚类的数目。
获得样本集X的概率公式，也即X的样本数据的联合概率函数为
式中，表示待估计的混合分布的参数，m是不确定性参数的样本量，即样本集合X中的数据个数，如表1中m＝30。
引入条件概率并对式(1)取对数，得到样本期望
式(2)中表示似然对数函数，βk(X(r))表示条件概率，即第r个样本点出现的概率值，X(r),r＝1,2,…,m表示m个不确定性参数样本量，其中，X(r)表示第r个样本点。
构造样本期望，利用期望最大化算法，根据小样本集合X和聚类数目K，通过迭代求解每一聚类集合的参数，即确定每一个聚类集合的概率，均值和协方差矩阵，以此实现对不确定性参数的类别划分。
利用期望最大化算法对不确定性参数进行划分，即求期望的最大化，即对式(2)分别对概率，均值和协方差矩阵求偏导，得到
ε是迭代终止条件，大小涉及到计算的精度，越小精度越高，一般取0.00001，指的是三个参数的最小二乘误差小于给定的误差限ε，也就是两次相邻的迭代结果的误差小于这个误差限的话，此时的结果作为最优的参数结果。
步骤4：对每个聚集进行椭球不确定性建模；
通过步骤3的聚类分析，可实现不确定参数的划分，最终利用高斯模型的等高椭球特性对每个聚集进行椭球不确定性建模，可获得不确定性参数的聚类多椭球模型为
针对本发明具体实例分析，车驾驶室结构不确定性参数样本如表1所示。车驾驶室有限元模型如图2所示，车驾驶室结构不确定性参数的聚类-多椭球模型为
式中x1和x2分别为纵梁支柱板厚度和前支柱外板厚度。对应的聚类-多椭球模型如图3所示。
步骤5：对不确定性参数进行空间转换，将原先的聚类多椭球空间转换成聚类单位球空间；
当进行空间转换后，原聚类-多椭球空间的每个聚集将被转换成对应的单位球空间，如图4所示。
式中δ表示单位球空间中的不确定性参数，Eδ表示单位球空间。
步骤6：对步骤2所构建的车辆安全性吸能结构的响应面方程在每个聚类单位球空间中获取响应极值点。
由椭球极值理论可知，结构响应的极值点出现不确定性椭球外壳上。显然，在单位球空间中极值点依然出现在求外壳上，且在极值点附近，结构吸能响应梯度矢量方向和极值点的矢量方向共线。因此可构造如下优化问题获取响应极值
s.t.δTδ＝1
式中▽g(δ)表示响应梯度，γ两矢量方向的夹角。按照如下的迭代更新搜索方向
直到迭代收敛便可获取响应极值点。
步骤7：将单位球空间中的极值点转换到原空间计算得到车辆安全性吸能结构响应的准确边界，并由此对车身结构尺寸进行再设计过程。
将极值点转换到原空间求解得到车驾驶室结构吸能响应准确上下界。
针对本发明具体实例，车驾驶室结构吸能响应准确上下界及极值点如表2所示
表2车驾驶室结构吸能响应极值点及边界值
综上所述，本发明实施例提供了一种车辆安全性设计的聚类-多椭球不确定性度量方法，有益效果如下：
(1)在少样本量的情况下，对车身结构不确定性参数进行聚类-多椭球建模，即可获得合理紧凑的参数不确定性域。
(2)聚类-多椭球模型描述了不确定性参数的多分布形貌，能够获得合理的不确定性参数域，使得传播分析更为准确。
(3)根据车辆安全性结构吸能响应的准确边界值，在不确定性参数椭球域内选取结构参数值以对车身结构进行再设计，提高了车身结构的安全性。
本领域技术人员可以理解，实现上述实施例方法的全部或部分流程，可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于计算机可读存储介质中。其中，所述计算机可读存储介质为磁盘、光盘、只读存储记忆体或随机存储记忆体等。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。

Claims (6)

1.一种车辆安全性设计的聚类-多椭球不确定性度量方法，其特征在于，包括以下步骤：
步骤1：测量车身部分结构的尺寸参数样本；进行车辆碰撞实验；测量相应部位碰撞后的形变量参数样本；
步骤2：建立所述车辆碰撞系统的有限元模型，对车辆碰撞前后获得的参数样本进行仿真，获得车辆碰撞安全性吸能结构的响应面方程；
步骤3：以混合高斯模型为基础，对步骤1获得的车身部分结构的尺寸参数样本进行聚类分析，将具有相似属性的样本点进行聚集；
步骤4：对每个聚集进行椭球不确定性建模：利用高斯模型的等高椭球特性对每个聚集进行椭球不确定性建模，获得不确定性参数的聚类多椭球模型；
步骤5：对不确定性参数进行空间转换，将原先的聚类多椭球空间转换成聚类单位球空间；
步骤6：对步骤2所构建的车辆安全性吸能结构的响应面方程在每个聚类单位球空间中获取响应极值点；
步骤7：将聚类单位球空间中的极值点转换到聚类多椭球空间计算得到车辆安全性吸能结构响应的准确边界，对车身结构尺寸进行再设计。
2.根据权利要求1所述的车辆安全性设计的聚类-多椭球不确定性度量方法，其特征在于，步骤1中所述尺寸参数样本进一步包括：
观测m个样本X(m)＝(x1,x2,…,xp)满足混合高斯分布，该分布由K个(K≥1)成分构成，每个成分代表一个不同的类别，K个成分即K个类别、也即K个独立的高斯分布，其中每个X(m)是p维的(p≥1)。
3.根据权利要求1所述的车辆安全性设计的聚类-多椭球不确定性度量方法，其特征在于，所述步骤3进一步包括：
构造样本期望，利用期望最大化算法，根据小样本集合X和聚类数目K，通过迭代求解每一聚类集合的参数，即确定每一个聚类集合的概率，均值和协方差矩阵，实现对不确定性参数的类别划分。
4.根据权利要求1所述的车辆安全性设计的聚类-多椭球不确定性度量方法，其特征在于，所述步骤5进一步包括：
当进行空间转换后，原聚类-多椭球空间的每个聚集被转换成对应的单位球空间。
5.根据权利要求1所述的车辆安全性设计的聚类-多椭球不确定性度量方法，其特征在于，所述步骤6进一步包括：
构造如下优化问题获取响应极值，
其中，表示响应梯度，γ两矢量方向的夹角，δ表示单位球空间中的不确定性参数，这里不确定性参数不泛指，而是计算过程中的不确定性参数；
按照如下的迭代更新搜索方向对搜索方向进行迭代更新直至迭代收敛，获取响应极值点。
6.根据权利要求2所述的车辆安全性设计的聚类-多椭球不确定性度量方法，其特征在于，所述步骤3进一步包括：
利用不确定性参数样本的联合概率函数构造样本期望。
CN201610412294.8A 2016-06-13 2016-06-13 一种车辆安全性设计的聚类-多椭球不确定性度量方法 Active CN106096136B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610412294.8A CN106096136B (zh) 2016-06-13 2016-06-13 一种车辆安全性设计的聚类-多椭球不确定性度量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610412294.8A CN106096136B (zh) 2016-06-13 2016-06-13 一种车辆安全性设计的聚类-多椭球不确定性度量方法

Publications (2)

Publication Number Publication Date
CN106096136A CN106096136A (zh) 2016-11-09
CN106096136B true CN106096136B (zh) 2019-03-26

Family

ID=57846283

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610412294.8A Active CN106096136B (zh) 2016-06-13 2016-06-13 一种车辆安全性设计的聚类-多椭球不确定性度量方法

Country Status (1)

Country Link
CN (1) CN106096136B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110895639B (zh) * 2019-11-27 2024-03-01 河北工业大学 一种基于高斯多椭球模型的机器人系统可靠性分析方法
CN113188444B (zh) * 2021-05-06 2023-01-13 上海航天测控通信研究所 激光测量系统不确定度测试实验与计算方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1758255A (zh) * 2005-11-03 2006-04-12 上海交通大学 基于响应面法的轿车车身零件轻量化方法
WO2010069614A1 (en) * 2008-12-17 2010-06-24 Telefonaktiebolaget Lm Ericsson (Publ) Reporting of positioning data
CN102495923A (zh) * 2011-11-23 2012-06-13 湖南大学 一种汽车碰撞安全的混合可靠性评估方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1758255A (zh) * 2005-11-03 2006-04-12 上海交通大学 基于响应面法的轿车车身零件轻量化方法
WO2010069614A1 (en) * 2008-12-17 2010-06-24 Telefonaktiebolaget Lm Ericsson (Publ) Reporting of positioning data
CN102495923A (zh) * 2011-11-23 2012-06-13 湖南大学 一种汽车碰撞安全的混合可靠性评估方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
基于凸模型和伪概率分布的不确定性结构响应分析;刘浩等;《固体力学学报》;20151230;第537-543页
基于凸模型的结构非概率可靠性优化;亢战;《力学学报》;20061130;第808-815页
基于响应面法的汽车吸能部件优化问题研究;王琥等;《系统仿真学报》;20070820;第3824-3829页
基于概率—凸集混合模型的汽车碰撞安全可靠性分析;张乾;《机械强度》;20121215;第916-920页
车辆正碰安全性的混合可靠度分析;李文学;《计算机仿真》;20120215;第358-361页

Also Published As

Publication number Publication date
CN106096136A (zh) 2016-11-09

Similar Documents

Publication Publication Date Title
CN103336869B (zh) 一种基于高斯过程联立mimo模型的多目标优化方法
Grigoriu et al. Spherical harmonic-based random fields for aggregates used in concrete
Feather et al. A numerical study into element type and mesh resolution for crystal plasticity finite element modeling of explicit grain structures
CN110458174A (zh) 一种无序点云关键特征点精确提取方法
Nie et al. A three-dimensional particle roundness evaluation method
CN101241521A (zh) 基于支持向量机的轿车车身综合性能指标建模方法
CN106096136B (zh) 一种车辆安全性设计的聚类-多椭球不确定性度量方法
Chung et al. Reconstruction of random two-phase polycrystalline solids using low-order probability functions and evaluation of mechanical behavior
CN105895089A (zh) 一种语音识别方法及装置
CN106644868A (zh) 一种二维非凸形随机骨料周围界面浓度的测定方法
CN104795063A (zh) 一种基于声学空间非线性流形结构的声学模型构建方法
Kleinert et al. Modeling of large scale granular systems using the discrete element method and the non-smooth contact dynamics method: A comparison
Iriawan et al. On the modeling of the average value of high school national examination in West Java using Bayesian hierarchical mixture normal approach
Jiang et al. Efficient and unbiased safety test for autonomous driving systems
CN110334450A (zh) 一种多块结构网格生成中物面投影错误的修复方法
EP3754534A1 (en) A fast method for computer-based simulation
Dietrich et al. Fast and flexible uncertainty quantification through a data-driven surrogate model
Wei et al. Simulation of rolling noise based on the mixed Lagrangian–Eulerian method
Xi et al. Validation metric for dynamic system responses under uncertainty
CN106682347B (zh) 基于岩体结构面随机强度的非连续变形分析方法
Fang et al. Multiobjective crashworthiness optimization with radial basis functions
Allen et al. Comparison of uncertainty propagation/response surface techniques for two aeroelastic systems
Wang A robust design optimization method considering correlated intervals
Choi et al. Confidence-Based Uncertainty Quantification and Reliability Assessment Using Limited Numbers of Input and Output Test Data and Validated Simulation Model
Cho Natural element hierarchic models for laminated composite plates

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant