en
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。
作者简介:

周博(1972-),男,教授,博士,研究方向为智能材料与结构力学、微尺度力学、石油工程力学等。E-mail:zhoubo@upc.edu.cn。

通讯作者:

林英松(1964-),女,教授,博士,研究方向为油气井工程、石油工程力学等。E-mail:linyingsong@upc.edu.cn。

中图分类号:O341

文献标识码:A

文章编号:1673-5005(2021)01-0151-07

DOI:10.3969/j.issn.1673-5005.2021.01.018

目录contents

    摘要

    基于修正偶应力理论推导任意截面形状 Bernoulli-Euler 微梁的变形能、弯曲刚度、外力功、动能等基本变量的偶应力理论表达式,进而通过 Hamilton 原理建立 Bernoulli-Euler 微梁的偶应力理论动力微分方程。 根据偶应力理论弯曲刚度表达式,研究任意截面形状 Bernoulli-Euler 微梁弯曲刚度的尺寸效应,分析泊松系数和截面形状对 Bernoulli-Euler 微梁弯曲刚度及其尺寸效应的影响。 基于偶应力理论动力微分方程,求解 Bernoulli-Euler 简支微梁的偶应力理论固有频率,据此研究任意面形状 Bernoulli-Euler 微梁固有频率的尺寸效应,分析泊松系数和截面形状对 Bernoulli-Euler 微梁固有频率及其尺寸效应的影响。 结果表明,建立的由基本变量偶应力理论表达式和偶应力理论动力微分方程构成的动力学模型能有效地描述任意截面形状 Bernoulli-Euler 微梁动力学特性的尺寸效应。

    Abstract

    Based on the modified couple stress theory, the couple stress theoretical expressions of the basic variables such as deformation energy, bending stiffness, external force work and kinetic energy of Bernoulli-Euler microbeam with arbitrary cross-section shape were derived. Then the couple stress theoretical dynamic differential equation of Bernoulli-Euler microbeam was established according to Hamilton principle. The scale effect of bending stiffness of Bernoulli-Euler microbeam with arbitrary section shape was assessed using the couple stress theoretical expression of the bending stiffness. And the influences of Poisson̍s coefficient and section shape on the bending stiffness and its scale effect of Bernoulli-Euler microbeam were analyzed. Based on the couple stress theoretical dynamic differential equation, the natural frequency of Bernoulli-Euler simply supported beam was obtained. The scale effect of the natural frequency of Bernoulli-Euler beams with arbitrary crosssection shape was assessed. And the influences of Poisson’s coefficient and cross-section shape on the natural frequency and its scale effect of Bernoulli-Euler beams were analyzed. The results show that the developed dynamic model, which consists of the couple stress theoretical expressions of the aforementioned basic variables and the couple stress theoretical dynamic differential equation, can effectively describe the scale effect of the dynamic characteristics of the Bernoulli-Euler microbeam.

  • 随着国际上微机电系统的快速发展,对微结构的力学性能分析成为了国际研究的热点。微梁在微机电系统中被广泛应用,比如微传感器、微驱动器和微开关。当金属或复合材料构件的尺寸减小到微纳米范围时,其力学性能随着几何尺寸与特征尺寸的改变而表现出较强的尺寸效应[1-2]。经典弹性理论将物体抽象为均匀连续介质模型,忽略了材料内部微结构对其力学性能的影响,无法描述微结构的变形特点及尺寸效应。 Toupin [3]提出的修正偶应力理论,在弹性体的本构关系和应变能密度中引入材料尺度参数,实现了对微结构尺寸效应的描述。 Park等[4]基于修正偶应力理论,通过材料尺度参数捕捉边界层效应,解析求解了一个剪切问题。 Armagan等[5]基于修正偶应力理论对功能梯度夹层微梁的弯曲、屈曲和自由振动进行了分析。 Nguyen等[6] 考虑尺寸效应和速率效应,建立了应变梯度黏塑性本构模型,描述了微米尺度金属的时效行为。贺丹等[7]考虑弯曲刚度的尺寸效应,基于修正偶应力理论建立了适用于各种边界条件的微尺度Euler梁模型,用于求解微尺度梁的弯曲和稳定问题。周博等[8]研究了Timoshenko微梁的挠度、转角和应力的尺寸效应,分析了Poisson比对Timoshenko微梁力学行为及其尺寸效应的影响。郑雪瑶等[9] 基于修正偶应力理论推导了梁的变形能计算公式,得到了反映尺寸效应的量纲归一化刚度公式并利用最小势能原理研究了Euler-Bernoulli微梁,得到了反映尺寸效应的量纲归一化挠曲线方程和量纲归一化最大应力方程。康泽天等[10]建立了微板的动力学模型, 研究了正交各向异性功能梯度四边简支微板的自由振动行为。 Zhou等[11]建立了黏弹性功能梯度(FG) Timoshenko微梁的尺寸依赖连续模型。 Ke等[12-13] 求解了功能梯度板的动力学问题。 Karimipour等[14] 研究了微环面板的尺寸依赖行为。 Liang等[15]提出一种考虑泊松效应和弱非局部应变梯度弹性效应的Bernoulli-Euler梁模型。 Attia等[16] 研究了小尺寸对自由振动的影响。 Dai等[17]提出一种非线性模型描述悬臂微梁的共振特性,并对其线性和非线性弹簧结构的定向动态响应问题进行研究。张大千等[18] 对微板进行了热稳定分析。刘建林等[19]探究了微纳米尺度下线弹性悬臂梁受集中载荷作用下的大变形问题。这些关于微梁力学特性尺寸效应方面的研究,主要是针对特定截面形状微梁结构进行的,且大多忽略了泊松系数对微梁力学行为及其尺寸效应的影响。笔者基于修正偶应力理论,推导任意截面形状Bernoulli-Euler微梁的弯曲刚度、变形能等基本变量的偶应力理论表达式,利用Hamilton变分原理建立任意截面形状Bernoulli-Euler微梁的偶应力理论动力微分方程。研究任意截面形状Bernoulli-Euler微梁弯曲刚度的尺寸效应。

  • 1 修正偶应力理论

  • Yang等[5]提出了修正偶应力理论。物体内的任一物质点有6 个自由度,包括3 个平动自由度和3 个转动自由度。 3 个平动自由度用位移矢量ui表示,3 个转动自由度用转动矢量 θi表示。转动矢量和位移矢量的关系为

  • θi=12eijkuk,j
    (1)
  • 式中,eijk 为置换张量。

  • 描述物体自由度和变形间关系的几何方程表示为

  • εij=12ui,j+uj,i
    (2)
  • χij=12θi,j+θj,i
    (3)
  • 式中,εijχij分别为应变张量和曲率张量。

  • 在小变形情况下各向同性弹性材料的本构方程为

  • σij=λεkkδij+2Gεij
    (4)
  • mij=2l2Gij
    (5)
  • 其中

  • λ=Eμ(1+μ)(1-2μ),G=E2(1+μ).

  • 式中,σijmij 分别为应力张量和偏斜偶应力张量; λG 为拉梅系数;l 为反映尺寸效应的材料参数, 称为材料特征尺寸;Eμ 分别为弹性模量和泊松系数。

  • 弹性体的变形能可表示为

  • U=12Ω σijεij+mijXijdV
    (6)
  • 式中,Ω 为弹性体体域;dV 为体积微元。

  • 2 动力学模型建立

  • 2.1 基本变量描述

  • Bernoulli-Euler梁理论认为梁截面在变形后仍保持为平面,且垂直于变形后的梁轴线,只是绕截面内的某一轴线旋转了一个角度。在计算过程中略去剪切变形,即假设只有弯曲变形,导出的应力和变形的计算公式对于细长梁是正确的。

  • 计算公式对于细长梁是正确的。图1 为分布载荷作用下的Bernoulli-Euler梁,x轴与梁的轴线重合,梁的截面中性轴与y轴重合,分布载荷集度为 q(x,t)。

  • 图1 横向载荷作用下的Bernoulli-Euler梁

  • Fig.1 Bernoulli-Euler beam under transverse load

  • Bernoulli-Euler梁位移场可表示为

  • u=-zw(x,t)x,v=0,w=w(x,t)
    (7)
  • 式中,u、v、w分别为x、y、z方向的位移分量。

  • 将式(7)代入式(1),得到转动分量

  • θx=0,θy=-wx,θz=0
    (8)
  • 根据式(7)和(2),将Bernoulli-Euler梁的非零应变分量表示为

  • εxx=-z2wx2εxx=-z2wx2
    (9)
  • 根据式(8)和(3),将Bernoulli-Euler梁的非零曲率分量表示为

  • χxy=χyx=-12w2x2
    (10)
  • 将式(8)代入式(4),得到Bernoulli-Euler梁的非零应力分量,即

  • σxx=-(λ+2G)z2wx2,σzz=σyy=-λz2wx2
    (11)
  • 将式(10) 代入式(5),得到Bernoulli-Euler梁的非零偏斜偶应力分量,即

  • mxy=myx=-Gl22wx2
    (12)
  • 将式(9)~(12)代入式(6),通过积分运算,得到Bernoulli-Euler梁的变形能,即

  • U=12L K2wx22dx
    (13)
  • 其中

  • K=(λ+2G)Iy+GAl2,Iy=A z2dA

  • 式中,K 为Bernoulli-Euler梁的偶应力理论弯曲刚度;A 为梁截面面积; Iy 为梁截面对其中性轴y轴的惯性矩。

  • Bernoulli-Euler梁的外载荷做功表示为

  • W=L q(x,t)w(x,t)dx
    (14)
  • Bernoulli-Euler梁的动能可近似表示为

  • T=12L ρAwt2dx
    (15)
  • 式中,ρ 为梁密度。

  • 2.2 动力微分方程

  • 利用Hamilton原理推导Bernoulli-Euler微梁动力微分方程。 Hamilton原理可描述为具有理想约束的完整系统在任一时间间隔[t 0,t 1 ]内,从初始位置到终止位置的所有可能运动中,真实运动的哈密顿作用量的变分等于零,它是一种求解动力学问题的变分原理。对于Bernoulli-Euler梁,Hamilton原理表示为

  • t0t1 (δT-δU+δW)dt=0
    (16)
  • 式中,δTδU δW 分别为Bernoulli-Euler梁的动能变分、变形能变分和外力功的变分。

  • 对Bernoulli-Euler梁的变形能方程式(13)进行变分运算,得到

  • δU=L K4wx4δwdx
    (17)
  • 对外载荷做功方程式(14)进行变分运算,得到

  • δW=L q(x,t)δwdx
    (18)
  • 对Bernoulli-Euler梁的动能方程式(15)进行变分运算,得到

  • δT=-L ρA2wt2δwdx
    (19)
  • 根据式(17)~(19),得到

  • δT-δU+δW=L -ρA2wt2-K4wx4+q(x,t)δwdx
    (20)
  • 由式(20)可知,若使Hamilton原理即式(16)对任一时间间隔[t 0,t 1 ]都成立,应满足

  • ρA2wt2+K4wx4-q(x,t)=0
    (21)
  • 式(21)为描述Bernoulli-Euler微梁的挠度、载荷集度间微分关系的动力微分方程。

  • 3 尺寸效应

  • 3.1 弯曲刚度比

  • 在材料特征尺寸 l=0 情况下,根据式(13)中描述的Bernoulli-Euler梁的偶应力理论弯曲刚度 K, 得到Bernoulli-Euler梁的经典弹性理论弯曲刚度,表示为

  • K0=(λ+2G)Iy
    (22)
  • 基于式(22) 将偶应力理论弯曲刚度和经典弹性理论弯曲刚度的比值表示为

  • KK0=1+1-2μ2-2μlry2
    (23)
  • 其中

  • ry=Iy/A

  • 式中,ry 为截面对中性轴惯性半径;l 为材料特征尺寸。

  • 惯性半径 ry 是Bernoulli-Euler梁的重要几何特征尺寸,在截面面积相同的情况下,截面惯性半径越大,其惯性矩越大,弯曲刚度越大。

  • 3.2 固有频率比

  • 令载荷集度函数 q( x, t)=0 并将其代入式(21),得到Bernoulli-Euler微梁的自由振动方程,表示为

  • ρA2wt2+K4wx4=0
    (24)
  • 根据分离变量法,设自由振动方程式(24) 的解为

  • w(x,t)=[Bcos(ωt)+Csin(ωt)]Φ(x)
    (25)
  • 式中,ω 为固有频率;Φ(x)为振型函数;BC 为待定系数。

  • 将式(25)代入式(24),经简化整理得

  • 4Φx4-β4Φ=0
    (26)
  • 其中

  • β4=ρAω2K.

  • 微分方程式(26)的通解为

  • Φ(x)=C1sin(βx)+C2cos(βx)+C3sinh(βx)+C4cosh(βx)
    (27)
  • 式中,C1C2C3C4 为待定系数。

  • 以简支梁为例,其左右两端挠度为零,弯矩为零,边界条件可描述为

  • Φx=0=Φx=L=0,d2Φdx2x=0=d2Φdx2x=L=0.
    (28)
  • 式中,L 为简支梁长度;Φ 为简支梁的振型函数。

  • 根据式(27)和(28),可得

  • C2=C3=C4=0,C1sin(βL)=0
    (29)
  • 根据式(27)和(29)可知,若 Φ(x)存在非零解, 必有 C1≠0,因此

  • sin(βL)=0
    (30)
  • 求解式(30),得

  • β=iπL,i=1,2,3,
    (31)
  • 根据式(26)和式(31),得简支梁的偶应力理论固有频率

  • ω=(iπ)2KρAL4,i=1,2,3,
    (32)
  • 令材料特征尺寸l=0,根据式( 32)、( 22) 和(23),可以得到Bernoulli-Euler简支梁的经典弹性理论固有频率,即

  • ω0=(iπ)2K0ρAL4,i=1,2,3,
    (33)
  • 根据式(33)、(22)和(23),可以将简支梁的偶应力理论固有频率和经典弹性理论固有频率的比值表示为

  • ωω0=1+1-2μ2-2μlry2
    (34)
  • 3.3 尺寸效应分析

  • 图2 为不同泊松系数时,Bernoulli-Euler梁的偶应力理论弯曲刚度与截面惯性半径间的关系。其纵坐标为Bernoulli-Euler梁的偶应力理论弯曲刚度和经典弹性理论弯曲刚度的比值,横坐标为截面惯性半径与材料特征尺寸的比值。由图2 可知:当惯性半径与材料特征尺寸比值较小时,随着惯性半径与材料特征尺寸比值增加,偶应力理论弯曲刚度与经典弹性理论弯曲刚度比值逐渐减小,Bernoulli-Euler微梁的弯曲刚度特性具有明显的尺寸效应;由图2 可以看出,当惯性半径与材料特征尺寸比值小于2 时,随着横坐标增加,Bernoulli-Euler梁的偶应力理论弯曲刚度与经典弹性理论弯曲刚度极速下降,之后缓慢下降,并最后趋于平缓,尺寸效应变得不明显、可以忽略;Bernoulli-Euler梁的泊松系数越小,惯性半径对弯曲刚度的影响越大、弯曲刚度特性的尺寸效应越明显。

  • 图2 不同泊松系数下弯曲刚度与惯性半径间关系

  • Fig.2 Relation between bending stiffness and inertial radius with different Poisson̍s ratio

  • 图3 为不同惯性半径与材料特征尺寸比值时, Bernoulli-Euler微梁的偶应力理论弯曲刚度和泊松系数间的关系曲线,其纵坐标为偶应力理论弯曲刚度和经典弹性理论弯曲刚度比值。由图3 可知: Bernoulli-Euler梁的偶应力理论弯曲刚度与经典弹性理论弯曲刚度比值随着泊松系数增大而减小;惯性半径与材料特征尺寸比值越大,泊松系数对弯曲刚度及其尺寸效应的影响越小。

  • 图3 不同惯性半径下弯曲刚度与泊松系数间关系

  • Fig.3 Relation between bending stiffness and Poisson̍s ratio with different inertial radius

  • 图4 为泊松系数等于0.45、截面面积相同,截面形状分别为矩形(高宽比2)、箱形(外高宽比2, 内外高比0.8,内外宽比0.8)、圆形、圆环(内外径比0.8)时,Bernoulli-Euler微梁的偶应力理论弯曲刚度和截面面积间的关系曲线,其纵坐标为偶应力理论弯曲刚度和经典弹性理论弯曲刚度比值,横坐标为截面面积的平方根与材料特征尺寸比值。由图4 可知:Bernoulli-Euler微梁的偶应力理论弯曲刚度与经典弹性理论弯曲刚度的比值随着截面面积增大而减小且逐渐趋近于1;截面形状对Bernoulli-Euler微梁的弯曲刚度及其尺寸效应有较大影响。

  • 图4 不同截面形状下弯曲刚度与截面面积间关系

  • Fig.4 Relational curves between bending stiffness and cross section area with different section shape

  • 图5 为内外径比为0.5 的圆环截面BernoulliEuler微梁,在不同泊松系数时,偶应力理论弯曲刚度与截面面积间的关系曲线,其纵坐标为偶应力理论弯曲刚度与经典弹性理论弯曲刚度比值,横坐标为横截面面积的平方根与材料特征尺寸比值。由图5 可知:Bernoulli-Euler梁的偶应力理论弯曲刚度与经典弹性理论弯曲刚度比值随着截面面积增大而减小且逐渐趋近于1;Bernoulli-Euler微梁的泊松系数越小,截面面积对偶应力理论弯曲刚度的影响越大、尺寸效应越明显。

  • 图5 不同泊松系数下弯曲刚度与截面面积关系

  • Fig.5 Relation between bending stiffness and cross section area with different Poisson̍s ratio

  • 图6 为不同泊松系数时,Bernoulli-Euler微梁的固有频率和截面惯性半径间的关系曲线,其纵坐标为Bernoulli-Euler梁的偶应力理论固有频率和经典弹性理论固有频率比值,横坐标为截面惯性半径和材料特征尺寸比值。由图6 可知:当惯性半径与材料特征尺寸比值较小时,随着惯性半径与材料特征尺寸比值增加,偶应力理论固有频率与经典弹性理论固有频率比值逐渐减小,Bernoulli-Euler梁的固有频率特性具有明显的尺寸效应;由图6 可以看出,当惯性半径与材料特征尺寸比值小于2 时,随着截面惯性半径和材料特征尺寸比值增加,BernoulliEuler梁的偶应力理论固有频率与经典弹性理论固有频率比值极速下降,之后缓慢下降,最后趋于平缓,尺寸效应变得不明显,可以忽略;Bernoulli-Euler梁的泊松系数越小,惯性半径对固有频率的影响越大,固有频率特性的尺寸效应越明显。

  • 图6 不同泊松系数下固有频率与惯性半径关系

  • Fig.6 Relation between natural frequency and inertial radius with different Poisson̍s ratio

  • 图7 为不同惯性半径与材料特征尺寸比值时, Bernoulli-Euler微梁的偶应力理论固有频率和泊松系数间的关系曲线,其纵坐标为偶应力理论固有频率和经典弹性理论固有频率比值。由图7 可知: Bernoulli-Euler梁的偶应力理论固有频率与经典弹性理论固有频率的比值,随着泊松系数的增加而减小;惯性半径与材料特征尺寸比值越大,泊松系数对固有频率的影响越小。

  • 图7 不同惯性半径下固有频率与泊松系数关系

  • Fig.7 Relation between natural frequency and Poisson̍s ratio with different inertial radius

  • 图8 为泊松系数等于0.2、截面面积相同,截面形状分别为矩形(高宽比2)、箱形(外高宽比2,内外高比0.8,内外宽比0.8)、圆形、圆环(内外径比0.8)时,Bernoulli-Euler微梁偶应力理论固有频率和截面面积间的关系曲线,其纵坐标为偶应力理论固有频率和经典弹性理论固有频率的比值,横坐标为截面面积的平方根与材料特征尺寸比值。由图8 可知:Bernoulli-Euler微梁的偶应力理论固有频率与经典弹性理论固有频率比值随着截面面积增大而减小且逐渐趋近于1;截面形状对Bernoulli-Euler微梁偶应力理论固有频率及其尺寸效应有较大影响。

  • 图8 不同截面形状下固有频率和截面面积间关系

  • Fig.8 Relation between natural frequency and cross section area with different section shape

  • 图9 为圆形截面Bernoulli-Euler微梁在不同泊松系数时偶应力理论固有频率与截面面积间的关系曲线,其纵坐标为偶应力理论固有频率与经典弹性理论固有频率比值,横坐标为截面面积的平方根与材料特征尺寸比值。由图9 可知:Bernoulli-Euler梁的偶应力理论固有频率与经典弹性理论固有频率比值随着截面面积增大而减小且逐渐趋于1;Bernoulli-Euler微梁的泊松系数越小,其截面面积对偶应力理论固有频率的影响越大,尺寸效应越明显。

  • 图9 不同泊松系数下固有频率和截面面积关系

  • Fig.9 Relation between natural frequency and cross section area with different Poisson̍s ratio

  • 4 结论

  • (1)当Bernoulli-Euler微梁截面惯性半径与材料特征尺寸比值较小时,Bernoulli-Euler微梁的偶应力理论弯曲刚度与经典弹性理论弯曲刚度的比值,随着惯性半径增加而减小,表现出明显的尺寸效应;当惯性半径与材料特征尺寸比值小于2 时, Bernoulli-Euler梁的偶应力理论弯曲刚度与经典弹性理论弯曲刚度呈现极速下降趋势,之后缓慢下降,并最后趋于平缓,尺寸效应变得不明显并可以忽略。

  • (2)泊松系数是影响Bernoulli-Euler微梁动力学特性及尺寸效应的重要因素,偶应力理论弯曲刚度和偶应力理论固有频率均随着泊松系数增大而减小;泊松系数越小,Bernoulli-Euler偶应力理论弯曲刚度和偶应力理论固有频率的尺寸效应越明显。

  • (3)截面形状对Bernoulli-Euler微梁的偶应力理论弯曲刚度、偶应力理论固有频率、弯曲刚度的尺寸效应和固有频率的尺寸效应均有明显影响。

  • 参考文献

  • 参考文献

  • 版权所有 中国石油大学学报(自然科学版)编辑部 Copyright©2008 All Rights Reserved
    主管单位:中华人民共和国教育部 主办单位:中国石油大学
    地址: 青岛市黄岛区长江西路66号中国石油大学期刊社 邮编:266580 电话:0532-86983496 E-mail: journal@upc.edu.cn
    本系统由:北京勤云科技发展有限公司设计