CN107315895A - 一种石墨烯纳米质量传感器振动特性分析的半解析方法 - Google Patents

一种石墨烯纳米质量传感器振动特性分析的半解析方法 Download PDF

Info

Publication number
CN107315895A
CN107315895A CN201610246796.8A CN201610246796A CN107315895A CN 107315895 A CN107315895 A CN 107315895A CN 201610246796 A CN201610246796 A CN 201610246796A CN 107315895 A CN107315895 A CN 107315895A
Authority
CN
China
Prior art keywords
mrow
msup
mtd
msub
mfrac
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
Application number
CN201610246796.8A
Other languages
English (en)
Other versions
CN107315895B (zh
Inventor
申志彬
姜人伟
唐国金
李飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201610246796.8A priority Critical patent/CN107315895B/zh
Publication of CN107315895A publication Critical patent/CN107315895A/zh
Application granted granted Critical
Publication of CN107315895B publication Critical patent/CN107315895B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00—Subject matter not provided for in other main groups of this subclass

Landscapes

  • Force Measurement Appropriate To Specific Purposes (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明公开了一种石墨烯纳米质量传感器振动特性分析的半解析方法，步骤包括：将石墨烯纳米质量传感器在x方向上沿附加质量的质心划分两个子结构，建立各个子结构的控制方程；将每一个子结构采用NE+1条结线划分为NE个矩形区域，每一个矩形区域为一个条形单元，根据每一个子结构的控制方程推导出该子结构中各个条形单元的运动微分方程，针对各个子结构组装总体运动微分方程并处理边界条件，求解总体运动微分方程得到石墨烯纳米质量传感器自由振动的固有频率ω。本发明能够处理复杂几何形状问题、解的形式统一、易于利用计算机编程，可用于相对较复杂的几何区域以及边界条件问题的求解，具有求解精度高、计算过程数据存储量少、计算效率高的优点。

Description

一种石墨烯纳米质量传感器振动特性分析的半解析方法
技术领域
本发明涉及弹性力学领域，具体涉及一种石墨烯纳米质量传感器振动特性分析的半解析方法。
背景技术
在实际工程问题中，只有极少数力学问题能够给出解析解，大多数问题需要通过数值方法来求解，含附加质量的石墨烯传感器振动特性的求解亦如此。在众多数值方法中，有限元法最常见，也是最实用的一种方法。有限元法最大的优点是不受求解区域、边界条件以及材料属性的限制，可以分析具有复杂几何形状的弹性力学问题。但是，石墨烯纳米质量传感器几何尺寸一般在纳米量级，要考虑尺度效应，因此需要引入微分型非局部本构模型，况且石墨烯纳米质量传感器还需要考虑附加集中质量的惯性力，这使有限元方法的求解效率以及求解精度大大降低。
条形传递函数方法是一种求解二维弹性力学问题的半解析数值方法。这种方法的思想类似于有限条法，也是将求解区域分为若干个条形区域，称为条形单元，在条形单元内利用多项式和连续函数近似横向和纵向位移，从而得到基于条形单元的整体微分方程，最后利用传递函数方法求解微分方程，得到半解析解。该方法的一个显著优点是，其既具有有限元方法的灵活性，可以分析复杂形状的几何区域，同时又能给出封闭形式的高精度半解析解，而且在求解含附加质量问题时候无需在控制方程中考虑附加质量项，从而给求解带来极大的便利性。然而，传统的条形传递函数方法基于Hamilton原理，需要先给出待求问题对应的能量泛函。而问题是，并非所有问题都可以很容易地给出其相应的能量泛函，如考虑尺度效应的微分型非局部本构模型的薄板弯曲问题，这使得条形传递函数方法的应用受到了限制。然而，Galerkin方法不需要先写出待研究问题的能量泛函，可以直接对微分方程进行近似求解。
发明内容
本发明要解决的技术问题：针对现有技术的上述问题，提供一种能够处理复杂几何形状问题、解的形式统一、易于利用计算机编程，可用于相对较复杂的几何区域以及边界条件问题的求解，求解精度高、计算过程数据存储量少、计算效率高的石墨烯纳米质量传感器振动特性分析的半解析方法。
为了解决上述技术问题，本发明采用的技术方案为：
一种石墨烯纳米质量传感器振动特性分析的半解析方法，步骤包括：
1)将包含附加质量的矩形结构的石墨烯纳米质量传感器在x方向上沿附加质量的质心划分为两个子区域Ω1和Ω2得到对应的两个子结构，建立各个子结构的控制方程；
2)将每一个子结构采用NE+1条结线划分为NE个矩形区域，每一个矩形区域为一个条形单元，第j个条形单元包含第j条结线和第j+1条结线以及4个结点；
3)根据每一个子结构的控制方程推导出该子结构中各个条形单元的运动微分方程；
4)针对各个子结构，将所有条形单元的运动微分方程进行组装总体运动微分方程，并针对所述总体运动微分方程处理边界条件；
5)对处理边界条件后的总体运动微分方程计算传递矩阵和边界矩阵，根据传递矩阵和边界矩阵求解处理边界条件后总体运动微分方程，得到石墨烯纳米质量传感器自由振动的固有频率ω并作为石墨烯纳米质量传感器的振动特性分析结果输出。
优选地，所述步骤1)中建立的各个子结构的控制方程如式(1)所示；
式(1)中，Ω1表示子区域Ω1对应的子结构，Ω2表示子区域Ω2对应的子结构，Mxx1为子区域Ω1在x方向的弯矩，Mxy1为子区域Ω1的扭矩，Myy1为子区域Ω1在y方向的弯矩，Mxx2为子区域Ω2在x方向的弯矩，Mxy2为子区域Ω2的扭矩，Myy2为子区域Ω2在y方向的弯矩，ρ为石墨烯的密度，h为石墨烯纳米质量传感器的厚度，w1为子区域Ω1对应的子结构对应的横向位移，w2为子区域Ω2对应的子结构对应的横向位移，t为时间，x为横坐标、y为纵坐标。
优选地，所述步骤1)中将每一个子结构采用NE+1条结线划分为NE个矩形区域后，每一个条形单元结线位移函数向量如式(2)所示、内部的横向位移函数如式(3)所示；
φ(x,t)＝{wj θj wj+1 θj+1}T (2)
式(2)中，φ(x,t)为第j个条形单元的结线位移函数向量，wj为第j条结线的横向位移，θj为第j条结线的转角，wj+1为第j+1条结线的横向位移，θj+1为第j+1条结线的转角；
w(x,y,t)＝N(y)φ(x,t) (3)
式(3)中，w(x,y,t)为第j个条形单元内部的横向位移函数，N(y)为选定的形函数矩阵，φ(x,t)为第j个条形单元的结线位移函数向量。
优选地，所述步骤3)推导出各个条形单元的运动微分方程如式(4)所示；
式(4)中，和为条形单元的刚度矩阵，和为条形单元的质量矩阵，φi为条形单元的结线位移函数向量，t为时间，x为横坐标、y为纵坐标，D为石墨烯纳米质量传感器弯曲刚度，N为选定的形函数矩阵N(y)，T为矩阵转置符号，ν为泊松比，ρ为石墨烯的密度，h为石墨烯纳米质量传感器的厚度，l为条形单元的宽度，μ为非局部参数。
优选地，所述选定的形函数矩阵N(y)为标准Euler梁单元的形函数。
优选地，所述步骤4)的详细步骤包括：
4.1)针对各个子结构，采用有限元法将各个条形单元的刚度矩阵和分别组装成总体刚度矩阵K(4)、K(2)和K(0)，将各个条形单元的质量矩阵和分别组装成总体质量矩阵M(2)和M(0)；
4.2)确定各个子结构的总体位移向量如式(5)所示；
Φi(x,t)＝{w1(x,t),θ1(x,t),w2(x,t),θ2(x,t),…,wNE(x,t),θNE(x,t)}T,i＝1,2 (5)
式(5)中，Φi(x,t)为第i个子结构的总体位移向量，w1(x,t)为第i个子结构中第1个条形单元的结线横向位移，θ1(x,t)为第i个子结构中第一个条形单元的结线转角，w2(x,t)为第i个子结构中第2个条形单元的结线横向位移，θ2(x,t)为第i个子结构中第2个条形单元的结线转角，wNE(x,t)为第i个子结构中第NE个条形单元的结线横向位移，θNE(x,t)为第i个子结构中第NE个条形单元的结线转角；
4.3)根据各个子结构的总体位移向量组装总体运动微分方程如式(6)所示；
式(6)中，Φi为第i个子结构的总体位移向量，K(4)、K(2)和K(0)为总体刚度矩阵，M(2)和M(0)为总体质量矩阵，t为时间，x为横坐标；
4.4)针对所述总体运动微分方程处理边界条件，得到总体运动微分方程如式(7)所示；
式(7)中，为第i个子结构处理边界条件后的总体位移向量，和为处理边界条件后的总体刚度矩阵，和为处理边界条件后的总体质量矩阵。
优选地，所述步骤5)的详细步骤包括：
5.1)将各个子结构处理边界条件后的总体运动微分方程取时间的Laplace变换，得到各个子结构的Laplace总体运动微分方程如式(8)所示；
式(8)中，为第i个子结构处理边界条件后的总体位移向量的Laplace变换，‘^’为Laplace变换符号，s为Laplace变换系数，和为处理边界条件后的总体刚度矩阵，为处理边界条件后的总体刚度矩阵，D(s)、D2、A均为中间矩阵；
5.2)基于定义的状态向量η(x,s)，将Laplace总体运动微分方程转换为式(9)所示：
式(9)中，η(x,s)为定义的状态向量，且定义的状态向量η(x,s)的函数表达式如式(10)所示，F(s)为传递矩阵，传递矩阵F(s)的函数表达式如式(11)所示；
式(10)中，η(x,s)为定义的状态向量，η1(x,s)和η2(x,s)均为中间向量，为第1个子结构处理边界条件后的总体位移向量的Laplace变换，为第2个子结构处理边界条件后的总体位移向量的Laplace变换；
式(11)中，F(s)为传递矩阵，为中间矩阵，为4×N1阶0元素方阵，N1为每个子结构的总体自由度数，D(s)、D2为式(8)中所示的中间矩阵，I为单位矩阵；
5.3)建立两个子结构连接处的边界条件如式(12)所示；
Mb(s)η(-0.5a,s)+Nb(s)η(0.5a,s)+Rb(s)η(x0,s)＝0 (12)
式(12)中，Mb(s)、Nb(s)和Rb(s)为边界矩阵，η为状态向量，a为石墨烯质量传感器在x方向上的长度，s为Laplace变换系数，x0为附加质量的质心的横坐标；其中边界矩阵Rb(s)如式(13)所示；
式(13)中，Rb(s)为边界矩阵，为4×N1阶0元素方阵，N1为每个子结构的总体自由度数，rb1和rb2均为中间矩阵，m0为附加质量的大小，I为单位矩阵，为处理边界条件后的总体刚度矩阵；
5.4)生成将式(9)所示转换后的Laplace总体运动微分方程、式(12)所示两个子结构连接处的边界条件的解如式(14)所示，且所述解的特征方程如式(15)所示；
式(14)中，Mb(s)、Nb(s)和Rb(s)为边界矩阵，a为石墨烯质量传感器在x方向上的长度，F(s)为传递矩阵，x0为附加质量的质心的横坐标，x为横坐标，s为Laplace变换系数；
式(15)中，det为求行列式符号，Mb(s)、Nb(s)和Rb(s)为边界矩阵，a为石墨烯质量传感器在x方向上的长度，F(s)为传递矩阵，x0为附加质量的质心的横坐标，x为横坐标，s为Laplace变换系数；
5.5)针对式(15)所示解的特征方程，令Laplace变换系数s为虚数单位i、石墨烯纳米质量传感器自由振动的固有频率ω两者之间的乘积，求解得到石墨烯纳米质量传感器自由振动的固有频率ω并作为石墨烯纳米质量传感器的振动特性分析结果输出。
本发明石墨烯纳米质量传感器振动特性分析的半解析方法具有下述优点：
1、本发明将包含附加质量的矩形结构的石墨烯纳米质量传感器在x方向上沿附加质量的质心划分为两个子区域得到对应的两个子结构，建立各个子结构的控制方程，因此无需在控制方程中考虑附加质量项，从而使问题容易求解、提高求解效率，具有求解精度高、计算过程数据存储量少、计算效率高的优点。
2、相对于其他解析或半解析方法，本发明石墨烯纳米质量传感器振动特性分析的半解析方法适用于多种边界条件，并且能够处理复杂几何形状问题；
3、本发明石墨烯纳米质量传感器振动特性分析的半解析方法解的形式统一，易于利用计算机编程，并且可以用于相对较复杂的几何区域以及边界条件问题的求解。
附图说明
图1为本发明实施例方法的基本流程示意图。
图2为本发明实施例中划分子结构的示意图。
图3为本发明实施例中将子结构划分为NE个矩形区域的示意图。
图4为本发明实施例中划分得到的一个条形单元的结构示意图。
具体实施方式
如图1所示，本实施例石墨烯纳米质量传感器振动特性分析的半解析方法的步骤包括：
1)将包含附加质量的矩形结构的石墨烯纳米质量传感器在x方向上沿附加质量的质心划分为两个子区域Ω1和Ω2得到对应的两个子结构(参见图2)，建立各个子结构的控制方程；
2)将每一个子结构采用NE+1条结线划分为NE个矩形区域，每一个矩形区域为一个条形单元，第j个条形单元包含第j条结线和第j+1条结线以及4个结点，如图3和图4所示，第j个条形单元的宽度为l，Oxy为条形单元局部坐标系；
3)根据每一个子结构的控制方程推导出该子结构中各个条形单元的运动微分方程；
4)针对各个子结构，将所有条形单元的运动微分方程进行组装总体运动微分方程，并针对总体运动微分方程处理边界条件；
5)对处理边界条件后的总体运动微分方程计算传递矩阵和边界矩阵，根据传递矩阵和边界矩阵求解处理边界条件后总体运动微分方程，得到石墨烯纳米质量传感器自由振动的固有频率ω并作为石墨烯纳米质量传感器的振动特性分析结果输出。
本实施例中，步骤1)中建立的各个子结构的控制方程如式(1)所示；
式(1)中，Ω1表示子区域Ω1对应的子结构，Ω2表示子区域Ω2对应的子结构，Mxx1为子区域Ω1在x方向的弯矩，Mxy1为子区域Ω1的扭矩，Myy1为子区域Ω1在y方向的弯矩，Mxx2为子区域Ω2在x方向的弯矩，Mxy2为子区域Ω2的扭矩，Myy2为子区域Ω2在y方向的弯矩，ρ为石墨烯的密度，h为石墨烯纳米质量传感器的厚度，w1为子区域Ω1对应的子结构对应的横向位移，w2为子区域Ω2对应的子结构对应的横向位移，t为时间，x为横坐标、y为纵坐标。
本实施例中，步骤1)中将每一个子结构采用NE+1条结线划分为NE个矩形区域后，每一个条形单元结线位移函数向量如式(2)所示、内部的横向位移函数如式(3)所示；
φ(x,t)＝{wj θj wj+1 θj+1}T (2)
式(2)中，φ(x,t)为第j个条形单元的结线位移函数向量，wj为第j条结线的横向位移，θj为第j条结线的转角，wj+1为第j+1条结线的横向位移，θj+1为第j+1条结线的转角；
w(x,y,t)＝N(y)φ(x,t) (3)
式(3)中，w(x,y,t)为第j个条形单元内部的横向位移函数，N(y)为选定的形函数矩阵，φ(x,t)为第j个条形单元的结线位移函数向量。本实施例中，选定的形函数矩阵N(y)为标准Euler梁单元的形函数N＝[N1 N2 N3 N4]。
本实施例中，步骤3)推导出各个条形单元的运动微分方程如式(4)所示；
式(4)中，和为条形单元的刚度矩阵，和为条形单元的质量矩阵，φi为条形单元的结线位移函数向量，t为时间，x为横坐标、y为纵坐标，D为石墨烯纳米质量传感器弯曲刚度，N为选定的形函数矩阵N(y)，T为矩阵转置符号，ν为泊松比，ρ为石墨烯的密度，h为石墨烯纳米质量传感器的厚度，l为条形单元的宽度，μ为非局部参数。
对于式(1)所示各个子结构的控制方程，其弯矩和扭矩的具体形式如式(4-1)所示；
式(4-1)中，μ为非局部参数，D为弹性矩阵，为Laplace算子，Mxxi为第i个子区域Ωi在x方向的弯矩，Mxyi为第i个子区域Ωi的扭矩，Myyi为第i个子区域Ωi在y方向的弯矩，wi为第i个子区域Ωi对应的子结构对应的横向位移。
根据式(4-1)，可以得到式(1)所示各个子结构的控制方程在y方向上的等效积分“弱”形式为式(4-2)所示；
式(4-2)中，l为条形单元的宽度，D为石墨烯纳米质量传感器弯曲刚度，为权函数，wi为第i个子区域Ωi对应的子结构对应的横向位移，ν为泊松比，x为横坐标、y为纵坐标，ρ为石墨烯的密度，h为石墨烯纳米质量传感器的厚度。
将式(3)所示内部的横向位移函数代入式(4-2)所示各个子结构的控制方程在y方向上的等效积分“弱”形式，并令得，其中，N＝[N1 N2 N3 N4]为形函数，取标准Euler梁单元的形函数，即可得到式(4-3)所示函数表达式，即可得到各个条形单元的运动微分方程如式(4)所示；
式(4-3)中，l为条形单元的宽度，D为石墨烯纳米质量传感器弯曲刚度，N为选定的形函数矩阵N(y)，φi为条形单元的结线位移函数向量，ν为泊松比，x为横坐标、y为纵坐标，ρ为石墨烯的密度，h为石墨烯纳米质量传感器的厚度。本实施例中，选定的形函数矩阵N(y)为标准Euler梁单元的形函数。
本实施例中，步骤4)的详细步骤包括：
4.1)针对各个子结构，采用有限元法将各个条形单元的刚度矩阵和分别组装成总体刚度矩阵K(4)、K(2)和K(0)，将各个条形单元的质量矩阵和分别组装成总体质量矩阵M(2)和M(0)；
4.2)确定各个子结构的总体位移向量如式(5)所示；
Φi(x,t)＝{w1(x,t),θ1(x,t),w2(x,t),θ2(x,t),…，wNE(x,t),θNE(x,t)}T,i＝1,2 (5)
式(5)中，Φi(x,t)为第i个子结构的总体位移向量，w1(x,t)为第i个子结构中第1个条形单元的结线横向位移，θ1(x,t)为第i个子结构中第一个条形单元的结线转角，w2(x,t)为第i个子结构中第2个条形单元的结线横向位移，θ2(x,t)为第i个子结构中第2个条形单元的结线转角，wNE(x,t)为第i个子结构中第NE个条形单元的结线横向位移，θNE(x,t)为第i个子结构中第NE个条形单元的结线转角；
4.3)根据各个子结构的总体位移向量组装总体运动微分方程如式(6)所示；
式(6)中，Φi为第i个子结构的总体位移向量，K(4)、K(2)和K(0)为总体刚度矩阵，M(2)和M(0)为总体质量矩阵，t为时间，x为横坐标；
4.4)针对总体运动微分方程处理边界条件，得到总体运动微分方程如式(7)所示；
式(7)中，为第i个子结构处理边界条件后的总体位移向量，和为处理边界条件后的总体刚度矩阵，和为处理边界条件后的总体质量矩阵。
本实施例中，步骤5)的详细步骤包括：
5.1)将各个子结构处理边界条件后的总体运动微分方程取时间的Laplace变换，得到各个子结构的Laplace总体运动微分方程如式(8)所示；
式(8)中，为第i个子结构处理边界条件后的总体位移向量的Laplace变换，‘^’为Laplace变换符号，s为Laplace变换系数，和为处理边界条件后的总体刚度矩阵，为处理边界条件后的总体刚度矩阵，D(s)、D2、A均为中间矩阵；
5.2)基于定义的状态向量η(x,s)，将Laplace总体运动微分方程转换为式(9)所示：
式(9)中，η(x,s)为定义的状态向量，且定义的状态向量η(x,s)的函数表达式如式(10)所示，F(s)为传递矩阵，传递矩阵F(s)的函数表达式如式(11)所示；
式(10)中，η(x,s)为定义的状态向量，η1(x,s)和η2(x,s)均为中间向量，为第1个子结构处理边界条件后的总体位移向量的Laplace变换，为第2个子结构处理边界条件后的总体位移向量的Laplace变换；
式(11)中，F(s)为传递矩阵，为中间矩阵，为4×N1阶0元素方阵，N1为每个子结构的总体自由度数，D(s)、D2为式(8)中所示的中间矩阵，I为单位矩阵；
5.3)建立两个子结构连接处的边界条件如式(12)所示；
Mb(s)η(-0.5a,s)+Nb(s)η(0.5a,s)+Rb(s)η(x0,s)＝0 (12)
式(12)中，Mb(s)、Nb(s)和Rb(s)为边界矩阵，η为状态向量，a为石墨烯质量传感器在x方向上的长度，s为Laplace变换系数，x0为附加质量的质心的横坐标；其中边界矩阵Rb(s)如式(13)所示；
式(13)中，Rb(s)为边界矩阵，为4×N1阶0元素方阵，N1为每个子结构的总体自由度数，rb1和rb2均为中间矩阵，m0为附加质量的大小，I为单位矩阵，为处理边界条件后的总体刚度矩阵；
5.4)生成将式(9)所示转换后的Laplace总体运动微分方程、式(12)所示两个子结构连接处的边界条件的解如式(14)所示，且解的特征方程如式(15)所示；
式(14)中，Mb(s)、Nb(s)和Rb(s)为边界矩阵，a为石墨烯质量传感器在x方向上的长度，F(s)为传递矩阵，x0为附加质量的质心的横坐标，x为横坐标，s为Laplace变换系数；
式(15)中，det为求行列式符号，Mb(s)、Nb(s)和Rb(s)为边界矩阵，a为石墨烯质量传感器在x方向上的长度，F(s)为传递矩阵，x0为附加质量的质心的横坐标，x为横坐标，s为Laplace变换系数；
5.5)针对式(15)所示解的特征方程，令Laplace变换系数s为虚数单位i、石墨烯纳米质量传感器自由振动的固有频率ω两者之间的乘积，求解得到石墨烯纳米质量传感器自由振动的固有频率ω并作为石墨烯纳米质量传感器的振动特性分析结果输出。
以上所述仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (7)

1.一种石墨烯纳米质量传感器振动特性分析的半解析方法，其特征在于步骤包括：
1)将包含附加质量的矩形结构的石墨烯纳米质量传感器在x方向上沿附加质量的质心划分为两个子区域Ω1和Ω2得到对应的两个子结构，建立各个子结构的控制方程；
2)将每一个子结构采用NE+1条结线划分为NE个矩形区域，每一个矩形区域为一个条形单元，第j个条形单元包含第j条结线和第j+1条结线以及4个结点；
3)根据每一个子结构的控制方程推导出该子结构中各个条形单元的运动微分方程；
4)针对各个子结构，将所有条形单元的运动微分方程进行组装总体运动微分方程，并针对所述总体运动微分方程处理边界条件；
5)对处理边界条件后的总体运动微分方程计算传递矩阵和边界矩阵，根据传递矩阵和边界矩阵求解处理边界条件后总体运动微分方程，得到石墨烯纳米质量传感器自由振动的固有频率ω并作为石墨烯纳米质量传感器的振动特性分析结果输出。
2.根据权利要求1所述的石墨烯纳米质量传感器振动特性分析的半解析方法，其特征在于，所述步骤1)中建立的各个子结构的控制方程如式(1)所示；
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&amp;Omega;</mi> <mn>1</mn> </msub> <mo>:</mo> </mrow> </mtd> <mtd> <mrow> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>M</mi> <mrow> <mi>x</mi> <mi>x</mi> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <mn>2</mn> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>M</mi> <mrow> <mi>x</mi> <mi>y</mi> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>M</mi> <mrow> <mi>y</mi> <mi>y</mi> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>-</mo> <mi>&amp;rho;</mi> <mi>h</mi> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>w</mi> <mn>1</mn> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>t</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;Omega;</mi> <mn>2</mn> </msub> <mo>:</mo> </mrow> </mtd> <mtd> <mrow> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>M</mi> <mrow> <mi>x</mi> <mi>x</mi> <mn>2</mn> </mrow> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <mn>2</mn> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>M</mi> <mrow> <mi>x</mi> <mi>y</mi> <mn>2</mn> </mrow> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>M</mi> <mrow> <mi>y</mi> <mi>y</mi> <mn>2</mn> </mrow> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>-</mo> <mi>&amp;rho;</mi> <mi>h</mi> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>w</mi> <mn>2</mn> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>t</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
式(1)中，Ω1表示子区域Ω1对应的子结构，Ω2表示子区域Ω2对应的子结构，Mxx1为子区域Ω1在x方向的弯矩，Mxy1为子区域Ω1的扭矩，Myy1为子区域Ω1在y方向的弯矩，Mxx2为子区域Ω2在x方向的弯矩，Mxy2为子区域Ω2的扭矩，Myy2为子区域Ω2在y方向的弯矩，ρ为石墨烯的密度，h为石墨烯纳米质量传感器的厚度，w1为子区域Ω1对应的子结构对应的横向位移，w2为子区域Ω2对应的子结构对应的横向位移，t为时间，x为横坐标、y为纵坐标。
3.根据权利要求2所述的石墨烯纳米质量传感器振动特性分析的半解析方法，其特征在于，所述步骤1)中将每一个子结构采用NE+1条结线划分为NE个矩形区域后，每一个条形单元结线位移函数向量如式(2)所示、内部的横向位移函数如式(3)所示；
φ(x,t)＝{wj θj wj+1 θj+1}T (2)
式(2)中，φ(x,t)为第j个条形单元的结线位移函数向量，wj为第j条结线的横向位移，θj为第j条结线的转角，wj+1为第j+1条结线的横向位移，θj+1为第j+1条结线的转角；
w(x,y,t)＝N(y)φ(x,t) (3)
式(3)中，w(x,y,t)为第j个条形单元内部的横向位移函数，N(y)为选定的形函数矩阵，φ(x,t)为第j个条形单元的结线位移函数向量。
4.根据权利要求3所述的石墨烯纳米质量传感器振动特性分析的半解析方法，其特征在于，所述步骤3)推导出各个条形单元的运动微分方程如式(4)所示；
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>k</mi> <mi>e</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </msubsup> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>4</mn> </msup> <msub> <mi>&amp;phi;</mi> <mi>i</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>x</mi> <mn>4</mn> </msup> </mrow> </mfrac> <mo>+</mo> <msubsup> <mi>k</mi> <mi>e</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msubsup> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>&amp;phi;</mi> <mi>i</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <msubsup> <mi>k</mi> <mi>e</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msubsup> <msub> <mi>&amp;phi;</mi> <mi>i</mi> </msub> <mo>+</mo> <msubsup> <mi>m</mi> <mi>e</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msubsup> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>&amp;phi;</mi> <mi>i</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>t</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <msubsup> <mi>m</mi> <mi>e</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msubsup> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>4</mn> </msup> <msub> <mi>&amp;phi;</mi> <mi>i</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> <mo>&amp;part;</mo> <msup> <mi>t</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>k</mi> <mi>e</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </msubsup> <mo>=</mo> <mrow> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>l</mi> </msubsup> <mrow> <msup> <mi>DN</mi> <mi>T</mi> </msup> <mi>N</mi> <mi>d</mi> <mi>y</mi> </mrow> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>k</mi> <mi>e</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msubsup> <mo>=</mo> <mrow> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>l</mi> </msubsup> <mrow> <mi>D</mi> <mrow> <mo>&amp;lsqb;</mo> <mrow> <msup> <mi>&amp;nu;N</mi> <mi>T</mi> </msup> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mi>N</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>-</mo> <mn>2</mn> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <mi>&amp;nu;</mi> </mrow> <mo>)</mo> </mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <msup> <mi>N</mi> <mi>T</mi> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>N</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>+</mo> <mi>&amp;nu;</mi> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msup> <mi>N</mi> <mi>T</mi> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mi>N</mi> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mi>d</mi> <mi>y</mi> </mrow> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>k</mi> <mi>e</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msubsup> <mo>=</mo> <mrow> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>l</mi> </msubsup> <mrow> <mi>D</mi> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msup> <mi>N</mi> <mi>T</mi> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mi>N</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mi>d</mi> <mi>y</mi> </mrow> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>m</mi> <mi>e</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msubsup> <mo>=</mo> <mo>-</mo> <mrow> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>l</mi> </msubsup> <mrow> <msup> <mi>&amp;rho;h&amp;mu;N</mi> <mi>T</mi> </msup> <mi>N</mi> <mi>d</mi> <mi>y</mi> </mrow> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>m</mi> <mi>e</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msubsup> <mo>=</mo> <mrow> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>l</mi> </msubsup> <mrow> <msup> <mi>&amp;rho;hN</mi> <mi>T</mi> </msup> <mi>N</mi> <mi>d</mi> <mi>y</mi> </mrow> </mrow> <mo>-</mo> <mrow> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>l</mi> </msubsup> <mrow> <msup> <mi>&amp;rho;h&amp;mu;N</mi> <mi>T</mi> </msup> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mi>N</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mi>d</mi> <mi>y</mi> </mrow> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
式(4)中，和为条形单元的刚度矩阵，和为条形单元的质量矩阵，φi为条形单元的结线位移函数向量，t为时间，x为横坐标、y为纵坐标，D为石墨烯纳米质量传感器弯曲刚度，N为选定的形函数矩阵N(y)，T为矩阵转置符号，ν为泊松比，ρ为石墨烯的密度，h为石墨烯纳米质量传感器的厚度，l为条形单元的宽度，μ为非局部参数。
5.根据权利要求4所述的石墨烯纳米质量传感器振动特性分析的半解析方法，其特征在于，所述选定的形函数矩阵N(y)为标准Euler梁单元的形函数。
6.根据权利要求5所述的石墨烯纳米质量传感器振动特性分析的半解析方法，其特征在于，所述步骤4)的详细步骤包括：
4.1)针对各个子结构，采用有限元法将各个条形单元的刚度矩阵和分别组装成总体刚度矩阵K(4)、K(2)和K(0)，将各个条形单元的质量矩阵和分别组装成总体质量矩阵M(2)和M(0)；
4.2)确定各个子结构的总体位移向量如式(5)所示；
Φi(x,t)＝{w1(x,t),θ1(x,t),w2(x,t),θ2(x,t),…,wNE(x,t),θNE(x,t)}T,i＝1,2 (5)
式(5)中，Φi(x,t)为第i个子结构的总体位移向量，w1(x,t)为第i个子结构中第1个条形单元的结线横向位移，θ1(x,t)为第i个子结构中第一个条形单元的结线转角，w2(x,t)为第i个子结构中第2个条形单元的结线横向位移，θ2(x,t)为第i个子结构中第2个条形单元的结线转角，wNE(x,t)为第i个子结构中第NE个条形单元的结线横向位移，θNE(x,t)为第i个子结构中第NE个条形单元的结线转角；
4.3)根据各个子结构的总体位移向量组装总体运动微分方程如式(6)所示；
<mrow> <msup> <mi>K</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </msup> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>4</mn> </msup> <msub> <mi>&amp;Phi;</mi> <mi>i</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>x</mi> <mn>4</mn> </msup> </mrow> </mfrac> <mo>+</mo> <msup> <mi>K</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msup> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>&amp;Phi;</mi> <mi>i</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <msup> <mi>K</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msup> <msub> <mi>&amp;Phi;</mi> <mi>i</mi> </msub> <mo>+</mo> <msup> <mi>M</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msup> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>&amp;Phi;</mi> <mi>i</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>t</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <msup> <mi>M</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msup> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>4</mn> </msup> <msub> <mi>&amp;Phi;</mi> <mi>i</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> <mo>&amp;part;</mo> <msup> <mi>t</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
式(6)中，Φi为第i个子结构的总体位移向量，K(4)、K(2)和K(0)为总体刚度矩阵，M(2)和M(0)为总体质量矩阵，t为时间，x为横坐标；
4.4)针对所述总体运动微分方程处理边界条件，得到总体运动微分方程如式(7)所示；
<mrow> <msup> <mover> <mi>K</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </msup> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>4</mn> </msup> <msub> <mover> <mi>&amp;Phi;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>x</mi> <mn>4</mn> </msup> </mrow> </mfrac> <mo>+</mo> <msup> <mover> <mi>K</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msup> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mover> <mi>&amp;Phi;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <msup> <mover> <mi>K</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msup> <msub> <mover> <mi>&amp;Phi;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>+</mo> <msup> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msup> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mover> <mi>&amp;Phi;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>t</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <msup> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msup> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>4</mn> </msup> <msub> <mover> <mi>&amp;Phi;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> <mo>&amp;part;</mo> <msup> <mi>t</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> 2
式(7)中，为第i个子结构处理边界条件后的总体位移向量，和为处理边界条件后的总体刚度矩阵，和为处理边界条件后的总体质量矩阵。
7.根据权利要求6所述的石墨烯纳米质量传感器振动特性分析的半解析方法，其特征在于，所述步骤5)的详细步骤包括：
5.1)将各个子结构处理边界条件后的总体运动微分方程取时间的Laplace变换，得到各个子结构的Laplace总体运动微分方程如式(8)所示；
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mfrac> <msup> <mo>&amp;part;</mo> <mn>4</mn> </msup> <mrow> <mo>&amp;part;</mo> <msup> <mi>x</mi> <mn>4</mn> </msup> </mrow> </mfrac> <msub> <mover> <mi>&amp;Phi;</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>D</mi> <mn>2</mn> </msub> <mfrac> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mrow> <mo>&amp;part;</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <mi>D</mi> <mo>(</mo> <mi>s</mi> <mo>)</mo> <mo>)</mo> </mrow> <msub> <mover> <mi>&amp;Phi;</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <mi>A</mi> <mrow> <mo>(</mo> <msup> <mover> <mi>K</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msup> <mo>+</mo> <msup> <mi>s</mi> <mn>2</mn> </msup> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>D</mi> <mn>2</mn> </msub> <mo>=</mo> <mo>-</mo> <mi>A</mi> <msup> <mover> <mi>K</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>A</mi> <mo>=</mo> <msup> <mrow> <mo>(</mo> <msup> <mover> <mi>K</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </msup> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
式(8)中，为第i个子结构处理边界条件后的总体位移向量的Laplace变换，‘^’为Laplace变换符号，s为Laplace变换系数，和为处理边界条件后的总体刚度矩阵，为处理边界条件后的总体刚度矩阵，D(s)、D2、A均为中间矩阵；
5.2)基于定义的状态向量η(x,s)，将Laplace总体运动微分方程转换为式(9)所示：
<mrow> <mfrac> <mo>&amp;part;</mo> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mi>&amp;eta;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>F</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mi>&amp;eta;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
式(9)中，η(x,s)为定义的状态向量，且定义的状态向量η(x,s)的函数表达式如式(10)所示，F(s)为传递矩阵，传递矩阵F(s)的函数表达式如式(11)所示；
式(10)中，η(x,s)为定义的状态向量，η1(x,s)和η2(x,s)均为中间向量，为第1个子结构处理边界条件后的总体位移向量的Laplace变换，为第2个子结构处理边界条件后的总体位移向量的Laplace变换；
<mrow> <mi>F</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mover> <mi>F</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <msub> <mn>0</mn> <mrow> <mn>4</mn> <msub> <mi>N</mi> <mn>1</mn> </msub> <mo>&amp;times;</mo> <mn>4</mn> <msub> <mi>N</mi> <mn>1</mn> </msub> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mn>0</mn> <mrow> <mn>4</mn> <msub> <mi>N</mi> <mn>1</mn> </msub> <mo>&amp;times;</mo> <mn>4</mn> <msub> <mi>N</mi> <mn>1</mn> </msub> </mrow> </msub> </mtd> <mtd> <mrow> <mover> <mi>F</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mover> <mi>F</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mi>I</mi> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>I</mi> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>I</mi> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>D</mi> <mn>2</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
式(11)中，F(s)为传递矩阵，为中间矩阵，为4×N1阶0元素方阵，N1为每个子结构的总体自由度数，D(s)、D2为式(8)中所示的中间矩阵，I为单位矩阵；
5.3)建立两个子结构连接处的边界条件如式(12)所示；
Mb(s)η(-0.5a,s)+Nb(s)η(0.5a,s)+Rb(s)η(x0,s)＝0 (12)
式(12)中，Mb(s)、Nb(s)和Rb(s)为边界矩阵，η为状态向量，a为石墨烯质量传感器在x方向上的长度，s为Laplace变换系数，x0为附加质量的质心的横坐标；其中边界矩阵Rb(s)如式(13)所示；
<mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mn>0</mn> <mrow> <mn>4</mn> <msub> <mi>N</mi> <mn>1</mn> </msub> <mo>&amp;times;</mo> <mn>4</mn> <msub> <mi>N</mi> <mn>1</mn> </msub> </mrow> </msub> </mtd> <mtd> <msub> <mn>0</mn> <mrow> <mn>4</mn> <msub> <mi>N</mi> <mn>1</mn> </msub> <mo>&amp;times;</mo> <mn>4</mn> <msub> <mi>N</mi> <mn>1</mn> </msub> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>r</mi> <mrow> <mi>b</mi> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>r</mi> <mrow> <mi>b</mi> <mn>2</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <msub> <mi>r</mi> <mrow> <mi>b</mi> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>I</mi> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mi>I</mi> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>I</mi> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>0</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msup> <mover> <mi>K</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <msub> <mi>r</mi> <mrow> <mi>b</mi> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <mi>I</mi> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>I</mi> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>I</mi> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <msup> <mover> <mi>K</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </msup> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
式(13)中，Rb(s)为边界矩阵，为4×N1阶0元素方阵，N1为每个子结构的总体自由度数，rb1和rb2均为中间矩阵，m0为附加质量的大小，I为单位矩阵，为处理边界条件后的总体刚度矩阵；
5.4)生成将式(9)所示转换后的Laplace总体运动微分方程、式(12)所示两个子结构连接处的边界条件的解如式(14)所示，且所述解的特征方程如式(15)所示；
<mrow> <mo>&amp;lsqb;</mo> <msub> <mi>M</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mn>0.5</mn> <mi>a</mi> <mi>F</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <msub> <mi>N</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mn>0.5</mn> <mi>a</mi> <mi>F</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <msub> <mi>x</mi> <mn>0</mn> </msub> <mi>F</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> </mrow> </msup> <mo>&amp;rsqb;</mo> <mi>&amp;eta;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
式(14)中，Mb(s)、Nb(s)和Rb(s)为边界矩阵，a为石墨烯质量传感器在x方向上的长度，F(s)为传递矩阵，x0为附加质量的质心的横坐标，x为横坐标，s为Laplace变换系数；
<mrow> <mi>det</mi> <mrow> <mo>(</mo> <msub> <mi>M</mi> <mi>b</mi> </msub> <mo>(</mo> <mi>s</mi> <mo>)</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mn>0.5</mn> <mi>a</mi> <mi>F</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <msub> <mi>N</mi> <mi>b</mi> </msub> <mo>(</mo> <mi>s</mi> <mo>)</mo> <msup> <mi>e</mi> <mrow> <mn>0.5</mn> <mi>a</mi> <mi>F</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <msub> <mi>x</mi> <mn>0</mn> </msub> <mi>F</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> </mrow> </msup> <mo>)</mo> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
式(15)中，det为求行列式符号，Mb(s)、Nb(s)和Rb(s)为边界矩阵，a为石墨烯质量传感器在x方向上的长度，F(s)为传递矩阵，x0为附加质量的质心的横坐标，x为横坐标，s为Laplace变换系数；
5.5)针对式(15)所示解的特征方程，令Laplace变换系数s为虚数单位i、石墨烯纳米质量传感器自由振动的固有频率ω两者之间的乘积，求解得到石墨烯纳米质量传感器自由振动的固有频率ω并作为石墨烯纳米质量传感器的振动特性分析结果输出。
CN201610246796.8A 2016-04-20 2016-04-20 一种石墨烯纳米质量传感器振动特性分析的半解析方法 Active CN107315895B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610246796.8A CN107315895B (zh) 2016-04-20 2016-04-20 一种石墨烯纳米质量传感器振动特性分析的半解析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610246796.8A CN107315895B (zh) 2016-04-20 2016-04-20 一种石墨烯纳米质量传感器振动特性分析的半解析方法

Publications (2)

Publication Number Publication Date
CN107315895A true CN107315895A (zh) 2017-11-03
CN107315895B CN107315895B (zh) 2019-12-20

Family

ID=60184231

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610246796.8A Active CN107315895B (zh) 2016-04-20 2016-04-20 一种石墨烯纳米质量传感器振动特性分析的半解析方法

Country Status (1)

Country Link
CN (1) CN107315895B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102004823A (zh) * 2010-11-11 2011-04-06 浙江中科电声研发中心 一种扬声器振动和声学特性的数值模拟方法
CN104897508A (zh) * 2014-03-06 2015-09-09 中国科学院上海微系统与信息技术研究所 一种功能材料热力学参数的测试方法
CN105426343A (zh) * 2015-11-02 2016-03-23 株洲时代新材料科技股份有限公司 一种基于傅里叶级数的复杂结构振动解析分析方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102004823A (zh) * 2010-11-11 2011-04-06 浙江中科电声研发中心 一种扬声器振动和声学特性的数值模拟方法
CN104897508A (zh) * 2014-03-06 2015-09-09 中国科学院上海微系统与信息技术研究所 一种功能材料热力学参数的测试方法
CN105426343A (zh) * 2015-11-02 2016-03-23 株洲时代新材料科技股份有限公司 一种基于傅里叶级数的复杂结构振动解析分析方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ZHI-BIN SHEN等: "Transverse vibration of nanotube-based micro-mass sensor via nonlocal Timoshenko beam theory", 《COMPUTATIONAL MATERIALS SCIENCE》 *
唐海利: "基于非局部弹性理论的碳纳米管质量传感器振动特性研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 *
申志彬: "基于非局部弹性理论的微纳米质量传感器振动特性研究", 《中国博士学位论文全文数据库 信息科技辑》 *

Also Published As

Publication number Publication date
CN107315895B (zh) 2019-12-20

Similar Documents

Publication Publication Date Title
Lin A mesh-based partition of unity method for discontinuity modeling
Olshevskiy et al. Three-dimensional solid brick element using slopes in the absolute nodal coordinate formulation
Carrera et al. Accuracy of refined finite elements for laminated plate analysis
Jeon et al. The MITC3+ shell element in geometric nonlinear analysis
Bediz A spectral-Tchebychev solution technique for determining vibrational behavior of thick plates having arbitrary geometry
Liu et al. A point interpolation mesh free method for static and frequency analysis of two-dimensional piezoelectric structures
Hassani et al. A review of homogenization and topology optimization III—topology optimization using optimality criteria
Schwab et al. Comparison of three-dimensional flexible beam elements for dynamic analysis: classical finite element formulation and absolute nodal coordinate formulation
CN108986220B (zh) 一种加速有限元求解物体网格模型弹性变形的方法
Vieira et al. A higher order thin-walled beam model including warping and shear modes
CN101944144A (zh) 一种基于无网格的布类仿真方法
Lu et al. A computationally efficient modeling method for the vibration analyses of two-dimensional system structures using reduced transfer matrix method for multibody system
Rezaiee-Pajand et al. Mixing dynamic relaxation method with load factor and displacement increments
Ruiz et al. Eigenvector sensitivity when tracking modes with repeated eigenvalues
Jain et al. Effect of self-weight on topological optimization of static loading structures
Dyniewicz Space–time finite element approach to general description of a moving inertial load
Wang Optimization of support positions to minimize the maximal deflection of structures
Yarasca et al. Best Theory Diagrams for cross-ply composite plates using polynomial, trigonometric and exponential thickness expansions
Eugster Numerical analysis of nonlinear wave propagation in a pantographic sheet
Puckett et al. Compound strip method for analysis of plate systems
Singh et al. Nonlinear free vibration analysis of composite plates with material uncertainties: A Monte Carlo simulation approach
Dai et al. The global nonlinear galerkin method for the analysis of elastic large deflections of plates under combined loads: A scalar homotopy method for the direct solution of nonlinear algebraic equations
CN107315895A (zh) 一种石墨烯纳米质量传感器振动特性分析的半解析方法
CN106909699A (zh) 基于Galerkin条形传递函数的薄板振动特性分析方法
CN107092730A (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