疏水多级微结构的纤维丛拓扑优化方法与流程

[0001]
本发明属于疏水多级微结构技术领域，本发明涉及疏水固体表面上周期性多级微结构的几何构型拓扑优化技术领域，特别涉及一种实现多级结构逆向设计的纤维丛拓扑优化方法。


背景技术：

[0002]
润湿是表面化学的一个重要研究方向。固体表面的润湿能力可以按疏水性，亲水性，疏油性，亲油性，双疏性和双亲性进行分类。在具有特殊润湿能力的人工表面微结构研究中，最新成果主要是针对疏水微结构。
[0003]
对于一个粗糙平面，润湿取决于固体表面自由能和表面微结构构型。粗糙平面上的润湿现象具有两种不同的模态，分别是wenzel模态和cassie-baxter模态。在wenzel 模态中，液体完全填充粗糙平面上的微结构。在cassie-baxter模态中，液汽界面和微结构之间存在汽囊，且汽囊会被束缚在微结构中；此时，粗糙表面相当于一个固体和气体的复合表面，因而粗糙表面疏水性得到有效的增强。当增大施加于液体的压强时，固体和液体的接触模态可以从cassie-baxter模态转变为wenzel模态。该转变过程中，液体会填满粗糙表面的微结构，进而导致疏水性的降低。因此，固体表面上微结构几何构型的合理性决定了微结构抵抗cassie-baxter模态向wenzel模态转变的能力。该抵抗能力可通过cassie-baxter模态中液汽界面凸起的体积衡量；并且体积越小，微结构抵抗模态转变的能力越高。
[0004]
固体表面微结构主要包括层状和多级微结构两类。层状微结构可视为对多级微结构的近似，且多级微结构能够支持更多处于亚稳态的cassie-baxter模态。因此，多级微结构的合理设计受到了更广泛的关注。多级微结构具有几何上的纤维丛构型。目前，多级微结构的设计普遍采用仿生方法，该方法通过模拟自然界植物表面上的纹理微结构获取相应的几何构型。该仿生学设计方法缺乏严格的数学模型，并且设计效率较低。因此，发展严格、高效的多级微结构逆向设计方法具有很强的必要性。
[0005]
在已有的结构逆向设计方法中，拓扑优化被认为是最为有效的。目前，密度法和水平集法是两种普遍采用的拓扑优化实现方法。相比于水平集法，密度法在初值依赖性、收敛效率和多约束处理等方面具有显著的优势。因此，以下采用密度法给出并实现疏水多级微结构的纤维丛拓扑优化方法。


技术实现要素：

[0006]
要解决的技术问题
[0007]
为了解决现有技术中表面微结构设计凭借研究人员直觉和仿生手段的局限性，本发明提出一种通过对固体表面多级微结构支撑的液汽界面建模，计算与该界面凸起体积极小值对应的微结构几何构型，实现疏水多级微结构的纤维丛拓扑优化方法。
[0008]
技术方案
[0009]
一种疏水多级微结构的纤维丛拓扑优化方法，其特征在于步骤如下：
[0010]
步骤1：设计变量定义，包括主机构设计变量和二级结构设计变量：
[0011]
所述的主结构设计变量：定义取值于[0,1]的设计变量z
m
，用以表达主结构表面σ；为获得光滑的主结构表面，对z
m
进行以下滤波操作：
[0012][0013]
对滤波处理后的主结构设计变量z
f
做进一步的滤波，以获得主结构上点的竖坐标z
s
并控制主结构高度：
[0014][0015]
其中，ω是z
m
的定义域，对应于归一化的固体表面周期性剖分单元；是二维平面内的梯度算子；i和j分别是x和y轴的方向向量；r
m
是滤波半径，其值为常数；b
z
是控制主结构高度的参数，该参数取非负值；主结构表面σ上的法向量为主结构表面σ称为主结构流形，该流形为二级结构的可变设计区域；
[0016]
所述的二级结构设计变量：定义取值于[0,1]的变量γ为主结构上二级结构的设计变量；二级结构分布于主结构表面σ上；为控制二级结构的特征尺寸和移除其中的灰度区域，将γ进行如下的滤波和投影处理：
[0017][0018][0019]
其中，是σ上的切向梯度算子；γ
p
为投影后的二级结构设计变量，将其命名为材质密度；r
f
是滤波半径，其值为常数；ξ和β是投影参数，其值由数值实验获得；σ上的切向梯度算子与ω上的梯度算子之间的关系为
[0020][0021]
因为依赖于z
s
，该算子对z
s
的一阶变分满足
[0022][0023]
其中，是z
s
的试函数；
[0024]
步骤2：基于以上的设计变量定义，所给出的拓扑优化方法将通过同时演化两套设计变量实现主结构和二级结构的最佳匹配；几何上，最佳匹配后的两套设计变量构成纤维丛(σ
×
γ
p
(σ),σ,proj1,γ
p
(σ))，其中σ为该纤维丛的底流形，γ
p
:σ
→
[0,1]为纤维， proj1为满足以下条件的自然映射：
[0025][0026]
因此，所给出的是一种纤维丛拓扑优化方法，并且该方法通过最佳匹配所定义的主结构和二级结构设计变量实现；
[0027]
当疏水多级微结构上支撑的液汽界面施加的压强逐渐增大时，处于的 cassie-baxter模态的液汽界面曲率随之逐渐增大，cassie-baxter模态的液汽界面最终完全由二级结构支撑，此时对应的亚稳态称为cassie-baxter模态的终态；此后， cassie-baxter模态开始向wenzel模态转变；因此，疏水多级微结构的拓扑优化目标是提高cassie-baxter模态的终态稳定性；
[0028]
基于不相溶两相流体的界面自由能极小原理，多级微结构支撑的cassie-baxter终态液汽界面是常平均曲率曲面；该曲面由young-laplace方程描述，其物理意义在于液汽界面上毛细压和表面张力之间的平衡；在归一化尺度下，无量纲young-laplace 方程为
[0029][0030]
其中，是液汽界面相对于σ的无量纲化位移；d0是液汽界面原始位移的数量级；κ是对应于σ平均曲率的静压强，将其命名为底流形压强；是无量纲化的表面张力；σ和p分别为表面张力和两相界面处的静压强；常量是无量纲化的液体表面张力；为保证无量纲young-laplace方程解的唯一性，液汽界面的边界设置为由滤波方程2可得进而底流形压强κ可变换为如下形式：
[0031][0032]
在疏水多级微结构拓扑优化中，对无量纲化表面张力进行基于材质密度的插值，并通过材质密度对底流形压强κ进行惩罚，具体如下：
[0033][0034][0035]
其中，是液固界面的无量纲化表面张力；q是用于调整材质插值和惩罚方程凹凸性的参数；p
κ
是底流形压强κ的惩罚因子，其最大和最小值分别为p
κ,max
和p
κ,min
；在表面张力的材质插值方程中，的理论值为正无穷大；在数值执行过程中，该参数取为足够大的正数，以同时保证数值计算的收敛性和液固界面的逼近精度；惩罚因子 p
κ
的作用是消除二级结构上的底流形压强，并在液汽界面上保持该压强；
[0036]
基于上述底流形压强变换、表面张力材质插值和底流形压强惩罚，方程8中的无量纲young-laplace方程可变换为
[0037][0038]
对cassie-baxter模态的终态稳定性采用如下最小二乘形式进行度量
[0039][0040]
其中，|σ|是底流形σ的面积，具体表达式为
[0041]
|σ|＝∫
σ
1dσ
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(14)
[0042]
方程13中的终态稳定性度量等价于单位底流形面积上液汽界面凸起体积的平方；因此，疏水多级微结构的拓扑优化目标设置为极小化方程13中的终态稳定性度量；
[0043]
步骤3：建立如下拓扑优化问题，以实现主结构流形和二级结构版图的最佳匹配，进而实现纤维丛结构拓扑优化：
[0044][0045]
for fiber bundle(σ
×
γ
p
(σ),σ,proj1,γ
p
(σ)),
[0046]
to minimizewithconstrained by
[0047][0048]
在方程15中，f
d
是二级结构于主结构流形上的占空比，其表达式为
[0049][0050]
f0是设计者给定的占空比；
[0051]
为求解方程15中的拓扑优化问题，采用基于梯度信息的迭代算法，其中的液汽界面终态稳定性度量和二级结构占空比的梯度信息可通过伴随分析法获得；基于伴随分析法，液汽界面终态稳定性度量的敏度为：
[0052][0053]
其中δ为一阶变分算子，δ(j|σ|2)为j|σ|2的伴随敏度：
[0054][0055]
其中，δz
m
和δγ分别是z
m
和γ的一阶变分；z
fa
和γ
fa
分别是z
f
和γ
f
的伴随变量；和分别是和的对偶空间；和分别是σ和ω上无穷光滑函数空间和的闭空间；方程18中的伴随变量由以下伴随方程求得：
[0056]
findwithonsatisfying
[0057][0058]
findsatisfying
[0059][0060]
findwithγ
fa
＝0 onsatisfying
[0061][0062]
findwith z
sa
＝0 onsatisfying
[0063][0064]
findwith z
fa
＝0 onsatisfying
[0065][0066]
其中，κ
a
和z
sa
分别是κ和z
s
的伴随变量；和分别是κ
a
，γ
a
，z
sa
和z
fa
的试函数；方程17中，|σ|的伴随敏度为：
[0067][0068]
其中的伴随变量z
fa
由以下伴随方程求得：
[0069]
findwith z
sa
＝0 onsatisfying
[0070][0071]
findwith z
fa
＝0 onsatisfying
[0072][0073]
二级结构占空比的敏度为：
[0074][0075]
其中，δ|σ|由方程24，25和26求得；f
d
|σ|的伴随敏度为
[0076][0077]
方程28中的伴随变量γ
fa
和z
fa
由以下方程求得：
[0078]
findwithγ
fa
＝0 onsatisfying
[0079][0080]
findwith z
sa
＝0 onsatisfying
[0081][0082]
findwith z
fa
＝0 onsatisfying
[0083][0084]
在进行伴随分析后，采用以下迭代步骤求解上述变分问题：
[0085]
(a)求解滤波方程1，2和3，并用方程4投影滤波后的二级结构设计变量；
[0086]
(b)由方程16计算当前设计变量对应的二级结构占空比；
[0087]
(c)求解无量纲young-laplace方程12，并由方程13和14计算液汽界面终态稳定性度量；
[0088]
(d)计算方程18中的伴随敏度，其中的伴随变量由方程19，20，21，22，23求得；
[0089]
(e)计算方程24中的伴随敏度，其中的伴随变量由方程25，26求得；
[0090]
(f)由方程17计算液汽界面终态稳定性度量的伴随敏度；
[0091]
(g)计算方程28中的伴随敏度，其中的伴随变量由方程29，30，31求得；
[0092]
(h)计算方程27中的二级结构占空比伴随敏度；
[0093]
(i)演化主结构和二级结构的设计变量；
[0094]
(j)判断是否满足收敛条件，若满足则终止迭代；若不满足则返回(a)。
[0095]
本发明技术方案更进一步的说：步骤2中在底流形压强的惩罚因子中，p
κ,max
和
p
κ,min
的取值分别为1和0。
[0096]
本发明技术方案更进一步的说：步骤2中f0的取值范围为(0,1)。
[0097]
本发明技术方案更进一步的说：步骤2中终止迭代的条件为：(1)迭代次数达到最大值315；(2)连续5次迭代的设计目标值变化量与5次迭代目标平均值的比值和占空比与f0的偏差都小于10-3
。
[0098]
本发明技术方案更进一步的说：步骤2中迭代过程中，β的初始值取为1，之后经每30次迭代其值翻倍一次；ξ取值0.5。
[0099]
有益效果
[0100]
本发明提出的一种疏水固体表面上的多级微结构几何构型的拓扑优化方法，该方法克服了目前表面微结构设计凭借研究人员直觉和仿生手段的局限性，提高了设计方法的适用性、灵活性和效率。由于该方法的主结构设计变量定义了二级结构的可变设计区域，本方法是一种具有可变设计区域的拓扑优化方法。由于多级微结构通过最佳匹配主结构流形和二级结构版图获得，且主结构流形和二级结构版图构成几何上的纤维丛，本方法是一种纤维丛拓扑优化方法。
附图说明
[0101]
图1为疏水多级微结构图：1a，1b分别为疏水多级微结构的立体和截面示意图。
[0102]
图2为归一化尺度上多级微结构底流形和二级结构版图，及对应的纤维丛沉降和缩放操作示意图。
[0103]
图3为主结构设计变量的两次滤波操作示意图。
[0104]
图4为二级结构设计变量的滤波和投影操作示意图。
[0105]
图5为疏水多级微结构上cassie-baxter模态的受压演化过程示意图。
[0106]
图6为正三角形相关图：图6a，6b1，6b2，6c1，6c2分别为固体表面的正三角形周期性划分，划分单元上获得的疏水多级微结构纤维丛构型的立体图、俯视图，沉降操作后的纤维丛构型立体图、俯视图。
[0107]
图7为正四边形相关图：图7a，7b1，7b2，7c1，7c2分别为固体表面的正四边形周期性划分，划分单元上获得的疏水多级微结构纤维丛构型立体图、俯视图，沉降操作后的纤维丛构型立体图、俯视图。
[0108]
图8为正六边形相关图：图8a1，8c1，8c2，8d1，8d2分别为固体表面的正六边形周期性划分，划分单元上获得的疏水多级微结构纤维丛构型立体图、俯视图，沉降操作后的纤维丛构型立体图、俯视图。
[0109]
图9为30
°
顶角等腰三角形相关图：图9a，9b1，9b2，9c1，9c2分别为固体表面30
°
顶角等腰三角形的周向对称性划分，划分单元上获得的疏水多级微结构纤维丛构型立体图、俯视图，沉降操作后的纤维丛构型立体图、俯视图。
[0110]
图10为15
°
顶角等腰三角形相关图：图10a，10b1，10b2，10c1，10c2分别为固体表面15
°
顶角等腰三角形的周向对称性划分，划分单元上获得的疏水多级微结构纤维丛构型立体图、俯视图，沉降操作后的纤维丛构型立体图、俯视图。
[0111]
图11为等边四边形相关图：图11a，11b1，11b2，11c1，11c2分别为固体表面的等边四边形手性划分，划分单元上获得的疏水多级微结构纤维丛构型立体图、俯视图，沉降操作
后的纤维丛构型立体图、俯视图。
[0112]
图12为等边九边形相关图：图12a，12b1，12b2，12c1，12c2分别为固体表面的等边九边形手性划分，划分单元上获得的疏水多级微结构纤维丛构型立体图、俯视图，沉降操作后的纤维丛构型立体图、俯视图。
[0113]
图13为菱形相关图：图13a，13b1，13b2，13c1，13c2分别为基于两个菱形单元的固体表面准周期性划分，划分单元上获得的疏水多级微结构纤维丛构型立体图、俯视图，沉降操作后的纤维丛构型立体图、俯视图。
具体实施方式
[0114]
现结合实施例、附图对本发明作进一步描述：
[0115]
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及具体实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，而不构成对本发明的限制。
[0116]
本发明疏水多级微结构的纤维丛拓扑优化方法包括以下步骤：
[0117]
固体表面微结构通常都是周期性分布的，并且通过周期性排布一个结构单元可获得疏水固体表面上的周期性微结构。因此，多级微结构的拓扑优化针对结构单元进行 (图1a)。多级微结构的结构单元由主结构和二级结构构成。几何上，这种结构对应着由底流形和纤维构成的纤维丛(图1b)。
[0118]
多级微结构上的液汽界面采用young-laplace描述，其解具有比例缩放的性质，因而以下在归一化尺度上进行多级微结构的拓扑优化，以保证所得结构在一定尺度范围内的适用性和数值求解精度；进而，多级微结构可通过沉降和缩放拓扑优化所得的归一化纤维丛构型获得(图2)。
[0119]
定义如下的多级微结构拓扑优化设计变量：
[0120]
主结构设计变量：定义取值于[0,1]的设计变量z
m
，用以表达主结构表面σ。为获得光滑的主结构表面，对z
m
进行以下滤波操作：
[0121][0122]
对滤波处理后的主结构设计变量z
f
做进一步的滤波，以获得主结构上点的竖坐标z
s
并控制主结构高度：
[0123][0124]
其中，ω是z
m
的定义域，对应于归一化的固体表面周期性剖分单元；是二维平面内的梯度算子；i和j分别是x和y轴的方向向量；r
m
是滤波半径，其值为常数；b
z
是控制主结构高度的参数，该参数取非负值；主结构表面σ上的法向量为
主结构表面σ称为主结构流形，该流形为二级结构的可变设计区域。上述主结构设计变量经两次滤波后获得主结构形貌的过程如图3所示。
[0125]
二级结构设计变量：定义取值于[0,1]的变量γ为主结构上二级结构的设计变量。二级结构分布于主结构表面σ上。为控制二级结构的特征尺寸和移除其中的灰度区域，将γ进行如下的滤波和投影处理：
[0126][0127][0128]
其中，是σ上的切向梯度算子；γ
p
为投影后的二级结构设计变量，将其命名为材质密度；r
f
是滤波半径，其值为常数；ξ和β是投影参数，其值由数值实验获得。σ上的切向梯度算子与ω上的梯度算子之间的关系为
[0129][0130]
因为依赖于z
s
，该算子对z
s
的一阶变分满足
[0131][0132]
其中，是z
s
的试函数。上述二级结构设计变量经滤波和投影后获得二级结构版图的过程如图4所示。
[0133]
基于以上的设计变量定义，所给出的拓扑优化方法将通过同时演化两套设计变量实现主结构和二级结构的最佳匹配。几何上，最佳匹配后的两套设计变量构成纤维丛 (σ
×
γ
p
(σ),σ,proj1,γ
p
(σ))，其中σ为该纤维丛的底流形，γ
p
:σ
→
[0,1]为纤维，proj1为满足以下条件的自然映射：
[0134][0135]
因此，所给出的是一种纤维丛拓扑优化方法，并且该方法通过最佳匹配所定义的主结构和二级结构设计变量实现。
[0136]
当如图1所示的疏水多级微结构上支撑的液汽界面施加的压强逐渐增大时，处于的cassie-baxter模态的液汽界面曲率随之逐渐增大。如图5所示，cassie-baxter模态的液汽界面最终完全由二级结构支撑，此时对应的亚稳态称为cassie-baxter模态的终态。此
后，cassie-baxter模态开始向wenzel模态转变。因此，疏水多级微结构的拓扑优化目标是提高cassie-baxter模态的终态稳定性。
[0137]
基于不相溶两相流体的界面自由能极小原理，多级微结构支撑的cassie-baxter终态液汽界面是常平均曲率曲面。该曲面由young-laplace方程描述，其物理意义在于液汽界面上毛细压和表面张力之间的平衡。在归一化尺度下，无量纲young-laplace 方程为
[0138][0139]
其中，是液汽界面相对于σ的无量纲化位移；d0是液汽界面原始位移的数量级；κ是对应于σ平均曲率的静压强，将其命名为底流形压强；是无量纲化的表面张力；σ和p分别为表面张力和两相界面处的静压强；常量是无量纲化的液体表面张力；为保证无量纲young-laplace方程解的唯一性，液汽界面的边界设置为由滤波方程2可得进而底流形压强κ可变换为如下形式：
[0140][0141]
在疏水多级微结构拓扑优化中，对无量纲化表面张力进行基于材质密度的插值，并通过材质密度对底流形压强κ进行惩罚，具体如下：
[0142][0143][0144]
其中，是液固界面的无量纲化表面张力；q是用于调整材质插值和惩罚方程凹凸性的参数；p
κ
是底流形压强κ的惩罚因子，其最大和最小值分别为p
κ,max
和p
κ,min
。在表面张力的材质插值方程中，的理论值为正无穷大；在数值执行过程中，该参数取为足够大的正数，以同时保证数值计算的收敛性和液固界面的逼近精度。在底流形压强的惩罚因子中，p
κ,max
和p
κ,min
的取值分别为1和0。惩罚因子p
κ
的作用是消除二级结构上的底流形压强，并在液汽界面上保持该压强。
[0145]
基于上述底流形压强变换、表面张力材质插值和底流形压强惩罚，方程8中的无量纲young-laplace方程可变换为
[0146][0147]
对图5所示的cassie-baxter模态的终态稳定性采用如下最小二乘形式进行度量
[0148][0149]
其中，|σ|是底流形σ的面积，具体表达式为
[0150]
|σ|＝∫
σ
1dσ
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(14) 方程13中的终态稳定性度量等价于单位底流形面积上液汽界面凸起体积的平方。因此，疏水多级微结构的拓扑优化目标设置为极小化方程13中的终态稳定性度量。
[0151]
综上所述，可建立如下拓扑优化问题，以实现主结构流形和二级结构版图的最佳匹配，进而实现纤维丛结构拓扑优化：
[0152][0153]
for fiber bundle(σ
×
γ
p
(σ),σ,proj1,γ
p
(σ)),
[0154]
to minimizewithconstrained by
[0155][0156]
在方程15中，f
d
是二级结构于主结构流形上的占空比，其表达式为
[0157][0158]
f0是设计者给定的占空比，其取值范围为(0,1)。
[0159]
为求解方程15中的拓扑优化问题，采用基于梯度信息的迭代算法，其中的液汽界面终态稳定性度量和二级结构占空比的梯度信息可通过伴随分析法获得。基于伴随分析法，液汽界面终态稳定性度量的敏度为：
[0160][0161]
其中δ为一阶变分算子。δ(j|σ|2)为j|σ|2的伴随敏度：
[0162][0163]
其中，δz
m
和δγ分别是z
m
和γ的一阶变分；z
fa
和γ
fa
分别是z
f
和γ
f
的伴随变量；和分别是和的对偶空间；和分别是σ和ω上无穷光滑函数空间和的闭空间。方程18中的伴随变量由以下伴随方程求得：
[0164]
findwithonsatisfying
[0165][0166]
findsatisfying
[0167][0168]
findwithγ
fa
＝0 onsatisfying
[0169][0170]
findwith z
sa
＝0 onsatisfying
[0171][0172]
findwith z
fa
＝0 onsatisfying
[0173][0174]
其中，κ
a
和z
sa
分别是κ和z
s
的伴随变量；和分别是κ
a
，γ
a
，z
sa
和z
fa
的试函数。方程17中，|σ|的伴随敏度为：
[0175][0176]
其中的伴随变量z
fa
由以下伴随方程求得：
[0177]
findwith z
sa
＝0 onsatisfying
[0178][0179]
findwith z
fa
＝0 onsatisfying
[0180][0181]
二级结构占空比的敏度为：
[0182][0183]
其中，δ|σ|由方程24，25和26求得；f
d
|σ|的伴随敏度为
[0184][0185]
方程28中的伴随变量γ
fa
和z
fa
由以下方程求得：
[0186]
findwithγ
fa
＝0 onsatisfying
[0187][0188]
findwith z
sa
＝0 onsatisfying
[0189][0190]
findwith z
fa
＝0 onsatisfying
[0191][0192]
在进行伴随分析后，采用以下迭代步骤求解上述变分问题：
[0193]
(a)求解滤波方程1，2和3，并用方程4投影滤波后的二级结构设计变量；
[0194]
(b)由方程16计算当前设计变量对应的二级结构占空比；
[0195]
(c)求解无量纲young-laplace方程12，并由方程13和14计算液汽界面终态稳定性度量；
[0196]
(d)计算方程18中的伴随敏度，其中的伴随变量由方程19，20，21，22，23 求得；
[0197]
(e)计算方程24中的伴随敏度，其中的伴随变量由方程25，26求得；
[0198]
(f)由方程17计算液汽界面终态稳定性度量的伴随敏度；
[0199]
(g)计算方程28中的伴随敏度，其中的伴随变量由方程29，30，31求得；
[0200]
(h)计算方程27中的二级结构占空比伴随敏度；
[0201]
(i)演化主结构和二级结构的设计变量；
[0202]
(j)判断是否满足收敛条件，若满足则终止迭代；若不满足则返回(a)，其中，终止迭代的条件为：(1)迭代次数达到最大值315；(2)连续5次迭代的设计目标值变化量与5次迭
代目标平均值的比值和占空比与f0的偏差都小于10-3
。
[0203]
在上述迭代过程中，β的初始值取为1，之后经每30次迭代其值翻倍一次，ξ取值0.5。
[0204]
在一个实施例中，采用上述疏水多级微结构的纤维丛拓扑优化方法，将相关参数设置为表1中的值，可分别得到如图6，7，8所示基于正三角形、正四边形和正六边形的三种固体表面轴对称周期性划分单元上的微结构纤维丛构型，如图9，10所示基于顶角分别为30
°
和15
°
等腰三角形的两种固体表面周向对称周期性划分单元上的微结构纤维丛构型，如图11，12所示基于等边四边形和等边九边形的两种固体表面手性划分单元上的微结构纤维丛构型，如图13所示基于两个菱形的固体表面准周期性划分单元上的微结构纤维丛构型。
[0205][0206]
表1
[0207]
基于所介绍的纤维丛拓扑优化方法，通过表1中的参数，求解方程15中的拓扑优化问题，可得如图6，7，8，9，10，11，12，13所示分别具有轴对称性、周向对称性、手性和准周期性的固体表面划分的单元多级微结构几何构型，而且它们可以使用双光子光刻等微纳三维打印类加工工艺制造出来。在缩放操作中，周期性单元的特征尺寸是其缩放系数；在沉降操作中，拉伸距离应该比液汽界面的深度要大，以避免液汽界面接触到微结构的底部而导致的cassie-baxter模态崩溃。如图6，7，8，9，10，11， 12，13所示，多级微结构单元的纤维丛构型具有峰状的主结构和密布于主结构顶端的二级结构。
[0208]
以上所述本发明的具体实施方式，并不构成对本发明保护范围的限定。任何根据本发明的技术构思所作出的各种其他相应的改变与变形，均应包含在本发明权利要求的保护范围内。