横截面参数化网格变形算法及耐撞性优化应用

薄壁梁结构是汽车等运载工具的主要承载构件,提高该类结构的耐撞性对乘员安全具有重要意义。然而,形状优化设计要求多组有限元模型与仿真分析,因此需要特定的建模技术或人工交互。本文提出了一种基于横截面形状的参数化网格变形方法,以实现已有有限元模型的有效重用。以给定有限元模型为输入,采用基于各向异性径向基函数网格变形方法,并结合骨架内嵌空间,可快速生成适用于仿真分析的有限元模型变体。以S形梁轴向冲击耐撞性设计为例,采用所提方法改变构件塑性铰区域的横截面形状,可快速(低于4 s)获取100组局部变形有限元模型,并采用代理模型技术和多目标遗传算法优化结构耐撞性。数值结果显示,构件耐撞性获显著提高,验证了所提参数化变形方法的有效性,展示了与一般形状优化框架的可集成性。     

基于塑性铰理论的车身概念设计正碰分析与优化方法

在概念设计阶段,车身碰撞安全性能评价是一个难点问题,需要详细的结构模型,本文基于塑性铰理论提出采用梁单元简化模型对框架车身进行概念设计阶段的耐撞性评估和优化设计方法。首先,介绍了关于箱型截面薄壁梁弯曲特性研究的理论模型与计算过程,接着赋予梁单元塑性铰的特性,模拟薄壁梁变形,再对框架车身进行了碰撞仿真。将仿真结果与详细模型对比,以分析简化模型的精度及可靠性。最后,以此为基础对框架车身进行耐撞性优化。结果表明,该简化模型易于创建,且有较高的精度,可用于概念设计阶段梁结构的设计工作。     

具有恒定冲击载荷的梯度泡沫金属材料设计

多胞材料可通过大变形大量地吸收冲击能量，引入密度梯度可进一步提高其耐撞性。梯度多胞材料的宏观力学响应对材料密度分布极为敏感，不同类型的细观构型的影响也极为不同。已有的研究工作主要局限在对给定的密度梯度分析其动态响应，较少对耐撞性设计方法进行研究。本文针对梯度闭孔泡沫金属材料，基于非线性塑性冲击波模型发展了耐撞性反向设计方法，以维持冲击物受载恒定为目标，运用级数法获得了简化模型和渐近解。利用变胞元尺寸法构建了连续梯度变化的三维Voronoi细观有限元模型，并利用ABAQUS/Explicit有限元软件对理论设计进行数值验证。结果表明，反向设计理论简化模型的渐近解对于梯度闭孔泡沫金属材料的耐撞性设计是有效的，所提出的耐撞性设计方法在控制冲击吸能过程和冲击物受载方面具有指导意义。     

基于响应面模型的薄壁构件耐撞性优化设计

建立了薄壁构件的结构耐撞性优化模型，通过试验设计在设计空间中选择少量样本获得设计的响应特性，建立响应面模型，并应用Pareto遗传算法进行结构耐撞性的优化设计。     

某型飞机机翼前缘抗鸟撞结构设计与试验验证

针对某型飞机机翼前缘抗鸟撞性能不满足适航要求的问题，采用"数值仿真-试验验证-再仿真"的研究思路，对该结构进行了抗鸟撞优化设计。首先，通过有限元数值仿真，分别对具有三角板结构和前墙结构的两种新型前缘结构抗鸟撞能力进行了分析。仿真结果表明:具有前墙结构的机翼前缘抗鸟撞能力优于原始结构和带三角板结构的机翼前缘；在鸟撞过程中，这种具有前墙结构的机翼前缘通过利用破损蒙皮继续变形吸能的方式提高了结构的抗鸟撞性能。基于此，对带前墙结构的机翼前缘进行了抗鸟撞试验，一方面验证了数值模拟方法的准确性，另一方面验证了前墙结构的抗鸟撞效果。最后，采用数值仿真方法对带前墙结构的机翼前缘结构进行了参数分析，得到了前缘蒙皮厚度和前墙厚度与结构抗鸟撞性能的关系，并基于结构承载和抗鸟撞能力的综合要求，确定了最终结构参数。分析表明，优化后的机翼前缘结构不仅满足抗鸟撞要求，而且实现结构减重30%。     

基于物理的变形曲线在结构形状优化中的应用

用各向异性模型定义NURBS曲线的变形能,基于所提出的变形能模型,用有限元法对NURBS曲线表达的设计边界进行模态计算,然后,将设计边界用模态向量的线性组合参数化表示。将这种基于边界特征向量的几何形状表示方法应用于优化参数定义．提出了一种适用于结构形状优化的自适应几何精化方法,它有效地将户型有限元分析、优化方法和设计边界的形状表达集成在一起。     

基于元胞自动机的变厚度薄壁梁侧向耐撞性优化设计方法

薄壁结构是汽车等运载工具的重要防护装置,除了其轴向防撞能力外,侧向耐撞性能分析与提升方式也非常重要。研究基于薄壁结构厚度合理分布的侧向耐撞性能提升方式和建立基于元胞自动机的变厚度薄壁梁侧向耐撞性优化方法。以汽车B柱受力环境和性能要求为设计需求,首先利用所建立的方法给出了连续变厚度的薄壁梁厚度分布设计,其性能较常规的等厚度薄壁梁最大侵入位移大幅下降(下降82%),验证了变厚度设计的有效性;然后,考虑单向变厚度便于柔性轧制工艺制成TRB,给出了轴向连续变厚度薄壁梁的厚度分布设计,该设计较等厚度梁最大侵入位移下降73%;与连续变厚度梁相比,在侵入位移降低量略小的情况下,实现了可制造性。设计实例表明本文提出的连续变厚度设计能够有效提高侧向耐撞性能,所建立的方法能够获得合理的厚度分布设计,是有效的耐撞性优化设计方法。     

一种有限元模型修正中的参数选择方法

有限元模型修正面临的首要问题是待修正参数的选择.目前主要是基于灵敏度分析的结果进行参数选择.本文从结构固有频率分析的能量法出发,论证了某阶固有频率主要受到同阶振型中发生弹性变形部位刚度的影响,据此提出一种待修正参数的选择方法,以振型中弹性变形为依据,选择弹性变形较大部位的参数作为待修正参数.用一个梁模型仿真验证了该方法的正确性,并将其应用于系杆拱桥的有限元模型修正中,修正后有限元计算结果与试验结果的最大误差缩小至3.3%.     

复合材料波纹梁盒段耐撞性的数值模拟

运用MSC/DYTRAN有限元软件,结合参数等效的方法对复合材料波纹梁盒段的碰撞过程进行了数值模拟,获得了波纹梁盒段碰撞时的破坏过程、载荷-时间曲线、能量吸收能力等数据,将计算结果与M.A.McCarthy的试验结果进行对比,可以发现两者比较吻合。从而说明采用参数等效的方法可以得到复合材料结构耐坠撞性设计中需要的数据。     

船舶碰撞动力学过程的数值仿真研究

在分析显式非线性有限元基本理论和关键技术的基础上 ,探讨了船舶动力学过程的数值仿真方法 ,并对碰撞过程中船体的力学性能进行了分析 ,获得并讨论了碰撞力、能量吸收和损伤变形的时序结果 ,所得结论具有一般性。碰撞仿真分析中 ,被撞船的舷侧受撞区域和撞击船的首部区域作为可变形结构处理 ,而其余区域则视为刚体。被撞船周围的水划分成Euler有限体积网格 ,船体结构则划分成Lagrange有限元网格 ,并采用耦合技术将它们联系在一起。撞击船周围的水的影响采用附连水质量进行处理。     

Implementation of Convergence in Adaptive Global Digital Image Correlation

In FE based global digital image correlation (DIC) a finite element mesh is used to describe the deformation of the region of interest (ROI). However, the identification of an optimal mesh is a difficult problem and is often obtained by using “mechanical” pre-knowledge of the solution. In Finite Element Analysis (FEA) an optimal mesh can be found without any pre-knowledge of the solution by using mesh adaptivity, where an initial (non optimal) mesh is refined until the optimal solution is obtained. Refinement of the mesh can be based on error and/or convergence estimators. Despite the fundamental differences between FEA and DIC, in the present article the convergence procedure is successfully used in a recently published global FE based DIC method. In the used global DIC method elements can receive higher order shape functions, also known as p-elements. Using the aforementioned algorithm, also called p-DIC, refinement to a non-uniform higher order mesh is possible. Using the non-uniform mesh, an optimal mesh can be obtained for each section of the ROI. The presented study shows that a convergence scheme can be used to automatically control the mesh refinement in a global DIC approach. The convergence boundary, in percentage, is a more intuitive boundary than the absolute error boundary used in the original p-DIC approach. The procedure is validated using numerical examples and the robustness to experimental variables is investigated. Finally, the complete procedure is tested against a wide range of practical examples.     

AERODYNAMIC OPTIMIZATION OF HIGH-SPEED TRAIN BASED ON RBF MESH DEFORMATION

基于局部型函数三维参数化方法、改进的蚁群算法和改进的克里金（Kriging）代理模型,开展了列车头型的三维气动减阻优化设计研究。为了避免复杂几何外形大变形情况下千万量级网格的重复生成,提高高速列车头型优化设计的效率,引入了缩减控制点的径向基函数网格变形技术。优化结果表明：径向基函数网格变形技术在不降低网格质量的情况下可以有效缩短网格变形的时间消耗,能够用于复杂几何外形的气动优化设计;在给定的设计空间内,控制鼻锥外形的6个关键设计参数对列车气动阻力的影响呈单调递增关系;优化后,在满足约束条件的情况下,简化外形列车的整车气动阻力减小5.41%,头尾车减阻效果明显,中间车气动阻力基本不变。     

Crashworthiness analysis and optimization of fourier varying section tubes

Thin-walled structures are widely used as energy absorption devices for their proven advantages on lightweight and crashworthiness. However, a majority of studies have being focus on exploring separately the crashworthiness of the thin-walled structure with a specific geometric section, such as circular, square, hexagon, octagon etc., and little research has investigated the relationship of crashworthiness among thin-walled structures with different sections systematically. This paper utilizes Fourier series expansion to generate a series of novel sectional configurations, namely Fourier varying sectional tubes (FVSTs), to look into their advantages of crashworthiness, thereby developing some FVSTs with highest possible energy absorption capacity. Based on the validated finite element (FE) models, parametric analysis is conducted to investigate the effects of cross-sectional configuration, perimeter and thickness of FVSTs on collapse mode and energy absorption. The results showed that the collapse modes of FVSTs are fairly sensitive to cross-sectional configuration, perimeter and wall thickness. Of these FVSTs generated, the highest specific energy absorption (SEA) increases 77.54% by increasing perimeter and 69.73% by decreasing wall thickness. Finally, a discrete optimization based on the orthogonal arrays is conducted to obtain the optimal FVST for maximizing SEA under the constraint of the initial peak crushing force (IPCF). The optimized FVSTs are of superior crashworthiness and great potential as an energy absorber.     

附加自由阻尼结构系统的有限元分析和拓扑优化设计方法研究

给出了一种适用于车身等复杂附加阻尼结构系统的有限元建模和动态特性分析方法,包括附加自由阻尼薄壁结构和附加阻尼材料粘弹性特性的有限元建模方法以及附加阻尼结构动态特性的有限元分析方法。在此基础上,又以强迫振动响应最小为优化设计目标,给出了一个附加阻尼结构的拓扑优化设计方法,包括优化设计问题的列式和适用的求解算法。通过对一个简单薄壁构件和一个车身地板上的附加阻尼结构的拓扑优化设计,验证了提出的有限元分析和拓扑优化设计方法的有效性。     

仿生多胞薄壁管耐撞性分析及优化

为提高薄壁管结构耐撞性，以雀尾螳螂虾螯为仿生原型，结合仿生学设计方法，设计一种含正弦胞元的多胞薄壁管结构。以初始峰值载荷、比吸能和碰撞力效率为耐撞性指标，通过有限元数值模拟分析了不同碰撞角度（0o、10o、20o和30o）条件下，仿生胞元数对薄壁管耐撞性的影响，通过多目标的复杂比例评估法获取仿生薄壁管的最优胞元数。基于不同碰撞角度权重因子组合，设置了4种单一角度工况和3种多角度工况，采用多目标粒子群优化方法获取了不同工况下薄壁管结构最优胞元高宽比和壁厚。复杂比例评估结果表明，胞元数为4的薄壁管为最优晶胞数仿生薄壁管。优化结果表明，单一角度工况下，最优结构参数高宽比的范围为0.88～1.50，壁厚的范围为0.36～0.60 mm，碰撞角度为0o和10o的最优高宽比明显小于碰撞角度为20o和30o的；多角度工况下，最优高宽比范围为1.01～1.10，壁厚范围为0.49～0.57 mm。     

考虑车辆高速和低速耐撞性的多目标优化设计

<正>面高速耐撞性设计,需确保车内乘员的人身安全,使车内乘员的人身伤害降到最低;正面低速耐撞性设计,要求尽量减少车辆在碰撞中的损伤,使车辆具有良好的碰撞损伤修复经济性。本文根据上述特点,将车辆的正面高速耐撞性与正面低速耐撞性相结合,并且考虑车身部件的轻量化,提出了一种优化设计方法。该方法将保险杠、吸能盒内外板以及前纵梁内外板五个部件厚度作为优化设计变量,在正面低速碰撞中,以前纵梁吸收的碰撞总能量作为其是否发生较大变形的依据,利用代理模型和遗传算法进行了优化计算。最后,将该方法应用在某一车型的耐撞性设计中。     

Elastic structural response of prismatic metal sandwich tubes to internal moving pressure loading

The analytical and numerical modeling of the structural response of a prismatic metal sandwich tube subjected to internal moving pressure loading is investigated in this paper. The prismatic core is equivalent to homogeneous and cylindrical orthotropic solids via homogenization procedure. The sandwich tube with the “effective” homogenized core is modeled using multi-layer sandwich theory considering the effects of transverse shear deformation and compressibility of the core; moreover, the solutions are obtained by using the precise integration method. Several dynamic elastic finite element (FE) simulations are carried out to obtain the structural response of the tube to shock loading moving at different velocities. The comparison between analytic solutions and FE simulations demonstrates that the transient analytical model, based on the proposed sandwich model, is capable of predicting the critical velocity and the dynamic structural response of the sandwich tube with the “effective” homogenized core with a high degree of accuracy. In addition, the critical velocity predicted using FE simulations of the complete model is not in agreement with that of the effective model. However, the structural response and the maximum amplification factors obtained using FE simulations of the complete model are nearly similar to that of the effective model, when the shock loading moves at the critical velocity. The influences of the relative density on the structural response are studied, and the capabilities of load bearing for sandwich tubes with different cores are compared with each other and with the monolithic tube. The results indicate that Kagome and triangle-6 are preferred among five topologies.     

Shape Optimization for 3D Turbulent Flows Using Automatic Differentiation

Abstract

We present through this paper some new results with our approach for optimal shape design based on a CAD-free framework for shape and unstructured mesh deformations, automatic differentiation by program for the gradient computation and mesh adaption by metric control. Automatic differentiation allows for an easy, reliable and fast discrete adjoint computation. We managed to get the Jacobian of our Navier-Stokcs solver including a k-e turbulence model and wall-laws. The CAD-free framework is shown to be particularly convenient for optimization when the mesh connectivities and the shape discretization arc variable during optimization. Using these ingredients constrained optimization for turbulent transonic flows has been investigated in both 2 and 3D.