CN107301295A - 适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法 - Google Patents

适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法 Download PDF

Info

Publication number
CN107301295A
CN107301295A CN201710487046.4A CN201710487046A CN107301295A CN 107301295 A CN107301295 A CN 107301295A CN 201710487046 A CN201710487046 A CN 201710487046A CN 107301295 A CN107301295 A CN 107301295A
Authority
CN
China
Prior art keywords
mrow
msubsup
attribute
auxetic
function
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
Application number
CN201710487046.4A
Other languages
English (en)
Other versions
CN107301295B (zh
Inventor
李好
高亮
付君健
肖蜜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN201710487046.4A priority Critical patent/CN107301295B/zh
Publication of CN107301295A publication Critical patent/CN107301295A/zh
Application granted granted Critical
Publication of CN107301295B publication Critical patent/CN107301295B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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]

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)
  • Laminated Bodies (AREA)

Abstract

本发明属于超材料设计相关技术领域，其公开了一种适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法，其包括以下步骤：(1)采用加权法建立每层的、同时考虑微观结构拉胀属性及宏观结构刚度属性的多目标优化模型；(2)控制所述超材料结构整体在一个方向上的拉胀属性及刚度属性呈梯度形式分布；(3)计算每层设计域上的微观位移、宏观位移、以及每层微结构的等效弹性属性；(4)基于计算获得的敏度，采用优化准则法更新设计变量，进而判断所述目标函数是否收敛。上述拓扑优化方法采用加权法建立同时考虑微观结构拉胀属性及宏观结构刚度属性的多目标优化模型继而进行优化，使得优化后的超材料同时具备拉胀属性及足够的刚度。

Description

适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法
技术领域
本发明属于超材料设计相关技术领域，更具体地，涉及一种适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法。
背景技术
超材料是指具有人工设计的微结构且呈现出天然材料不具备或者较少具备的超常物理属性的复合材料。拉胀材料是一类具有负泊松比的新型力学超材料，它具有受拉时垂直方向膨胀(受压时垂直方向收缩)的非常规力学特性，以此增强材料的剪切模量、吸能能力、抗裂性和断裂韧性，拥有正泊松比材料所不具备的功能特点。
基于拉胀超材料的工程产品应具备多功能性，例如舰艇的防护结构和车辆防撞结构，更加强调特殊功能性和承载性的综合：一方面，产品在受冲击时，拉胀效应(负泊松比效应)使得材料瞬间朝受冲击处聚集，从而吸收更多能量，防止结构被破坏；另一方面，产品必须具备足够的刚度以保障结构的稳定性及静、动态性能。通常来讲，拉胀超材料的细观多孔微结构必须满足疏松，以实现其中铰链机构的收缩与伸展，进而呈现出拉胀特性，然而这也导致基于拉胀超材料的工程结构在承载时缺乏刚度。相应地，本领域存在着发展一种能够同时考虑拉胀属性及足够刚度的超材料拓扑优化方法的技术需求。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提供了一种适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法，其同时考虑微结构构型在空间上的变化、不同性能指标的综合以及宏观载荷与边界条件对为微结构的影响等多个关键因素，解决了基于仿真、试验或者实验的设计方法无法实现或者实现成本过高的问题，使得经拓扑优化后的超材材料结构同时具有拉胀属性及足够的刚度。
为实现上述目的，本发明提供了一种适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法，其包括以下步骤：
(1)将待优化的超材料结构沿着所述超材料结构的梯度方向进行分层，并采用加权法建立每层的、同时考虑微观结构拉胀属性及宏观结构刚度属性的多目标优化模型；
(2)通过每层的微结构的体积比来控制所述超材料结构整体在一个方向上的拉胀属性及刚度属性呈梯度形式分布；
(3)计算每层设计域上的微观位移、宏观位移、以及每层微结构的等效弹性属性，计算所述多目标优化模型的目标函数及体积约束函数的敏度；
(4)基于计算获得的敏度，采用优化准则法更新设计变量，进而判断所述目标函数是否收敛，若收敛，则输出最优宏微观结构构型及分布，否则转至步骤(3)。
进一步地，采用加权法建立每层的、同时考虑微观结构拉胀属性及宏观结构刚度属性的多目标优化模型包括以下子步骤：
(a)分别建立柔度优化模型及拉胀属性优化模型，其中所述柔度优化模型的宏观目标函数由每层微结构的柔度值之和来表示，所述拉胀属性优化模型的微观目标函数由弹性材料矩阵中分量的线性组合来表示；
(b)采用归一化指数加权准则对所述柔度优化模型及所述拉胀属性优化模型进行加权处理；
(c)将所述宏观目标函数及所述微观目标函数进行线性加权，以得到所述多目标优化模型。
进一步地，所述宏观目标函数fMA如以下公式所示：
式中，m＝1,2,…,M，表示M层微结构中的第m层；上标MA表示宏观量；上标MI表示微观量；表示第m层的宏观位移；表示第m层的微观位移；表示第m层的宏观应变；表示微结构等效弹性矩阵；表示微观水平集函数；ΩMA表示宏观域。
进一步地，所述微观目标函数如以下公式所示：
式中，分别对应材料弹性矩阵中的分量，通过控制优化过程中出现铰链机构形式的负泊松比结构，通过增加结构的剪切性能。
进一步地，所述多目标优化模型Jm由以下公式表示：
式中， exp(·)表示指数函数；加权因子wm,a和wm,b分别表示宏观目标函数和微观目标函数在m层的权重，wm,a和wm,b之和等于1；分别代表宏观、微观目标函数的最小值、最大值；q＝2是归一化指数加权准则的指数常数。
进一步地，所述超材料结构的每层由多个单元设计域组成，每个单元设计域上设置有连接域，相邻两层的所述单元设计域通过对应的所述连接域相连接。
进一步地，相邻两个所述单元设计域的连接域的位置相对应，所有所述连接域沿所述超材料结构的梯度方向设置，且相邻两个单元设计域上的连接域的尺寸一致。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，本发明提供的适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法主要具有以下有益效果：
(1)所述拓扑优化方法将待优化的超材料结构沿着所述超材料结构的梯度方向进行分层，使得分层后的超材料结构优化问题被分解为层优化问题，且多个层优化并行进行，增加了宏观结构整体刚度；
(2)所述拓扑优化方法采用加权法建立同时考虑微观结构拉胀属性及宏观结构刚度属性的多目标优化模型继而进行优化，使得优化后的超材料同时具备拉胀属性及足够的刚度；
(3)所述超材料结构的每层由多个单元设计域组成，每个单元设计域上设置连接域，相邻两层的单元设计域通过对应的所述连接域相连接，保证了微结构之间的连续性；
(4)所有连接域沿超材料结构的梯度方向设置，且相邻两个单元设计域上的连接域的尺寸一致，保证了微结构的数值精度和周期性。
附图说明
图1是本发明较佳实施方式提供的适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法的流程图；
图2是图1中的适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法涉及的分层设计示意图；
图3是图1中的适用于具有功能梯度和拉胀属性的超材料的拓扑优化方法涉及的超材料的不同单元设计域的连接示意图；
图4是采用采用图1中的适用于具有功能梯度和拉胀属性的超材料的拓扑优化方法进行拓扑优化的矩形梁初始设计域的示意图；
图5是图4中的矩形梁经过优化后的超材料示意图；
图6是图4中的矩形梁在一个梯度方向的放大图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
请参阅图1至图3，本发明较佳实施方式提供的适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法，所述拓扑优化方法是一系统的优化方法，其同时考虑微结构构型在空间上的变化、不同性能指标的综合以及宏观载荷与边界条件对为微结构的影响等多个关键因素，解决了基于仿真、试验或者实验的设计方法无法实现或者实现成本过高的问题，使得经拓扑优化后的超材料同时具有拉胀属性及足够的刚度。
本实施方式中，所述的适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法主要具有以下有益效果：
步骤一，将待优化的超材料结构沿着所述超材料结构的梯度方向进行分层，使得分层后的超材料结构优化被分解为层优化，并采用加权法建立同时考虑微观结构拉胀属性及宏观结构刚度属性的多目标优化模型。
具体地，将待优化的超材料结构进行分层，将超材料结构优化问题分为多个子优化问题，对每一层建立一种同时考虑微观结构拉胀属性和宏观结构刚度属性的多目标优化模型。
对于一个被分为M层的超材料结构，宏观目标函数由每一层微结构的柔度值之和表示：
式中，m＝1,2,…,M，表示M层微结构中的第m层；上标MA表示宏观量；上标MI表示微观量；表示第m层的宏观位移；表示第m层的微观位移；表示第m层的宏观应变；表示微结构等效弹性矩阵；表示微观水平集函数；ΩMA表示宏观域。
考虑各向同性和各向正交异性材料的微结构设计，第m层微结构的拉胀属性可表示为：
式中，分别对应材料弹性矩阵中的分量，通过控制优化过程中出现铰链机构形式的负泊松比结构，通过增加结构的剪切性能。
为建立多种性能属性的多目标优化模型，首先采用归一化指数加权准则对各单目标优化模型(柔度优化模型及拉胀属性优化模型)进行处理。
式中，exp(·)表示指数函数；加权因子wm,a和wm,b分别表示宏观目标函数和微观目标函数在m层的权重，wm,a和wm,b之和等于1； 分别代表宏观、微观目标函数的最小值、最大值；q＝2是归一化指数加权准则的指数常数。
同时考虑微观结构拉胀属性和宏观结构刚度属性的多目标优化模型通过对两种目标函数的线性加权得到，所述多目标优化模型由以下公式表示：
步骤二，通过每层的微结构的体积比来控制所述超材料整体结构在某一方向上的拉胀属性及刚度属性呈梯度形式分布。具体地，超材料结构的每层由周期性微结构组成，通过微结构的体积比控制拉胀材料整体在某一方向上的拉胀属性和刚度属性呈梯度形式分布；超材料结构的第m层微结构的最大体积分数定义如下：
式中，[V1,VM]表示体积约束范围，ξ＝1时体积分数呈线性变化，否则体积分数呈非线性变化。
步骤三，对每层微结构分别进行微观有限元分析及宏观有限分析以获得设计域上的宏观位移及微观位移，同时计算每层微结构的等效弹性属性，进而计算所述多目标优化模型的目标函数及体积约束函数的敏度。
具体地，首先，对每层微结构进行微观有限元分析，以获得设计域上的微观位移。微观有限元模型的双线性能量形式和线性载荷形式分别表示为：
式中，表示微观平衡方程的弱形式；和分别为定义在微观设计域内的周期性位移；Dpqrs为微观实体材料的弹性矩阵；和分别为微观扰动应变；为单位测试应变；H为Heaviside函数。
之后，通过数值均匀化方法计算每一层微结构的等效弹性属性。在水平集方法的框架下，用数值均匀化方法来近似计算每一层微结构的等效弹性属性，二维结构中的等效弹性量计算公式如下所示：
式中，i,j,k,l＝1,2；Dpqrs为微观实体材料的弹性矩阵；上标MI表示微观量；ΩMI微观结构设计域；|ΩMI|表示微观结构的面积；ΦMI是微观水平集函数；是单位测试应变(对于二维问题，包括水平单位应变、竖直单位应变和剪切单位应变)。
接着，采用计算获得的等效弹性模量对每层微结构进行宏观有限元分析以获得设计域上的宏观位移及柔度值。宏观有限元模型的双线性能量形式和线性载荷形式分别表示为：
式中，aMA(uMA,vMA,uMI,ΦMI)＝lMA(vMA)表示宏观平衡方程的弱形式；p表示体积力；τ是应用在边界ΓMA上的牵引力；uMA和vMA分别是定义在ΩMA域上的非周期性位移。
最后，基于计算获得的所述微观位移、宏观位移及等效弹性属性来计算所述多目标优化模型的目标函数及体积约束函数的敏度。所述目标函数的敏度计算公式如下所示：
式中，为参数化水平集方法中径向基插值函数的扩展系数。
所述体积约束函数的敏度计算公式如下所示：
步骤四，基于计算获得的敏度，采用优化准则法更新设计变量(水平基函数)以达到微结构优化的目的，进而判断所述目标函数是否收敛，若收敛，则输出最优宏微观结构构型及分布，否则转至步骤三。
为了增加宏观结构的整体刚度，所有层的优化都并行进行。每层由多个单元设计域组成，每个单元设计域上设置连接域D，相邻两层的单元设计域通过对应的所述连接域D相连接，以保证微结构之间的连续性。本实施方式中，相邻两个单元设计域的连接域D的位置相对应，所有连接域D沿超材料结构的梯度方向F设置，且相邻两个单元设计域上的连接域D的尺寸一致，以保证微结构的数值精度和周期性。预先设定的连接域D在优化过程中被设置为非设计域。
请参阅图4至图6，以下以一个具体实例来对本发明进行进一步的详细说明。以一个长L＝30cm宽H＝10cm的矩形梁为例，下边界中点处有一集中载荷F＝50N，上边界固定约束，结构分为10层，梯度方向为Y轴方向，体积分数沿Y轴负方向逐渐增大。本实例中仅考虑各向正交异性的微结构，宏观结构离散为30×10＝300个四节点单元。每一层上的微结构离散为40×40＝1600个四点节单元。
为实现功能梯度属性，每一层的微结构体积分数在梯度方向上从0.25到0.55呈非线性变化，取ξ＝2，10层微结构的体积分数约束分别为0.25、0.2537、0.2648、0.2833、0.3039、0.3406、0.3833、0.4315、0.4870、0.55。同时考虑微观结构拉胀属性和宏观结构刚度属性的多目标优化模型中取wm,a＝0.4和wm,b＝0.6
优化后得到的具有拉胀属性的超材料如图5所示，图6所示为一放大的矩形梁的超材料体积单元示意图，总柔度值fMA＝157.9代表体积单元由一系列典型的具有拉胀属性的超材料组成；正交各向异性微结构的泊松比v12＝-1.147，v21＝-0.2608；每一层的微结构之间有较好的连接性，由于采用的基于参数化水平集的拓扑优化方法，微结构具有光滑的边界。
本发明提供的适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法，其将待优化的超材料结构沿着所述超材料结构的梯度方向进行分层，使得分层后的超材料结构优化问题被分解为多个子优化问题，且多个子化问题并行优化，增加了宏观结构整体刚度。同时，所述拓扑优化方法采用加权法建立同时考虑微观结构拉胀属性及宏观结构刚度属性的多目标优化模型继而进行优化，使得优化后的超材料同时具备拉胀属性及足够的刚度。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法，其特征在于，该拓扑优化方法包括以下步骤：
(1)将待优化的超材料结构沿着所述超材料结构的梯度方向进行分层，并采用加权法建立每层的、同时考虑微观结构拉胀属性及宏观结构刚度属性的多目标优化模型；
(2)通过每层的微结构的体积比来控制所述超材料结构整体在一个方向上的拉胀属性及刚度属性呈梯度形式分布；
(3)计算每层设计域上的微观位移、宏观位移、以及每层微结构的等效弹性属性，进而计算所述多目标优化模型的目标函数及体积约束函数的敏度；
(4)基于计算获得的敏度，采用优化准则法更新设计变量，进而判断所述目标函数是否收敛，若收敛，则输出最优宏微观结构构型及分布，否则转至步骤(3)。
2.如权利要求1所述的适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法，其特征在于：采用加权法建立每层的、同时考虑微观结构拉胀属性及宏观结构刚度属性的多目标优化模型包括以下子步骤：
(a)分别建立柔度优化模型及拉胀属性优化模型，其中所述柔度优化模型的宏观目标函数由每层微结构的柔度值之和来表示，所述拉胀属性优化模型的微观目标函数由弹性材料矩阵中分量的线性组合来表示；
(b)采用归一化指数加权准则对所述柔度优化模型及所述拉胀属性优化模型进行加权处理；
(c)将所述宏观目标函数及所述微观目标函数进行线性加权，以得到所述多目标优化模型。
3.如权利要求2所述的适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法，其特征在于：所述宏观目标函数如以下公式所示：
<mrow> <msup> <mi>f</mi> <mrow> <mi>M</mi> <mi>A</mi> </mrow> </msup> <mo>(</mo> <mrow> <msup> <mi>u</mi> <mrow> <mi>M</mi> <mi>A</mi> </mrow> </msup> <mo>,</mo> <msup> <mi>u</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msup> <mo>,</mo> <msup> <mi>&amp;Phi;</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msup> </mrow> <mo>)</mo> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msubsup> <mi>f</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>A</mi> </mrow> </msubsup> <mo>(</mo> <mrow> <msubsup> <mi>u</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>A</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>u</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>&amp;Phi;</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> </mrow> <mo>)</mo> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mo>&amp;Integral;</mo> <msubsup> <mi>&amp;Omega;</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>A</mi> </mrow> </msubsup> </msub> <msubsup> <mi>&amp;epsiv;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>u</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>A</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <msubsup> <mi>D</mi> <mrow> <mi>i</mi> <mi>j</mi> <mi>k</mi> <mi>l</mi> </mrow> <mi>H</mi> </msubsup> <mo>(</mo> <mrow> <msubsup> <mi>u</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>&amp;Phi;</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> </mrow> <mo>)</mo> <msub> <mi>&amp;epsiv;</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> <mo>(</mo> <msubsup> <mi>u</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>A</mi> </mrow> </msubsup> <mo>)</mo> <msubsup> <mi>d&amp;Omega;</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>A</mi> </mrow> </msubsup> </mrow>
式中，m＝1,2,…,M，表示M层微结构中的第m层；上标MA表示宏观量；上标MI表示微观量；表示第m层的宏观位移；表示第m层的微观位移；表示第m层的宏观应变的转置；表示微结构等效弹性矩阵；表示微观水平集函数；ΩMA表示宏观域。
4.如权利要求3所述的适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法，其特征在于：所述微观目标函数如以下公式所示：
<mrow> <msubsup> <mi>f</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>u</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>&amp;Phi;</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mi>D</mi> <mn>1122</mn> <mi>H</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>u</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>&amp;Phi;</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>D</mi> <mn>2211</mn> <mi>H</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>u</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>&amp;Phi;</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>-</mo> <mn>2</mn> <msubsup> <mi>D</mi> <mn>1212</mn> <mi>H</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>u</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>&amp;Phi;</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow>
式中，分别对应材料弹性矩阵中的分量，通过控制优化过程中出现铰链机构形式的负泊松比结构，通过增加结构的剪切性能。
5.如权利要求4所述的适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法，其特征在于：所述多目标优化模型Jm由以下公式表示：
<mrow> <msub> <mi>J</mi> <mi>m</mi> </msub> <mrow> <mo>(</mo> <msubsup> <mi>u</mi> <mi>m</mi> <mrow> <mi>M</mi> <mn>4</mn> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>u</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>&amp;Phi;</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mi>J</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>A</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>u</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>A</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>u</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>&amp;Phi;</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>J</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>u</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>&amp;Phi;</mi> <mi>m</mi> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> 1
式中， exp(·)表示指数函数；加权因子wm,a和wm,b分别表示宏观目标函数和微观目标函数在m层的权重，wm,a和wm,b之和等于1；分别代表宏观、微观目标函数的最小值、最大值；q＝2是归一化指数加权准则的指数常数。
6.如权利要求1所述的适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法，其特征在于：所述超材料结构的每层由多个单元设计域组成，每个单元设计域上设置有连接域，相邻两层的所述单元设计域通过对应的所述连接域相连接。
7.如权利要求6所述的适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法，其特征在于：相邻两个所述单元设计域的连接域的位置相对应，所有所述连接域沿所述超材料结构的梯度方向设置，且相邻两个单元设计域上的连接域的尺寸一致。
CN201710487046.4A 2017-06-23 2017-06-23 适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法 Active CN107301295B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710487046.4A CN107301295B (zh) 2017-06-23 2017-06-23 适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710487046.4A CN107301295B (zh) 2017-06-23 2017-06-23 适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法

Publications (2)

Publication Number Publication Date
CN107301295A true CN107301295A (zh) 2017-10-27
CN107301295B CN107301295B (zh) 2018-06-12

Family

ID=60135049

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710487046.4A Active CN107301295B (zh) 2017-06-23 2017-06-23 适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法

Country Status (1)

Country Link
CN (1) CN107301295B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107563010A (zh) * 2017-08-08 2018-01-09 西北工业大学 基于形状特征的多尺度结构材料一体化设计方法
CN109472056A (zh) * 2018-10-15 2019-03-15 上海交通大学 任意泊松比超材料的拓扑优化形成方法
CN109657284A (zh) * 2018-11-27 2019-04-19 华中科技大学 一种面向超材料的等几何拓扑优化方法
CN109871574A (zh) * 2018-12-28 2019-06-11 华中科技大学 一种基于代理模型的多尺度拓扑优化方法
CN110955938A (zh) * 2019-11-06 2020-04-03 华中科技大学 一种具有梯度多孔夹芯的夹层结构拓扑优化方法
CN116432330A (zh) * 2022-12-23 2023-07-14 华中科技大学 功能梯度拉胀超材料填充的多尺度壳体设计方法及设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080183436A1 (en) * 2007-01-30 2008-07-31 Industrial Technology Research Institute Method of evolutionary optimization algorithm for structure design
CN104584326A (zh) * 2012-05-09 2015-04-29 杜克大学 超材料设备及使用该超材料设备的方法
CN105005647A (zh) * 2015-07-02 2015-10-28 哈尔滨工程大学 一种基于阻尼最大化的复合三明治板壳结构的拓扑优化方法
CN106457748A (zh) * 2014-01-24 2017-02-22 墨尔本皇家理工大学 结构化多孔超材料

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080183436A1 (en) * 2007-01-30 2008-07-31 Industrial Technology Research Institute Method of evolutionary optimization algorithm for structure design
CN104584326A (zh) * 2012-05-09 2015-04-29 杜克大学 超材料设备及使用该超材料设备的方法
CN106457748A (zh) * 2014-01-24 2017-02-22 墨尔本皇家理工大学 结构化多孔超材料
CN105005647A (zh) * 2015-07-02 2015-10-28 哈尔滨工程大学 一种基于阻尼最大化的复合三明治板壳结构的拓扑优化方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
WEIHONG ZHANG等: "《Scale-related topology optimization of cellular materials and structures》", 《INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING》 *
李好: "《改进的参数化水平集拓扑优化方法与应用研究》", 《中国学位论文全文数据库》 *
王昱等: "《负泊松比负热膨胀超材料微结构拓扑优化设计》", 《上海理工大学学报》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107563010A (zh) * 2017-08-08 2018-01-09 西北工业大学 基于形状特征的多尺度结构材料一体化设计方法
CN107563010B (zh) * 2017-08-08 2020-09-25 西北工业大学 基于形状特征的多尺度结构材料一体化设计方法
CN109472056A (zh) * 2018-10-15 2019-03-15 上海交通大学 任意泊松比超材料的拓扑优化形成方法
CN109657284A (zh) * 2018-11-27 2019-04-19 华中科技大学 一种面向超材料的等几何拓扑优化方法
CN109657284B (zh) * 2018-11-27 2020-12-29 华中科技大学 一种面向超材料的等几何拓扑优化方法
CN109871574A (zh) * 2018-12-28 2019-06-11 华中科技大学 一种基于代理模型的多尺度拓扑优化方法
CN110955938A (zh) * 2019-11-06 2020-04-03 华中科技大学 一种具有梯度多孔夹芯的夹层结构拓扑优化方法
CN116432330A (zh) * 2022-12-23 2023-07-14 华中科技大学 功能梯度拉胀超材料填充的多尺度壳体设计方法及设备
CN116432330B (zh) * 2022-12-23 2024-03-19 华中科技大学 功能梯度拉胀超材料填充的多尺度壳体设计方法及设备

Also Published As

Publication number Publication date
CN107301295B (zh) 2018-06-12

Similar Documents

Publication Publication Date Title
CN107301295B (zh) 适用于具有功能梯度及拉胀属性的超材料的拓扑优化方法
Liu et al. General regression neural network for prediction of sound absorption coefficients of sandwich structure nonwoven absorbers
Mir et al. Review of mechanics and applications of auxetic structures
Chen et al. Free vibration analysis of composite sandwich plates with different truss cores
Hu et al. Topology optimization of lightweight lattice structural composites inspired by cuttlefish bone
Li et al. Optimization design of lightweight structure inspired by glass sponges (Porifera, Hexacinellida) and its mechanical properties
Biswas et al. A micromorphic computational homogenization framework for auxetic tetra-chiral structures
Wheeler et al. Soft origami: Classification, constraint, and actuation of highly compliant origami structures
Tamijani et al. Vibration analysis of curvilinearly-stiffened functionally graded plate using element free Galerkin method
Yarasca et al. Best Theory Diagrams for cross-ply composite plates using polynomial, trigonometric and exponential thickness expansions
Laurien et al. Nonlocal wrinkling instabilities in bilayered systems using peridynamics
Pham et al. Nonlocal free vibration of functionally graded porous nanoplates using higher-order isogeometric analysis and ANN prediction
Liu et al. Mechanical investigation on failure modes of reticular porous metal foams under different loadings in engineering applications
CN109002614A (zh) 一种稳定成孔的改进水平集拓扑优化方法
Arshid et al. Aerodynamic stability and free vibration of FGP-Reinforced nano-fillers annular sector microplates exposed to supersonic flow
Nascimento et al. Numerical analysis of honeycomb-shaped polymeric foams using the FEM and the RPIM
Talatahari et al. Effect of the electric field on the sound transmission loss of double-walled electro-rheological fluid sandwich plates with functionally graded carbon nanotube reinforced composite facesheets
Lim Refined shear correction factor for very thick simply supported and uniformly loaded isosceles right triangular auxetic plates
Sahoo Laminated composite stiffened cylindrical shell panels with cutouts under free vibration
Tiwari et al. Dynamic analysis of composite cylinders using 3-D degenerated shell elements
Li et al. Architectural synthesis and analysis of dual-cellular fluidic flexible matrix composites for multifunctional metastructures
Yang Ship Vibration and Noise Reduction with Metamaterial Structures
Tornabene et al. Higher-Order Weak Formulation for Arbitrarily Shaped Doubly-Curved Shells
Li et al. Non-linear buckling and postbuckling behavior of 3d braided cylindrical panels under axial compression in thermal environments
Lu et al. Radial Basis Function Neural Network in Vibration Control of Civil Engineering Structure

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