CN103218488A - 一种二维液-液声子晶体拓扑优化方法 - Google Patents
一种二维液-液声子晶体拓扑优化方法 Download PDFInfo
- Publication number
- CN103218488A CN103218488A CN2013101261945A CN201310126194A CN103218488A CN 103218488 A CN103218488 A CN 103218488A CN 2013101261945 A CN2013101261945 A CN 2013101261945A CN 201310126194 A CN201310126194 A CN 201310126194A CN 103218488 A CN103218488 A CN 103218488A
- Authority
- CN
- China
- Prior art keywords
- liquid
- centerdot
- phonon crystal
- pixel
- unit cell
- 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.)
- Pending
Links
Images
Landscapes
- Liquid Crystal (AREA)
Abstract
一种二维液-液声子晶体拓扑优化方法,将声子晶体原胞分割成N×N正方形像素型结构;然后,根据二维液-液声子晶体色散关系所满足的声波波动方程,开发出任意材料布局的二维液-液像素型声子晶体能带快速计算的平面波展开法程序,计算其带隙;最后利用遗传优化算法,根据对带隙的要求,搜寻二维液-液声子晶体原胞最优的材料布局。过程包括:输入待求解的参数;初始化;计算个体适应度;依次执行选择、交叉和变异遗传操作生成下一代种群,使种群向前进化;检验种群是否满足停止条件。这种拓扑优化方法使声子晶体的可设计性变强,得到具有最优带隙特性的新颖的声子晶体结构;同时减少计算时间,提高计算效率,使所设计的声子晶体达到最好的技术经济性能。
Description
技术领域
本发明属声带隙材料设计领域,涉及一种声带隙材料拓扑优化设计方法,特别涉及一种基于遗传算法和快速平面波展开法(PWE),对二维液-液声子晶体原胞材料分布进行拓扑优化设计,以获得所需要的声带隙特征。
背景技术
声子晶体是指具有不同声学特性的材料按周期性复合在一起、具有声波带隙的周期复合材料。当声波在受到材料常数的周期性调制时,可能会产生声带隙,即一定频率范围的声波的传播被抑制或禁止。声子晶体的这种特性具有极大的理论价值,在无源隔音、声滤波器等新型声学功能材料等方面具有广泛的应用前景。目前,声子晶体带隙特性的研究已经成为一个新的热点。
传统的声子晶体设计思路一般是:在特定的晶格类型条件下,采用有限的几种对称形状(如二维问题中,采用圆柱、正方柱、椭圆柱等)作为原胞的散射体,通过调整这些散射体的几何参数来打开声子晶体带隙或改善带隙的特性,进而经验性地提炼可能的设计规律,然后指导性地寻找更好的声子晶体新结构。
然而,在计算声子晶体能带特性之前,必须事先设定声子晶体的原胞结构。但对何种形式的声子晶体原胞结构具有所期望的最优带隙,则是一个难以预测的问题,这在很大程度上束缚了人们对声子晶体的应用。
发明内容
本发明的目的在于,为了克服现有的声子晶体设计分析方法的不足,采用拓扑优化技术,根据对二维液-液声子晶体带隙的要求,自动寻找对应的声子晶体材料最优布局,得到具有最优带隙特性的新颖的声子晶体结构,使其达到最好的技术经济性能。
本发明解决其技术问题所采用的技术方案是:用改进的快速平面波展开法计算二维液-液声子晶体的色散关系,获取相应的带隙值;然后应用遗传优化算法,根据带隙所要达到目标,搜索声子晶体最优材料拓扑布局。其特点是包括以下步骤:
1)输入待求解的问题参数:如材料常数、平面波个数、像素尺度N、种群大小、交叉率及变异率等。
2)初始化:随机生成初始种群,用N×N矩阵二进制数对原胞进行编码,构成染色体(个体)。
3)计算个体适应度:采用改进的快速平面波展开法,计算每个遗传个体对应的声子晶体原胞的色散关系,得到对应的禁带数值。根据优化的目标,用声子晶体禁带构造目标函数,然以目标函数值来度量遗传个体的适应度。
4)依次执行选择、交叉和变异遗传操作生成下一代种群,使种群向前进化,不断更新种群。其中,选择是根据遗传个体适应度值的大小,采用保留精英的轮盘赌机制;交叉操作采用平均交叉方式,随机选择染色体实施行交叉或列交叉,以增强算法的全局搜索能力;个体变异采用位变异机制。由于“优胜劣汰,适者生存”自然选择和遗传进化机理,最优秀个体的拓扑结构就越来越逼近目标要求的声子晶体原胞最优拓扑布局。
5)检验种群是否满足停止条件(例如固定的进化代数或种群是否稳定)。如果满足则输出最优结果,退出;否则返回步骤3)。
6)所述步骤3)中遗传个体对应的声子晶体原胞为正方晶格,晶格常数为a。声子晶体原胞被分割为N×N(N为偶数)正方形像素型结构。将每个像素随机配置所选定的液体声学材料(仅考虑两种材料,故采用二进制字符串形式表示)。声子晶体原胞设定为关于x、y轴对称和关于z轴90°旋转对称模型(如图2所示)。此时,声子晶体的结构设计问题等价于像素中声学材料的选择问题。显然,当N足够大时原则上可以表达任意形状的结构,材料边界也更光滑,但此时搜索规模往往呈指数形式增长,而且像素粒度过小也给工艺制备带来了难度。本发明综合考虑运算能力和制备工艺的限制,研究像素粒度适中的声子晶体结构拓扑优化设计问题。
7)所述步骤3)中二维液-液像素型声子晶体色散关系所应满足的本构方程为:
其中,ρ为材料密度;C11=ρcl 2是液体中声波弹性常数,cl是液体中声波声速。求解该方程组,即可得到二维液-液像素型声子晶体带隙特性。
8)所述步骤3)中用改进的快速平面波展开法求解方程组(1)时,组元材料常数g(ρ-1、C11 -1)的傅里叶变换g(G)至关重要。而对于像素型声子晶体而言,可以根据位移特性来计算任意分布的组元材料常数的傅里叶变换g(G)。
设Pa为原胞中填充介质ga的像素集合。首先计算中心点像素P0∈Pa(图3)材料常数g0的傅里叶变换g0(G)为:
根据位移特性,任一填充ga介质的像素Pr∈Pa的傅里叶变换为:
gr(G)=g0(G)eiG·r (3)
由(3)式可以计算像素型声子晶体材料常数分布的傅里叶变换g(G)为:
设 原胞结构的像素材料常数分布可表示为δ(r)=[δ(r1),δ(r2),…δ(rk),…,δ(rL)],rk表示原胞中第k个像素的中心点位置,L表示声子晶体原胞分割后总的像素数。则(4)式等价于
g(G)=g0(G)e(G)·δ(r) (5)
由于e(G)可以离线计算,(5)式表明,给定任意原胞结构的δ(r)后,既可以采用加法计算g(G)。将g(G)离线计算并存储,供每次带隙计算调用,从而大大降低运算量,提高运算速度。
本发明的有益效果是:可以根据对二维液-液声子晶体带隙的要求,自动快速寻找具有全局最优的液-液声子晶体拓扑形状,得到新颖的声子晶体结构;摆脱了传统的经验设计思路,实现了根据带隙需要主动设计声子晶体结构的目标,使声子晶体的可设计性变强;同时减少了计算时间,提高了计算效率,并使所设计的声子晶体达到最好的技术经济性能。
附图说明
图1是本发明基于遗传算法和快速平面波展开法对二维液-液声子晶体结构拓扑优化设计流程框图;
图2是本发明中关于x、y轴对称和关于z轴90°旋转对称的声子晶体原胞10×10像素结构示意图和基因表达图;
图3是本发明任一填充ga材料像素单元的傅里叶变换计算时的位移变换图;
图4是本发明实例最优声子晶体3×3原胞图;
图5是本发明实例最优声子晶体能带图;
图6是本发明实例遗传算法进化过程中的种群平均适应度及最优个体适应度变化图;
图7是本发明实例遗传算法进化过程中的最优个体第三能级和第二能级间的带隙上下限及带隙宽度变化图。
具体实施方式
以下结合说明书附图对本发明的具体实施例加以说明:
本实施例研究由两种液体材料构成的二维正方晶格液-液声子晶体,晶格常数a=1m。材料A为水,其密度ρA=1025kg/m3、声波波速clA=1531m/s;材料B为水银,其密度ρB=13500kg/m3、声波波速为clB=1450m/s。
本实施例包括以下步骤:
1)输入待求解的问题参数:平面波个数为221个平面波、像素尺度N=20、种群规模为60、交叉率为0.5、变异率为0.1~0.8。
2)初始化:随机生成初始种群,用20×20矩阵二进制数对原胞进行编码,构成染色体(个体)。
3)计算个体适应度:声子晶体原胞被分割为20×20正方形像素型结构。将每个像素随机配置所选定的材料(水或水银),采用二进制字符串形式表示,1和0分别表示材料水和材料水银。本发明综合考虑运算能力和制备工艺的限制,研究声子晶体原胞关于x、y轴对称和关于z轴90°旋转对称条件下,声子晶体结构拓扑优化设计问题,遗传个体二进制字符串长度为
采用改进的快速平面波展开法,计算每个遗传个体对应的声子晶体原胞的色散关系,得到对应的禁带数值。拓扑优化的目标设为第三能级和第二能级间的带隙相对值最大,即:
minω3(k)为第三能级本征频率的最小值、maxω2(k)为第二能级本征频率的最大值。然后以此目标函数值来度量遗传个体的适应度。
4)依次执行选择、交叉和变异遗传操作生成下一代种群,使种群向前进化,不断更新种群。
5)检验种群是否满足停止条件(种群是否稳定)。如果满足则输出最优结果,退出;否则返回步骤3)。优化流程见图1。
经过100次循环迭代,得到二维水-水银声子晶体最优拓扑结构,其3×3原胞见图4。该最优声子晶体结构为复式格子的正方晶格结构,散射体材料为水(图4中的黑色部分),基体材料为水银(图4中的白色部分),散射体为两种不同尺寸的具有倒角的正方形,散射体填充率为34%;最优声子晶体的能带图见图5,图中横轴坐标为波矢,纵轴坐标为频率(单位:kHz),ΓXM三角形为不可约布里渊区,求色散曲线时,波矢沿Γ→X→M→Γ方向扫描。由图5可见,第三能级和第二能级间的完全带隙相对宽度为最大,此带隙的上边缘最小值为1589.3Hz,带隙的下边缘最大值为594.2Hz,带隙宽度为995.1Hz;种群的平均适应度和最优个体的适应度随优化过程而变化(图6),图中横轴坐标为遗传优化进化次数(单位:次),纵轴坐标为第三和第二能级间带隙相对值表示的适应度值。种群的平均适应度为种群每个个体适应度的平均值,适应度最大的个体为该种群的最优个体。由图6可以发现,遗传算法优化进化前期,进化速度很快,而在后期,特别是接近最优解时,速度很慢,最优解所对应的个体适应度为0.9115;进化过程中声子晶体第三能级和第二能级间的完全带隙上下限(分别指带隙上边缘最小值和带隙下边缘最大值)及带隙宽度(带隙上限与带隙下限的差值)变化见图7,图中横轴坐标为遗传优化进化次数(单位:次),纵轴坐标为频率(单位:kHz)。由图7可以发现,随着优化进化的进行,第三能级和第二能级间的完全带隙的上边缘最小值在变大,下边缘最大值在变小,带隙的宽度在变大。
由结果表明:随着迭代的进行,种群向前进化,最后得到了新颖的最优二维水-水银声子晶体,其第三能级和第二能级间带隙完全打开,相对带隙值达到最大,实现了优化的目标;这种全局拓扑优化方法超越了经验设计的局限,实现了液-液声子晶体的主动设计。
最后应说明的是:以上实施例仅用以说明本发明而并非限制本发明所描述的技术方案;因此,尽管本说明书参照上述的各个实施例对本发明已进行了详细的说明,但是,本领域的普通技术人员应当理解,仍然可以对本发明进行修改或等同替换;而一切不脱离发明的精神和范围的技术方案及其改进,其均应涵盖在本发明的权利要求范围当中。
Claims (3)
1.一种二维液-液声子晶体拓扑优化方法,用改进的快速平面波展开法(PWE)计算二维液-液声子晶体的色散关系,获取相应的带隙值;然后应用遗传优化算法,根据带隙所要达到目标,搜索声子晶体最优材料拓扑布局;其特征在于:包括以下步骤:
1)输入待求解的参数:包括,材料常数、平面波个数、像素尺度N、种群大小、交叉率及变异率;
2)初始化:随机生成初始种群,用N×N矩阵二进制数对原胞进行编码,构成染色体(个体);
3)计算个体适应度:二维液-液像素型声子晶体的色散关系所满足的本构方程为:
其中,ρ为材料密度,C11=ρcl 2是液体中声波弹性常数,cl是液体中声波声速;采用改进的快速平面波展开法求解方程组(1)时,任意分布的组元材料常数的傅里叶变换g(G)可以根据位移特性来计算;设Pa为原胞中填充介质ga的像素集合,首先计算中心点像素P0∈Pa材料常数g0的傅里叶变换g0(G)为:
gr(G)=g0(G)eiG·r (3)
由(3)式可以计算像素型声子晶体材料常数分布的傅里叶变换g(G)为:
设 原胞结构的像素材料常数分布可表示为δ(r)=[δ(r1),δ(r2),…δ(rk),…,δ(rL)],rk表示原胞中第k个像素的中心点位置,L表示声子晶体原胞分割后总的像素数;则(4)式等价于
g(G)=g0(G)e(G)·δ(r) (5);
将g(G)离线存储,供每次带隙计算调用;采用改进的快速平面波展开法,计算每个遗传个体对应的声子晶体原胞的色散关系,得到对应的禁带数值;根据优化的目标,用二维液-液像素型声子晶体禁带构造目标函数,然以目标函数值来度量遗传个体的适应度;
4)依次执行选择、交叉和变异遗传操作生成下一代种群,使种群向前进化;
5)检验种群是否满足停止条件,满足则输出最优结果,退出;否则返回步骤3)。
2.根据权利要求1所述的一种二维液-液声子晶体拓扑优化方法,其特征在于:所述的遗传个体对应的声子晶体原胞为正方晶格;声子晶体原胞被分割为偶数的N×N正方形像素型结构;将每个像素随机配置所选定的声学材料;声子晶体原胞设定为关于x、y轴对称和关于z轴90°旋转对称模型。
3.根据权利要求1所述的一种二维液-液声子晶体拓扑优化方法,其特征在于:拓扑优化的目标为带隙相对值最大,即:
minωn+1(k)为第n+1能级本征频率的最小值,maxωn(k)为第n能级本征频率的最大值,然后以此目标函数值来度量遗传个体的适应度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013101261945A CN103218488A (zh) | 2013-04-12 | 2013-04-12 | 一种二维液-液声子晶体拓扑优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013101261945A CN103218488A (zh) | 2013-04-12 | 2013-04-12 | 一种二维液-液声子晶体拓扑优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103218488A true CN103218488A (zh) | 2013-07-24 |
Family
ID=48816269
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2013101261945A Pending CN103218488A (zh) | 2013-04-12 | 2013-04-12 | 一种二维液-液声子晶体拓扑优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103218488A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103793555A (zh) * | 2014-01-02 | 2014-05-14 | 南通大学 | 用于木工板件零件的智能排样优化方法 |
CN104464715A (zh) * | 2014-11-24 | 2015-03-25 | 广东工业大学 | 一种声子晶体分束器 |
CN104485098A (zh) * | 2014-11-24 | 2015-04-01 | 广东工业大学 | 一种调节声子晶体分束效率的方法 |
CN108229090A (zh) * | 2017-12-29 | 2018-06-29 | 西北工业大学 | 一种晶体结构的预测方法 |
GB2561685A (en) * | 2017-02-24 | 2018-10-24 | Ford Global Tech Llc | Design and manufacture of vibration damping structures |
CN117133385A (zh) * | 2023-07-27 | 2023-11-28 | 中南大学 | 一种弹性超材料非梯度多目标拓扑优化方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101118611A (zh) * | 2007-09-07 | 2008-02-06 | 北京航空航天大学 | 基于遗传算法的业务过程模型资源配置优化方法 |
US20090271341A1 (en) * | 2004-12-01 | 2009-10-29 | Matsushita Electric Industrial Co., Ltd. | Optimization processing method using a distributed genetic algorithm |
CN101777346A (zh) * | 2010-01-21 | 2010-07-14 | 北京工业大学 | 二维阿基米德格子声带隙材料 |
CN101794573A (zh) * | 2010-01-21 | 2010-08-04 | 北京工业大学 | 二维串联结构隔音声带隙材料 |
CN102314971A (zh) * | 2010-07-05 | 2012-01-11 | 田多贤 | 基于声子物理属性的高温超导材料制作方法及该方法制作的高温超导材料 |
-
2013
- 2013-04-12 CN CN2013101261945A patent/CN103218488A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090271341A1 (en) * | 2004-12-01 | 2009-10-29 | Matsushita Electric Industrial Co., Ltd. | Optimization processing method using a distributed genetic algorithm |
CN101118611A (zh) * | 2007-09-07 | 2008-02-06 | 北京航空航天大学 | 基于遗传算法的业务过程模型资源配置优化方法 |
CN101777346A (zh) * | 2010-01-21 | 2010-07-14 | 北京工业大学 | 二维阿基米德格子声带隙材料 |
CN101794573A (zh) * | 2010-01-21 | 2010-08-04 | 北京工业大学 | 二维串联结构隔音声带隙材料 |
CN102314971A (zh) * | 2010-07-05 | 2012-01-11 | 田多贤 | 基于声子物理属性的高温超导材料制作方法及该方法制作的高温超导材料 |
Non-Patent Citations (8)
Title |
---|
BIN WU等: "BAND STRUCTURES OF 2-D SONIC CRYSTALS WITH SQUARE LATTICES BASED FINITE ELEMENT METHOD", 《PIEZOELECTRICITY,ACOUSTIC WAVES AND DEVICE APPLICATIONS(SPAWDA),2011 SYMPOSIUM ON》 * |
BIN WU等: "Research on two-dimensional phononic crystal with magnetorheological material", 《ULTRASONICS SYMPOSIUM,2008.IUS 2008.IEEE 》 * |
ZHAO HUANYU等: "Study on the properties of sound isolation for acoustic band gaps materials with tandem structure", 《2010 INTERNATIONAL CONFERENCE ON ELECTRICAL AND CONTROL ENGINEERING》 * |
温激鸿: "声子晶体振动带隙及减振特性研究", 《中国博士学位论文全文数据库工程科技Ⅰ辑 》 * |
董浩文等: "基于有限元法和遗传算法的二维声子晶体的优化设计", 《北京力学会第19届学术年会论文集》 * |
郁殿龙: "基于声子晶体理论的梁板类周期结构振动带隙特性研究", 《中国博士学位论文全文数据库工程科技Ⅱ辑》 * |
钟会林 等: "遗传算法在二维声子晶体带隙优化中的应用", 《物理学报》 * |
龚春娟等: "快速遗传算法优化计算二维光子晶体", 《半导体学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103793555A (zh) * | 2014-01-02 | 2014-05-14 | 南通大学 | 用于木工板件零件的智能排样优化方法 |
CN104464715A (zh) * | 2014-11-24 | 2015-03-25 | 广东工业大学 | 一种声子晶体分束器 |
CN104485098A (zh) * | 2014-11-24 | 2015-04-01 | 广东工业大学 | 一种调节声子晶体分束效率的方法 |
CN104464715B (zh) * | 2014-11-24 | 2017-12-29 | 广东工业大学 | 一种声子晶体分束器 |
GB2561685A (en) * | 2017-02-24 | 2018-10-24 | Ford Global Tech Llc | Design and manufacture of vibration damping structures |
US10737478B2 (en) | 2017-02-24 | 2020-08-11 | Ford Global Technologies, Llc | Manufacture of vibration damping structures |
GB2561685B (en) * | 2017-02-24 | 2022-02-23 | Ford Global Tech Llc | Design and manufacture of vibration damping structures |
CN108229090A (zh) * | 2017-12-29 | 2018-06-29 | 西北工业大学 | 一种晶体结构的预测方法 |
CN108229090B (zh) * | 2017-12-29 | 2021-11-09 | 西北工业大学 | 一种晶体结构的预测方法 |
CN117133385A (zh) * | 2023-07-27 | 2023-11-28 | 中南大学 | 一种弹性超材料非梯度多目标拓扑优化方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103218488A (zh) | 一种二维液-液声子晶体拓扑优化方法 | |
CN103218529A (zh) | 一种二维固-固声子晶体xy模态拓扑优化方法 | |
CN103246767B (zh) | 基于遗传算法和有限元法的液-液声子晶体拓扑优化方法 | |
Dou et al. | Optimization of wind turbine yaw angles in a wind farm using a three-dimensional yawed wake model | |
Dong et al. | Topology optimization of two-dimensional asymmetrical phononic crystals | |
CN103760907B (zh) | 一种基于布谷鸟搜索算法的水下潜器三维路径规划方法 | |
du Pont | Dune morphodynamics | |
CN101599104A (zh) | 一种航空涡轮发动机叶片颤振边界的模拟方法 | |
Neshat et al. | New insights into position optimisation of wave energy converters using hybrid local search | |
CN103324786B (zh) | 基于遗传算法的二维固相声子晶体xy模带隙优化方法 | |
Bernard et al. | High-order h-adaptive discontinuous Galerkin methods for ocean modelling | |
CN103310046B (zh) | 基于超椭圆曲线的带辐板涡轮盘内腔形状优化设计方法 | |
Yoo et al. | Flow analysis and optimization of a vertical axis wind turbine blade with a dimple | |
CN103903607A (zh) | 一种调节声子晶体定向辐射频率的方法 | |
CN103823916A (zh) | 一种基于多维黎曼解的任意拉格朗日欧拉方法 | |
CN103246807A (zh) | 一种二维固-固声子晶体混合模态带隙优化方法 | |
Zeng et al. | Hydrodynamic interactions among wave energy converter array and a hierarchical genetic algorithm for layout optimization | |
Apsley et al. | Bed-load sediment transport on large slopes: model formulation and implementation within a RANS solver | |
He et al. | Optimization of a wave energy converter square array based on the differential evolution algorithm | |
Chen et al. | Numerical study of liquid mixing in microalgae-farming tanks with baffles | |
Wang et al. | Performance prediction of array arrangement on ducted composite material marine current turbines (CMMCTs) | |
CN104464715B (zh) | 一种声子晶体分束器 | |
Neshat et al. | A new insight into the position optimization of wave energy converters by a hybrid local search | |
CN103279594B (zh) | 一种二维固-固声子晶体z模态带隙优化方法 | |
CN105424800B (zh) | 基于格栅效应的室内周期矩形声扩散体散射系数预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20130724 |
|
RJ01 | Rejection of invention patent application after publication |