1第二章晶体三极管基础2.1双极型... 2.1.1BJT的工作原理1.BJT结构:...
黄佳 尹进 张昭 张洪武
摘要：基于声子晶体特性提出一种层叠式方柱型声子晶体单胞结构，通过冲击响应谱分析和随机振动分析考察其周期性结构减振特性，并与相应非周期结构对比.计算结果表明该声子晶体周期性结构对于冲击载荷和随机振动载荷均有较好的减振特性.在三组分材料参数中影响单胞带隙特性的主要因素是外层材料密度和中间层材料弹性模量.声子晶体结构的减振效果随周期数的增加而愈加明显.组分材料力学性能参数和周期型结构周期数是声子晶体结构减振设计的主要因素.
关键词：声子晶体； 单胞结构； 减振设计； 隔振材料； 低频带隙
中图分类号： TH113； O328
文献标志码：B
Abstract：A phononic crystal unit cell structure with stacked square columns is designed on the basis of the phononic crystal characteristics. The vibration reduction characteristics of the periodic structure are investigated through the shock response spectrum analysis and the random vibration analysis and compared with that of the corresponding non-periodic structure. The results show that the periodic structure of phononic crystal demonstrates better vibration reduction characteristics in the shock load and random vibration load problems. In the parameters of 3-components material， the density of the material in the outer layer and the elastic modulus of the material in the middle layer are the main factors that affect the band gap of unit cell. The vibration reduction of phononic crystal structure is obviously increased with the increase of the periodic numbers. The mechanical performance parameters of each component material and the periodic numbers of periodic structure are the main factors for the vibration reduction design of phononic crystal structure.
Key words：phononic crystal； unit cell structure； vibration reduction design； vibration isolation material； low frequency band gap
0引言
在航空航天、船舶等领域，冲击和随机振动是结构设计中需要考虑的重要力学环境，如火箭发动机点火产生的瞬态冲击[1-2]、水下武器对船舶攻击时的非接触爆炸[3]等，因此考虑冲击载荷和随机振动载荷作用的减振设计，对于结构服役安全具有重要意义.
声子晶体是一种空间周期性排布的多组分复合材料，具有阻挡某些频率范围内弹性波和声波传播的特性，对应的频率范围称为带隙或禁带.[4]该周期性复合材料的最小单元称为单胞.由声子晶体单胞构成的有限周期结构称为声子晶体结构.作为由天然晶体中电子带隙类推而来的周期结构，声子晶体可阻挡弹性波传播的特性，使得弹性波特别是声波在声子晶体中的传播特性研究受到越来越广泛的关注.[5-6]声子晶体的禁带特性使得其在减振隔振[7]、降噪等领域都有潜在、广阔的应用前景.值得注意的是，美国国防部高级计划局大力资助声子晶体在减振降噪方面的研究，使声子晶体的应用范围扩大到飞行器发动机底座、航空电子设备底盘等.[4]对于声子晶体的研究大多以带隙特性为重点，以梁板类结构[8-10]、理论分析和一维弯曲、扭转及二维面内面外振动带隙[11-15]为主.
目前，对声子晶体低频减振特性的三维振动特性研究相对较少.研究三维声子晶体的结构特性对其进一步工程化应用具有现实意义.声子晶体的尺寸和禁带的频率范围是需要重点考虑的因素，此外，声子晶体周期数对减振效果的影响以及声子晶体单胞材料参数变化对带隙特性的影响也不容忽视.鉴于此，建立一种具有隔振效果的层叠式方柱型声子晶体三维模型，考察其周期结构对冲击载荷和随机振动载荷的隔振效果，并与非周期结构进行对比.采用多重多级子结构方法[16-18]考察其三维模型的振动带隙特性和单胞材料参数对带隙特性的影响，对比分析周期结构和非周期结构传输特性，考察周期数对结构传输特性的影响.
1基于有限元的带隙计算方法
声子晶体研究中沿用晶格理论描述其空间结构的周期性.周期结构的最小重复单元称为原胞，理想周期结构是由原胞沿着3个不共面的基本平移矢量周期分布构成，这3个基本平移矢量称为基矢[19]，记为
2.2冲击响应谱振动分析
考察2种结构的冲击响应谱振动特性.分别将声子晶体结构和非周期结构的4根立柱底面设置为载荷输入面，沿z方向输入某固体火箭发动机点火时加速度冲击响应谱[2]，频率范围为10～370 Hz，输入功率谱最大值约100g（g为重力加速度），可计算得到整体结构响应.因上部结构钢板的响应为考察对象，故单独取钢板模型及其结果进行查看，见图4.钢板中部z方向加速度响应很小，约为0.23g，整个钢板最大加速度响应也很小，约为0.33g.
取非周期结构上部的钢板响应为考察对象，单独取出钢板模型及其结果进行查看，见图5.由此可知：钢板中部z方向加速度响应较大，约为32.4g；整个钢板z方向最大加速度响应位于钢板中心部位.可见，声子晶体结构上部钢板中心位置z方向加速度响应仅为一般结构的0.71%，钢板最大加速度仅为一般结构的1.02%.与非周期结构相比，声子晶体结构响应值更小，即弹性波穿过声子晶体结构后具有更小的能量.因此，本文提出的声子晶体结构具有更好的抗冲击能力，可在一定程度上减弱冲击载荷对其周围结构的不利影响.
2.3随机振动功率谱分析
考察声子晶体结构对随机振动载荷环境的减振效果.分别将声子晶体结构和非周期结构4根立柱底面设置为载荷输入面.随机振动功率谱输入选取某运载火箭贮箱增压系统随机振动工作条件[20]，载荷沿z方向.选取上部钢板结构中心点z方向加速度响应结果进行分析.声子晶体结构和非周期结构钢板中心位置z方向功率谱密度（Power Spectral Density， PSD）见图6.PSD描述随机振动条件下结构振动能量密度在频域上的分布特性.由图6可知，在给定的随机振动条件下，声子晶体结构和非周期结构上部钢板振动响应能量密度主要集中在20 Hz以上，其中非周期结构PSD最大值为1 641.33 Hz-1，最小值为1.66×10-6 Hz-1；声子晶体结构PSD最大值为7.03 Hz-1，最小值为7.64×10-20 Hz-1.在该频率范围内，声子晶体结构对应的PSD普遍显著小于非周期结构，即随机振动使声子晶体结构产生的振动在相应频率上分布的能量密度小于非周期结构.值得注意的是，在60 Hz以上的频率范围中，声子晶体结构产生随机振动的PSD甚至可以比非周期结构低20个数量级.
为使分析结果更直观，计算2种结构的PSD均方根值的积分（Cumulative Root-Mean-Square， CRMS）.PSD的CRMS表征能量概念.计算得到声子晶体结构和非周期结构钢板中心位置z方向PSD的CRMS结果见图7.作为对随机振动响应量能量大小的衡量标准，PSD的CRMS的大小可以反映结构对随机振动的响应情况.在给定的随机振动输入条件下，声子晶体结构上部钢板的响应能量远小于非周期结构的，其中在50 Hz以上的频率范围中，声子晶体结构的PSD的CRMS值一直小于非周期结构，且二者差值越来越大，这与PSD及其CRMS之间的积分关系一致.由2种结构PSD及其CRMS的结果对比可知，声子晶体结构能够显著减小随机振动对上部结构的影响，说明所提出的声子晶体结构具有更好的减振特性.
3单胞带隙特性计算与分析
3.1单胞带隙特性计算
为进一步对本文提出的声子晶体隔振特性进行研究，计算该声子晶体单胞带隙特性，并分析其影响因素.对前文中的单胞形式进行带隙特性计算，得到带隙特性，见图8.由此可知：在80～310 Hz出现宽度很大的禁带，说明层叠型方柱声子晶体单胞具有良好的带隙特性，可以使频率范围较宽且频率较低的弹性波大幅衰减.这可以解释该声子晶体的周期型结构能够对相应的冲击和随机振动问题产生明显减振效果的现象.
3.2材料参数对单胞带隙特性的影响
为探讨该型声子晶体单胞带隙的可调控性，研究单胞材料参数对单胞带隙特性的影响.当中心质量块与中间层材料不变时，考察外层金属材料参数对该声子晶体单胞带隙特性的影响.外层金属主要材料参数有弹性模量、密度和泊松比.分别计算单胞带隙上下边界和带隙中心频率随3个材料参数的变化规律，计算结果见图9a～9c.由此可以看出：外层金属的弹性模量和泊松比对带隙几乎没有影响，外层金属的密度只影响带隙的下边界；随着外层金属密度变大为1 600～31 600 kg/m3，带隙下边界先降低后趋于不变，为66～146 Hz，从而带隙中心频率也呈现先降低而后趋于不变的趋势，范围为230～190 Hz.
当中心质量块与外层金属材料不变时，考察中间层材料参数对该声子晶体单胞带隙特性的影响.
中间层主要材料参数包括弹性模量、密度和泊松比，分别计算单胞带隙上下边界和带隙中心频率随3个材料参数的变化规律，计算结果见图9d～9f.由此可知，中间层的弹性模量对带隙的上下边界均有影响，中间层弹性模量（0.5～3 MPa）越大，带隙上下边界越高，上边界为223～540 Hz，下边界为56～137 Hz，且上边界增高的程度随弹性模量增大而更显著.因此，中间层弹性模量越大，带隙宽度Δf越大，为167～403 Hz，带隙中心频率fc越高，为139.5～338.5 Hz，带隙宽度变化程度更为显著，即相对带隙rbg=Δf/fc越大；中间层泊松比变大，带隙中心频率稍有提高，但带隙宽度变化不大；中间层的材料密度对带隙几乎没有影响.可见，中间层材料对于该单胞带隙特性影响较为显著的参数为弹性模量.
4声子晶体结构传输特性计算分析
4.1传输特性计算
为验证单胞带隙特性并更直观地展示其衰减弹性波的作用，考察声子晶体周期结构和非周期结构的传输特性.分别将声子晶体结构和非周期结构4根立柱底面设置为载荷输入面，输入载荷为z方向单位加速度，频率范围为0～400 Hz.考察钢板中心点位置z方向加速度响应a/a0.计算得到2种结构传输特性曲线见图10.
参考结构（非周期结构）相对加速度对数值在0～400 Hz的大部分频率范围内都大于0，即参考结构在对应频率范围内会将载荷输入放大.声子晶体结构在本文考虑的频率范围内的绝大频域上产生的相对加速度对数值都小于0，即声子晶体结构在对应频率范围内会将载荷输入衰减.其中，在30～400 Hz频率范围内存在宽频率段衰减现象，衰减范围与上文计算得到的声子晶体单胞带隙特性一致，且有限周期声子晶体结构在比带隙范围更大的频域范围内也存在使载荷输入衰减的效果，因此声子晶体具有良好的减振作用.
4.2周期数对振动特性的影响
为考察声子晶体周期数量对周期结构隔振效果的影响，分别建立1周期和2周期声子晶体模型并计算其传输特性，有限元模型见图11.
分别对1，2和3周期时声子晶体结构传输特性进行计算，得到其传输特性曲线，见图12.3种周期数的声子晶体结构z方向加速度均在50 Hz处开始衰减，在50～310 Hz范围内有很显著的衰减效果，这
与上述计算得到的单胞带隙特性一致.3种声子晶体结构相对加速度最小值均出现在200～250 Hz范
围内，1周期结构加速度响应最小值较输入小6个数量级，2周期结构小10个数量级，3周期结构小
12个数量级.
由此可见，周期数分别为1，2和3时，声子晶体结构在带隙范围内均有减振效果，周期数越大，减振效果越显著.在对于减振要求不高的情况下，单胞形式的1周期声子晶体结构即可实现一定的减振功能.
5结论
基于声子晶体概念建立层叠式方柱型隔振材料单胞，具有能够衰减较宽的低频弹性波的特性.通过对该型单胞的带隙特性及其周期结构振动特性分析，得到以下结论.
1）层叠式方柱型声子晶体单胞的有限周期结构具有很好的隔振效果，可以大幅衰减振动频率在其单胞带隙范围内的弹性波，较一般非周期结构有明显的隔振效果.
2）层叠式方柱型声子晶体单胞具有良好的带隙特性，可以使相应频率范围内的弹性波大幅衰减.
3）声子晶体结构中影响带隙特性的主要设计参数为外层材料密度和中间层材料弹性模量.
4）周期声子晶体结构具有良好的减振特性，并且随着周期数的增加，减振效果愈加明显.
参考文献：
[1]易磊， 寇军强， 李卫鹏， 等. 战术固体火箭发动机点火冲击试验研究[J]. 弹箭与制导学报， 2010， 30（6）： 150-151.
YI Lei， KOU Junqiang， LI Weipeng， et al. Experiment research on ignition impulsion of solid rocket motor[J]. J Projectiles， Rockets， Missiles & Guidance， 2010， 30（6）： 150-151.
[2]王翠荣， 施广富， 郭军. 固体火箭发动机冲击信号响应谱分析[J]. 固体火箭技术， 2003， 26（2）： 57-60.
WANG Cuirong， SHI Guangfu， GUO Jun. Analysis of impact signal response spectrum of solid rocket motors[J]. J Solid Rocket Technol， 2003， 26（2）： 57-60.
[3]石菲. 船舶发动机排气消声器抗冲击性能研究[D]. 哈尔滨： 哈尔滨工程大学， 2006.
[4]郁殿龙. 基于声子晶体理论的梁板类周期结构振动带隙特性研究[D]. 长沙： 国防科学技术大学， 2006.
[5]LIU Z， ZHANG X， MAO Y， et al. Locally resonant sonic materials[J]. Scince， 2000（289）： 1734-6.
[6]MALDOVAN M. Sound and heat revolutions in phononics[J]. Scince， 2013（503）： 209-217.
[7]温激鸿， 王刚， 郁殿龙， 等. 声子晶体振动带隙及减振特性研究[J]. 中国科学： E辑： 技术科学， 2007， 37（9）： 1126-1139.
WEN Jihong， WANG Gang， YU Dianlong， et al. Vibration attenuation and band gap characteristics of phononic crystals[J]. Scientia Sinica： Ser E： Technol Sci， 2007， 37（9）： 1126-1139.
[8]宋卓斐， 王自东， 王艳林， 等. 一维杆状声子晶体的带隙特性[J]. 振动与冲击， 2010， 29（2）： 145-148.
SONG Zhuofei， WANG Zidong， WANG Yanlin， et al. Band gap property of a one dimension rod phononic crystal[J]. J Vibration & Shock， 2010， 29（2）： 145-148.
[9]卢天健， 高国钦， 马守林. 二维四方排列半圆铝管/空气声子晶体的禁带特性[J]. 中国科学： E辑： 技术科学， 2009， 39（1）： 57-64.
LU Tianjian， GAO Guoqin， MA Shoulin. Properties of band gaps in 2D half Aluminum rods of hollows-air phononic crystal with square lattice[J]. Scientia Sinica： Ser E： Technol Sci， 2009， 39（1）： 57-64.
[10]丁红星， 沈中华， 李加， 等. 复合兰姆波声子晶体中超宽部分禁带[J]. 物理学报， 2012， 61（19）： 394-401 .
DING Hongxing， SHEN Zhonghua， LI Jia， et al. Large partial band-gaps for Lamb waves in multiple phononic crystals thin plates[J]. Acta Phisica Sinica， 2012， 61（19）： 394-401.
[11]SHENG P， ZHANG X， LIU Z， et al. Locally resonant sonic materials[J]. Physica B： Condensed Mat， 2003， 338（1-4）： 201-205.
[12]文岐华， 左曙光， 魏欢. 多振子梁弯曲振动中的局域共振带隙[J]. 物理学报， 2012， 61（3）： 240-246.
WEN Qihua， ZUO Shuguang， WEI Huan. Locally resonant elastic wave band gaps in flexural vibration of multi-oscillators beam[J]. Acta Phisica Sinica， 2012， 61（3）： 240-246.
[13]舒海生， 刘少刚， 王威远， 等. 集中质量边界条件下声子晶体杆的纵向振动传递特性研究[J]. 振动与冲击， 2012， 31（19）： 113-117.
SHU Haisheng， LIU Shaogang， WANG Weiyuan， et al. Transmission characteristics of longitudinal vibration of a phononic crystal rod with concentrated mass boundary condition[J]. J Vibration & Shock， 2012， 31（19）： 113-117.
[14]YU D， WEN J， ZHAO H， et al. Vibration reduction by using the idea of phononic crystals in a pipe-conveying fluid[J]. J Sound & Vibration， 2008， 318（1-2）： 193-205.
[15]ZHANG S， WU J， HU Z. Low-frequency locally resonant band-gaps in phononic crystal plates with periodic spiral resonators[J]. J Appl Phys， 2013， 113（16）： 163511（1-8）.
[16]钟万勰. 一个多用途的结构分析程序JIGFEX（一）[J]. 大连工学院学报， 1977， 17（3）： 19-42.
ZHONG Wanxie. A general-purpose program for structure analysis JIGFEX（Ⅰ）[J]. J Dalian Univ Technol， 1977， 17（3）： 19-42.
[17]钟万勰. 一个多用途的结构分析程序JIGFEX（二）[J]. 大连工学院学报， 1977， 17（4）： 14-35.
ZHONG Wanxie. A general-purpose program for structure analysis JIGFEX（Ⅱ）[J]. J Dalian Univ Technol， 1977， 17（4）： 14-35.
[18]LU X， LIN J， ZHONG W. Subspace iteration method in multi-level substructure systems[J]. Comput & Structures， 1989， 33（2）： 459-462.
[19]温熙森， 温激鸿， 郁殿龙， 等. 声子晶体[M]. 北京： 国防工业出版社， 2009： 29-41.
[20]王翀， 梁国柱. 随机振动条件火箭贮箱增压系统工作过程仿真[J]. 北京航空航天大学学报， 2013， 39 （10）： 1314-1318.
WANG Chong， LIANG Guozhu. Simulation of the working process of launch vehicle pressurization system in random vibration condition[J]. J Beijing Univ Aeronautics & Astronautics， 2013， 39 （10）： 1314-1318.
（编辑武晓英）
