考虑稳定性约束的五零能模式材料结构设计方法研究

五零能模式材料是一种声学超材料,虽然它一般采用固体材料设计,并且单胞常应用微桁架构型,但理论上这种材料设计的结构是不稳定的。本文从分析设计结构整体出发,探讨一种通过施加预载荷,确保五零能模式材料结构稳定性的方法,并给出了保证结构稳定性的充分条件,且应用典型的五零能模式材料单胞构型加以验证。同时,为方便提出的稳定性条件在设计过程中的实施,将五零能模式材料结构的设计过程描述为五零能模式材料的变换过程,采用结构尺寸优化模型来描述。当材料单胞为微桁架构型时,这一结构尺寸优化模型简单,结构稳定性约束和五零能模式材料性质约束均为线性约束。     

两种五零能模式材料单胞构型及其性能分析

五零能模式材料是一种新型的人工超材料,虽属于弹性材料,但组成其单胞的特殊构型使其宏观静态表现为仅能承载一种受力状态,动态表现为仅能传播一种弹性波。本文首先构造了两种五零能模式材料的单胞构型,其具有不同的弹性特性,其中一种材料可传播弹性膨胀波,另一种可传播弹性剪切波。然后分别采用代表体元法和均匀化法分析这两种单胞的等效弹性模量。五零能模式材料的分析分为两步更直观,开始从单胞桁架模型入手,检验单胞构型是否满足五零能模式的定义,然后分析单胞实体模型,考察单胞构型的结构参数与其等效弹性模量的关系。研究表明对于这种低密度弹性材料的分析,代表体元法更适合。     

基于声学超材料的穿孔板缺陷态俘能特性研究

论文提出了一套适用于在低频段区间工作的声学超材料俘能装置.通过在含孔的声学超材料结构中制造一个局域共振缺陷态,进而将入射声波的弹性应变能集中在缺陷区域中,并利用压电晶片实现能量的转化.论文采用有限元的方法研究了元胞含孔的声学超材料在共振频率下俘获功率与电压的性能.进一步地,通过逐渐改变孔的尺寸,探究了孔径大小对俘能效果...     

基于拓扑描述函数的特定性能材料设计方法

研究了拓扑描述函数在材料设计中的应用,给出了一种基于拓扑描述函数的特定性能材料设计问题的提法和求解方法.将拓扑描述函数表示成含参数的基函数之和,将材料微结构拓扑优化问题转化为设计基函数描述参数的尺寸优化问题,使问题求解更方便.基于拓扑描述函数的方法可以准确确定设计域上任意点的材料分布,避免了变密度法常出现的棋盘效应、设计变量和有限元单元相关的缺点;与传统的水平集方法相比,其优化模型可以利用现有的优化方法求解,避免了差分法求偏微分方程.具有正泊松比和负泊松比的特定弹性性能材料的设计算例,说明了基于拓扑描述函数的材料设计方法的有效性.     

旋转壳的数值传递函数方法

应用数值传递函数方法建立一种用于分析旋转壳静力、动力响应的截锥壳单元,在本方法中,单元的位移在环向展开为Fourier级数的形式,应用薄壳理论可以得到解耦的微分方程,通过Laplace变换可以将方程转化为频域内的常微分方程,将其表示为状态空间形式后,可以应用数值传递函数方法求解,对复杂的系统可以应用与有限元类似的方法,划分多个单元组合求解,文中给出了几种旋转壳的动力、静力问题的数值算例,并与其它方法进行了比较,表明本文方法具有精度高,计算方便等特点。     

二阶非自伴两点边值问题Garlerkin有限元的超收敛算法

提出了基于改进位移模式的二阶非自伴两点边值问题Garlerkin有限元的超收敛算法. 用常规有限元解的位移模式与高阶有限元解的位移模式之和构造新的位移模式,基于Garlerkin 方法,采用积分形式推导了单元平衡方程. 对于线性单元,本文给出了有代表性的算例,结点和单元的位移、导数都达到了h⁴阶的超收敛精度.     

一类可压缩超弹性柱体中微孔的径向增长

研究了一类可压缩超弹性材料组成的柱体在给定的径向拉伸作用下预存微孔的增长问题.利用能量变分原理将此类问题的数学模型归结为一个二阶非线性常微分方程的边值问题,并通过变数变换求得了方程的参数型解析解,得到了微孔的增长与给定的拉伸之间的精确关系.通过相应的数值算例,解释了材料参数和结构参数对微孔增长的影响,分析了柱体中应力的分布及径向位移的变化情况.     

非线性动力学常微分方程组高精度数值积分方法

建立了一种求解非线性动力学常微分方程组初值问题的新方法．若非线性函数一阶导数存在，则给出解的积分方程表达式，计算得到按规定误差要求的高精度数值解．引入一般自治或非自治非线性系统的首次近似Jacobi矩阵，不作任何假设重构等价的非线性常微分方程组，简捷而有广泛的适应性，不改变方程的本质，但其主项构成线性化方程组，其它项则代表非线性函数高阶余项而不涉及Taylor级数展开计算，给出该方程组初值问题的Duhamel卷积分解析表达式，在时间步长内进行数值积分选代求解，在指定误差内快速收敛，逐步递推获得非线性常微分方程的瞬态响应和全时域高精度数值解．积分解连续满足微分方程组而不是在离散的步长端点上满足代数方程组，打破了传统用增量法在离散点上建立的代数方程组迭代求解，从而使传统Euler型逐步积分法的各种差分格式算法改变成真正的积分格式算法．数值计算中给出指数矩阵递增展开式，变矩阵乘法为乘积系数的加法，避免了大量矩阵自乘而大大提高计算效率．算法验证为无条件稳定，则保证对线性常微分方程而言，计算中舍入误差的传播不会扩散，不出现计算机字长有限而引起舍入误差导致计算不确定性问题．基于以上理论和数值方法，计算了线性非线性算例并进行了分析，验证了本方法简捷而有广泛的适应性，可以有足够的精确性．     

三维声学特征频率灵敏度分析的边界元法

声系统特征频率的灵敏度分析为其优化设计提供了基础,具有重要意义。边界元法在声学问题的求解中具有独特优势,但因其系统方程系数矩阵的频率相关性导致的非线性特征值问题给声学特征频率的灵敏度分析带来了很大困难。为此,本文首先对非线性特征值问题进行了线性化处理,利用围道积分投影方法将非线性特征方程转换为小规模广义特征方程,然后对其关于设计变量直接求导,并引入左特征向量和转换矩阵构造了一种适用于内外声场的三维声学单/重特征频率灵敏度分析的边界元法。数值算例验证了该方法的适用性,以及对单/重特征频率灵敏度的计算精度。     

低速冲击下压电层合板的非线性动力响应与疲劳损伤寿命分析

基于损伤力学,考虑压电效应、几何非线性和损伤演化的影响,建立了集中载荷作用下压电层合板的非线性运动微分方程,考虑撞击物与压电层合板之间的弹性接触效应,确定了层合板所承受的冲击力.对此非线性动力问题,采用有限差分法进行求解,且引入损伤层合板单元模型,采用有限元方法模拟了冲击荷载作用下压电层合板的损伤演化.算例中,讨论了撞击物的冲击速度、压电材料以及损伤对压电层合板非线性动力响应的影响,以及周期冲击荷载下冲击速度和压电材料对压电层合板疲劳损伤寿命的影响.     

A novel design scheme for acoustic cloaking of complex shape based on region partitioning and multi-origin coordinate transformation

Acoustic cloaking is an important application of acoustic metamaterials. This article proposes a novel design scheme for acoustic cloaking based on the region partitioning and multi-origin coordinate transformation. The cloaked region is partitioned into multiple narrow strips. For each strip, a local coordinate system is established with the local origin located at the strip center, and a coordinate transformation in the local coordinate system is conducted to squeeze the material along the strip length direction to form the cloaked region. To facilitate the implementation of the acoustic cloak, the multilayer effective medium is used to approximate the non-uniform anisotropic material parameters. The effectiveness of the proposed coordinate transformation method is verified by comparing the results from our method with those in the literature. Firstly, the results of a circular acoustic cloak in the literature are reproduced by using our finite element (FE) simulations for validation. Then, a comparison is made between the traditional coordinate transformation scheme and our new scheme for simulating an elliptical acoustic cloak. The results indicate that the proposed multi-origin coordinate transformation method has a better cloaking effect on the incident wave along the ellipse minor axis direction than the traditional method. This means that for the same object, an appropriate transformation scheme can be selected for different incident wave directions to achieve the optimal control effect. The validated scheme is further used to design an arch-shaped cloak composed of an upper semicircular area and a lower rectangular area, by combining the traditional single-centered coordinate transformation method for the semicircular area and the proposed multi-origin method for the rectangular area. The results show that the designed cloak can effectively control the wave propagation with significantly reduced acoustic pressure level. This work provides a flexible acoustic cloak design method applicable for arbitrary shapes and different wave incident directions, enriching the theory of acoustic cloaking based on coordinate transformation.     

SPACE-TIME FINITE ELEMENT METHOD FOR SCHRDINGER EQUATION AND ITS CONSERVATION

Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.     

SPACE-TIME FINITE ELEMENT METHOD FOR SCHRODINGER EQUATION AND ITS CONSERVATION

Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.     

敷设多孔吸声材料声腔特征值分析的径向积分边界元法

由于Helmholtz方程的基本解是频率的函数,因此传统边界元法在处理声场特征值问题时具有天生的缺陷。本文采用Laplace方程基本解生成积分方程,通过径向积分法将在此过程中产生的域积分项转化为边界积分。此方法克服了传统边界元法系数矩阵对频率的依赖,同时克服了特解积分法对特解的依赖,并通过对表面声导纳的多项式逼近,将敷设多孔吸声材料声腔特征值问题转化为矩阵多项式,从而避免了复杂的非线性求解。通过数值算例验证了算法的有效性。     

非规则区域中拉普拉斯方程的有限分析5点格式

建立了非规则区域的有限分析5点格式,增加了有限分析法对不规则边界的适应性。应用所提出的方法对水利工程中常见的有压和无压流动进行了计算,与实验和前人的计算结果相比较,本文的方法都能得到较为满意的结果。本文的计算格式也可以应用到其他非规则区域的计算中。     

Sliding frictional contact between a rigid punch and a laterally graded elastic medium

Analytical and computational methods are developed for contact mechanics analysis of functionally graded materials (FGMs) that possess elastic gradation in the lateral direction. In the analytical formulation, the problem of a laterally graded half-plane in sliding frictional contact with a rigid punch of an arbitrary profile is considered. The governing partial differential equations and the boundary conditions of the problem are satisfied through the use of Fourier transformation. The problem is then reduced to a singular integral equation of the second kind which is solved numerically by using an expansion–collocation technique. Computational studies of the sliding contact problems of laterally graded materials are conducted by means of the finite element method. In the finite element analyses, the laterally graded half-plane is discretized by quadratic finite elements for which the material parameters are specified at the centroids. Flat and triangular punch profiles are considered in the parametric analyses. The comparisons of the results generated by the analytical technique to those computed by the finite element method demonstrate the high level of accuracy attained by both methods. The presented numerical results illustrate the influences of the lateral nonhomogeneity and the coefficient of friction on the contact stresses.     

FINITE ANALYTIC NUMERICAL SOLUTION OF NAVIER-STOKES EQUATIONS

The finite analytic method is used in the present study to calculate the turbulent flow field described with Navier-Stokes equation in body-fitted curvilinear coordinate system The finite analytic method invokes the analytic solution of governing partial differential equation in formulating the algebraic equation that relates a nodal value in an element to its neighbour nodal values according to the direction and the magnitude of convection. It is shown that the finite analytic method has good numerical stability and accuracy. The turbulent flow fields through a single and a tandem cascades of airfoil are numerically simulated by using finite analytic method respectively in this paper. The k ? ε turbulence model and wall function are employed in. the present study. The agreement of numerical solution with experiment result is quite good.     

硬脑膜受压膨出挠度的力学分析

硬脑膜是一种粘弹性材料,为控制硬脑膜在脑压作用下的膨出度,对粘弹性薄膜受压膨出挠度作力学分析。以位移为未知量,从粘弹性材料的分型本构关系出发将Foepple薄膜大挠度理论从弹性推广到粘弹性膜,得到一组非线性积分偏微分方程。先在空间上运用Galerkin方法将积分偏微分方程组化为积分常微分方程组。然后,在时间域上运用数值积分和有限差分将方程离散为非线性代数方程组。本文对四周固定夹紧的圆形、椭圆形和矩形薄膜进行了求解,并将求解结果用于颅底缺损重建膜的膨出量计算,计算值与实验值吻合,为颅底外科提供一个理论分析方法。     

Design and analysis of nonlinear-transformation-based broadband cloaking for acoustic wave propagation

A coordinate-transformation method can be used to design invisibility cloaks for many types of waves, including acoustic waves. The traditional method for designing a cloak depends on a transformation from a virtual space to a physical space. Previous acoustic cloaks that are mainly designed with linear-transformation-based acoustics have drawbacks that acoustic wave trajectories in the cloaks cannot be controlled and tuned. This work uses a nonlinear mapping from a ray trajectory perspective to construct acoustic cloaks with tunable non-singular material properties. Use of a ray trajectory equation is a straightforward and alternate way to study propagation characteristics of different types of waves, which allows more flexibility in controlling the waves. A broadband cylindrical cloak for acoustic waves in an inviscid fluid is realized with layered non-singular, homogeneous, and isotropic materials based on a nonlinear transformation. Some advantages and improvements of the invisibility nonlinear-transformation cloak over a traditional linear-transformation cloak are analyzed. The invisibility capability of the nonlinear-transformation cloak can be tuned by adjusting a design parameter that is shown to have influence on the acoustic wave energy flowing into the region inside the cloak. Numerical examples show that the nonlinear-transformation cloak is more effective for making a domain undetectable by acoustic waves in an inviscid fluid and shielding acoustic waves from outside the cloak than the linear-transformation cloak in a broad frequency range. The methodology developed here can be used to design nonlinear-transformation cloaks for other types of waves.