CN113191044A - 一种单材料多孔结构的拓扑优化设计方法 - Google Patents
一种单材料多孔结构的拓扑优化设计方法 Download PDFInfo
- Publication number
- CN113191044A CN113191044A CN202110396411.7A CN202110396411A CN113191044A CN 113191044 A CN113191044 A CN 113191044A CN 202110396411 A CN202110396411 A CN 202110396411A CN 113191044 A CN113191044 A CN 113191044A
- Authority
- CN
- China
- Prior art keywords
- porous structure
- function
- density
- level set
- material porous
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 80
- 239000000463 material Substances 0.000 title claims abstract description 68
- 238000005457 optimization Methods 0.000 title claims abstract description 55
- 239000011159 matrix material Substances 0.000 claims abstract description 38
- 238000006073 displacement reaction Methods 0.000 claims abstract description 23
- 230000000903 blocking effect Effects 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 230000035945 sensitivity Effects 0.000 claims abstract description 4
- 238000001914 filtration Methods 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 7
- 230000010354 integration Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 230000006870 function Effects 0.000 description 74
- 238000004519 manufacturing process Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 3
- 239000000945 filler Substances 0.000 description 3
- 235000015842 Hesperis Nutrition 0.000 description 2
- 235000012633 Iberis amara Nutrition 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000010146 3D printing Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 210000003323 beak Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000000265 homogenisation Methods 0.000 description 1
- 238000011423 initialization method Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Complex Calculations (AREA)
Abstract
本发明属于单材料结构拓扑优化设计相关技术领域,其公开了一种单材料多孔结构的拓扑优化设计方法,包括以下步骤:(1)构建单材料多孔结构,进而通过扩展系数矩阵分块法进行插值以计算得到初始扩展系数;(2)计算单材料多孔结构中连续变化的单元伪密度、结构总体积及每个单元的局部伪密度,并求取局部伪密度最大值;(3)对单材料多孔结构的弹性模量进行插值以得到等效单元弹性模量;(4)在求解得到整体结构的位移场,进而计算单材料多孔结构最小柔度拓扑优化模型的目标函数;接着,计算目标函数、结构总体积及局部伪密度对设计变量的灵敏度,继而得到单材料多孔结构的结构拓扑。本发明提高了优化求解效率,降低了优化复杂性及计算量。
Description
技术领域
本发明属于单材料结构拓扑优化设计相关技术领域,更具体地,涉及一种单材料多孔结构的拓扑优化设计方法。
背景技术
增材制造(AM,也被称为3D打印)能够制造出结合复杂度远远超过传统制造技术的组件,这使得一些复杂的轻量化制造成为可能。自然界中有许多常见的轻型壳多孔结构,如植物茎、鸟喙和人类骨骼。多孔结构是指结构由一个固体壳与多孔的内部所组成,而不完全为固体。多孔结构因其具有高的比刚度、良好吸声性和减震性受到了广泛研究和工程应用。
目前,多孔结构设计可分为单尺度设计和多尺度设计两大类。单尺度拓扑优化设计是将壳体和填充物分别视为两种材料。然而,尽管填充体可以按照预先定义的微观结构进行设计,但单尺度的填充体结构拓扑优化限制了结构性能,使得设计与实际的最优值相差较大。因此,常采用多尺度设计方法对壳体填充结构进行设计,多尺度方法设计多孔结构是以数值均匀化方法将微观结构和宏观材料性能之间联系起来进行设计多孔结构,但从宏观和微观同时设计多孔结构,这无疑给多孔结构拓扑优化设计带来了计算上的挑战。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种单材料多孔结构的拓扑优化设计方法,所述设计方法是一种整体式系统结构优化设计方法,能够优化出满足一定要求的单材料多孔结构,可以满足航空航天、火箭、导弹等军工制造中的结构性能设计需求。
为实现上述目的,按照本发明的一个方面,提供了一种单材料多孔结构的拓扑优化设计方法,所述设计方法包括以下步骤:
(1)分别基于初始水平集函数及径向基函数构建单材料多孔结构,进而通过扩展系数矩阵分块法进行插值以计算得到初始扩展系数,即设计变量初值;
(2)计算单材料多孔结构中连续变化的单元伪密度,并计算结构总体积;接着计算结构的每个单元的局部伪密度,并使用p范数方法求取局部伪密度最大值;
(3)基于得到的单元伪密度对单材料多孔结构的弹性模量进行插值以得到等效单元弹性模量;
(4)基于等效单元弹性模量在结构设计域中求解得到整体结构的位移场,进而计算单材料多孔结构最小柔度拓扑优化模型的目标函数;接着,计算目标函数、结构总体积及局部伪密度对设计变量的灵敏度,并更新全局设计变量,继而得到单材料多孔结构的结构拓扑。
进一步地,分别采用两种初始化结构构建初始水平集函数及径向基函数以构建单材料多孔结构,并通过离散小波变换方法压缩初始水平集函数及径向基函数矩阵,基于压缩后的初始水平集函数和径向基函数通过扩展系数矩阵分块法计算得到初始扩展系数,即设计变量初值。
进一步地,通过隐式水平集函数构建单材料多孔结构,以插值隐式水平集函数的扩展系数为结构设计变量。
进一步地,采用全局径向基函数对水平集函数进行插值。
进一步地,采用高斯积分将原始单元的4个节点上2*2个水平集函数值插值为41*41个水平集函数值。
进一步地,基于Heaviside函数获得单材料多孔结构中连续变化的单元伪密度,并计算结构总体积;接着,通过圆形过滤法计算结构的每个单元的局部伪密度,并使用p范数方法求取局部伪密度最大值。
进一步地,步骤(2)包括以下子步骤:
首先,将设计域Ω均匀划分,并采用Heaviside函数将水平集函数值映射成为有限元模型中的单元伪密度,单元伪密度由以下公式求得:
ρi=∫ΩH(Φi)dΩ
式中,Φ为插值后的单元中节点的水平集函数值,Φi是在插值之前的单元四个节点的水平集函数值;Ω为设计域;H为Heaviside函数;
接着,通过单元伪密度计算结构的总体积;
之后,通过圆形过滤法求解局部伪密度,Ve为任意单元e为中心的局部伪密度:
其中,N表示在该单元的所有单元的数量,
N={i|||ρi-ρc||2≤R}
其中,R和ρc分别表示过滤半径和圆心单元;
接着,使用p范数求解单元伪密度的最大值。
进一步地,使用单元伪密度插值得到等效单元弹性模量的计算公式如下:
E=ρi*E0
其中,E为插值后的等效单元弹性模量;ρi为单元伪密度;E0为材料的弹性模量。
进一步地,单材料多孔结构最小柔度拓扑优化模型的表达式为:
FIND:α=[α1 α2 ... αN]
αi,min≤αi≤αi,max
式中,Vmax表示为整体结构体积约束;Ve和Vp分别表示结构局部伪密度及局部伪密度约束;u和v分别表示允许的位移空间U中的实位移场和虚位移场;u0表示狄利克雷边界上的位移;H(Φ)是Heaviside函数;α是设计变量,用来表示径向基函数插值后的扩展系数,且αi,max和αi,min表示设计变量的上下界限;a(u,v)=l(v)表示为弹性平衡方程弱形式。
进一步地,使用离散小波变换方法压缩优化过程中的整体刚度矩阵、扩展系数矩阵及水平集函数矩阵。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,本发明提供的单材料多孔结构的拓扑优化设计方法主要具有以下有益效果:
1.采用基于水平集方法的局部伪密度约束生成单材料多孔结构,这样不仅可以从宏观角度生成多孔结构,而且优化后的多孔结构具有光滑边界且不存在中间密度。
2.采用矩阵分块法来降低优化过程中所需的内存,可以在有限内存下实现大规模设计域的多孔结构拓扑优化设计工作。
3.引入全局径向基函数GSRBF对水平集函数进行插值,克服其在优化过程中需要重新初始化、速度扩展、迭代步长需满足逆风差分格式CFL条件及无法与许多成熟优化算法结合等数值缺陷,同时采用离散小波变换DWT方法压缩插值系数矩阵,从而形成一种极其稀疏的矩阵插值系统,在保证足够插值精度的前提下,进一步提高优化求解效率。
4.本发明提供的单材料多孔结构拓扑优化设计方法,与现有的技术方案相比,所研究结构是基于参数化水平集方法的单材料多孔结构,具有更加清晰的结构边界,可以满足航空航天、火箭、导弹等军工制造中的结构性能设计需求。
5.所述优化设计方法开发了参数化水平集方法与多孔结构问题相结合,可以得到具有较高稳定性且具有最佳材料相分布的多孔结构。
6.通过矩阵分块、矩阵范数及局部伪密度法的应用,可以实现仅在宏观尺度上对多孔结构的设计,降低优化复杂性及计算量。
附图说明
图1是本发明提供的单材料多孔结构的拓扑优化设计方法的流程示意图;
图2中的(a)、(b)及(c)分别是本发明实施例长悬臂梁初始设计域及初始孔洞位置的示意图;
图3中的(a)、(b)、(c)及(d)分别是本发明实施例均匀初始化下悬臂梁的优化结果示意图;(a)对应的过滤半径为10,体积约束为0.5;(b)对应的过滤半径为10,体积约束为0.6;(c)对应的过滤半径为15,体积约束为0.5;(d)对应的过滤半径为15,体积约束为0.6;
图4中的(a)、(b)、(c)及(d)分别是本发明实施例非均匀初始化下悬臂梁的优化结果示意图;(a)对应的过滤半径为10,体积约束为0.4;(b)对应的过滤半径为10,体积约束为0.5;(c)对应的过滤半径为15,体积约束为0.4;(d)对应的过滤半径为15,体积约束为0.5。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
请参阅图1,本发明提供的单材料多孔结构的拓扑优化设计方法,所述设计方法包括以下步骤:
步骤一,分别采用两种初始化结构构建初始水平集函数及径向基函数以构建单材料多孔结构,并通过离散小波变换方法压缩初始水平集函数及径向基函数矩阵,基于压缩后的初始水平集函数和径向基函数通过扩展系数矩阵分块法计算得到初始扩展系数,即设计变量初值。
具体地,通过隐式水平集函数构建单材料多孔结构,以插值隐式水平集函数的扩展系数为结构设计变量,在N个固定水平集节点处的结构隐式水平集函数为:
式中,x=x1,x2,...,xN表示所有插值节点坐标,即水平集节点;N表示节点总数;αn表示水平集函数在节点n处的扩展系数;Φ表示结构中表示单材料多孔结构的水平集函数;t为时间。该水平集函数由高斯径向基函数φn(x)插值;φn(x)表示高斯径向基函数,公式为:
其中,c是形状参数,等于水平集网格面积或者体积的倒数;xn表示水平集函数的第n个节点的坐标;||x-xn||是用于计算当前采样点x到xn节点距离的欧几里得范数。
为了提高优化效率,离散小波变换方法是减轻全插值矩阵造成的计算成本的关键,采用离散小波变换方法压缩插值系数矩阵,从而形成一种极其稀疏的矩阵插值系统,在保证足够插值精度的前提下,进一步提高优化求解效率。具体地,将原始插值矩阵A转换为其相同大小的小波形式采用小波基的矩阵可以很容易地区分其重要元素和冗余元素;因此,采用一个阈值方法来清除中适当数量的无用元素,并且重新构造一个更稀疏的插值矩阵最后,利用稀疏矩阵可以有效的计算水平集函数。
由于对水平集函数采用径向基函数插值后,在优化求解过程中,径向基函数矩阵A会随着设计域的增大而增大,这将会占用大部分的内存,导致在优化求解时出现矩阵大小超出计算机内存的问题。为了计算大规模的多孔结构,采用矩阵分块的方法对径向基函数矩阵进行分块计算,这种矩阵分块的方法要求Aii(i=1,2,…)的逆存在,以划分为2*2的矩阵块为例,求解扩展系数矩阵C=A-1·f。
扩展系数矩阵C为:
本实施方式中,初始化孔洞采用了两种形式对设计域进行初始化,一种初始化孔洞位置采用均匀密排圆孔结构,另一种为非均匀随机密排圆孔结构,通过这两种形式的初始化方法可以获得两种不同的多孔结构。
步骤二,基于Heaviside函数获得单材料多孔结构中连续变化的单元伪密度,并计算结构总体积;接着,通过圆形过滤法计算结构的每个单元的局部伪密度,并使用p范数方法求取局部伪密度最大值。
具体地,首先将设计域Ω均匀划分,并采用Heaviside函数将水平集函数值映射成为有限元模型中的单元伪密度,对应的公式为:
式中,ζ是一个非常小的正数,本实施方式取ζ=0.001,Δ约等于Heaviside函数带宽的一半。
单元伪密度可以由以下公式求得:
ρi=∫ΩH(Φi)dΩ。
其中,为了计算的更为准确,使用高斯积分将原始单元的4个节点上2*2个水平集函数值插值为41*41个水平集函数值,插值公式如下:
Ni=(1+ξ0)·(1+η0)/4
在上式中,
ξ0=ξi·ξ,η0=ηi·η
ξi={-1,1,1-1},ηi={-1,-1,1,1}
ξ∈[-1:0.05:1],η∈[-1:0.05:1];
Φ为插值后的单元中节点的水平集函数值;Φi是在插值之前的单元四个节点的水平集函数值。
接着,通过单元伪密度计算结构的总体积,公式为:V=∫ΩρidΩ。
之后,通过圆形过滤法求解局部伪密度,公式为:
Ve为任意单元e为中心的局部伪密度:
其中,N表示在该单元的所有单元的数量,
N={i|||ρi-ρc||2≤R}
其中,R和ρc分别表示过滤半径和圆心单元。
接着,使用p范数求解单元伪密度的最大值,公式如下:
当p趋于无穷大时,有:
maxVe=‖Ve‖p。
步骤三,基于连续变化单材料的单元伪密度,对单材料多孔结构弹性模量进行插值,以得到等效单元弹性模量。
具体地,定义材料的弹性模量为E0,泊松比为v,将单材料多孔结构的节点i的位移置为0,即
Uix=0,Uiy=0
其中,Uix和Uiy分别为节点i的x方向位移和y方向位移。
使用单元伪密度插值得到等效单元弹性模量的计算公式如下:
E=ρi*E0
其中,E为插值后的等效单元弹性模量。
步骤四,基于参数化水平集理论建立单材料多孔结构最小柔度拓扑优化模型,并基于得到的等效弹性模量在结构设计域中通过有限元分析求解整体结构的位移场,根据得到的位移场计算单材料多孔结构最小柔度拓扑优化模型的目标函数;接着,基于自伴随方法计算目标函数、结构总体积和局部伪密度对设计变量的灵敏度,并采用MMA移动渐近线算法更新全局设计变量,继而确定单材料多孔结构中材料的最优分布。
具体地,单材料多孔结构最小柔度拓扑优化模型的表达式为:
FIND:α=[α1 α2 ... αN]
αi,min≤αi≤αi,max
式中,Vmax表示为整体结构体积约束,Ve和Vp分别表示结构局部伪密度及局部伪密度约束,u和v分别表示允许的位移空间U中的实位移场和虚位移场,u0表示狄利克雷边界上的位移,H(Φ)是Heaviside函数,α是设计变量,用来表示径向基函数插值后的扩展系数,且αi,max和αi,min表示设计变量的上下界限。a(u,v)=l(v)表示为弹性平衡方程弱形式。
基于虚功原理,针对有限元平衡方程弱形式进行计算,对应的弱形式如下:
其中,a表示双线性能量式;l表示单线性荷载形式;dΩ为结构设计域的积分算子;H表示Heaviside函数,用于表征结构形式的特征函数;ε为应变场;T表示矩阵的转置;u表示结构场的位移;v表示在动力学上允许的位移空间U中的一个虚拟位移;τ表示应用在边界的部分边界上的牵引力;P表示结构设计域的体积力;δ表示Dirac函数,为Heaviside函数的一阶微分;表示差分算子。
基于水平集方法构建单材料多孔结构最小柔度拓扑优化模型包括以下步骤:
(4.1)初始化设计参数、整体体积约束Vmax和局部伪密度约束Vup。
(4.2)基于单元伪密度求得的等效弹性模量E计算结构单元刚度Ke,并组装得到整体刚度矩阵K;接着进行有限元分析求解结构位移场u。
(4.3)基于步骤(4.2)中求解的结构位移场u计算单材料多孔结构最小柔度拓扑优化模型的目标函数J:
其中,ε为应变场;T表示矩阵的转置;u表示结构场的位移;v表示在动力学上允许的位移空间U中的一个虚拟位移;E是通过单元伪密度插值后的等效弹性模量;dΩ为结构设计域的积分算子。
(4.4)基于自伴随方法求解目标函数与约束函数针对设计变量进行灵敏度分析,并采用MMA移动渐近线算法更新全局设计变量,判断模型是否满足收敛条件,若否,则返回到步骤(4.2),若是,则输出单材料多孔结构的最优拓扑结构。其中,通过计算本次与前一步的目标函数差值是否大于1x10-6来决定是否进行下一次优化求解,即判断两次扩展系数的差值是否小于阈值;若是,则优化结束,输出最新的多孔结构拓扑;否则,转入下一次优化求解。
具体地,根据自伴随方法计算目标函数与约束函数针对优化设计变量的一阶微分,计算如下:
其中,
其中,α表示单材料多孔结构的设计变量,是高斯径向基函数插值时的扩展系数,仅与时间变量相关;ε为应变场;u表示结构场的位移;E是单材料多孔结构的等效弹性模量;Φ表示多孔结构的水平集函数;τ表示应用在边界的部分边界上的牵引力;F表示结构设计域的体积力;H表示Heaviside函数;δ表示Dirac函数,为Heaviside函数的一阶微分;dΩ为结构设计域的积分算子;为散度算子;p为p范数的指数;∑(·)为求和符号;Ve为局部伪密度;表示第i个高斯径向基函数。
实施例
本实施例首先给出一个长悬臂梁结构作为算例说明本研究方法的有效性,如图2中的(a)所示为一个厚度为1的矩形区域,其结构左端被完全约束,在最右端中点处施加一个垂直向下的单位力,图2中的(b)为结构采用均匀孔洞进行初始化的结果,图2中的(c)为结构采用非均匀孔洞进行初始化的结果,比较两种初始化结果可以看出,均匀初始化可以明显看出孔洞位置及大小,而非均匀初始化由于孔洞的随机性,导致会有孔洞重叠而呈现出一种散点图的形式。在优化过程中,材料的弹性模量和泊松比分别为1和0.3,且结构被均匀离散成200×100的四节点方形单元。无量纲材料特性的设置是为了方便不同设计之间的比较,以下算例中均采用这种策略。
如图3所示给出了在不同的过滤半径R和体积约束下的均匀初始孔洞位置的优化结果,从优化结果中可以看出,通过调节过滤半径R的大小可以控制生成的多孔结构中分支结构的密度,从优化结果中还可以观察出,当初始孔洞为均匀分布时,生成的多孔结构在外边缘处会生成一圈类似壳状的填充,而分支结构是存在于壳内的。
如图4所示给出了在不同的过滤半径R和体积约束下的非均匀初始孔洞位置的优化结果。从优化结果中可以得出与均匀初始化优化结果相同的结论,即通过调节过滤半径R的大小可以控制生成的多孔结构中分支结构的密度。但与均匀初始化所得结论不同的是,当初始孔洞为非均匀分布时,没有像均匀初始孔洞优化时生成一圈壳体包罗,生成的分支结构是存在于整个结构中的。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种单材料多孔结构的拓扑优化设计方法,其特征在于,所述方法包括以下步骤:
(1)分别基于初始水平集函数及径向基函数构建单材料多孔结构,进而通过扩展系数矩阵分块法进行插值以计算得到初始扩展系数,即设计变量初值;
(2)计算单材料多孔结构中连续变化的单元伪密度,并计算结构总体积;接着计算结构的每个单元的局部伪密度,并使用p范数方法求取局部伪密度最大值;
(3)基于得到的单元伪密度对单材料多孔结构的弹性模量进行插值以得到等效单元弹性模量;
(4)基于等效单元弹性模量在结构设计域中求解得到整体结构的位移场,进而计算单材料多孔结构最小柔度拓扑优化模型的目标函数;接着,计算目标函数、结构总体积及局部伪密度对设计变量的灵敏度,并更新全局设计变量,继而得到单材料多孔结构的结构拓扑。
2.如权利要求1所述的单材料多孔结构的拓扑优化设计方法,其特征在于:分别采用两种初始化结构构建初始水平集函数及径向基函数以构建单材料多孔结构,并通过离散小波变换方法压缩初始水平集函数及径向基函数矩阵,基于压缩后的初始水平集函数和径向基函数通过扩展系数矩阵分块法计算得到初始扩展系数,即设计变量初值。
3.如权利要求2所述的单材料多孔结构的拓扑优化设计方法,其特征在于:通过隐式水平集函数构建单材料多孔结构,以插值隐式水平集函数的扩展系数为结构设计变量。
4.如权利要求1所述的单材料多孔结构的拓扑优化设计方法,其特征在于:采用全局径向基函数对水平集函数进行插值。
5.如权利要求4所述的单材料多孔结构的拓扑优化设计方法,其特征在于:采用高斯积分将原始单元的4个节点上2*2个水平集函数值插值为41*41个水平集函数值。
6.如权利要求1所述的单材料多孔结构的拓扑优化设计方法,其特征在于:基于Heaviside函数获得单材料多孔结构中连续变化的单元伪密度,并计算结构总体积;接着,通过圆形过滤法计算结构的每个单元的局部伪密度,并使用p范数方法求取局部伪密度最大值。
7.如权利要求6所述的单材料多孔结构的拓扑优化设计方法,其特征在于:步骤(2)包括以下子步骤:
首先,将设计域Ω均匀划分,并采用Heaviside函数将水平集函数值映射成为有限元模型中的单元伪密度,单元伪密度由以下公式求得:
ρi=∫ΩH(Φi)dΩ
式中,Φ为插值后的单元中节点的水平集函数值,Φi是在插值之前的单元四个节点的水平集函数值;Ω为设计域;H为Heaviside函数;
接着,通过单元伪密度计算结构的总体积;
之后,通过圆形过滤法求解局部伪密度,Ve为任意单元e为中心的局部伪密度:
其中,N表示在该单元的所有单元的数量,
N={i|||ρi-ρc||2≤R}
其中,R和ρc分别表示过滤半径和圆心单元;
接着,使用p范数求解单元伪密度的最大值。
8.如权利要求1-7任一项所述的单材料多孔结构的拓扑优化设计方法,其特征在于:使用单元伪密度插值得到等效单元弹性模量的计算公式如下:
E=ρi*E0
其中,E为插值后的等效单元弹性模量;ρi为单元伪密度;E0为材料的弹性模量。
9.如权利要求1-7任一项所述的单材料多孔结构的拓扑优化设计方法,其特征在于:单材料多孔结构最小柔度拓扑优化模型的表达式为:
FIND:α=[α1 α2 ... αN]
αi,min≤αi≤αi,max
10.如权利要求1-7任一项所述的单材料多孔结构的拓扑优化设计方法,其特征在于:使用离散小波变换方法压缩优化过程中的整体刚度矩阵、扩展系数矩阵及水平集函数矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110396411.7A CN113191044B (zh) | 2021-04-13 | 2021-04-13 | 一种单材料多孔结构的拓扑优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110396411.7A CN113191044B (zh) | 2021-04-13 | 2021-04-13 | 一种单材料多孔结构的拓扑优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113191044A true CN113191044A (zh) | 2021-07-30 |
CN113191044B CN113191044B (zh) | 2023-03-28 |
Family
ID=76975608
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110396411.7A Active CN113191044B (zh) | 2021-04-13 | 2021-04-13 | 一种单材料多孔结构的拓扑优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113191044B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113742970A (zh) * | 2021-08-30 | 2021-12-03 | 华中科技大学 | 一种点阵结构优化设计的基于水平集的变密度方法 |
CN115062510A (zh) * | 2022-06-14 | 2022-09-16 | 华南理工大学 | 一种结果逼近的多目标拓扑优化方法、系统、介质和设备 |
CN115631817A (zh) * | 2022-10-31 | 2023-01-20 | 华中科技大学 | 基于复杂物理场的多材料全尺度拓扑优化设计方法及应用 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108875125A (zh) * | 2018-04-27 | 2018-11-23 | 北京航空航天大学 | 一种位移和全局应力混合约束下的连续体双材料结构拓扑优化方法 |
CN109190233A (zh) * | 2018-08-24 | 2019-01-11 | 华南理工大学 | 一种结构拓扑优化方法 |
CN110795873A (zh) * | 2019-09-30 | 2020-02-14 | 北京擎靖天启科技服务有限公司 | 一种考虑尺寸控制的跨尺度拓扑优化方法 |
CN111460622A (zh) * | 2020-03-06 | 2020-07-28 | 华中科技大学 | 一种基于参数化彩色水平集的多组件布局拓扑优化方法 |
CN112100882A (zh) * | 2020-08-27 | 2020-12-18 | 华南理工大学 | 一种具有光滑边界的连续体结构密度演化拓扑优化方法 |
-
2021
- 2021-04-13 CN CN202110396411.7A patent/CN113191044B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108875125A (zh) * | 2018-04-27 | 2018-11-23 | 北京航空航天大学 | 一种位移和全局应力混合约束下的连续体双材料结构拓扑优化方法 |
CN109190233A (zh) * | 2018-08-24 | 2019-01-11 | 华南理工大学 | 一种结构拓扑优化方法 |
CN110795873A (zh) * | 2019-09-30 | 2020-02-14 | 北京擎靖天启科技服务有限公司 | 一种考虑尺寸控制的跨尺度拓扑优化方法 |
CN111460622A (zh) * | 2020-03-06 | 2020-07-28 | 华中科技大学 | 一种基于参数化彩色水平集的多组件布局拓扑优化方法 |
CN112100882A (zh) * | 2020-08-27 | 2020-12-18 | 华南理工大学 | 一种具有光滑边界的连续体结构密度演化拓扑优化方法 |
Non-Patent Citations (2)
Title |
---|
HAO LI ET AL.: "Spatial-varying multi-phase infill design using density-based topology optimization", 《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》 * |
李好: "改进的参数化水平集拓扑优化方法与应用研究", 《中国博士学位论文全文数据库》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113742970A (zh) * | 2021-08-30 | 2021-12-03 | 华中科技大学 | 一种点阵结构优化设计的基于水平集的变密度方法 |
CN113742970B (zh) * | 2021-08-30 | 2023-12-26 | 华中科技大学 | 一种点阵结构优化设计的基于水平集的变密度方法 |
CN115062510A (zh) * | 2022-06-14 | 2022-09-16 | 华南理工大学 | 一种结果逼近的多目标拓扑优化方法、系统、介质和设备 |
CN115631817A (zh) * | 2022-10-31 | 2023-01-20 | 华中科技大学 | 基于复杂物理场的多材料全尺度拓扑优化设计方法及应用 |
Also Published As
Publication number | Publication date |
---|---|
CN113191044B (zh) | 2023-03-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113191044B (zh) | 一种单材料多孔结构的拓扑优化设计方法 | |
CN112100882B (zh) | 一种具有光滑边界的连续体结构密度演化拓扑优化方法 | |
CN110069800B (zh) | 具有光滑边界表达的三维结构拓扑优化设计方法及设备 | |
WO2022000132A1 (zh) | 基于三周期极小曲面的三维多孔散热结构的设计与优化方法 | |
CN111737835A (zh) | 基于三周期极小曲面的三维多孔散热结构的设计与优化方法 | |
CN111709085B (zh) | 一种约束阻尼薄板结构拓扑优化设计方法 | |
CN112765856B (zh) | 一种功能梯度多孔结构拓扑优化的混合水平集方法 | |
CN114372438B (zh) | 基于轻量化深度强化学习的芯片宏单元布局方法及系统 | |
CN111489447B (zh) | 一种适用于格子Boltzmann方法的直角网格自适应建模方法 | |
CN108875132B (zh) | 一种多工况下应力约束类桁架连续体结构的拓扑优化方法 | |
CN111523270A (zh) | 一种改进的连续体结构拓扑优化后处理方法 | |
CN114239356A (zh) | 一种基于有限元网格的共形点阵材料设计方法 | |
CN114254408A (zh) | 一种基于代理模型的梯度点阵等几何拓扑优化方法 | |
CN113887095A (zh) | 一种基于等几何分析的渐进式结构拓扑优化方法 | |
CN115203997A (zh) | 一种基于多变量设计的点阵-实体复合结构拓扑优化方法 | |
CN108897956B (zh) | 一种多孔机械零部件优化设计方法 | |
CN115310209B (zh) | 基于vae的气动形状迁移优化方法及相关装置 | |
CN111079279A (zh) | 一种多构型晶格结构多尺度拓扑优化设计方法 | |
CN114241026A (zh) | 基于平坦度划分的点云精简算法和装置 | |
CN111047687B (zh) | 一种基于三维t样条的异质材料实体建模方法 | |
CN114972695A (zh) | 一种点云生成方法、装置、电子设备及存储介质 | |
CN114757002A (zh) | 一种约束阻尼板壳的结构-材料分层拓扑优化设计方法 | |
CN113987860B (zh) | 一种基于动态高斯核卷积滤波的结构优化方法和装置 | |
Alexa | Polycover: Shape approximating with discrete surface orientation | |
Yap et al. | Adaptive image restoration based on hierarchical neural networks |
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 |