CN108563905B - 汽车b柱加强板碳纤维增强复合材料优化设计方法 - Google Patents
汽车b柱加强板碳纤维增强复合材料优化设计方法 Download PDFInfo
- Publication number
- CN108563905B CN108563905B CN201810409682.XA CN201810409682A CN108563905B CN 108563905 B CN108563905 B CN 108563905B CN 201810409682 A CN201810409682 A CN 201810409682A CN 108563905 B CN108563905 B CN 108563905B
- Authority
- CN
- China
- Prior art keywords
- optimization
- reinforcing plate
- composite material
- column
- column reinforcing
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Body Structure For Vehicles (AREA)
Abstract
本发明公开了一种汽车B柱加强板碳纤维增强复合材料优化设计方法,为克服现有技术存在有限元整车碰撞模型计算时间过长及解决碳纤维增强复合材料铺层厚度、铺层角度与铺层顺序设计的问题,所述的优化设计方法的步骤为:1.从整车有限元模型中解耦B柱子结构;2.碳纤维增强复合材料B柱加强板优化:1)B柱加强板静态工况建立;2)碳纤维增强复合材料B柱加强板静态工况铺层优化:(1)拓扑优化;(2)尺寸优化;(3)顺序优化;3)B柱加强板优化结果于解耦子结构动态工况下验证计算:(1)改进B柱加强板;(2)对带有优化后碳纤维增强复合材料B柱加强板的模型进行模拟仿真与分析;(3)查看优化后B柱与原车型中的B柱的对比数据。
Description
技术领域
本发明涉及汽车被动安全性研究领域中的一种设计方法,更确切地说,本发明涉及一种汽车B柱子结构解耦以及汽车B柱加强板碳纤维增强复合材料优化设计方法。
背景技术
随着能源的紧缺以及新能源汽车的兴起,轻量化已经成为当前汽车行业中不可或缺的研究内容。现今主要乘用车轻量化设计方式包含:采用拓扑等方式对于车身部件进行结构设计、高强度钢的大量应用、轻质金属以及高分子材料、陶瓷材料和复合材料等非金属材料的使用。其中,碳纤维增强复合材料凭借其密度低、比强度高和工艺性好等优点,已经在复合材料车身结构中有初步使用,如BMW i7 B柱在金属结构的关键部位局部附加碳纤维增强复合材料,构成复合结构的做法,既能充分发挥碳纤维增强复合材料的优势,同时可降低成本,减少用量,是目前自主品牌汽车中对于碳纤维增强复合材料的主要应用发展方向。为适应当前社会发展中所需的更高的抗撞性和轻量化要求,并且随着碳纤维增强复合材料成型技术的不断发展,碳纤维增强复合材料将会在汽车车身结构中有更为广泛的应用。
然而当前于车身轻量化研究中大多只是对传统金属材料进行简单替换,对抗冲击性能不满足部位通常采用增加板厚等方式以满足要求。但这种方式没有充分利用复合材料可设计性能,并且轻量化效果不明显,有必要对碳纤维增强复合材料进行材料-结构一体化设计。
B柱结构为汽车车身在侧面碰撞情况中的主要吸能和承力部件,也是垂直方向构成传力路径的部件,它可以将力向上传递给车顶盖纵梁,力再通过车顶支撑横梁进行横向传递。其抗撞性能是在侧面碰撞工况下决定对人体造成伤害大小的最为关键的因素,B柱所受到的变形情况将会对整个车辆的侧面方向结构都会产生一定量的影响。侧面碰撞时人体的伤害主要是由于头、胸、脊椎以及骨盆部位的冲击所造成的不可逆转的伤害。若其不具有足够的刚度,则将会造成较为严重的伤亡情况,综上所述,在侧面碰撞情况中,对于B柱结构的研究,具有很重要的意义。碳纤维增强复合材料B柱加强板优化目的在于充分利用纤维增强复合材料的比强度好、比模量高、耐腐蚀、抗疲劳等优点,其铺层厚度、角度以及顺序对于材料性能都具有至关重要的关系,现将拓扑优化、尺寸优化以及顺序优化应用到碳纤维增强复合材料B柱加强板优化中以达到汽车轻量化效果。
另一方面,目前由于计算机性能的限制,整车碰撞模型计算时间过长,影响效率;而往往当前所存在的子结构解耦方法所实现的子结构模型与整车中的侵入量以及侵入速度曲线拟合程度低。本文所提出的解耦方法可以较为准确地将整车中的复杂碰撞情况进行简化,大大节约了计算时间,提高计算效率。
通过进行国内外相关文献检索,未发现有类似的解耦及碳纤维增强复合材料B柱加强板参数优化分析方法。
发明内容
本发明所要解决的技术问题是克服了现有技术存在的有限元整车碰撞模型计算时间过长及解决碳纤维增强复合材料铺层厚度、铺层角度和铺层顺序设计的问题,提供了一种汽车B柱加强板碳纤维增强复合材料优化设计方法。
为解决上述技术问题,本发明是采用如下技术方案实现的:所述的汽车B柱加强板碳纤维增强复合材料优化设计方法的步骤如下:
1)从整车有限元模型中解耦B柱子结构:
(1)根据原有整车侧面碰撞模型,保留与B柱相连接的部件来保证其变形情况的准确性,即门槛梁、侧围板、顶盖横梁、车门管状防撞梁,再加上B柱结构作为子结构模型;
(2)采用与B柱相连接部件中与B柱相连接位置接近的区域以及侧面碰撞台车缓冲蜂窝铝最前端表面作为碰撞参考点,进行侧面碰撞有限元计算,通过Ls-prepost后处理软件提取碰撞参考点的位移-时间曲线数据;
(3)将提取出的曲线数据通过MATLAB编程软件转化为Ls-dyna仿真计算软件可识别的数据信息,再通过文本编辑器导入简化后只保留B柱和其周围部件以及侧面碰撞台车前端的子工况模型的计算k文件,便准确模拟整车侧面碰撞情况;
2)碳纤维增强复合材料B柱加强板优化:
(1)B柱加强板静态工况建立;
(2)碳纤维增强复合材料B柱加强板静态工况铺层优化:
A.拓扑优化;
B.尺寸优化;
C.顺序优化;
(3)B柱加强板优化结果于解耦子结构动态工况下验证计算。
按照技术方案中所述的B柱加强板静态工况建立是指:
(1)结构简化
根据已有B柱加强板结构形状的几何尺寸,将B柱加强板厚度设定成统一初始厚度,并且保留整车模型中原金属B柱加强板中的形状结构参数,将简化结构中预留出螺栓固定孔结构,并进行其网格质量检查;
(2)碳纤维增强复合材料B柱加强板的建立
a.将即将优化的部件B柱加强板更改为碳纤维增强复合材料,根据经验进行铺层角度、铺层顺序、铺层厚度设定,使其于碰撞仿真中可以保持原抗撞性能,将此铺层信息作为接下来优化的初始参数;
b.由于B柱加强板皆为设计区域,故采用8号各向异性材料模型创建铺层形成层合板结构,属性定义针对设计区域设置选取PCOMPP,由于Optistruct优化软件只能对材料进行减薄而不具备增厚的功能,故需建立所谓的初始厚度较厚的超级层,该超级层的厚度通过初始优化参数进行设置;
(3)加载的设置
提取B柱模型的抗冲击性仿真中的B柱加强板最大相对位移来计算所需等效力的数值,将其加载在与碰撞台车相接触的B柱加强板结构优化模型中,并按照加载点的个数,将等效力进行平均分配加载;根据x、y、z轴的力的数量级进行合理忽略简化,所得的静态仿真模型可以与动态仿真模型获得较好的位移拟合结果,之后在此基础上进行静态优化仿真模拟。
按照技术方案中所述的拓扑优化是指:
a.定义设计变量,对于拓扑优化首先创建变量,定义质量和参考点位移两个响应,二者分别作为优化响应和约束响应;创建约束,将参考点位移响应作为约束,设置约束上限值,将质量最小定义为优化目标;
b.调整优化控制参数,包括最大允许迭代次数、目标容差值以及优化步长;
c.定义控制卡片,在此处需要在控制卡片部分定义OUTPUT(输出)卡片关键字为FSTOSZ,然后保存或输出文件;
d.最后进行检查,如果没有错误信息则可开始进行优化,优化结束后查看结果文件,在HyperView中打开优化结果文件检查单元厚度。
按照技术方案中所述的尺寸优化是指:
a.为保证铺层信息的正确传递,单个方向的厚度信息以铺层形状的形式进行输出,首先重新设置控制卡片,删除“OMIT”卡片,将一种类型关键字FSTOSZ改为另一种类型关键字SZTOSH;
b.编辑层合板Laminate参数,将层合板形式设置为Total,分别定义每层尺寸变量,确保修改每层厚度上限,将不能修改的规整为相近整数,依旧将参考点位移值作为约束条件,优化目标仍然是质量最小;
c.指定可加工的厚度的值,在优化设置后进行检查优化,并查看铺层厚度和数目。
按照技术方案中所述的顺序优化是指:对于复合材料来说,铺层顺序与其层合板的各向异性力学性能以及工艺可行性都有着紧密的联系,复合材料铺层顺序的变化会影响其结构的整体性能,在计算最优的层叠次序之前,首先对铺层顺序变量设置制造约束,定义0度、90度、45度和-45度四个方向铺层最多可以有2层相同铺层连续出现;+45度和-45度成对出现来保证对称性,最大限度避免层间应力的产生,优化约束和优化目标不变,最后进行铺层顺序优化;检查输出的结果文件,在浏览器中查看各个迭代步的铺层顺序。
按照技术方案中所述的B柱加强板优化结果于解耦子结构动态工况下验证计算是指:
(1)根据铺层结果进行规整分区来改进B柱加强板;
(2)对带有优化后碳纤维增强复合材料B柱加强板的模型进行解耦子工况侧撞台车碰撞的模拟仿真计算;
(3)查看优化后B柱与原车型中的B柱的对比数据,即通过其B柱加强板上对应人体胸部的参考点的侵入量和侵入速度最大值的对比来验证其性能的改善程度,同时获得了其轻量化效果。
与现有技术相比本发明的有益效果是:
1.参阅图3,本发明所述的汽车B柱加强板碳纤维增强复合材料优化设计方法有效地对于所提取出的整车相关碰撞参考点的位移时间曲线进行整合,数据处理后得到简化后的解耦B柱子结构碰撞模型,与现有的子工况解耦方式相比,该模型能够更为准确地与整车碰撞侵入速度及侵入量曲线拟合。
2.利用本发明所述的汽车B柱加强板碳纤维增强复合材料优化设计方法在车身抗撞性概念设计阶段可以准确并且快速地计算出所设计的碳纤维增强复合材料B柱的抗撞性能,通过对新型材料碳纤维铺层角度、厚度以及铺层顺序的优化得到最优B柱参数方案,该方法有效减少了计算机仿真时间,降低了开发成本,大大缩短了设计周期。
3.本发明所述的汽车B柱加强板碳纤维增强复合材料优化设计方法将拓扑优化、尺寸优化以及顺序优化依次应用到碳纤维增强复合材料B柱加强板的轻量化设计中。所提出的设计思路可应用于车用复合材料的材料-结构一体化设计中,实现对复合材料独有的材料设计,充分利用复合材料与金属材料相比可设计性的优势,快速高效设计出满足汽车性能要求的轻量化部件。
附图说明
下面结合附图对本发明作进一步的说明:
图1为本发明所述的汽车B柱加强板碳纤维增强复合材料优化设计方法的流程框图;
图2为本发明所述的汽车B柱加强板碳纤维增强复合材料优化设计方法的解耦B柱子结构模型;
图3-1为本发明所述的汽车B柱加强板碳纤维增强复合材料优化设计方法所获得的B柱子结构模型相对于整车模型侵入速度拟合图;
图3-2为本发明所述的汽车B柱加强板碳纤维增强复合材料优化设计方法所获得的B柱子结构模型相对于整车模型侵入量拟合图;
图4为本发明所述的汽车B柱加强板碳纤维增强复合材料优化设计方法中B柱加强板形状结构参数模型图;
图5为本发明所述的汽车B柱加强板碳纤维增强复合材料优化设计方法中采用施加等效静载力的方法所得B柱加强板加载及约束示意图;
图6为本发明所述的汽车B柱加强板碳纤维增强复合材料优化设计方法中B柱加强板的拓扑优化后单元厚度需要减薄处的示意图;
图7为本发明所述的汽车B柱加强板碳纤维增强复合材料优化设计方法中B柱加强板优化结果的分区图;
具体实施方式
下面结合附图对本发明作详细的描述:
本发明所述的汽车B柱加强板碳纤维增强复合材料优化设计方法分为两个阶段,首先将B柱子结构从整车有限元侧面碰撞模型中解耦,在实现与整车模型中位于人体胸部位置碰撞参考点的位移-时间曲线拟合良好的情况下进行后续计算,大大减少计算机运行时间,提高计算效率;然后在解耦的子结构模型基础上,根据碳纤维增强复合材料试验参数建模轻量化B柱加强板结构,获得优化初始数据,即初始纤维板铺层和厚度;采用施加等效静载力的方法进行碳纤维B柱加强板的铺层角度、厚度以及铺层顺序的优化;再根据优化结果对B柱加强板进行改进,并于最后将其代入动态仿真模型中计算来验证其轻量化效果及碰撞性能。
汽车B柱加强板碳纤维增强复合材料优化设计方法的步骤如下:
1.从整车有限元模型中解耦B柱子结构
1)参阅图2,根据原有整车侧面碰撞模型,保留与B柱相连接的部件来保证其变形情况的准确性,即门槛梁、侧围板、顶盖横梁、车门管状防撞梁,再加上B柱结构作为子结构模型。
2)采用与B柱相连接部件中与B柱相连接位置接近的区域以及侧面碰撞台车缓冲蜂窝铝最前端表面作为碰撞参考点,进行侧面碰撞有限元计算,通过Ls-prepost后处理软件提取碰撞参考点的位移-时间曲线数据。
3)将提取出的曲线数据通过MATLAB编程软件转化为Ls-dyna仿真计算软件可识别的数据信息,再通过文本编辑器导入简化后只保留B柱和其周围部件以及侧面碰撞台车前端的子工况模型的计算k文件,便可准确模拟整车侧面碰撞情况,并且大大减少了计算时间,提高工作效率。
2.碳纤维增强复合材料B柱加强板优化
本发明首先于前步骤中所得的B柱子结构碰撞模型中将金属B柱加强板更改为碳纤维增强复合材料,仿真参数通过复合材料力学性能试验所得,包括拉伸试验及剪切试验。通过有限元碰撞仿真计算保证所建立的碳纤维增强复合材料B柱加强板与原金属B柱加强板性能相同,采取的铺层顺序、铺层角度皆根据经验设定,此铺层设定为之后优化的初始参数。简化乘用车B柱加强板几何模型,采用施加等效静载力的方法,根据B柱模型的抗冲击性仿真中的B柱加强板最大相对位移来计算所需等效力的数值,将其加载在与碰撞台车相接触的B柱加强板结构优化模型处,并按照加载点的个数,将接触力进行平均分配加载。利用Optistruct优化软件对B柱加强板依次进行拓扑、尺寸及顺序优化来获得合理的铺层材料分布、铺层厚度、角度及顺序;最后对静态铺层优化结果进行分区解读设置后带入到Ls-Dyna仿真计算软件验证动态工况下是否满足要求。
所述的碳纤维增强复合材料B柱加强板优化步骤如下:
1)B柱加强板静态工况建立
(1)结构简化
根据已有B柱加强板结构形状的几何尺寸,将B柱加强板厚度设定成统一初始厚度,并且保留整车模型中原金属B柱加强板中的形状结构参数,将简化结构中预留出螺栓固定孔结构,并进行其网格质量检查;
(2)碳纤维增强复合材料B柱加强板的建立
a.首先将即将优化的部件B柱加强板更改为碳纤维增强复合材料,根据经验进行铺层角度、铺层顺序、铺层厚度设定,使其于碰撞仿真中可以保持原抗撞性能,将此铺层信息作为接下来优化的初始参数。
b.由于B柱加强板皆为设计区域,故采用8号各向异性材料模型创建铺层形成层合板结构,属性定义针对设计区域设置选取PCOMPP(一种材料属性形式),由于Optistruct优化软件只能对材料进行减薄而不具备增厚的功能,故需建立所谓的初始厚度较厚的超级层。该超级层的厚度通过初始优化参数进行设置。
(3)加载的设置
提取B柱模型的抗冲击性仿真中的B柱加强板最大相对位移来计算所需等效力的数值,将其加载在与碰撞台车相接触的B柱加强板结构优化模型中,并按照加载点的个数,将等效力进行平均分配加载;根据x、y、z轴的力的数量级进行合理忽略简化,所得的静态仿真模型可以与动态仿真模型获得较好的位移拟合结果,之后在此基础上进行静态优化仿真模拟;
2)碳纤维增强复合材料B柱加强板静态工况铺层优化
(1)拓扑优化
a.定义设计变量,对于拓扑优化首先创建变量,定义质量和参考点位移两个响应,二者分别作为优化响应和约束响应;创建约束,将参考点位移响应作为约束,设置约束上限值,将质量最小定义为优化目标;
b.调整优化控制参数,包括最大允许迭代次数、目标容差值以及优化步长;
c.定义控制卡片,在此处需要在控制卡片部分定义OUTPUT(输出)卡片关键字为FSTOSZ(一种类型关键字),然后保存或输出文件;
d.最后进行检查,如果没有错误信息则可开始进行优化,优化结束后查看结果文件,在HyperView中打开优化结果文件检查单元厚度;
(2)尺寸优化
a.为保证铺层信息的正确传递,单个方向的厚度信息以铺层形状的形式进行输出,首先重新设置控制卡片,删除“OMIT”卡片,将FSTOSZ(一种类型关键字)改为SZTOSH(一种类型关键字);
b.编辑Laminate(层合板)参数,将层合板形式设置为Total(总共),分别定义每层尺寸变量,确保修改每层厚度上限,将不能修改的规整为相近整数,依旧将参考点位移值作为约束条件,优化目标仍然是质量最小;
c.指定可加工的厚度的值,在优化设置后进行检查优化,并查看铺层厚度和数目;
(3)顺序优化
对于复合材料来说,铺层顺序与其层合板的各向异性力学性能以及工艺可行性都有着紧密的联系,复合材料铺层顺序的变化会影响其结构的整体性能。在计算最优的层叠次序之前,首先对铺层顺序变量设置制造约束,定义0度、90度、45度和-45度四个方向铺层最多可以有2层相同铺层连续出现;+45度和-45度成对出现来保证对称性,最大限度避免层间应力的产生,优化约束和优化目标不变,最后进行铺层顺序优化;检查输出的结果文件,在浏览器中查看各个迭代步的铺层顺序;
3)B柱加强板优化结果于解耦子结构动态工况下验证计算
(1)根据铺层结果进行规整分区来改进B柱加强板;
(2)对带有优化后碳纤维增强复合材料B柱加强板的模型进行解耦子工况侧撞台车碰撞的模拟仿真计算;
(3)查看优化后B柱与原车型中的B柱的对比数据,即通过其B柱加强板上对应人体胸部的参考点的侵入量和侵入速度最大值的对比来验证其性能的改善程度,同时获得了其轻量化效果。
实施例:
参阅图1,本发明接下来结合实例介绍利用所述的汽车B柱加强板碳纤维增强复合材料优化设计方法;
1.从整车有限元模型中解耦B柱子结构
1)参阅图2,根据原有整车侧面碰撞模型,保留与B柱相连接的部件来保证其变形情况的准确性,即门槛梁、侧围板、顶盖横梁、车门管状防撞梁,再加上B柱结构作为子结构模型。
2)采用与B柱相连接部件中与B柱相连接位置接近的区域以及侧面碰撞台车缓冲蜂窝铝最前端表面作为碰撞参考点,进行侧面碰撞有限元计算,通过Ls-prepost后处理软件提取碰撞参考点的位移-时间曲线数据。
3)将提取出的曲线数据通过MATLAB编程软件转化为Ls-dyna仿真计算软件可识别的数据信息,再通过文本编辑器导入简化后只保留B柱和其周围部件以及侧面碰撞台车前端的子工况模型的计算k文件,便可准确模拟整车侧面碰撞情况,并且大大减少了计算时间,提高工作效率。
2.碳纤维增强复合材料B柱加强板优化
1)B柱加强板静态工况建立
(1)结构简化
参阅图4,根据原有整车有限元模型,只保留其原有金属B柱加强板形状模型结构,包括其位于中部的B柱加强板有限元模型中螺栓孔的位置。将B柱加强板设定为统一初始优化厚度,为2.5mm。进行其几何表面的二维网格质量检查,要求有翘曲度大于5、纵横比大于5、单元长度小于7.5、雅克比小于0.7。
(2)碳纤维增强复合材料B柱加强板的建立
a.首先在Ls-dyna仿真计算软件中采用54-55号各向异性材料模型创建层合板结构。材料参数如表1所示。采取角度分别为0度、45度、-45度和90度进行铺层,单层厚度为0.25mm,总厚度为2.5mm。根据经验铺层顺序设置为[0/45/-45/90/45/45/90/-45/45/0]。进行子工况侧面碰撞仿真计算后可满足原金属B柱加强板性能,即保证了其位于人体胸部参考点的侵入速度以及侵入量不增大,且重量由1.034kg减少为0.2881kg,减少72.14%。
表1
RHO | EA | EB | μ21 | GAB | GBC |
1.75e-9 | 145200 | 8100 | 0.002 | 4400 | 1530 |
GCA | Xt | Xc | Yt | Yc | S |
1530 | 2455 | 2455 | 35 | 35 | 61 |
b.在Optistruct优化软件中,B柱加强板采用8号各向异性材料模型创建铺层形成层合板结构,材料参数见表1。由于该软件只能对铺层进行减薄,而不具备加厚的功能,故需设计厚度较厚的超级层。根据前步骤中所得满足碰撞性能的碳纤维增强复合材料铺层信息作为初始优化参数。属性定义针对设计区域设置为PCOMPP(一种材料属性形式),设定4个铺层方向,即0度、45度、-45度和90度,各层厚度分别为0.5mm、1mm、0.5mm、0.5mm,总厚度为2.5mm。
(3)加载的设置
参阅图5,采用施加等效静载力的方法,根据B柱模型的抗冲击性仿真中的B柱加强板最大相对位移计算所需等效力的数值,将其加载在与碰撞台车相接触的B柱加强板结构优化模型处,并按照加载点的个数,将接触力进行平均分配加载。由于X、Z两轴方向与Y轴的接触力峰值相比,数量级相差甚远,故X、Z两轴接触力可忽略不计,只提取出Y轴的峰值力进行等效分析。设置8份1N的集中力均匀地加载在B柱结构优化模型上与碰撞台车相接触的部位,作为均布力进行仿真模拟计算,基于已定义好的B柱加强板结构进行分析计算可实现与碰撞模型中相同的相对位移,为37.317mm。
2)碳纤维增强复合材料B柱加强板静态工况铺层优化
(1)拓扑优化
a.首先定义设计变量Topo(定义的变量名称),选中设计空间的铺层。定义质量响应M和位移响应D,分别作为优化目标响应和约束响应。创建约束,将位移响应D作为约束,设置上限值为37.35mm,将质量响应M最小作为优化目标。
b.调整优化控制参数,最大允许迭代次数为80,目标容差值使用最小值0.005以及优化步长使用0.5。
c.定义控制卡片,在此处需要在控制卡片部分定义OUTPUT(输出)卡片关键字为FSTOSZ(一种类型关键字),然后保存或输出文件。
d.最后进行检查,如果没有错误信息则进行优化。优化结束后查看结果文件,在HyperView中打开优化结果文件检查单元厚度,B柱加强板厚度需要主要减薄部分如图6所示。
(2)尺寸优化
a.首先设置卡片,将FSTOSZ(一种类型关键字)改为SZTOSH(一种类型关键字),删除OMIT(一种类型关键字),以保证铺层信息的正确传递,即单个方向的厚度信息以铺层形状的形式进行输出。
b.编辑Laminate(层合板)参数,将层合板形式设置为Total(总共)。定义尺寸设计变量Size,将每层厚度上限修改为单层厚度的二倍即0.5mm,各个层均需要更新,将不能修改规整的归整为相近的整数。优化目标仍然是质量最小。
c.指定可加工厚度为0.25mm,进行检查后若没有错误信息则进行优化。优化结束后查看结果文件,查看铺层厚度和数目,共有16个铺层项目,单层厚度为0.25mm。
(3)顺序优化
在对最优的层叠次序进行优化之前,首先设置制造约束,定义0度、90度、45度和-45度四个方向铺层最多可以有2层连续出现,+45度和-45度成对出现。优化约束和优化目标仍是上一步尺寸优化中的位移和质量最小目标。检查输出的结果文件,在浏览器中查看各个迭代步的铺层顺序。
3)B柱加强板优化结果于解耦子结构动态工况下验证计算
(1)根据铺层优化结果对于B柱加强板进行工程解读,整合分为四个区域,如图7所示。1区厚度为2.5mm,铺10层,铺层顺序为[45/-45/0/0/45/-45/0/0/90/90];2区厚度为0.5mm,铺2层,铺层顺序为[45/-45];3区厚度为1mm,铺4层,铺层顺序为[45/-45/90/90];4区厚度为1mm;铺4层,铺层顺序为[0/0/45/-45]。
(2)对带有优化后碳纤维增强复合材料B柱加强板的模型进行解耦子工况侧撞台车碰撞的模拟仿真计算;
(3)查看仿真计算结果数据,相比于原结构其性能有所改善,即B柱加强板上对应于人体胸部的参考点的相对侵入量有所降低,如表2所示。原金属B柱加强板质量为1.034千克,优化前碳纤维增强复合材料B柱加强板质量为0.2881千克,经过铺层优化后可将其质量将为0.2174千克,优化前后减重比例可达24.54%,与原金属材料相比最终可实现减重比例高达78.97%,减重效果明显。
表2
综上所述,本发明提出的B柱子结构解耦以及碳纤维增强复合材料B柱加强板优化分析方法,通过子结构解耦来结构简化进行有限元分析计算,在拟合程度优异的情况下大大节省了计算时间,提高了计算效率;然后在此基础上进行B柱加强板静态工况的铺层优化,包括拓扑优化、尺寸优化及顺序优化三部分;最后将改进碳纤维增强复合材料B柱加强板在动态工况下进行验证仿真。本发明充分利用了复合材料与金属材料相比的材料可设计性能,在保证B柱抗撞性能的基础上轻量化效果明显。
Claims (6)
1.一种汽车B柱加强板碳纤维增强复合材料优化设计方法,其特征在于,所述的汽车B柱加强板碳纤维增强复合材料优化设计方法的步骤如下:
1)从整车有限元模型中解耦B柱子结构:
(1)根据原有整车侧面碰撞模型,保留与B柱相连接的部件来保证其变形情况的准确性,即门槛梁、侧围板、顶盖横梁、车门管状防撞梁,再加上B柱结构作为子结构模型;
(2)采用与B柱相连接部件中与B柱相连接位置接近的区域以及侧面碰撞台车缓冲蜂窝铝最前端表面作为碰撞参考点,进行侧面碰撞有限元计算,通过Ls-prepost后处理软件提取碰撞参考点的位移-时间曲线数据;
(3)将提取出的曲线数据通过MATLAB编程软件转化为Ls-dyna仿真计算软件可识别的数据信息,再通过文本编辑器导入简化后只保留B柱和其周围部件以及侧面碰撞台车前端的子工况模型的计算k文件,便准确模拟整车侧面碰撞情况;
2)碳纤维增强复合材料B柱加强板优化:
(1)B柱加强板静态工况建立;
(2)碳纤维增强复合材料B柱加强板静态工况铺层优化:
A.拓扑优化;
B.尺寸优化;
C.顺序优化;
(3)B柱加强板优化结果于解耦子结构动态工况下验证计算。
2.按照权利要求1所述的汽车B柱加强板碳纤维增强复合材料优化设计方法,其特征在于,所述的B柱加强板静态工况建立是指:
(1)结构简化
根据已有B柱加强板结构形状的几何尺寸,将B柱加强板厚度设定成统一初始厚度,并且保留整车模型中原金属B柱加强板中的形状结构参数,将简化结构中预留出螺栓固定孔结构,并进行其网格质量检查;
(2)碳纤维增强复合材料B柱加强板的建立
a.将即将优化的部件B柱加强板更改为碳纤维增强复合材料,根据经验进行铺层角度、铺层顺序、铺层厚度设定,使其于碰撞仿真中可以保持原抗撞性能,将此铺层信息作为接下来优化的初始参数;
b.由于B柱加强板皆为设计区域,故采用8号各向异性材料模型创建铺层形成层合板结构,属性定义针对设计区域设置选取PCOMPP,由于Optistruct优化软件只能对材料进行减薄而不具备增厚的功能,故需建立所谓的初始厚度较厚的超级层,该超级层的厚度通过初始优化参数进行设置;
(3)加载的设置
提取B柱模型的抗冲击性仿真中的B柱加强板最大相对位移来计算所需等效力的数值,将其加载在与碰撞台车相接触的B柱加强板结构优化模型中,并按照加载点的个数,将等效力进行平均分配加载;根据x、y、z轴的力的数量级进行合理忽略简化,所得的静态仿真模型可以与动态仿真模型获得较好的位移拟合结果,之后在此基础上进行静态优化仿真模拟。
3.按照权利要求1所述的汽车B柱加强板碳纤维增强复合材料优化设计方法,其特征在于,所述的拓扑优化是指:
a.定义设计变量,对于拓扑优化首先创建变量,定义质量和参考点位移两个响应,二者分别作为优化响应和约束响应;创建约束,将参考点位移响应作为约束,设置约束上限值,将质量最小定义为优化目标;
b.调整优化控制参数,包括最大允许迭代次数、目标容差值以及优化步长;
c.定义控制卡片,在此处需要在控制卡片部分定义OUTPUT(输出)卡片关键字为FSTOSZ,然后保存或输出文件;
d.最后进行检查,如果没有错误信息则可开始进行优化,优化结束后查看结果文件,在HyperView中打开优化结果文件检查单元厚度。
4.按照权利要求1所述的汽车B柱加强板碳纤维增强复合材料优化设计方法,其特征在于,所述的尺寸优化是指:
a.为保证铺层信息的正确传递,单个方向的厚度信息以铺层形状的形式进行输出,首先重新设置控制卡片,删除“OMIT”卡片,将一种类型关键字FSTOSZ改为另一种类型关键字SZTOSH;
b.编辑层合板Laminate参数,将层合板形式设置为Total,分别定义每层尺寸变量,确保修改每层厚度上限,将不能修改的规整为相近整数,依旧将参考点位移值作为约束条件,优化目标仍然是质量最小;
c.指定可加工的厚度的值,在优化设置后进行检查优化,并查看铺层厚度和数目。
5.按照权利要求1所述的汽车B柱加强板碳纤维增强复合材料优化设计方法,其特征在于,所述的顺序优化是指:
对于复合材料来说,铺层顺序与其层合板的各向异性力学性能以及工艺可行性都有着紧密的联系,复合材料铺层顺序的变化会影响其结构的整体性能,在计算最优的层叠次序之前,首先对铺层顺序变量设置制造约束,定义0度、90度、45度和-45度四个方向铺层最多可以有2层相同铺层连续出现;+45度和-45度成对出现来保证对称性,最大限度避免层间应力的产生,优化约束和优化目标不变,最后进行铺层顺序优化;检查输出的结果文件,在浏览器中查看各个迭代步的铺层顺序。
6.按照权利要求1所述的汽车B柱加强板碳纤维增强复合材料优化设计方法,其特征在于,所述的B柱加强板优化结果于解耦子结构动态工况下验证计算是指:
(1)根据铺层结果进行规整分区来改进B柱加强板;
(2)对带有优化后碳纤维增强复合材料B柱加强板的模型进行解耦子工况侧撞台车碰撞的模拟仿真计算;
(3)查看优化后B柱与原车型中的B柱的对比数据,即通过其B柱加强板上对应人体胸部的参考点的侵入量和侵入速度最大值的对比来验证其性能的改善程度,同时获得了其轻量化效果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810409682.XA CN108563905B (zh) | 2018-05-02 | 2018-05-02 | 汽车b柱加强板碳纤维增强复合材料优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810409682.XA CN108563905B (zh) | 2018-05-02 | 2018-05-02 | 汽车b柱加强板碳纤维增强复合材料优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108563905A CN108563905A (zh) | 2018-09-21 |
CN108563905B true CN108563905B (zh) | 2021-09-21 |
Family
ID=63537641
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810409682.XA Expired - Fee Related CN108563905B (zh) | 2018-05-02 | 2018-05-02 | 汽车b柱加强板碳纤维增强复合材料优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108563905B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108984994B (zh) * | 2018-09-27 | 2023-05-30 | 东风小康汽车有限公司重庆分公司 | 一种基于有限元的汽车钣金件结构设计方法 |
CN110110480B (zh) * | 2019-05-21 | 2020-09-11 | 吉林大学 | 一种考虑铺层相容性的碳纤维保险杠防撞梁结构优化设计方法 |
CN110826222B (zh) * | 2019-11-05 | 2023-03-24 | 上海波客实业有限公司 | 一种汽车碳纤维增强复合材料覆盖件正向开发方法 |
CN111143946A (zh) * | 2019-12-28 | 2020-05-12 | 江苏亨睿碳纤维科技有限公司 | 一种复合材料汽车b柱的设计方法 |
CN113128079A (zh) * | 2019-12-30 | 2021-07-16 | 上海荣泰健康科技股份有限公司 | 按摩椅缓冲件的顶压测试方法 |
CN111169039A (zh) * | 2020-01-10 | 2020-05-19 | 广东亚太新材料科技有限公司 | 一种汽车用碳纤维复合材料b柱加强件的制备方法 |
CN111324980B (zh) * | 2020-01-21 | 2022-08-09 | 重庆长安汽车股份有限公司 | 一种汽车结构轻量化分级优化设计方法 |
CN113033051B (zh) * | 2021-03-24 | 2023-09-12 | 天津中德应用技术大学 | 碳纤维增强复合材料薄板的模压成型工艺优化方法 |
CN113312701B (zh) * | 2021-04-30 | 2024-03-19 | 中铝材料应用研究院有限公司 | 基于拓扑及尺寸优化的全铝客车车身门立柱结构设计方法 |
CN113536619B (zh) * | 2021-05-24 | 2023-05-16 | 天津中德应用技术大学 | 一种碳纤维增强复合材料模压成型工艺参数优化方法 |
CN113420376B (zh) * | 2021-06-17 | 2022-07-22 | 吉林大学 | 基于多尺度的碳纤维复合材料抗冲击力学性能仿真方法 |
CN114379655B (zh) * | 2022-01-21 | 2023-03-03 | 广东汇天航空航天科技有限公司 | 后车身框架和车辆 |
CN114741846B (zh) * | 2022-03-15 | 2023-12-01 | 中国民航大学 | 一种超低热膨胀复合材料支撑结构的设计方法 |
CN116522420B (zh) * | 2023-06-25 | 2023-08-29 | 山东石油化工学院 | 用于曲面柱型复合材料的力学性能仿真设计方法及*** |
CN117910323B (zh) * | 2024-03-15 | 2024-05-31 | 西北工业大学 | 一种Double-Double铺层复合材料C型梁变厚度优化方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102013214717A1 (de) * | 2013-07-29 | 2015-01-29 | Bayerische Motoren Werke Aktiengesellschaft | Biegeträger aus faserverstärktem Kunststoff einer Kraftfahrzeugkarosserie |
CN105653768A (zh) * | 2015-12-25 | 2016-06-08 | 上海交通大学 | 基于粒子群优化算法的轿车车身轻量化结构的实现方法 |
CN106484979A (zh) * | 2016-09-29 | 2017-03-08 | 吉林大学 | 基于独立评价工况的前防撞梁总成轻量化设计方法 |
CN107679343A (zh) * | 2017-10-31 | 2018-02-09 | 吉林大学 | 连续纤维增强热固性复合材料座椅骨架优化方法 |
-
2018
- 2018-05-02 CN CN201810409682.XA patent/CN108563905B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102013214717A1 (de) * | 2013-07-29 | 2015-01-29 | Bayerische Motoren Werke Aktiengesellschaft | Biegeträger aus faserverstärktem Kunststoff einer Kraftfahrzeugkarosserie |
CN105653768A (zh) * | 2015-12-25 | 2016-06-08 | 上海交通大学 | 基于粒子群优化算法的轿车车身轻量化结构的实现方法 |
CN106484979A (zh) * | 2016-09-29 | 2017-03-08 | 吉林大学 | 基于独立评价工况的前防撞梁总成轻量化设计方法 |
CN107679343A (zh) * | 2017-10-31 | 2018-02-09 | 吉林大学 | 连续纤维增强热固性复合材料座椅骨架优化方法 |
Non-Patent Citations (2)
Title |
---|
基于多工况的乘用车前防撞梁总成轻量化设计;张君媛 等;《汽车安全与节能学报》;20171231;第252-260页 * |
碳纤维复合材料汽车B柱简化结构参数分析;薛姣;《中国优秀硕士学位论文全文数据库》;20130815;C035-110 * |
Also Published As
Publication number | Publication date |
---|---|
CN108563905A (zh) | 2018-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108563905B (zh) | 汽车b柱加强板碳纤维增强复合材料优化设计方法 | |
CN107679343B (zh) | 连续纤维增强热固性复合材料座椅骨架优化方法 | |
Duan et al. | Parametric modeling and multiobjective crashworthiness design optimization of a new front longitudinal beam | |
Duan et al. | Multi-objective reliability-based design optimization for the VRB-VCS FLB under front-impact collision | |
Hou et al. | Multi-objective and multi-constraint design optimization for hat-shaped composite T-joints in automobiles | |
CN111475977B (zh) | 一种基于碰撞性能优化的汽车白车身轻量化设计方法 | |
Rajasekaran et al. | Multi-Objective Optimization of Material Layout for Body-In-White using Design of Experiments. | |
CN109117532B (zh) | 汽车轻量化优化方法 | |
Tisza | Recent development trends in sheet metal forming | |
Liu et al. | Structural design and global sensitivity analysis of the composite B-pillar with ply drop-off | |
CN111950080A (zh) | 一种车身轻量化设计方法 | |
CN103569236A (zh) | 一种汽车车身前部结构的设计方法及其汽车的设计方法 | |
Liu et al. | Multi-objective optimization of an auto panel drawing die face design by mesh morphing | |
CN109255141B (zh) | 一种汽车车身正向概念设计截面形状优化方法 | |
CN112464382B (zh) | 汽车仪表板横梁尺寸优化设计方法 | |
Baskin et al. | A case study in structural optimization of an automotive body-in-white design | |
CN105468826A (zh) | 复合材料的设计方法 | |
Meng et al. | Intelligent layout optimization of reconfigurable flexible fixture for assembling multiple aircraft panels | |
Botkin | Structural optimization of automotive body components based on parametric solid modeling | |
CN116882058A (zh) | 一种用于汽车地板的复合材料轻量化设计与优化方法 | |
Ma et al. | Parallel optimization of design and manufacture for carbon fiber reinforced plastic oil pan based on the thickness distribution | |
Botkin | Modelling and optimal design of a carbon fibre reinforced composite automotive roof | |
Zhang | Lightweight and performance of anti-collision strength of automobiles based on carbon fiber composites | |
Ma et al. | Optimal design of carbon fiber b-pillar structure based on equal stiffness replacement | |
Chen et al. | Carbody structural lightweighting based on implicit parameterized model |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210921 |