一种低碳约束下的结构拓扑优化方法

技术领域

本发明属于结构优化相关技术领域，具体为一种低碳约束下的结构拓扑优化方法。

背景技术

随着经济社会的发展，环境和社会问题(例如环境污染、全球变暖、能源危机等)正日益凸显。低碳发展作为一种新的发展模式，被认为是应对全球变暖和能源危机的有效模式。如何采取有效措施减缓温室气体排放已经成为目前产品设计发展的热点方向。

一般来说，零部件的结构布局越好，其生命周期所产生的碳排放就越少。而结构优化设计方法能够得到更好的产品结构布局，被越来越多地用于产品优化改进，具体可分为尺寸优化、形状优化和拓扑优化等。拓扑优化是指在给定的设计空间内，通过添加一些特定的约束条件，在满足工程人员期望的强度、刚度、稳定性等结构性能的前提下，依旧可以减少材料使用、降低制造成本的一种设计方法，对于产品低碳的绿色化生产具有重要意义。

目前，结构拓扑优化大多以结构体积分数、质量、应力、位移、固有频率等作为约束，对结构进行优化改进，但在制造时可能会产生更多的碳排放。如果在结构优化时，将降低碳排放作为约束条件添加入拓扑优化模型中，就能得到满足性能和力学要求，并且碳排放低的拓扑结构，目前专门针对这方面的研究还比较少。

发明内容：

本发明的目的是提供一种低碳约束下的结构拓扑优化方法，所要解决的问题：一是解决如何确定结构碳排放限值问题；二是解决怎么将碳排放限值作为低碳约束添加入结构拓扑优化中的问题。

为了实现上述发明，提出的一种低碳约束下的结构拓扑优化方法，以结构柔度最小化为目标，结构碳排放限值为约束，技术方案实施步骤如下：

步骤一：定义目标产品结构工作条件、设计域，设置边界条件，确定结构受力点以及计算所受载荷。

步骤二：采用变密度法以单元的相对密度为设计变量，以结构的柔度最小为优化目标，以结构碳排放限值为约束条件，建立拓扑优化模型。

步骤三：将初始设计域离散成nelx×nely个有限元单元，通过变密度法引入中间密度单元，设定惩罚因子、过滤半径，初始化设计变量。

步骤四：采用改进的固体各向同性材料惩罚模型(SIMP)对单元的弹性模量进行计算，得到单元的刚度矩阵，组装得到结构的总体刚度矩阵，在设定的边界条件和载荷下对结构进行有限元分析，计算得到单元节点的位移。

步骤五：累加每个有限元单元的碳排放，得到结构总体碳排放，求解目标柔度函数以及结构总体碳排放约束对于设计变量的灵敏度。

步骤六：构造目标模型的拉格朗日函数，使最优解满足库尔塔克(Kuhn-Tucker)条件，即K-T条件，运用灵敏度过滤方法调整目标柔度函数及低碳约束条件的灵敏度，采用优化准则法(OC法)得到设计变量的启发式更新格式，更新设计变量。

步骤七：判断结果是否满足优化收敛条件，若不满足则转至步骤四，并依次进行迭代计算；若满足，则终止拓扑优化进程，获得满足低碳约束的拓扑优化模型。

进一步地，步骤二中所述的拓扑优化模型为：

find：x

min：C(x

s.t.g

KU＝F

0＜x

其中，x

进一步地，步骤二中所述结构碳排放限值约束函数表达式如下：

其中，g

进一步地，步骤四中所述的改进的SIMP插值模型如下：

E

其中，E

进一步地，步骤五中所述的结构碳排放计算方法如下：

式中，x

进一步的，步骤五中所述结构碳排放约束对设计变量的灵敏度如下：

进一步地，步骤五中所述结构柔度目标函数对设计变量的灵敏度如下：

进一步地，步骤六中在低碳约束下构造的满足K-T条件的拉格朗日函数表达式如下：

进一步地，步骤六中求解的拉格朗日函数结果由不动点思想得到的迭代格式为：

进一步地，步骤六中由优化准则法(OC法)可得设计变量的启发式更新格式如下：

其中，k为迭代次数；ξ为阻尼因子，可以保证结果收敛，一般为0.5；m为正向移动极限，可以提高迭代过程的稳定性，一般取0.1～0.3。

进一步的，步骤七中所述判断优化收敛条件为：

其中，ε为收敛精度。

对于上述提出的发明方法，相比较于现有技术，本发明提供的一种低碳约束下的结构拓扑优化方法主要具有以下有益效果：

(1)本发明在结构拓扑优化的过程中，能够达到保证结构刚度的同时，降低结构碳排放，更加符合当代社会提出的“碳达峰”、“碳中和”的“双碳目标”；

(2)本发明提出的碳排放限值计算方法，在面对组成材料不止一种的产品结构时，可以进行函数拓展，在进行拓扑优化时能够得到更好的材料布局，有利于材料节约，因此本发明适用范围广泛；

(3)本发明采用拓扑优化领域广泛采用的优化准则法(OC法)进行优化求解计算，并且使用改进后的SIMP插值模型，迭代速度更快，有利于后期拓展并推广，适用范围广泛。

附图说明

图1是本发明一种低碳约束下的结构拓扑优化方法的流程示意图。

图2是本发明的实例：汽车后防撞梁的初始结构、边界条件以及载荷作用点示意图。

图3是基于本发明获得的优化后汽车后防撞梁的拓扑结构示意图。

图4是基于本发明获得的优化过程中汽车后防撞梁结构柔度随迭代次数变化的示意图。

具体实施方式

下面结合附图及具体实施对本发明做进一步详细说明：

如图1所示，本发明公开了一种低碳约束下的结构拓扑优化方法，包括以下步骤：

步骤一：定义目标产品结构工作条件、设计域，设置边界条件，确定结构受力点以及计算所受载荷。

步骤二：采用变密度法以单元的相对密度为设计变量，以结构的柔度最小为优化目标，以结构碳排放限值为约束条件，建立拓扑优化模型。

所述的低碳约束下的结构拓扑优化模型为：

其中，x

在结构的使用寿命内，通过温室气体减排率与初始碳排放值积累得到结构在使用寿命期间的碳排放总值，根据结构使用寿命期间的目标温室气体减排率设置碳排放最大值，以此作为拓扑优化的低碳约束条件，所述结构碳排放限值约束函数表达式如下：

其中，g

根据IPCC第五次评估第三工作组调查的大量温室气体减缓情景，以及对未来全球人为温室气体减排率与大气CO2浓度之间对应关系的预测，做出假设：碳排放限值与温室气体减排率成正比，并随着时间呈确定性和线性变化。依据假设，提出全球温室气体减排率的计算公式如下：

其中，Y为结构的设定使用年份。

步骤三：将初始设计域离散成nelx×nely个有限元单元，通过变密度法引入中间密度单元，设定惩罚因子、过滤半径，初始化设计变量。

步骤四：采用改进的固体各向同性材料惩罚模型(SIMP)对单元的弹性模量进行计算，得到单元的刚度矩阵，组装得到结构的总体刚度矩阵，在设定的边界条件和载荷下对结构进行有限元分析，计算得到单元节点的位移。

所述的改进的SIMP插值模型如下：

E

其中，E

步骤五：累加每个有限元单元的碳排放，得到结构总体碳排放，求解目标柔度函数以及结构总体碳排放约束对于设计变量的灵敏度。

所述结构碳排放计算公式为：

式中，x

所述拓扑优化模型的灵敏度求解如下：

对于荷载来说，因为荷载为外界施加，与x

结构的碳排放对x

对于结构的位移和刚度约束KU＝F，两边同时对x

结构的总体刚度矩阵是对称的，则有F

对式(8)左乘一个U

结构的总体柔度表达式如式(1)所示，两边同时对x

将式(10)、(11)代入式(12)中，根据改进的SIMP插值模型，得出结构柔度对x

步骤六：构造目标模型的拉格朗日函数，使最优解满足库尔塔克(Kuhn-Tucker)条件，即K-T条件，运用灵敏度过滤方法调整目标柔度函数及低碳约束条件的灵敏度，采用优化准则法(OC法)得到设计变量的启发式更新格式，更新设计变量。

所述拓扑优化模型的拉格朗日函数为：

其中，λ

所述拉格朗日函数在极值点

对于上式，进行分情况讨论如下：

(1)当x

(2)当x

(3)当x

由上述讨论，当材料的相对密度在取值范围内时，即x

将式(6)、式(8)及式(13)代入式(16)中，得到结果如下：

上式基于不动点迭代的思想，转化为显式的结果如下：

令

则有

所述的由优化准则法(OC法)可得设计变量的启发式更新格式如下：

其中，k为迭代次数；ξ为阻尼因子，可以保证结果收敛；m为正向移动极限，可以提高迭代过程的稳定性。

在迭代过程中，采用二分法计算低碳约束的Lagrange乘子λ

其中k为迭代次数；

令

步骤七：判断结果是否满足优化收敛条件，若不满足则转至步骤四，并依次进行迭代计算；若满足，则终止拓扑优化进程，获得满足低碳约束的拓扑优化模型。

所述判断优化收敛条件为：

参照图2-图4，以某款汽车后防撞梁为实例，进而对本发明做进一步说明。防撞梁是车身结构的重要组成部分，它可以直接保护驾驶员和乘客在碰撞中的安全，因此强度要求较高，而且耗材较多，在其制造过程中会产生大量碳排放。汽车后防撞梁的初始结构、边界条件以及载荷作用位置如图2所示，现对其最前面的梁体做结构优化，初始结构设计空间为120cm*12cm，其下侧有两处为固定端，在最上侧中点处受一个竖直向下的集中力，大小设为1kN。汽车后防撞梁整体结构采用铝合金，泊松比为0.3，弹性模量为70GPa，密度为2.7×10

图3为按约束条件优化后的防撞梁结构图，优化后的结构相比较原结构，可以看出，其结构更加迎合整体实际受力和固定约束，材料分布更加合理。

图4为在低碳约束下，迭代过程中防撞梁的结构柔度随迭代次数变化曲线。由图4可知，防撞梁结构柔度从2602.44N·m减小为最终的635.65N·m，降低约75.60％。原实体防撞梁的结构碳排放值为965.78kg CO

本发明在基于SIMP的结构拓扑优化的基础上，通过温室气体减排率设定了结构的碳排放限值，并将其作为低碳约束添加入拓扑优化方法中。通过灵敏度分析，结合优化准则法得到了基于低碳约束下的拓扑优化模型的迭代格式。采用MATLAB软件完成了对汽车后防撞梁实例的拓扑优化改进。

本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。