CN112800608B - 一种考虑多组件布局的多相材料拓扑优化设计方法 - Google Patents
一种考虑多组件布局的多相材料拓扑优化设计方法 Download PDFInfo
- Publication number
- CN112800608B CN112800608B CN202110123949.0A CN202110123949A CN112800608B CN 112800608 B CN112800608 B CN 112800608B CN 202110123949 A CN202110123949 A CN 202110123949A CN 112800608 B CN112800608 B CN 112800608B
- Authority
- CN
- China
- Prior art keywords
- component
- level set
- function
- set function
- matrix
- 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
- 239000000463 material Substances 0.000 title claims abstract description 192
- 238000005457 optimization Methods 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 43
- 239000011159 matrix material Substances 0.000 claims abstract description 85
- 230000014509 gene expression Effects 0.000 claims abstract description 4
- 238000006073 displacement reaction Methods 0.000 claims description 32
- 238000004458 analytical method Methods 0.000 claims description 6
- 230000004069 differentiation Effects 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 6
- 238000010206 sensitivity analysis Methods 0.000 claims description 5
- 239000000758 substrate Substances 0.000 claims description 3
- 239000002131 composite material Substances 0.000 description 10
- 239000011797 cavity material Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 230000009467 reduction Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 235000015842 Hesperis Nutrition 0.000 description 1
- 235000012633 Iberis amara Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Complex Calculations (AREA)
Abstract
本发明属于多相材料拓扑优化设计相关技术领域，其公开了一种考虑多组件布局的多相材料拓扑优化设计方法，该方法用显式水平集函数及隐式水平集函数分别构造多组件形状及基体结构；通过彩色水平集多相材料理论和显式水平集函数得到多个材料以及多组件的拓扑表达，基于Heaviside函数得到单元密度；接着用材料插值方案将不同相材料的弹性模量与多组件的弹性模量插值得到等效弹性模量；其次，基于参数化水平集拓扑优化方法和彩色水平集多相材料理论建立考虑多组件布局的多相材料拓扑优化设计模型，实现多组件布局与多相材料基体的并行优化设计及多相材料中多组件集成系统结构整体性能达到最优。本发明实现了考虑多组件布局的多相材料拓扑优化设计。
Description
技术领域
本发明属于多相材料拓扑优化设计相关技术领域，更具体地，涉及一种考虑多组件布局的多相材料拓扑优化设计方法。
背景技术
复合材料是由具有明显可区分特性的多种成分材料相制成，可以产生多种组合功能。复合多相材料是两种或者两种以上不同材料的不同组合，与单相材料相比，复合多相材料可以实现不同材料的性能。一般地，复合多相材料具有较高的比刚度和比强度，也就说明材料重量轻，而且强度和刚度较高。除此之外，复合多相材料具有较好的抗疲劳性能，其疲劳寿命比一般的金属长数倍。在减振降噪方面，复合多相材料同样优于单相材料，具有较好的减振隔音性能。伴随着3D技术的兴起，多相材料具备优异的性能且成本较低，在工程中的应用变得越来越重要。
结合复合多相材料性能，特别符合航空航天领域中结构设计特点，不仅可以实现飞行器结构的轻量化设计，同时可使整体结构具有较强的刚度和优异的减振性能，且设备零件具有较高的安全性能。多组件集成布局优化设计是航空航天、汽车、军工等领域中经常遇到的技术问题，有时需要将一个或者几个固定形状的几何体嵌入到设计结构中来满足我们需要的特定功能。所以就要同时考虑多组件在结构中的布局设计以及连接组件的基体的结构优化设计。这就要求在优化组件位置分布与角度的同时，优化基体材料的分布，以实现结构性能最优。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提供了一种考虑多组件布局的多相材料拓扑优化设计方法，所述设计方法是一种整体式系统结构优化设计方法，能够同时实现优化多组件布局及多相材料基体的拓扑结构，且基体材料是复合多相材料，具有更加优越的材料性能和组合功能，可以满足航空航天、火箭、导弹等军工制造中的结构性能设计需求。
为实现上述目的，按照本发明的一个方面，提供了一种考虑多组件布局的多相材料拓扑优化设计方法，所述方法包括以下步骤：
(1)采用两套设计变量构造多组件形状和不同材料相的基体结构，其中，采用扩展系数插指隐式水平集函数构建不同材料相的基体结构，同时采用位置坐标和方向角度几何参数建立显式水平集函数，并采用显式水平集函数分别构造各个组件形状；接着，采用max函数将所有组件的水平集函数集成到一个水平集函数中，以得到集成有所有组件的多组件水平集函数；
(2)基于彩色水平集多相理论和多组件水平集函数构造基体材料、多组件及空洞的材料相，并通过组件和基体的水平集函数的不同组合得到不同材料相的拓扑模型；接着，基于Heaviside函数获得不同材料相连续变化的单元密度；
(3)基于不同材料相连续变化的单元密度，将基体多相材料的弹性模量与嵌入设计域中的各个组件的弹性模量相插值，以得到多组件结构系统的等效弹性模量；
(4)基于参数化水平集和彩色水平集多相材料理论建立考虑多组件集成系统布局优化的多相材料最小柔度拓扑优化模型，并基于得到的等效弹性模量在结构设计域中通过有限元分析求解整体结构的位移场，根据得到的位移场计算多相材料最小柔度拓扑优化模型的目标函数；接着，基于形状微分与伴随变量法对基体和组件的设计变量进行灵敏度分析，并采用MMA移动渐近线算法更新全局设计变量，继而确定各个组件的最优布局和不同材料相基体的最优分布。
进一步地，基体结构的隐式水平集函数为：
其中，x＝x1,x2,...,xN表示所有插值节点坐标，即水平集节点；N表示节点总数；t为伪时间变量；αl,n表示在第l个水平集函数节点n处的扩展系数；Φl(l＝1,2)表示基体结构材料1和基体结构材料2的水平集函数，该水平集函数由高斯径向基函数φl,n(x)插值；φl,n(x)表示第l个高斯径向基函数。
进一步地，
其中，c是形状参数，等于水平集网格面积或者体积的倒数；xn表示水平集函数的第n个节点的坐标；||x-xn||是用于计算当前采样点x到xn节点距离的欧几里得范数。
进一步地，组件的水平集函数为：
其中，Φci(x,y)表示第i个组件的水平集函数；xci表示第i个组件的横坐标；yci表示第i个组件的纵坐标；rci表示第i个圆形组件的半径；(x,y)为采样点坐标。
进一步地，多组件水平集函数为：
Φc(x,y)＝max(Φc1,Φc2,...,Φcnc)
其中，nc为所有组件的数量。
进一步地，不同材料相的拓扑模型如下：
其中，Φ1(x)表示基体材料1的水平集函数；Φ2(x)表示基体材料2的水平集函数。
进一步地，材料相的特征函数χr为：
其中，Φ1表示基体材料1的水平集函数；Φ2表示基体材料2的水平集函数；r＝1,2,3,4；ΦC是集成所有组件的水平集函数；H表示Heaviside函数：
进一步地，多相材料最小柔度拓扑优化模型的表达式为：
Find:αl,n＝[αl,1,αl,2,...,αl,N]；SCi＝{xC1,yC1,θC1,xC2,yC2,θC2,...,xCnc,yCnc,θCnc}
其中，α1,n表示基体材料1的设计变量，是高斯径向基函数插值时的扩展系数，仅与时间变量相关；α1,n,max和α1,n,min分别表示基体材料1设计变量的上限和下限；α2,n表示基体材料2的设计变量，是高斯径向基函数插值时的扩展系数，仅与时间变量相关；α2,n,max和α2,n,min分别表示基体材料2设计变量的上限和下限；SCi为组件的设计变量，即为xCi,yCi,θCi，其中xCi表示第i个组件的横坐标，yCi表示第i个组件的纵坐标，θCi表示第i个组件与平轴之间逆时针旋转的角度；nc是组件的总数量；xCi,max和xCi,min分别是组件质心横坐标移动的上限和下限；yCi,max和yCi,min分别是组件质心纵坐标移动的上限和下限；ε为应变场；T表示矩阵的转置；u表示结构场的位移；v表示在动力学上允许的位移空间U中的一个虚拟位移；a表示双线性能量式；l表示单线性负载式；E表示多组件集成系统的等效弹性模量；Φ1表示基体材料1的水平集函数；Φ2表示基体材料2的水平集函数；ΦC是集成所有组件的水平集函数；dΩ为结构设计域的积分算子；H表示Heaviside函数，用于表征结构形式的特征函数；g1表示基体材料1的体积约束函数；g2表示基体材料2的体积约束函数；fv1是基体材料1约束的体积分数；fv2是基体材料2约束的体积分数；V0表示结构设计域体积。
进一步地，步骤(4)包括以下子步骤：
(4.1)初始化设计参数和基体不同材料的水平集函数以及组件水平集函数；
(4.2)基于求得的多组件集成系统的等效弹性模量E计算结构单元刚度Ke，并组装得到整体刚度矩阵K；接着进行有限元分析求解结构位移场u；
(4.3)基于步骤(4.2)中求解的结构位移场u计算考虑多组件集成系统布局优化的多相材料最小柔度拓扑优化模型的目标函数J：
其中，ε为应变场；T表示矩阵的转置；u表示结构场的位移；v表示在动力学上允许的位移空间U中的一个虚拟位移；E是多组件集成系统的等效弹性模量；dΩ为结构设计域的积分算子；
(4.4)基于形状微分与伴随变量法，就模型的目标函数与约束函数对基体设计变量进行灵敏度分析，并采用MMA移动渐近线算法更新全局设计变量，判断模型是否满足收敛条件；若否，则返回到步骤(4.2)，若是，则输出基体材料的最优拓扑结构和多组件的最优布局。
进一步地，多相材料最小柔度拓扑优化模型的目标函数与约束函数针对组件设计变量的一阶微分为：
其中，u表示结构场的位移；Φ1表示基体材料1的水平集函数；Φ2表示基体材料2的水平集函数；Φci(x,y)表示第i个组件的水平集函数；ΦC是集成所有组件的水平集函数；gl表示基体材料l的体积约束函数，l＝1,2SCi为组件的设计变量，即为xCi,yCi,θCi，其中xCi表示第i个组件的横坐标，yCi表示第i个组件的纵坐标，θCi表示第i个组件与平轴之间逆时针旋转的角度；表示集成所有多组件的水平集函数对组件设计变量的一阶微分。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，本发明提供的考虑多组件布局的多相材料拓扑优化设计方法主要具有以下有益效果：
1.本发明提供的考虑多组件布局的多相材料拓扑优化设计方法，与现有的技术方案相比，基体材料是复合多相材料，具有更加优越的材料性能和组合功能，可以满足航空航天、火箭、导弹等军工制造中的结构性能设计需求。
2.所述优化设计方法开发了两个及两个以上的材料相复合结构的通用并行设计方式，可以实现多相材料的基体与不同材料的组件同时优化设计，在获得组件最优布局的同时，得到基体结构不同相材料的最佳分布。
3.所述优化设计方法采用显式和隐式水平集函数分别表示多组件和基体结构，采用两套设计变量，即插值基体隐式水平集函数的扩展系数和构造组件显式水平集函数的位置坐标及方向角度几何参数，极大地降低了设计变量数量。
4.所述优化设计方法可以方便地处理不同相材料的多个体积约束，精确控制多个组件的几何形状，同时可以获得清晰、独特的不同材料相的拓扑结构边界。
附图说明
图1是本发明提供的考虑多组件布局的多相材料拓扑优化设计方法的流程示意图；
图2是本发明实施例的初始结构设计域的示意图；
图3是图2中的初始结构设计域嵌入两个圆形组件后的示意图；
图4是基体材料1、基体材料2及空洞材料中嵌入两个圆形组件的初始布局示意图；
图5是拓扑优结构柔度J＝522.0177时的拓扑结构示意图；
图6是基体材料1的体积分数迭代曲线示意图；
图7是基体材料2的体积分数迭代曲线示意图；
图8是目标函数的迭代曲线示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
请参阅图1，本发明提供的考虑多组件布局的多相材料拓扑优化设计方法，所述多相材料拓扑优化设计方法主要包括以下步骤：
(1)采用两套设计变量构造多组件形状和不同材料相的基体结构，其中，采用扩展系数插指隐式水平集函数构建不同材料相的基体结构，同时采用位置坐标和方向角度几何参数建立显式水平集函数，并采用显式水平集函数分别构造各个组件形状；接着，采用max函数将所有组件的水平集函数集成到一个水平集函数中，以得到集成有所有组件的多组件水平集函数。
具体地，通过隐式水平集函数构建不同材料相的基体结构，以插值隐式水平集函数的扩展系数为基体设计变量，在N个固定水平集结点处的基体结构隐式水平集函数为：
其中，x＝x1,x2,...,xN表示所有插值节点坐标，即水平集节点；N表示节点总数；t为伪时间变量；αl,n表示在第l个水平集函数节点n处的扩展系数；Φl(l＝1,2)表示基体结构材料1和基体结构材料2的水平集函数，该水平集函数由高斯径向基函数φl,n(x)插值；φl,n(x)表示第l个高斯径向基函数，公式为：
其中，c是形状参数，等于水平集网格面积或者体积的倒数；xn表示水平集函数的第n个节点的坐标；||x-xn||是用于计算当前采样点x到xn节点距离的欧几里得范数。
以组件位置坐标和方向角度为设计变量构造显式水平集函数，使用max函数将所有组件集成到一个水平集函数中。以下以圆形组件为例，组件的水平集函数为：
其中，Φci(x,y)表示第i个组件的水平集函数；xci表示第i个组件的横坐标；yci表示第i个组件的纵坐标；rci表示第i个圆形组件的半径；(x,y)为采样点坐标。
矩形组件的水平集函数为：
其中，Φci(x,y)表示第i个组件的水平集函数；Q是一个相对较大的偶数，这里取值Q＝6；xci表示第i个组件的横坐标；yci表示第i个组件的纵坐标；θi表示第i个组件与平轴之间逆时针旋转的角度；Li为第i个矩形组件长度的一半；ti为第i个矩形组件宽度的一半；x'为嵌入组件的局部横坐标；y'为嵌入组件的局部纵坐标。
接着，使用max函数将所有函数集成到一个水平集函数Φc(x,y)中，此水平集函数集成了所有组件的几何参数：
Φc(x,y)＝max(Φc1,Φc2,...,Φcnc)
其中，nc为所有组件的数量。
(2)基于彩色水平集多相理论和多组件水平集函数构造基体材料、多组件及空洞的材料相，并通过组件和基体的水平集函数的不同组合得到不同材料相的拓扑模型；接着，基于Heaviside函数获得不同材料相连续变化的单元密度。
具体地，通过彩色水平集多相材料理论和显式水平集函数构造基体材料1、基体材料2等基体材料和多组件、空洞等材料相。本实施方式中以基体三相材料为例说明，但本发明所提供的方法不限于三相材料。通过组件和基体的水平集函数的不同组合，分别得到不同材料相的拓扑模型如下：
其中，Φ1(x)表示基体材料1的水平集函数，Φ2(x)表示基体材料2的水平集函数。
接着；利用Heaviside函数获得不同材料相连续变化的单元密度，不同材料相的特征函数χr(r＝1,2,3,4)为：
其中，Φ1表示基体材料1的水平集函数，Φ2表示基体材料2的水平集函数；H表示Heaviside函数：
(3)基于不同材料相连续变化的单元密度，将基体多相材料的弹性模量与嵌入设计域中的各个组件的弹性模量相插值，以得到多组件结构系统的等效弹性模量。
具体地，将基体各相材料的弹性模量与各个组件的弹性模量相插值得到多组件集成系统的等效弹性模量E:
特别地，当嵌入的所有组件材料属性相同时：
其中，Er(r＝1,2,3,4)表示第r种材料相的弹性模量；Eci表示第i个组件材料的弹性模量；Ec表示当所有组件材料相同时的弹性模量；H表示Heaviside函数；ΦC是集成所有组件的水平集函数。
(4)基于参数化水平集和彩色水平集多相材料理论建立考虑多组件集成系统布局优化的多相材料最小柔度拓扑优化模型，并基于得到的等效弹性模量在结构设计域中通过有限元分析求解整体结构的位移场，根据得到的位移场计算多相材料最小柔度拓扑优化模型的目标函数；接着，基于形状微分与伴随变量法对基体和组件的设计变量进行灵敏度分析，并采用MMA移动渐近线算法更新全局设计变量，继而确定各个组件的最优布局和不同材料相基体的最优分布。
具体地，所述的基于参数化水平集和彩色水平集多相材料理论建立的考虑多组件集成系统布局优化的多相材料最小柔度拓扑优化模型的表达式为：
Find:αl,n＝[αl,1,αl,2,...,αl,N]；SCi＝{xC1,yC1,θC1,xC2,yC2,θC2,...,xCnc,yCnc,θCnc}
其中，α1,n表示基体材料1的设计变量，是高斯径向基函数插值时的扩展系数，仅与时间变量相关；α1,n,max和α1,n,min分别表示基体材料1设计变量的上限和下限；α2,n表示基体材料2的设计变量，是高斯径向基函数插值时的扩展系数，仅与时间变量相关；α2,n,max和α2,n,min分别表示基体材料2设计变量的上限和下限；SCi为组件的设计变量，即为xCi,yCi,θCi，其中xCi表示第i个组件的横坐标，yCi表示第i个组件的纵坐标，θCi表示第i个组件与平轴之间逆时针旋转的角度；nc是组件的总数量；xCi,max和xCi,min分别是组件质心横坐标移动的上限和下限；yCi,max和yCi,min分别是组件质心纵坐标移动的上限和下限；ε为应变场；T表示矩阵的转置；u表示结构场的位移；v表示在动力学上允许的位移空间U中的一个虚拟位移；a表示双线性能量式；l表示单线性负载式；E表示多组件集成系统的等效弹性模量；Φ1表示基体材料1的水平集函数；Φ2表示基体材料2的水平集函数；ΦC是集成所有组件的水平集函数；dΩ为结构设计域的积分算子；H表示Heaviside函数，用于表征结构形式的特征函数；g1表示基体材料1的体积约束函数；g2表示基体材料2的体积约束函数；fv1是基体材料1约束的体积分数；fv2是基体材料2约束的体积分数；V0表示结构设计域体积。
基于虚功原理，针对有限单元平衡方程进行计算，对应的弱变分形式如下：
a(u,v,Φ1,Φ2,ΦC)＝∫DεT(u)Eε(v)dΩ
l(u,v,Φ1,Φ2,ΦC)＝
∫DpvH(Φ1)H(Φ2)(1-H(ΦC))dΩ+∫DpvH(Φ1)(1-H(Φ2))(1-H(ΦC))dΩ+
∫DpvH(ΦC)dΩ+∫Dτvδ(Φ1)H(Φ2)(1-H(ΦC))▽Φ1|dΩ-
∫Dτvδ(Φ2)H(Φ1)(1-H(ΦC))dΩ+∫Dτvδ(ΦC)dΩ
其中，a表示双线性能量式；l表示单线性负载式；Φ1表示基体材料1的水平集函数；Φ2表示基体材料2的水平集函数；ΦC是集成所有组件的水平集函数；dΩ为结构设计域的积分算子；H表示Heaviside函数，用于表征结构形式的特征函数；ε为应变场；T表示矩阵的转置；u表示结构场的位移；v表示在动力学上允许的位移空间U中的一个虚拟位移；τ表示应用在边界的部分边界上的牵引力；P表示结构设计域的体积力；δ表示Dirac函数，为Heaviside函数的一阶微分；▽表示差分算子。
构建考虑多组件集成系统布局优化的多相材料最小柔度拓扑优化模型具体包括以下步骤：
(4.1)初始化设计参数和基体不同材料的水平集函数以及组件水平集函数。具体还包括基体和组件材料属性参数、设计域长宽；
(4.2)基于求得的多组件集成系统的等效弹性模量E计算结构单元刚度Ke，并组装得到整体刚度矩阵K；接着进行有限元分析求解结构位移场u；
(4.3)基于步骤(4.2)中求解的结构位移场u计算考虑多组件集成系统布局优化的多相材料最小柔度拓扑优化模型的目标函数J：
其中，ε为应变场；T表示矩阵的转置；u表示结构场的位移；v表示在动力学上允许的位移空间U中的一个虚拟位移；E是多组件集成系统的等效弹性模量；dΩ为结构设计域的积分算子。
(4.4)基于形状微分与伴随变量法求解目标函数与约束函数针对基体设计变量进行灵敏度分析，并采用MMA移动渐近线算法更新全局设计变量，判断模型是否满足收敛条件，若否，则返回到步骤(4.2)，若是，则输出基体材料的最优拓扑结构和多组件的最优布局。
具体地，根据链式求导法则计算目标函数与约束函数针对组件设计变量的一阶微分，计算如下：
目标函数与约束函数针对基体设计变量αl(l＝1,2)的一阶微分：
其中：
其中，αl,n表示基体材料l(l＝1,2)的设计变量，是高斯径向基函数插值时的扩展系数，仅与时间变量相关；ε为应变场；T表示矩阵的转置；u表示结构场的位移；v表示在动力学上允许的位移空间U中的一个虚拟位移；E是多组件集成系统的等效弹性模量；Φ1表示基体材料1的水平集函数；Φ2表示基体材料2的水平集函数；ΦC是集成所有组件的水平集函数；Er表示第r种材料相的弹性模量，r＝1,2,3,4；τ表示应用在边界的部分边界上的牵引力；P表示结构设计域的体积力；H表示Heaviside函数；δ表示Dirac函数，为Heaviside函数的一阶微分；χr表示第r种材料相的特征函数；φl,n(x)表示第l个高斯径向基函数；g1表示基体材料1的体积约束函数；g2表示基体材料2的体积约束函数；dΩ为结构设计域的积分算子；dV表示体积积分。
目标函数与约束函数针对组件设计变量的一阶微分为：
其中，u表示结构场的位移；Φ1表示基体材料1的水平集函数；Φ2表示基体材料2的水平集函数；Φci(x,y)表示第i个组件的水平集函数；ΦC是集成所有组件的水平集函数；gl表示基体材料l的体积约束函数，l＝1,2SCi为组件的设计变量，即为xCi,yCi,θCi，其中xCi表示第i个组件的横坐标，yCi表示第i个组件的纵坐标，θCi表示第i个组件与平轴之间逆时针旋转的角度；表示集成所有多组件的水平集函数对组件设计变量的一阶微分。
根据链式规则可以得如下：
使用伴随法可以得到目标函数和体积约束函数对集成所有组件的水平集函数的一阶微分和如下：
其中，ε为应变场；T表示矩阵的转置；u表示结构场的位移；E是多组件集成系统的等效弹性模量；Φ1表示基体材料1的水平集函数；Φ2表示基体材料2的水平集函数；ΦC是集成所有组件的水平集函数；χr表示第r种材料相的特征函数，r＝1,2,3,4；gl(l＝1,2)表示基体材料l的体积约束函数，表示多组件集成系统的等效弹性模量对集成所有组件的水平集函数的一阶微分；表示第r种材料相的特征函数对集成所有组件的水平集函数的一阶微分；dΩ为结构设计域的积分算子；dV表示体积积分。
结合步骤(2)得到的多组件集成系统的等效弹性模量E可以得到：
其中，Er表示第r种材料相的弹性模量；χr表示第r种材料相的特征函数，ΦC是集成所有组件的水平集函数。
进一步地，可以计算特征函数的形状导数为：
其中，H表示Heaviside函数；δ表示Dirac函数，为Heaviside函数的一阶微分；Φ1表示基体材料1的水平集函数；Φ2表示基体材料2的水平集函数；ΦC是集成所有组件的水平集函数。
注意到，当组件材料属性不同时可以直接求
其中，Φci(x,y)表示第i个组件的水平集函数；δ表示Dirac函数，为Heaviside函数的一阶微分。
上面已完成和的求解过程，可以由Heaviside函数求导法则轻松获得。
更新设计变量并判断优化模型是否满足收敛条件，若否，则返回到步骤(4.2)；若是，则输出基体材料的最优拓扑结构和多组件的最优布局。本步骤关键在于的求解，由于ΦC＝max(ΦC1,ΦC2,...,ΦCnc)中包含max算子，所以在求解微分时，可以采用K-S方程来近似max算子，从而可以计算的解析求解。K-S方程方程如下：
其中，当A取一个较大的正数时(例如A＝50)，γ≈max(γ1,γ2,...,γn)，所以可以由K-S方程近似求得。
构造多组件的显式水平集函数对设计变量的一阶微分可根据显式水平集函数ΦCi对几何参数Si(xi,yi,θi)进行求导获得，所以到此目标函数和体积约束函数对多组件设计变量的敏度分析完毕。
实施例
本发明实施例提供了一种考虑多组件布局的多相材料拓扑优化设计方法，且本实施例中以基体三相材料为例说明本方法，三相基体材料即为基体材料1、基体材料2和空洞材料。实施例中，2D结构设计域如图2所示，结构尺寸定义为120×60，有限单元网格划分为120×60，定义基体材料1的弹性模量为E1＝1，基体材料2的弹性模量为E2＝5，空洞的弹性模量定义为E3＝0.001，基体所有材料的泊松比为μ＝0.3。定义基体材料1的体积分数为fv1＝0.25，基体材料2的体积分数为fv2＝0.25，空洞的体积分数为0.5。在本例中嵌入相同材料属性的多组件，所有组件的弹性模量为EC＝10。如图2所示定义设计域的左边固定约束，右边界中点施加有F＝-5N的集中载荷。将两个圆形组件嵌入到结构设计域中的初始布局如图3所示，圆形组件的半径r＝5，组件的初始质心位置为(54.5，16.3)、(54.5，43.7)。图4展示了基体材料1、基体材料2和空洞材料中嵌入两个圆形组件初始布局示意图。在优化中，设置优化过程最大迭代步数为200步，设置迭代终止的条件为目标函数迭代前后的目标函数值的差小于等于10-7时，终止迭代。
经过106步迭代，考虑多组件布局的多相材料拓扑优化设计的最优拓扑结构如图5所示，其柔度值为J＝522.0177，可以看到圆形组件移动到结构的关键承载路径上，基体材料1和基体材料2的边界区分清晰、独特，并且整体结构边界光滑、连接性较好。如图6是基体材料1的实际体积分数与约束体积分数差值的迭代曲线，可以看到开始时由于基体材料1的实际体积分数远小于约束值，所以优化前期波动较大，差值为负值，迭代20步后慢慢趋于稳定。图7是基体材料2的实际体积分数与约束体积分数差值的迭代曲线，可以看到开始时由于实际体积分数大于约束值，所以开始时出现剧烈下降，随后趋于稳定，可以看到该实施例在优化中可以方便地处理不同规模的多个体积约束。如图8所示是考虑多组件布局的多相材料拓扑优化设计模型的目标函数迭代曲线，在开始时由于全局体积分数没有达到材料体积的约束值，所以在前几步出现波动，当基体材料1和基体材料2的体积分数达到约束值时，目标函数开始稳定变化，在迭代10步以后目标函数逐渐收敛。结果表明所提出的方法具有较高的优化效率。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种考虑多组件布局的多相材料拓扑优化设计方法，其特征在于，该方法包括以下步骤：
(1)采用两套设计变量构造多组件形状和不同材料相的基体结构，其中，采用扩展系数插指隐式水平集函数构建不同材料相的基体结构，同时采用位置坐标和方向角度几何参数建立显式水平集函数，并采用显式水平集函数分别构造各个组件形状；接着，采用max函数将所有组件的水平集函数集成到一个水平集函数中，以得到集成有所有组件的多组件水平集函数；
(2)基于彩色水平集多相理论和多组件水平集函数构造基体材料、多组件及空洞的材料相，并通过组件和基体的水平集函数的不同组合得到不同材料相的拓扑模型；接着，基于Heaviside函数获得不同材料相连续变化的单元密度；
(3)基于不同材料相连续变化的单元密度，将基体多相材料的弹性模量与嵌入设计域中的各个组件的弹性模量相插值，以得到多组件结构系统的等效弹性模量；
(4)基于参数化水平集和彩色水平集多相材料理论建立考虑多组件集成系统布局优化的多相材料最小柔度拓扑优化模型，并基于得到的等效弹性模量在结构设计域中通过有限元分析求解整体结构的位移场，根据得到的位移场计算多相材料最小柔度拓扑优化模型的目标函数；接着，基于形状微分与伴随变量法对基体和组件的设计变量进行灵敏度分析，并采用MMA移动渐近线算法更新全局设计变量，继而确定各个组件的最优布局和不同材料相基体的最优分布。
2.如权利要求1所述的考虑多组件布局的多相材料拓扑优化设计方法，其特征在于：基体结构的隐式水平集函数为：
其中，x＝x1,x2,...,xN表示所有插值节点坐标，即水平集节点；N表示节点总数；t为伪时间变量；αl,n表示在第l个水平集函数节点n处的扩展系数；Φl(l＝1,2)表示基体结构材料1和基体结构材料2的水平集函数，该水平集函数由高斯径向基函数φl,n(x)插值；φl,n(x)表示第l个高斯径向基函数。
3.如权利要求2所述的考虑多组件布局的多相材料拓扑优化设计方法，其特征在于：
其中，c是形状参数，等于水平集网格面积或者体积的倒数；xn表示水平集函数的第n个节点的坐标；||x-xn||是用于计算当前采样点x到xn节点距离的欧几里得范数。
4.如权利要求2所述的考虑多组件布局的多相材料拓扑优化设计方法，其特征在于：组件的水平集函数为：
其中，Φci(x,y)表示第i个组件的水平集函数；xci表示第i个组件的横坐标；yci表示第i个组件的纵坐标；rci表示第i个圆形组件的半径；(x,y)为采样点坐标。
5.如权利要求4所述的考虑多组件布局的多相材料拓扑优化设计方法，其特征在于：多组件水平集函数为：
Φc(x,y)＝max(Φc1,Φc2,...,Φcnc)
其中，nc为所有组件的数量。
6.如权利要求4所述的考虑多组件布局的多相材料拓扑优化设计方法，其特征在于：不同材料相的拓扑模型如下：
其中，Φ1(x)表示基体材料1的水平集函数；Φ2(x)表示基体材料2的水平集函数。
7.如权利要求6所述的考虑多组件布局的多相材料拓扑优化设计方法，其特征在于：材料相的特征函数χr为：
其中，Φ1表示基体材料1的水平集函数；Φ2表示基体材料2的水平集函数；r＝1,2,3,4；ΦC是集成所有组件的水平集函数；H表示Heaviside函数：
8.如权利要求7所述的考虑多组件布局的多相材料拓扑优化设计方法，其特征在于：多相材料最小柔度拓扑优化模型的表达式为：
Find:αl,n＝[αl,1,αl,2,...,αl,N]；SCi＝{xC1,yC1,θC1,xC2,yC2,θC2,...,xCnc,yCnc,θCnc}
Min:
其中，α1,n表示基体材料1的设计变量，是高斯径向基函数插值时的扩展系数，仅与时间变量相关；α1,n,max和α1,n,min分别表示基体材料1设计变量的上限和下限；α2,n表示基体材料2的设计变量，是高斯径向基函数插值时的扩展系数，仅与时间变量相关；α2,n,max和α2,n,min分别表示基体材料2设计变量的上限和下限；SCi为组件的设计变量，即为xCi,yCi,θCi，其中xCi表示第i个组件的横坐标，yCi表示第i个组件的纵坐标，θCi表示第i个组件与平轴之间逆时针旋转的角度；nc是组件的总数量；xCi,max和xCi,min分别是组件质心横坐标移动的上限和下限；yCi,max和yCi,min分别是组件质心纵坐标移动的上限和下限；ε为应变场；T表示矩阵的转置；u表示结构场的位移；v表示在动力学上允许的位移空间U中的一个虚拟位移；a表示双线性能量式；l表示单线性负载式；E表示多组件集成系统的等效弹性模量；Φ1表示基体材料1的水平集函数；Φ2表示基体材料2的水平集函数；ΦC是集成所有组件的水平集函数；dΩ为结构设计域的积分算子；H表示Heaviside函数，用于表征结构形式的特征函数；g1表示基体材料1的体积约束函数；g2表示基体材料2的体积约束函数；fv1是基体材料1约束的体积分数；fv2是基体材料2约束的体积分数；V0表示结构设计域体积。
9.如权利要求1-8任一项所述的考虑多组件布局的多相材料拓扑优化设计方法，其特征在于：步骤(4)包括以下子步骤：
(4.1)初始化设计参数和基体不同材料的水平集函数以及组件水平集函数；
(4.2)基于求得的多组件集成系统的等效弹性模量E计算结构单元刚度Ke，并组装得到整体刚度矩阵K；接着进行有限元分析求解结构位移场u；
(4.3)基于步骤(4.2)中求解的结构位移场u计算考虑多组件集成系统布局优化的多相材料最小柔度拓扑优化模型的目标函数J：
其中，ε为应变场；T表示矩阵的转置；u表示结构场的位移；v表示在动力学上允许的位移空间U中的一个虚拟位移；E是多组件集成系统的等效弹性模量；dΩ为结构设计域的积分算子；
(4.4)基于形状微分与伴随变量法，就模型的目标函数与约束函数对基体设计变量进行灵敏度分析，并采用MMA移动渐近线算法更新全局设计变量，判断模型是否满足收敛条件；若否，则返回到步骤(4.2)，若是，则输出基体材料的最优拓扑结构和多组件的最优布局。
10.如权利要求9所述的考虑多组件布局的多相材料拓扑优化设计方法，其特征在于：多相材料最小柔度拓扑优化模型的目标函数与约束函数针对组件设计变量的一阶微分为：
其中，u表示结构场的位移；Φ1表示基体材料1的水平集函数；Φ2表示基体材料2的水平集函数；Φci(x,y)表示第i个组件的水平集函数；ΦC是集成所有组件的水平集函数；gl表示基体材料l的体积约束函数，l＝1,2 SCi为组件的设计变量，即为xCi,yCi,θCi，其中xCi表示第i个组件的横坐标，yCi表示第i个组件的纵坐标，θCi表示第i个组件与平轴之间逆时针旋转的角度；表示集成所有多组件的水平集函数对组件设计变量的一阶微分。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110123949.0A CN112800608B (zh) | 2021-01-29 | 2021-01-29 | 一种考虑多组件布局的多相材料拓扑优化设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110123949.0A CN112800608B (zh) | 2021-01-29 | 2021-01-29 | 一种考虑多组件布局的多相材料拓扑优化设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112800608A CN112800608A (zh) | 2021-05-14 |
| CN112800608B true CN112800608B (zh) | 2024-04-19 |
Family
ID=75812742
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110123949.0A Active CN112800608B (zh) | 2021-01-29 | 2021-01-29 | 一种考虑多组件布局的多相材料拓扑优化设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112800608B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114330053B (zh) * | 2021-12-16 | 2022-12-02 | 北京应用物理与计算数学研究所 | 模拟金属表面氢化点腐蚀形貌演化的方法及装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111460622A (zh) * | 2020-03-06 | 2020-07-28 | 华中科技大学 | 一种基于参数化彩色水平集的多组件布局拓扑优化方法 |
| CN111950149A (zh) * | 2020-08-13 | 2020-11-17 | 北京航空航天大学 | 基于参数化水平集法的连续体结构非概率拓扑优化方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106096158B (zh) * | 2016-06-16 | 2019-04-09 | 华南理工大学 | 一种柔性铰链的拓扑优化设计方法 |
-
2021
- 2021-01-29 CN CN202110123949.0A patent/CN112800608B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111460622A (zh) * | 2020-03-06 | 2020-07-28 | 华中科技大学 | 一种基于参数化彩色水平集的多组件布局拓扑优化方法 |
| CN111950149A (zh) * | 2020-08-13 | 2020-11-17 | 北京航空航天大学 | 基于参数化水平集法的连续体结构非概率拓扑优化方法 |
Non-Patent Citations (1)
| Title |
|---|
| 基于水平集方法的连续体结构拓扑优化;任晓辉;封建湖;;建筑科学与工程学报;20070330(01);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112800608A (zh) | 2021-05-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107844676B (zh) | 一种基于多性能约束的结构拓扑优化设计方法 | |
| CN111737835B (zh) | 基于三周期极小曲面的三维多孔散热结构的设计与优化方法 | |
| CN110795873B (zh) | 一种考虑尺寸控制的跨尺度拓扑优化方法 | |
| CN112765856B (zh) | 一种功能梯度多孔结构拓扑优化的混合水平集方法 | |
| CN109670200A (zh) | 一种等几何材料密度场结构拓扑优化方法 | |
| CN110222356B (zh) | 综合考虑稳定性与振动特性的板/壳结构轻量化拓扑优化设计方法 | |
| WO2022000132A1 (zh) | 基于三周期极小曲面的三维多孔散热结构的设计与优化方法 | |
| CN111460622B (zh) | 一种基于参数化彩色水平集的多组件布局拓扑优化方法 | |
| CN111950149A (zh) | 基于参数化水平集法的连续体结构非概率拓扑优化方法 | |
| CN113821887B (zh) | 基于无网格efgm和plsm的各向异性结构热力耦合拓扑优化方法 | |
| CN109657284A (zh) | 一种面向超材料的等几何拓扑优化方法 | |
| CN112800608B (zh) | 一种考虑多组件布局的多相材料拓扑优化设计方法 | |
| CN113204906B (zh) | 一种考虑结构稳定性的多相材料拓扑优化设计方法和系统 | |
| CN109002614B (zh) | 一种稳定成孔的改进水平集拓扑优化方法 | |
| CN115295097A (zh) | 一种考虑增材制造的纤维增强复合材料多材料多尺度变刚度优化设计方法 | |
| CN117373579B (zh) | 一种时域动载荷作用下多微结构多尺度并行拓扑优化方法及系统 | |
| CN117556672A (zh) | 面向智能制造中结构应力最小化设计的高效拓扑优化方法 | |
| CN112084697A (zh) | 一种可实现变形能传递解耦的柔顺机构设计方法 | |
| CN116362079A (zh) | 一种基于新型插值模型的多材料结构拓扑优化方法 | |
| CN111127491A (zh) | 一种针对蜂窝结构拓扑优化的多变量水平分割方法及设备 | |
| CN111651923B (zh) | 一种多自由度被动隔振系统优化设计方法 | |
| CN110941879A (zh) | 一种考虑连接性的三明治阻尼复合结构拓扑优化设计方法 | |
| CN116341160B (zh) | 针对筒状耐压结构内向增稳的拓扑优化方法及系统 | |
| Kim et al. | Parametric Aeroservoelastic-Gust Reduced-Order Models with State Consistence Enforcement | |
| CN113722947B (zh) | 一种基于多目标粒子群算法的重型机床基础设计方法 |
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 |