CN112765884A - 一种基于扩展等几何分析法的声子晶体材料结构设计方法 - Google Patents
一种基于扩展等几何分析法的声子晶体材料结构设计方法 Download PDFInfo
- Publication number
- CN112765884A CN112765884A CN202110061986.3A CN202110061986A CN112765884A CN 112765884 A CN112765884 A CN 112765884A CN 202110061986 A CN202110061986 A CN 202110061986A CN 112765884 A CN112765884 A CN 112765884A
- Authority
- CN
- China
- Prior art keywords
- scatterer
- phononic crystal
- extended
- crystal material
- nurbs
- 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
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/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- 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
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- 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)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于扩展等几何分析法的声子晶体材料结构设计方法。通过改变散射体形状,设计隔声频带宽和隔声频率低的声子晶体材料。本发明的扩展等几何分析法的计算网格独立于结构内部的几何,因此散射体形状变化时无需重构计算网格;采用基于LR‑NURBS的扩展等几何分析法分析声子晶体材料的能带结构特性时,在散射体界面附近采用小尺度网格,其他区域采用大尺度网格,既可提高精度,又能节约时间。散射体形状采用NURBS样条函数描述,这样可以构造任意复杂形状的散射体。本方法得到的散射体形状光滑,可以直接用于设计;可以精确地描述复杂形状的声子晶体周期胞元,实现CAD与CAE无缝结合地设计声子晶体材料的结构。
Description
技术领域
本发明属于声学超材料的结构设计领域,特别是涉及一种基于扩展等几何分析法的声子晶体材料结构设计方法。
背景技术
声子晶体是一种具有周期性结构并具有声/弹性波带隙的人工新型声学功能材料。弹性波在传播过程中会受到弹性常数的周期性调制,可能会产生声子带隙,即一定频率范围的弹性波传播被抑制。正是由于声子晶体的这种声学特性:奇特的带隙、局域、负折射特征,使得其在振动隔离、噪音控制、声学波导器件及声学透镜等方面具有广阔的应用前景。
声子晶体材料的结构设计包括散射体形状和其在基体中的分布。自二十世纪九十年代以来,随着声子晶体的发现和发展,人们开始关注应用声子晶体技术来设计具有声学禁带特性的隔声结构。在声子晶体中,可以通过增加基体和散射体的密度比、声速度比、物理参数,设置缺陷、增加组分、加设吸收层方法来获得一定范围的声学禁带或带隙特性,但上述方法都会不同程度的增加禁带宽度。在工程领域实际应用中,主要根据不同应用场合和对象,有选择地采用以上方法来实现声学禁带特性。但由于所获得的禁带区间往往不在人耳敏感的频率范围,迄今为止,大多数文献报道的声子晶体的禁带频率均处于超声范围或高声范围,都不能达到理想的隔声频带宽和隔声频率低的状态。
综上所述,寻找一种能够达到理想的隔声频带宽和隔声频率低的状态的声子晶体材料结构设计方法,成为研究人员关注的问题。
发明内容
为了解决上述背景技术提出的技术问题,本发明旨在提供一种基于扩展等几何分析法的声子晶体材料结构设计方法,通过改变散射体形状,设计隔声频带宽和隔声频率低的声子晶体材料。
为了实现上述技术目的,本发明的技术方案为:
一种基于扩展等几何分析法的声子晶体材料结构设计方法,具体包括以下步骤:
S1、基于周期性边界条件和Bloch定理,将声子晶体周期结构的离散特征方程化归到声子晶体周期胞元内的复系数特征方程;并将所述复系数特征方程分为实部和虚部两组方程,得到实数特征方程和虚数特征方程;
S2、采用NURBS样条集函数来构建声子晶体的散射体几何模型,即散射体边界NURBS曲线模型;
S3、对所述散射体几何模型进行优化,并构建声子晶体的隔声频带最宽和隔声频率最低的双目标优化模型;
S4、基于LR-NURBS的扩展等几何分析法,建立声子晶体周期胞元分析模型;
S5、采用粒子群算法求解所述双目标优化模型问题,获得声子晶体周期胞元系统的能带结构;
S6、对所述粒子群算法进行改进,并基于改进后的粒子群算法来求解所述双目标优化模型问题,得到多个描述散射体形状的控制点;并根据所述控制点,得到NURBS曲线,即得到所述声子晶体的散射体几何形状。
优选地,所述散射体边界NURBS曲线模型是由散射体边界曲线上一系列控制点构成,所述控制点为P1,P2,…,Pn,其中n为控制点数。
优选地,所述散射体边界曲线的表达式为:
优选地,所述步骤S3中的优化是对所述一系列控制点进行优化设计变量。
优选地,所述步骤S4具体为:
S41、将声子晶体周期胞元的整个区域均匀分成若干个计算网络;
S42、基于符号距离函数和NURBS曲线函数,得到声子晶体任意形状的水平集函数;
S43、依据控制点水平集函数值,在所述声子晶体周期胞元的整个区域中找出散射体区域,在所述散射体区域内进行局部网格细化;
S44、基于所述水平集函数,确定加强的控制点;
S45、依据弱不连续问题,计算加强函数,得到扩展等几何分析的位移逼近。
优选地,所述LR-NURBS的扩展等几何分析法中的LR-NURBS样条基函数为:
优选地,所述步骤S5具体为:优化分析的每次迭代正分析时,在边界上应用Bloch定理,采用扩展等几何分析法求解所述实数特征方程,获得系统的能带结构。
优选地,所述步骤S6中改进的方法为:在所述粒子群算法中引入多样性反馈的权重、混沌映射和混合策略的变异行为。
与现有技术相比,本发明的有益效果在于:
(1)本发明使用扩展等几何分析法求解系统的能带结构,每次优化迭代无需重划网格,相较常规有限元法和等几何分析法节约了大量的计算时间。另外,可以精确地描述复杂形状的声子晶体周期胞元,实现CAD与CAE无缝结合地设计声子晶体材料的结构;
(2)本发明采用NURBS样条函数描述散射体几何,可以描述任意复杂形状的散射体;
(3)本发明采用基于LR-NURBS的扩展等几何分析法计算系统的能带结构,可以在散射体附近采用小尺度网格,其他区域采用大尺度网格,提高精度和节省计算量。相对LR-B样条而言,LR-NURBS可以描述复杂形状的声子晶体周期胞元。相对NURBS样条而言,LR-NURBS可以进行局部细化;
(4)本发明使用了多样性反馈的权重及混沌映射等机制,并嵌入混合策略的变异行为,提高了优化算法的效率;
(5)本发明可以得到不同精英解的描述散射体形状的控制点,设计者可以根据需要设计出相应的散射体形状;
(6)本发明得到的散射体形状光滑,可以直接用于设计。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明方法流程图;
图2为本发明实施例声子晶体周期胞元模型示意图;
图3本发明实施例声子晶体周期胞元模型的初始散射体示意图;
图4为本发明实施例声子晶体周期胞元模型的扩展等几何分析网格示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1
参照图1所示,本发明提供一种基于扩展等几何分析法的声子晶体材料结构设计方法,具体包括以下步骤:(为说明本发明的方法流程,以如图2所示的声子晶体周期胞元为例,该胞元含有1个散射体。)
S1、利用周期性边界条件和Bloch定理,将声子晶体周期结构的离散特征方程化归到声子晶体一个周期胞元内的复系数特征方程,将复系数特征方程分为实部和虚部两组方程;
S2、采用NURBS样条基函数来描述声子晶体的散射体几何模型(散射体NURBS边界曲线),该曲线模型由一系列控制点构成,控制点为P1,P2,…,Pn(n为控制点数),NURBS样条基函数为R1,R2,…,Rn,参照图3所示。(采用NURBS样条函数描述散射体几何,可以构造任意复杂形状的散射体,并不仅限本实施例这一种散射体形状。)
其中,散射体边界曲线C的表达式为:
步骤3、以控制点P1,P2,…,P6为优化设计变量,得到的散射体几何是光滑的,可以直接用于设计,建立隔声频带最宽和隔声频率最低的双目标优化模型。
步骤4、建立声子晶体周期胞元的基于LR-NURBS的扩展等几何分析模型。不需要考虑声子晶体的散射体几何,先在整个区域生成均匀的计算网格,然后在散射体可能出现的区域进行网格局部细化,既能提高精度又能节省时间,参照图4所示。
扩展等几何分析中的LR-NURBS样条基函数为:
基于水平集函数值,确定加强的控制点;
依据弱不连续问题,计算加强函数,得到扩展等几何分析的位移逼近。
步骤5、采用粒子群算法求解双目标优化问题。优化分析的每次迭代正分析时,在边界上应用Bloch定理,采用扩展等几何分析法求解实数特征方程,获得系统的能带结构。
步骤6、为了提高优化算法的效率,在粒子群算法中引入多样性反馈的权重、混沌映射和混合策略的变异行为。采用改进的粒子群算法求解多目标优化问题,得到不同精英解的描述散射体形状的控制点,从而求得相应的NURBS曲线,即散射体形状(设计者可以根据具体需要确定散射体形状)。
综上,扩展等几何分析法的计算网格独立于结构内部的几何,因此散射体形状变化时无需重构计算网格;基于LR-NURBS的扩展等几何分析法还具有几何精确、高精度、高阶连续和局部细化等特性。采用基于LR-NURBS的扩展等几何分析法分析声子晶体材料的能带结构特性时,在散射体界面附近采用小尺度网格,其他区域采用大尺度网格,既可提高精度,又能节约时间。散射体形状采用NURBS样条函数描述,这样可以构造任意复杂形状的散射体。以描述散射体形状的控制点为优化变量,建立隔声频带最宽和隔声频率最低的双目标优化模型。采用改进的粒子群算法求解多目标优化问题,得到不同精英解的描述散射体形状的控制点,从而求得相应的NURBS曲线,即散射体形状。本方法得到的散射体形状光滑,可以直接用于设计;可以精确地描述复杂形状的声子晶体周期胞元,实现CAD与CAE无缝结合地设计声子晶体材料的结构。
以上所述的实施例仅是对本发明的优选方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
Claims (8)
1.一种基于扩展等几何分析法的声子晶体材料结构设计方法,其特征在于,具体包括如下步骤:
S1、基于周期性边界条件和Bloch定理,将声子晶体周期结构的离散特征方程化归到声子晶体周期胞元内的复系数特征方程;并将所述复系数特征方程分为实部和虚部两组方程,得到实数特征方程和虚数特征方程;
S2、采用NURBS样条集函数来构建声子晶体的散射体几何模型,即散射体边界NURBS曲线模型;
S3、对所述散射体几何模型进行优化,并构建声子晶体的隔声频带最宽和隔声频率最低的双目标优化模型;
S4、基于LR-NURBS的扩展等几何分析法,建立声子晶体周期胞元的力学分析模型;
S5、采用基于LR-NURBS的扩展等几何分析法求解实数特征方程,获得声子晶体周期胞元系统的能带结构;
S6、对所述粒子群算法进行改进,并基于改进后的粒子群算法来求解所述双目标优化模型问题,得到多个描述散射体形状的控制点;并根据所述控制点,得到NURBS曲线,即得到所述声子晶体的散射体几何形状。
2.根据权利要求1所述的基于扩展等几何分析法的声子晶体材料结构设计方法,其特征在于:所述散射体边界NURBS曲线模型是由散射体边界曲线上一系列控制点构成,所述控制点为P1,P2,…,Pn,其中n为控制点数。
4.根据权利要求2所述的基于扩展等几何分析法的声子晶体材料结构设计方法,其特征在于:所述步骤S3中的优化是对所述一系列控制点进行优化设计变量。
5.根据权利要求1所述的基于扩展等几何分析法的声子晶体材料结构设计方法,其特征在于:所述步骤S4具体为:
S41、将声子晶体周期胞元的整个区域均匀分成若干个计算网络;
S42、基于符号距离函数和NURBS曲线函数,得到声子晶体任意形状的水平集函数;
S43、依据控制点水平集函数值,在所述声子晶体周期胞元的整个区域中找出散射体区域,在所述散射体区域内进行局部网格细化;
S44、基于所述水平集函数,确定加强的控制点;
S45、依据弱不连续问题,计算加强函数,得到扩展等几何分析的位移逼近。
7.根据权利要求1所述的基于扩展等几何分析法的声子晶体材料结构设计方法,其特征在于:所述步骤S5具体为:优化分析的每次迭代正分析时,在边界上应用Bloch定理,采用扩展等几何分析法求解所述实数特征方程,获得所述声子晶体周期胞元系统的能带结构。
8.根据权利要求1所述的基于扩展等几何分析法的声子晶体材料结构设计方法,其特征在于:所述步骤S6中改进的方法为:在所述粒子群算法中引入多样性反馈的权重、混沌映射和混合策略的变异行为。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110061986.3A CN112765884B (zh) | 2021-01-18 | 2021-01-18 | 一种基于扩展等几何分析法的声子晶体材料结构设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110061986.3A CN112765884B (zh) | 2021-01-18 | 2021-01-18 | 一种基于扩展等几何分析法的声子晶体材料结构设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112765884A true CN112765884A (zh) | 2021-05-07 |
CN112765884B CN112765884B (zh) | 2021-09-21 |
Family
ID=75702772
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110061986.3A Active CN112765884B (zh) | 2021-01-18 | 2021-01-18 | 一种基于扩展等几何分析法的声子晶体材料结构设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112765884B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114491769A (zh) * | 2022-02-17 | 2022-05-13 | 河海大学 | 一种基于等几何分析法自由曲面结构一体化形态创构方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103176272A (zh) * | 2011-12-21 | 2013-06-26 | 北京邮电大学 | 二维光子晶体最大绝对带隙结构优化方法 |
US20130185035A1 (en) * | 2012-01-13 | 2013-07-18 | California Institute Of Technology | Systems and Methods of Analysis of Granular Elements |
CN110598275A (zh) * | 2019-08-23 | 2019-12-20 | 南京理工大学 | 基于响应面建模和改进粒子群算法的车轮型面优化方法 |
CN111402851A (zh) * | 2020-03-13 | 2020-07-10 | 中国农业大学 | 一种仿生声子晶体及其制作方法 |
CN112001004A (zh) * | 2020-08-21 | 2020-11-27 | 西安交通大学 | 一种解析中高频振动结构能量密度场的nurbs等几何分析方法 |
-
2021
- 2021-01-18 CN CN202110061986.3A patent/CN112765884B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103176272A (zh) * | 2011-12-21 | 2013-06-26 | 北京邮电大学 | 二维光子晶体最大绝对带隙结构优化方法 |
US20130185035A1 (en) * | 2012-01-13 | 2013-07-18 | California Institute Of Technology | Systems and Methods of Analysis of Granular Elements |
CN110598275A (zh) * | 2019-08-23 | 2019-12-20 | 南京理工大学 | 基于响应面建模和改进粒子群算法的车轮型面优化方法 |
CN111402851A (zh) * | 2020-03-13 | 2020-07-10 | 中国农业大学 | 一种仿生声子晶体及其制作方法 |
CN112001004A (zh) * | 2020-08-21 | 2020-11-27 | 西安交通大学 | 一种解析中高频振动结构能量密度场的nurbs等几何分析方法 |
Non-Patent Citations (1)
Title |
---|
CHRISTOPHER ZIMMERMANN,ROGER A. SAUER: "Adaptive local surface refinement based on LR NURBS", 《COMPUTATIONAL MECHANICS》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114491769A (zh) * | 2022-02-17 | 2022-05-13 | 河海大学 | 一种基于等几何分析法自由曲面结构一体化形态创构方法 |
CN114491769B (zh) * | 2022-02-17 | 2022-12-09 | 河海大学 | 一种基于等几何分析法自由曲面结构一体化形态创构方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112765884B (zh) | 2021-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Acoustic topology optimization of sound absorbing materials directly from subdivision surfaces with isogeometric boundary element methods | |
Liu et al. | Shape optimization of sound barrier using an isogeometric fast multipole boundary element method in two dimensions | |
CN112989680B (zh) | 减少网格使用量的fvfd远场积分边界条件计算方法 | |
Chen et al. | Subdivision surfaces—boundary element accelerated by fast multipole for the structural acoustic problem | |
CN113158527B (zh) | 一种基于隐式fvfd计算频域电磁场的方法 | |
CN112765884B (zh) | 一种基于扩展等几何分析法的声子晶体材料结构设计方法 | |
Zheng et al. | An efficient method to reduce the numerical dispersion in the ADI-FDTD | |
CN112446163B (zh) | 基于参数化水平集的能量有限元拓扑优化方法 | |
Xie et al. | Isogeometric topology optimization based on energy penalization for symmetric structure | |
CN112836416B (zh) | 一种用于抑制弹性波传播的声子晶体结构优化设计方法 | |
CN111079326A (zh) | 二维各向异性网格单元度量张量场光滑化方法 | |
CN113378440A (zh) | 一种流固耦合数值模拟计算方法、装置及设备 | |
Wang et al. | Accurate and efficient hydrodynamic analysis of structures with sharp edges by the Extended Finite Element Method (XFEM): 2D studies | |
Yang et al. | Isogeometric double-objective shape optimization of free-form surface structures with Kirchhoff–Love shell theory | |
Li et al. | Material design with topology optimization based on the neural network | |
CN116401498A (zh) | 一种基于mb-crwg基函数的电磁散射分析方法、设备及介质 | |
CN115035965A (zh) | 一种基于等几何形状优化的单相声子晶体板带隙优化方法 | |
Ruan et al. | Band gap characteristics of bionic acoustic metamaterials based on spider web | |
CN111046542B (zh) | 一种评估二十面体的le网格计算特性方法及离散方法 | |
Cheng et al. | Re-analysis method for inversion of block matrix based on change threshold | |
CN112364559A (zh) | 一种波浪能发电装置布局优化方法及装置 | |
Liu et al. | Heterogeneous CPU-GPU Accelerated Parallel Subgridding FDTD Algorithm | |
CN114996992B (zh) | 一种旋转对称光学器件的本征模式求解方法及系统 | |
Hwang et al. | Locally implicit hybrid algorithm for steady and unsteady viscous flows | |
CN117932754B (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 |