基于Matlab求解时谐复系数弹性波方程的超弱变分方法  

AN ULTRA WEAK VARIATIONAL FORMULATION FOR SOLVING TIME-HARMONIC ELASTIC WAVE EQUATIONS WITH COMPLEX WAVE COEFFICIENTS BASED ON MATLAB

作　　者：袁龙 赵茂先 YUAN Long;ZHAO Maoxian(College of Mathematics and Systems Science,Shandong University of Science and Technology,Qingdao 266590,Shandong,China)

机构地区：[1]山东科技大学数学与系统科学学院,山东青岛266590

出　　处：《力学与实践》2021年第2期289-293,共5页Mechanics in Engineering

基　　金：国家自然科学基金(11501529);山东省自然科学基金面上项目(ZR2020MA046);青岛市应用基础研究项目(17-1-1-9-jch)资助。

摘　　要：众所周知,弹性波方程被用来描述工程应用中的弹性波传输问题。高效数值求解弹性波方程是一个重要的研究课题。提出了基于平面波离散的间断有限元超弱变分法。首先,选取满足剖分单元上齐次偏微分方程和其对偶方程的解析解分别作为检验空间和测试空间,得到等价的间断有限元变分形式;然后分别定义二维空间下的矢量平面波检验离散空间和测试离散空间,进而得到离散变分形式。基于Matlab软件的数值实验验证了新算法具有较高的精度。As is known,some engineering problems,such as the elastic wave transmission,can be described by the elastic wave equation.The efficient numerical algorithms for solving elastic wave equations are important.An ultra weak variational formulation based on the plane wave discretizations is proposed.First,by choosing the analytic solutions satisfying the homogeneous elastic wave equations and its adjoint equations to form the trial space and the test space,respectively,we obtain the equivalent Trefftz variational formulation;then we define the two-dimensional vector plane wave discretized trial space and test space to discretize the proposed continuous variational formulation.Numerical results based on Matlab have validated the proposed method.

关 键 词：弹性波方程 超弱变分方法 平面波 

分 类 号：G642[文化科学—高等教育学] O24[文化科学—教育学]