基于人工神经网络的声子晶体逆向设计

声子晶体是一种人工周期性复合材料, 其带隙特性使其在减振、隔声、滤波和声学功能器件等领域具有潜在的应用价值. 如何准确操纵声波和机械波是声子晶体设计的主要挑战. 现有设计方法是基于对结构几何参数与材料参数的分析调整使其匹配特定的应用特性, 设计效率不高且无法达到最佳性能. 为此, 本文以一维层状声子晶体为例, 提出了一种基于Softmax逻辑回归和多任务学习的人工神经网络声子晶体逆向设计方法, 其中, Softmax逻辑回归实现分层结构各区域材料种类的选择, 通过多任务学习确定各区域材料的分布, 从而, 将声子晶体逆向设计问题转化为对单位胞元拓扑结构多组分材料的分类问题. 首先, 随机生成大量声子晶体拓扑结构样本; 然后, 采用有限元法进行并行计算得到所有样本的带隙分布; 接着, 通过神经网络建立带隙分布和拓扑结构之间的映射关系; 最后, 利用训练好的神经网络设计具有目标带隙特性的声子晶体, 即以目标带隙作为神经网络的输入, 网络将直接输出对应的声子晶体单元胞元拓扑结构. 算例表明本方法可根据应用需求快速高效地得到具有目标带隙的一维声子晶体. 该方法为声子晶体的逆向设计提供了一种新颖思路.      

基于材料相变的声子晶体带隙可调结构设计与性能分析

提出一种基于材料相变的穿孔型带隙可调声子晶体结构.其结构形式为含缝隙的形状记忆合金和环氧树脂的组合体,通过温度变化诱发相变引起的形状记忆合金材料性质的变化,实现声子晶体的带隙变化;通过合理布置缝隙与形状记忆合金相变材料的位置,实现声子晶体带隙性质的可调设计.基于有限元方法,建立了可调声子晶体的分析模型,分析了形状记忆合金的填充分数以及相变等对带隙性能的影响规律.分析结果表明,通过合理设计微结构形式,材料相变可实现带隙位置和宽度的调节,同时可实现特定频段内带隙的有无.     

包裹层结构对轻质声子晶体低频振动带隙的影响

基于局域共振机理提出新型轻质声子晶体包裹层结构设计方法。借助有限元法,计算新型声子晶体能带结构、本征模态,分析包裹层总缺口度数一定时不同包裹层缺口数量及布置位置对第一完全带隙截止频率的影响;设计包裹层与散射体连接形式为线连接与点连接两种新型声子晶体模型,分别得到起始频率为37.4Hz及19.0Hz的第一完全带隙,分析第一完全带隙起始频率处散射体本征模态,揭示新型声子晶体极低第一完全带隙起始频率产生机理;进而,与传统面连接型声子晶体通过增加散射体质量降低第一完全带隙起始频率的方法对比。研究结果表明,包裹层总缺口度数一定时,采用缺口数量更多且缺口位置距离连接短板更大的包裹层布置形式能得到更宽的第一完全带隙;提出的包裹层与散射体线连接与点连接型声子晶体在获得极低第一完全带隙起始频率的同时,显著降低了声子晶体质量,突破了传统声子晶体通过增加散射体质量降低第一完全带隙起始频率的限制,为轻质声子晶体获得极低局域共振带隙起始频率的研究设计提供了参考。     

最优二维固/固声子晶体带隙特性研究

声子晶体是一种具有弹性波带隙特性的周期性复合材料,其带隙特性受单胞拓扑形状的影响.通过拓扑优化技术,能够突破传统设计方法的局限,实现对声子晶体的主动设计.论文基于遗传算法和改进的平面波展开法,通过两阶段的优化过程,得到具有最大相对带隙的二维铜/环氧树脂声子晶体结构,并进一步研究不同带隙最优声子晶体单胞拓扑形状及其带隙特性.结果表明,利用所开发的带隙优化程序,能够得到满足带隙要求的具有全局最优的声子晶体结构;声子晶体最低带隙所对应的单胞是简单晶格结构,其散射体形状简单,而且对应的带隙频率最低,相对带隙最大,对于隔音减振最有实用价值.     

SH波在局部非均匀声子晶体中的传播

使用传输矩阵法,研究了水平剪切波(SH波)在一维局部材料非均匀声子晶体中的传播。首先使用皮亚诺方程给出了SH波在特定非均匀梯度介质声场的传输矩阵。然后使用传输矩阵和布洛赫定律给出了该声子晶体的禁带结构和传输系数。通过具体计算,发现这种非均匀性能明显影响声子晶体的频散曲线和禁带结构。然后比较了幂指数取不同值时,禁带结构与归一化频率和填充比的关系。发现幂指数不同,声子晶体的禁带不同。     

声光子晶体带隙特性与声光耦合作用研究综述

声光子晶体是一种同时具有光子和声子带隙的周期性结构.声光子晶体为同时控制电磁波和弹性波的传播提供了一个系统的平台,并在光学、声学及声光多功能器件、腔光力学等领域展现了十分广阔的应用前景.论文首先介绍了声光子晶体的基本概念,包括声光子晶体的材料、空间周期分类、能带结构的计算方法等;阐述了不同体系下声光子晶体双重带隙的特性,以及拓扑优化方法在声光子晶体带隙优化方面的应用;然后简要介绍了腔光力学,以及计算声光耦合作用的准静态方法和光力耦合系数方法,并针对当前各种声光子晶体结构中声光耦合作用的研究进行了阐述;还进一步介绍了声光子晶体波导和传感器的相关研究;最后,基于当前声光子晶体的研究进展对未来的研究方向进行了展望,其中涉及到增强声光子晶体谐振腔的声光相互作用、三维声光子晶体的研究、声光超构材料的设计、声光子晶体器件设计与应用等.     

随机声子晶体不确定性分析的直接概率积分法

由于制造工艺存在大量不确定因素，声子晶体材料属性不可避免地具有随机不确定性，使得声子晶体的物理响应呈现随机性，进而对声子晶体的减振降噪性能造成不利影响。如果采用传统的蒙特卡洛方法对随机声子晶体的物理响应进行不确定性量化，则计算代价昂贵。为此，本文基于高效的直接概率积分法对含随机材料参数的声子晶体开展不确定性量化研究。首先，在直接概率积分法框架下，对随机声子晶体能带结构的上下边界频率、带隙宽度和频率响应函数进行了不确定性量化，考察了随机参数大变异性对声子晶体带隙宽度的影响。然后，建立了声子晶体减振降噪可靠度计算公式，对考虑随机不确定性影响下的声子晶体减振降噪性能进行了定量评估。通过与蒙特卡洛方法比较，两个算例验证了直接概率积分法在随机声子晶体不确定性传播和减振降噪可靠性评估中的准确性和高效性。最后，基于直接概率积分法对局域共振型随机声子晶体进行了可靠度分析。结果表明，橡胶材料的随机性对局域共振型声子晶体减振降噪性能有较大影响。     

镜像对称准周期固/液声子晶体的传输特性

以准周期Fibonacci序列为基础,构建了准周期镜像对称结构的一维固液介质声子晶体。利用转移矩阵法,研究了弹性波在所构建的声子晶体中的透射特性。结果表明:弹性波正入射时,转移矩阵与变介质的特性阻抗有关,由于声子晶体结构的局域失谐,在透射谱的中心频率附近对称出现谐振模;对于Q7序列的声子晶体,带隙内谐振模带宽小于10Hz,能实现超窄带声滤波;弹性波斜入射时,带隙向高频区移动;入射角为5°时,谐振模也向高频移动;入射角为10°时,由于谐振的增强,声子晶体结构对频率的选择性提高,使带隙内的谐振模消失。     

考虑可制造性约束的声子晶体多目标拓扑优化

对声子晶体进行拓扑优化可得到具有目标带隙特性的声子晶体结构, 在减振、隔声等领域具有潜在应用价值. 然而, 优化结果常会出现孤立材料单元而导致的制造困难问题. 针对该问题, 本文提出协同考虑带隙性能和可制造性约束的二维多相声子晶体多目标拓扑优化方法. 以在特定频率段带隙最宽和结构质量最小为优化目标, 在对微结构进行连通性分析的基础上, 引入考虑可制造性因素的附加约束, 并利用有限元法和具有精英选择策略的非支配排序遗传算法对该优化模型进行数值求解. 通过数值算例验证了本文模型及方法的合理性和有效性. 结果表明, 附加可制造性约束的拓扑优化模型可有效避免二维声子晶体构型中出现孤立材料单元的情况, 优化结果在满足带隙性能预期指标的同时也可兼顾到实际可制造性要求. 与仅仅考虑带隙性能的单目标优化结果相比, 本文提出的同时兼顾带隙性能和可制造因素的多目标优化模型可以针对实际应用和制造条件, 实现不同目标间的平衡, 具有显著优势和良好的应用前景.      

一维声子晶体带隙的非局部辛分析

周期性弹性复合结构(声子晶体)中传播的弹性波存在特殊的色散关系:弹性波只能在某段频率范围内无损耗的传播,该频率范围称为通带.一维声子晶体的色散问题可以看作分层介质中弹性波的传播问题,利用二维弹性理论予以分析.为了研究非局部效应对声子晶体带隙特性的影响,将Eringen的二维非局部弹性理论引入到Hamilton体系下,利用精细积分与扩展的Wittrick Williams算法可获取任意频率范围内的本征解.通过对不同算例的数值计算,分析和对比了非局部理论方法与传统局部理论方法的差别.并进一步指出了该套算法的适用性和优势所在.     

An analytical study on the nonlinear vibration of functionally graded beams

Nonlinear vibration of beams made of functionally graded materials (FGMs) is studied in this paper based on Euler-Bernoulli beam theory and von Kármán geometric nonlinearity. It is assumed that material properties follow either exponential or power law distributions through thickness direction. Galerkin procedure is used to obtain a second order nonlinear ordinary equation with quadratic and cubic nonlinear terms. The direct numerical integration method and Runge-Kutta method are employed to find the nonlinear vibration response of FGM beams with different end supports. The effects of material property distribution and end supports on the nonlinear dynamic behavior of FGM beams are discussed. It is found that unlike homogeneous beams, FGM beams show different vibration behavior at positive and negative amplitudes due to the presence of quadratic nonlinear term arising from bending-stretching coupling effect.     

Stochastic analysis of a thermoelastic problem in functionally graded plates with uncertain material properties

The statistics (i.e., mean and variance) of temperature and thermal stress are analytically obtained in functionally graded material (FGM) plates with uncertainties in the thermal conductivity and coefficient of linear thermal expansion. These FGM plates are assumed to have arbitrary nonhomogeneous thermal and mechanical properties through the entire thickness of plate and are subjected to deterministic convective heating. The stochastic temperature and thermal stress fields are analysed by assuming the FGM plate is multilayered with distinct, random thermal conductivity and coefficient of linear thermal expansion in each layer. Vodicka’s method, which is a type of integral transform method, and a perturbation method are employed to obtain the analytical solutions for the statistics. The autocorrelation coefficients of each random property and cross-correlation coefficients between different random properties are expressed in exponential function forms as a non-homogeneous Markov random field of discrete space. Numerical calculations are performed for FGM plates composed of partially stabilised zirconia (PSZ) and austenitic stainless steel (SUS304), which have the largest dispersion of the random properties at the place where the volume fractions of the two constituent materials are both 0.5. The effects of the spatial change in material composition, thermal boundary condition and correlation coefficients on the standard deviations of the temperature and thermal stress are discussed.     

Mechanical behavior of functionally graded material plates under transverse load—Part I: Analysis

An elastic, rectangular, and simply supported, functionally graded material (FGM) plate of medium thickness subjected to transverse loading has been investigated. The Poisson’s ratios of the FGM plates are assumed to be constant, but their Young’s moduli vary continuously throughout the thickness direction according to the volume fraction of constituents defined by power-law, sigmoid, or exponential function. Based on the classical plate theory and Fourier series expansion, the series solutions of power-law FGM (simply called P-FGM), sigmoid FGM (S-FGM), and exponential FGM (E-FGM) plates are obtained. The analytical solutions of P-, S- and E-FGM plates are proved by the numerical results of finite element method. The closed-form solutions illustrated by Fourier series expression are given in Part I of this paper. The closed-form and finite element solutions are compared and discussed in Part II of this paper. Results reveal that the formulations of the solutions of FGM plates and homogeneous plates are similar, except the bending stiffness of plates. The bending stiffness of a homogeneous plate is Eh³/12(1 − ν²), while the expressions of the bending stiffness of FGM plates are more complicated combination of material properties.     

Magnetothermoelastic stress and perturbation of magnetic field vector in a functionally graded hollow sphere

In this article, a closed-form solution for one-dimensional magnetothermoelastic problem in a functionally graded material (FGM) hollow sphere placed in uniform magnetic and temperature fields subjected to an internal pressure is obtained using the infinitesimal theory of magnetothermoelasticity. Hyper-geometric functions are employed to solve the governing equation. The material properties through the graded direction are assumed to be nonlinear with an exponential distribution. The nonhomogeneity of the material in the radial direction is assumed to be exponential. The temperature, displacement and stress fields and the perturbation of magnetic field vector are determined and compared with those of the homogeneous case. Hence, the effect of inhomogeneity on the stresses and the perturbation of magnetic field vector distribution are demonstrated. The results of this study are applicable for designing optimum FGM hollow spheres.     

Effects of material gradients on transient dynamic mode-III stress intensity factors in a FGM

This paper presents a transient dynamic crack analysis for a functionally graded material (FGM) by using a hypersingular time-domain boundary integral equation method. The spatial variations of the material parameters of the FGM are described by an exponential law. A numerical solution procedure is developed for solving the hypersingular time-domain traction BIE. To avoid the use of time-dependent Green’s functions which are not available for general FGM, a convolution quadrature formula is adopted for approximating the temporal convolution, while a Galerkin method is applied for the spatial discretization of the hypersingular time-domain traction BIE. Numerical results for the transient dynamic stress intensity factors for a finite crack in an infinite and linear elastic FGM subjected to an impact anti-plane crack-face loading are presented and discussed. The effects of the material gradients of the FGM on the transient dynamic stress intensity factors and their dynamic overshoot over the corresponding static stress intensity factors are analyzed.     

Elasticity solution for an FGM cylindrical panel integrated with piezoelectric layers

The three-dimensional elasticity solution for static analysis of a functionally graded material (FGM) cylindrical panel with simply supported edges is developed. The modulus of elasticity varies continuously throughout the thickness direction in the form of an exponential function. The panel is bonded with piezoelectric layers. Using Fourier series expansions in the axial and circumferential directions, the state equations are derived. The stress, displacement and electric potential distributions are obtained by solving these state equations. The influences of the material gradient index, applied voltage, and radius to thickness ratio on the static behavior of FGM shell are also studied.     

The mixed mode crack problem in an FGM layer bonded to a homogeneous half-plane

In this paper, the plane elasticity problem of an arbitrarily oriented crack in a FGM layer bonded to a homogeneous half-plane is considered. The problem is modeled by assuming that the elastic properties of the FGM layer are exponential functions of the thickness coordinate and are continuous at the interface of the FGM layer and the half-plane.The Fourier transform technique is used to reduce the problem to the solution of a system of Cauchy-type singular integral equations, which are solved numerically. The stress intensity factors are computed for various crack orientations, crack locations and material parameters. The results show that crack length, crack orientation and the non-homogeneity parameter of the strip material have significant effect on the fracture of the FGM layer.     

Moving crack for functionally grated material in an infinite length strip under antiplane shear

Considered is a Yoffe crack in an infinite strip of functionally grated material (FGM) subjected to antiplane shear. The shear moduli in two directions of FGM are assumed to be of exponential form. The dynamic stress intensity factor and strain energy density factor at the crack tip are obtained by using integral transforms and dual-integral equations. The numerical results show that the decrease of the strain energy density factor varies with the shear moduli gradient, and the increase of the strain energy density factor varies with the increase of the moving crack speed. The ratio of shear moduli in material vertical orientation has a great influence on the strain energy density factor.     

功能梯度材料结构分析的半解析数值方法研究

一种典型的半解析数值方法——线法被引入功能梯度材料的结构分析。首先推导了功能梯度材料位移形式的平衡方程和边界条件,然后阐述了线法功能梯度材料结构分析的基本步骤和数值原理。该方法的基本思想是通过有限差分将问题的控制方程半离散为定义在沿梯度方向离散节线上的常微分方程组,然后应用B样条函数Gauss配点法求解该常微分方程组得到问题的解答。为演示线法在功能梯度材料结构分析中的应用,给出了线性梯度和指数梯度功能梯度材料板分别受恒定位移、均匀拉伸载荷和弯曲载荷作用的数值算例。与相应问题解析解和其他数值方法的比较表明,线法的计算结果具有很高的精度,而且不需要任何特殊的考虑就能够有效模拟材料内部物性参数的连续变化,也无需事先选取满足特定条件的待定场函数,是一种非常适合功能梯度材料结构形式和材料特点的半解析数值方法。     

Three-dimensional stress and free vibration analyses of functionally graded plates with circular holes by the use of the graded finite element method

Static and free vibration analyses of plates with circular holes are performed based on the three-dimensional theory of elasticity. The plates are made of a functionally graded material (FGM), and the volume fractions of the constituent materials vary continuously across the plate. The effective properties of the FGM plate are estimated by using the Mori–Tanaka homogenization method. A graded finite element method based on the Rayleigh–Ritz energy formulation is used to solve the problem. Effects of different volume fractions of the materials and hole sizes on the behavior of FGM plates under uniaxial tension are investigated. Natural frequencies of a fully clamped FGM plate with a circular cutout are derived. The results obtained are compared with available experimental data.