更全的杂志信息网

带边角裂纹悬臂Mindlin板的振动特性研究*

更新时间:2016-07-05

引言

裂纹广泛存在于工程结构中,如果出现裂纹,结构的阻尼和刚度将会发生改变,从而导致结构的振动特性发生变化.梁、板元件是工业中应用最广泛的结构,特别是在航天、航空领域,由于板出现裂纹经常导致重大安全事故的发生.因此,研究裂纹对板振动特性的影响有重要的理论意义和广泛的工程应用前景.

对于含边角穿透裂纹的振动特性研究,Yuan,Dickinson[1]应用区域分解法,并在内边界人为地附加虚拟弹簧,研究了带边角裂纹的四边简支矩形板的弯曲振动问题.Yuan和Young[2]也用了区域分解法研究了完全自由的带边角和内部裂纹环形板的振动问题.Liew[3]等人用同样的方法得到了各种边界条件下裂纹矩形薄板的固有频率.他们将裂纹板假设成可以用适当方程表达的子区域的组合体,并得到了裂纹薄板振动问题的特征值方程.Leissa,Huang[4,5]假设裂纹板的模态函数由两部分构成,一个是代数多项式,一个是附加的转角函数,并结合Ritz法研究了带V型裂纹或边角裂纹矩形的自由振动.Huang[6,7]用同样的方法研究了四边简支或悬臂的矩形Mindlin板和功能材料Reddy板.本课题组陈丽华、孙悦[8,9]等人研究了带边角裂纹悬臂矩形薄板的振动特性.在此基础上本文考虑剪切变形的影响对带边角裂纹悬臂矩形Mindlin板的振动特性进行研究.

本文用梁函数组合法得到理想完整板的模态函数,这样由完整板的模态函数和角函数组成的特殊模态函数能更好的描述板的振动特性并且有了更好的物理意义.最后,本文研究了裂纹长度、裂纹高度、裂纹角度对Mindlin悬臂板振动特性的影响,并利用Ansys有限元软件进行验证.

1 控制方程

基于一阶剪切变形理论,Mindlin板自由振动的动力学方程为:

(1)

(2)

(3)

矩形悬臂板的边界条件为:

固定端:y=0:w=0,φy=0,φx=0

(4a)

自由端: y=b:Mx=0, Mxy, Qx=0

对于处在成长期的青少年,冬季运动更有促进身体成长的作用,不过相比其他季节,冬季运动也更应做好“预备功课”。

(4b)

x=0,a:My=0, Mxy, Qy=0

(4c)

式中,w表示板的法向位移;J=h3/12;h为板的厚度;ρ为板的质量密度;w表示中性面上的横向位移;φxφy是中性面法线转过去的角度.这里:

活动引导式教学模式改变了传统的“灌输式”教学模式,这种教学模式能够很好地启发学生思维,突显以学生为中心的教学理念,对现阶段我国全方位进行的教学改革具有重要意义。在倡导素质教育、创新性人才培养大形势下,应该给以足够的重视,并加以推广。

(5a)

(5b)

(5c)

(5d)

式中,D为板的抗弯刚度,D=Eh3/12(1-v2);G为剪切模量,G=E/(2(1-v));κτ=π2/12为剪切修正系数;E,v分别为材料的弹性模量和Poisson比.

2 位移函数

裂纹板自由振动解的形式可设为:

每个实验重复3次,结果表示为x±s,利用Statistix 8.1软件进行数据统计分析,P<0.05表示差异显著。SigmaPlot12.5软件作图。

(6a)

φx(x,y,t)=Φx(x,y)φ(t)

(6b)

φy(x,y,t)=Φy(x,y)φ(t)

(6c)

W(x,y),φx(x,y)和φy(x,y)为裂纹板的振型函数.

以产学研的发展思路为指导思想,以建立优良的电子商务专业教学实训创新创业基地为基础,采用现代科学的教育教学方法为核心,为把学生培养成为社会的栋梁以及学校的电子商务专业的建设成为国家级特色重点专业而做出努力!

本文采用双向梁函数组合的级数形式来逼近完整板振动的真实振型.对于完整矩形Mindlin板的振型可以设为:

(7a)

(7b)

(7c)

其中Xi(x),Yi(y)分别为板x,y方向与边界条件相对应的梁位移函数,Φi(x),Φj(y)分别为与板x,y边界条件相应之梁转角函数,Aij,Bij,Cij为待定的振型系数.

下面根据梁函数组合法来求完整板的模态函数.

如图1, x方向(自由-自由)的各阶梁函数为:

X1=1

(8a)

(8b)

(i=3,4,5,…)

(8c)

式中,ki是与梁频率有关的系数,为梁函数系数,此处

1403 自制腹膜反折悬吊装置在腹侧肾肿瘤后腹腔镜下肾部分切除术中的应用 王 磊,储传敏,干思舜,叶剑青,曲发军,杨 炜,田毅君,潘秀武,杨启维,崔心刚

y方向(固定-自由)的各阶梁函数为:

(9)

(3) 裂纹角度对悬臂Mindlin板频率的影响

由文献 [10]可以得出转角函数:

(10a)

(10b)

Xi‴表示Xi(x)对x的三阶偏导,Yj‴表示Yj(y)对y的三阶偏导.

图1 矩形板的尺寸以及(r,θ)坐标 Fig. 1 Dimension and coordinate of a rectangular plate with a side-crack

将式(8a~8c)和(9)分别代入式(10a)和(10b)中便可得到转角的Φi(x),Φj(y)的振型函数,再代入到(7a~7c)便可以得到完整矩形Mindlin板的振型函数.

考虑到带边角裂纹矩形板的振动,只用式(7a~7c)来描述其振型函数显然是不合理的,因为不满足裂纹处位移和转角不连续的条件,因此需要在完整板的基础上附加上一个角函数来描述裂纹附近的性质.针对裂纹板,其振型函数由两部分构成,一部分是用梁函数组合法得到的理想完整板Mindlin板的振型,另一部分是利用裂纹的尖端奇异性理论来构造描述裂纹附近位移和转角的角函数,即运用半角三角函数的性质,在裂纹两边构造挠度和转角不连续性质的表达式.所以带裂纹矩形Mindlin板的振型函数可设为:

(11a)

(11b)

(11c)

式中ψkc(r,θ)(k=1,2,3)为描述裂纹的角函数.

附加描述裂纹性质的角函数可设为:

图表中给出了带边角裂纹Mindlin板的前三阶无量纲频率板的长宽比为b/a=2;厚度和宽度的比为:h/a=0.03;裂纹的角度α=15°或α=-15°.

的极小值可以求得板的振动频率.将裂纹板的模态函数(11a~11c)代入式(16a)和(16b)中,其中是待定系数,也相当于独立的广义坐标,于是,方程(17)可以简化为多元函数的极值条件:Umax-Tmax对待定系数中任意一个的偏导数均为零,即:

(12)

这里是待求的振型系数,r是板上任意一点到裂纹尖端的距离,θ是板上任意一点和裂纹尖端连线与裂纹的夹角(图1).由于角函数表达式中含有的独立变量rθ,需要与完整板中采用的直角坐标系统一,还要进行坐标变换.需要注意的是,随着坐标的变化,角坐标θ与直角坐标x,y的关系式在板的不同区域内是不同的,如图1所示,在坐标变换过程中需要将板分成四个区域.

(r,θ)坐标的原点选在裂纹尖端P处,且有-πθπ.θ在裂纹两边分别为-ππ,利用三角函数sinθ处不连续的性质来描述裂纹两边位移不连续的现象,而利用三角函数cosθ求一阶导数以后在处不连续性质来描述转角不连续的现象.如图1所示,在裂纹板的4个区域内,有相同的极半径r表达式:

w(x,y,t)=W(x,y)φ(t)

(13)

但是θ的表达式在各个区域的形式有所不同,具体表达式如下:

(14)

式中,x0=a-dcosα,y0=c-dsinα,如图1所示,α是裂纹顺时针方向与水平线的夹角,ab分别是矩形板的长和宽,cd分别是裂纹的高度和长度.

将式(13)中r的表达式和(14)中不同区域的θ表达式代入(12)中就得到了直角坐标下的角函数表达式.

电视面板占一台成品电视的七成左右成本,虽然我国面板产业已经逐步走上正轨,但进口面板依旧占据我国电视行业的很大占比,面板价格仍然会受到外来因素的一部分影响。2017年-2018年面板价格走势呈现先高后低趋势,2018第三季度65吋、55吋UHD面板同比下滑36.8%,18.7%。面板年初至8月一直保持下滑趋势,9月有所回涨。

在(12)式中gk(x,y)是满足各种几何边界条件的简单函数.对于y=0处固支的悬臂板其表达式可表示为:

gk(x,y)=y (k=1,2,3)

(15)

3 Ritz

本文研究的是带裂纹矩形Mindlin板,如图1所示.基于Mindlin板理论,用Ritz法可以求出裂纹板的固有频率.在Ritz法中,带裂纹矩形Mindlin板的最大势能(Umax)和最大动能(Tmax)分别为:

(16a)

(16b)

式中,ω是带裂纹矩形Mindlin板的固有频率;W, ΦxΦy分别是式(11a~11c)中的模态函数;这些函数下标里逗号后面的变量表示函数对相应变量的偏导.通过求能量方程:

∏=Vmax-Tmax

以高中语文课本中的文言文《陈情表》为例,在教学这篇课文时,应当细分目标。首先要让学生扎实掌握语言基本功,对课文要能够做到准确翻译。对古汉语素材的吸收是实现语言构建的重要机制,而且有助于使用优美的语言。其次,体会本文中的强烈感情,感受文中的语言美、形象美、亲情美。再次,开掘其中的传统文化教育要素,让学生从中获得深刻感染。值得注意的是,教学目标的分割有别于文本的分割,在教学中要注重文本的整体性,引导学生熟读精思,获得语文核心素养的全面提高。

(17)

(l=0,1,2,…,nn=1,2,3,…)

(18)

经过计算,便可得到3ij+3N(N+3)个关于 这些系数的线性代数方程式,令方程组的系数行列式为零,解这个方程可得到固有频率ω,将求出的各阶ω的值代入线性方程组中求出系数,便得到各阶频率相对应的模态函数.本文计算悬臂Mindlin板的固有频率和模态时分别取(i=4,j=8,N=3)得到一个142阶的矩阵.

淮北平原范围为淮河干流以北到沙颍河以南地区,地处我国南北气候过渡地带,属暖温带半湿润季风气候区,四季分明,气温和降雨自东南向西北递减。年内年际降水丰枯变化悬殊,6—8月雨水较多,能够满足一般树种的生长需要。年平均气温14~15℃,无霜期 200~220天,光热水等条件较好,适于农林业的综合发展。

4 结果分析

4.1 带裂纹悬臂Mindlin板的频率分析

(1)裂纹长度对悬臂Mindlin板频率的影响

表1考察相同裂纹角度和裂纹位置下不同裂纹长度对板固有频率的影响.裂纹位置c/b=1/2;裂纹的角度α=15°或α=-15°;裂纹长度的变化为d/a=0.1,0.2,0.3,…,0.6.

1 不同裂纹长度悬臂板的频率参数 Table 1 Frequency parameters ωa2 of the cantilevered square plates with different crack lengths

Modeαc/bd/a12315°0.503.44714.64121.1660.13.43614.50821.0070.23.40214.27120.6640.33.37413.89120.1850.43.32013.34019.5900.53.20412.68418.9340.63.12811.96418.283-15°0.50.13.43614.62920.1240.23.39914.41220.0540.33.38114.03519.6730.43.34913.51519.0170.53.30712.93118.2120.63.25012.31217.374

在裂纹角度和位置相同的情况下,改变裂纹的长度对固有频率有一定的影响,从表1可以看出,不论对于第几阶固有频率,均是裂纹长度越长,板固有频率就越小,这是由于局部刚度减小的原因造成的.当裂纹长度大于或等于板的宽度的一半时,固有频率减小得更快.

(2)裂纹位置对悬臂Mindlin板频率的影响

图2考查了相同裂纹角度和裂纹长度下不同裂纹位置对板的前三阶固有频率的影响,裂纹位置的变化为c/b=0.1,0.15,0.2,0.25,0.3,…,0.75.

在裂纹长度和角度相同的情况下,改变裂纹的位置对固有频率也有一定的影响,从图2可以得出:当α±15°时,随着裂纹位置由固定端向自由端移动,裂纹板的第一阶固有频率逐渐增大,第二阶固有频率先减小后增大,即当裂纹处在板中间位置时对固有频率的减小比在两端时有更显著的影响.在α为15°时,对于第三阶固有频率,由于第三阶振动模态有两条节线,所以固有频率的变化趋势是先增大后减小再增大,在越靠近节线位置附近频率减小的越显著.当α为-15°时,固有频率的趋势是先减小后增大.

图2 d/a=0.4时不同裂纹位置悬臂板的前三阶频率参数 Fig. 2 Firstthree order Frequency parameters ωa2 of the cantilevered square plate with different crack locations (d/a=0.4)

此处

图3(a)~3(c)考查了相同裂纹位置和裂纹长度下不同裂纹角度对板的前三阶固有频率的影响,横坐标的单位为(°).

在裂纹长度和位置相同的情况下,改变裂纹的角度对固有频率有一定的影响,当裂纹角度与x轴的夹角从正方向增大时,第一阶固有频率逐渐减小,第二阶和第三阶固有频率均是一直呈增大的趋势;当裂纹角度与x轴的夹角从负方向增大时,前两阶阶固有频率均呈现增大的趋势,第三阶固有频率先减小后增大.

图3 d/a=0.4,c/b=0.5时α不同裂纹角度的前三阶频率参数 Fig. 3 First three order Frequency parameters ωa2 of different crack angles at crack angles at d/a=0.4 and c/b=0.5

4.2 带裂纹悬臂Mindlin板的模态分析

得到各阶固有频率之后,本文计算了前三阶的振型,针对Mindlin板分别给出了挠度W和转角ΦxΦy的模态图.

采购管理模块为产品制造提供配套物资保证,采购管理主要功能需求包括:采购需求、采购计划、价格审批管理、采购订单、采购到货检验、采购发票管理等。

图4 裂纹板在d/a=0.4,c/b=0.5,α=15°时 有关W的前三阶模态函数图 Fig. 4 First three order model of W for the crack plate at d/a=0.4,c/b=0.5,α=15°

图5 裂纹板在d/a=0.4,c/b=0.5,α=15°时 有关Φx的前三阶模态函数图 Fig. 5 First three order model of Φx for the crack plate at d/a=0.4,c/b=0.5,α=15°

图6 裂纹板在d/a=0.4,c/b=0.5,α=15°时 有关Φy的前三阶模态函数图 Fig. 6 First three order model of Φy for the crack plate at d/a=0.4,c/b=0.5,α=15°

从图4~6中裂纹板的模态函数图可以看出,裂纹的存在使得各阶的模态函数图中出现了位移和转角不连续的现象.

5 有限元仿真

本文应用Ansys有限元分析软件进行建模和求解,由于本文研究的是悬臂裂纹Mindlin板,因此选取壳单元进行建模.图7为Ansys中截取的一张有限元模型图.

图7 Ansys中裂纹板模型图以及局部放大图 Fig. 7 The FE model of crack plate and its partial enlarged detail in ANSYS

有限元中材料和尺寸比例与前面的理论计算保持一致.其中计算结果与理论计算结果对比和误差分析如表2和表3.表2对不同裂纹长度下仿真与理论计算的结果进行对比.表3对不同裂纹位置下仿真与理论计算的结果进行对比.

图8为Ansys有限元仿真得到的裂纹板前三阶模态函数图.

将图8与上一节理论计算中的图4对比可以看出,同样的裂纹参数下,其振动的模态函数图形式是一致的,可以验证理论计算求解模态函数的正确性.

2 不同裂纹长度悬臂板的前三阶频率参数 Table 2 The first three order frequency parameters ωa2 of the cantilevered square plate with different crack lengths (c/b=0.5,α=15°)

d/aTheoreticalresultsAnsysError12312312303.44714.64121.1663.44214.71121.3730.15%0.48%0.97%0.13.43614.50821.0073.43114.61521.2650.15%0.73%1.21%0.23.40214.27120.6643.41014.35520.8680.23%0.59%0.98%0.33.37413.89120.1853.37113.90620.2950.09%0.11%0.54%0.43.32013.34019.5903.31013.2819.5980.30%0.45%0.04%0.53.20412.68418.9343.22112.49218.8370.53%1.54%0.51%0.63.12811.96418.2833.09811.57918.0870.97%3.32%1.08%

3 不同裂纹位置悬臂板的前三阶频率参数 Table 3 The first three order frequency parameters ωa2 of the cantilevered square plate with different crack location (d/a=0.4,α=15°)

cTheoreticalresultsAnsysError123123123Intactplate3.44714.6421.163.44214.7121.370.15%0.48%0.97%0.15b3.15113.6820.423.03312.9520.303.89%5.64%0.59%0.20b3.15313.5320.803.05812.9320.643.11%4.64%0.78%0.25b3.18113.3421.043.10112.920.832.58%3.41%1.01%0.30b3.20813.1921.003.14712.9120.811.94%2.17%0.91%0.35b3.24713.1420.833.19312.9620.591.69%1.39%1.17%0.40b3.28513.1720.513.23613.0420.261.51%1.00%1.23%0.45b3.32013.2320.063.27513.1519.901.37%0.61%0.80%0.50b3.34513.3419.593.3113.2819.601.06%0.45%0.05%0.55b3.36013.4719.303.34113.4319.400.57%0.30%0.52%0.60b3.37513.6119.153.36713.5919.340.24%0.15%0.98%0.65b3.39213.6719.123.38813.7619.410.12%0.65%1.49%0.70b3.40013.9419.173.40413.9319.620.12%0.07%2.29%0.75b3.14814.1019.333.41714.0919.927.87%0.07%2.96%

图8 时有限元仿真的前三阶模态函数图 Fig. 8 First three order model of crack plate at in ANSYS

将仿真的数据与理论计算结果进行对比,发现当裂纹参数变化时,两种结果均呈现了相同的变化规律,并且有限元仿真与理论结果的误差可以控制在一定的范围内.通过验证可以说明理论推导的正确性.

6 结论

本文利用Ritz法求得了带边角裂纹悬臂板的固有频率和模态函数,进而研究了不同裂纹参数对板固有频率的影响.并应用Ansys有限元分析软件进行验证.通过本文的研究发现:

(1)本文推导的附加角函数的模态函数可以表示裂纹处位移及转角的不连续性.

(2)通过对本文理论分析得到的固有频率和模态与有限元软件Ansys得到的结果进行对比验证了本文理论推导和计算的正确性.

(3)通过分析裂纹参数对固有频率的影响表明:①裂纹的存在使得局部刚度减小,不论对于第几阶固有频率,均是裂纹长度越长,板固有频率就越小.②裂纹靠近节线位置时频率减小的越显著.③随着角度从负向正增大时,第一阶固有频率逐渐减小,第二阶和第三阶固有频率均先减小后增大.

近年来,网络技术发展也非常的迅速,而在智能电网中的继电保护装置,其实本质上也是一个计算机装置,作用就是对系统进行全面的监控工作,对整个网络的运行进行指导,让智能电网的系统运行能够高效和安全。同时,继电保护装置在对电力信息的数据和故障信息进行获取的时候,会及时反馈信息到系统的网络控制中心。另外,随着我国智能电网的发展越来越迅速,在自动化程度方面,电力系统也做得越来越好,这样就能够保障了电力系统的安全和高效。

本文对裂纹板的研究在损伤探测方面具有重要的应用价值,可利用本文的研究成果来探测裂纹的位置和长度,并且为以后裂纹板的非线性振动研究奠定了理论基础.

参 考 文 献

1Yuan J, Dickinson S M. The flexural vibration of rectangular plate systems approached by using articial springs in the Rayleigh-Ritz method. Journal of Sound and Vibration, 1992,159:39~55

2Yuan J, Young P G, Dickinson S M. Natural frequencies of circular and annular plates with radial or circumferential cracks. Computer & Structures, 1994,53:327~34

3Liew K M, Hung K C, Lim M K. Solution method for analysis of cracked plates under vibration. Engineering Fracture Mechanics, 1994,48:393~404

4Huang C S, Leissa A W, Liao S C. Vibration analysis of rectangular plates with edge V-notches. International Journal of Mechanical Sciences, 2008,50:1255~1262

5Huang C S, Leissa A W. Vibration analysis of rectangular plates with side cracks via the Ritz method. Journal of Sound and Vibration, 2009,323:974~988

6Huang C S, Leissa A W, Li R S. Accurate vibration analysis of thick, cracked rectangular plates. Journal of Sound and Vibration, 2011,330:2079~2093

两组学生中都有大约三分之二的学生热衷使用在线词典,在线词典虽然在用法例句上与纸质词典差别不大,但从另一个侧面却反映了学生在学习上的惰性,表面上在线词典省去了学生翻阅书页的时间,但是他们同样也要输入单词。调查中还发现有不少学生在使用手机词典,而手机在线词典每屏所显示的信息量太少,要不停地来来回回翻页,才可以把一个词的用法读完,而纸质词典就没有这方面的问题,学生可以一览无余,反复阅读,以增强记忆。数据分析显示,目前在线词典的使用效果不如纸质词典,这和学生对电子词典的使用方法和电子词典的质量有关。

封国生介绍,当前,通过对病种指标统计、病种现状绘制、病种缺陷查找、病种发展规划的递进式分析,朝阳医院形成了以病种为核心的精细化分析报告。分析系统还向纵深建立了床位管理模型、术前住院日管理模型、医保费用管理模型、大型医疗设备管理模型、高值耗材管理模型、平均住院日分析等多个模型。如平均住院日模型将术前住院日、住院天数、病种例数和主治医师等因素进行分析;医师执业能力评价模型则从医师管理病种的难度、质量、效率、效益、缺陷等维度对医师进行综合评价。

7Huang C S, McGee O G, Chang M J. Vibrations of cracked rectangular FGM thick plates. Composite Structures 2011,93:1747~1764

8陈丽华,孙玥,张伟. 三阶剪切变形板的振动特性研究. 动力学与控制学报, 2013,11 (Chen L H Sun Y, Zhang W. The study of the vibration characteristic of plate using the third order shear deformation theory. Journal of Dynamics and Control, 2013,11 (in Chinese))

一些英语典故性成语不能直译其字面含义,只能翻译其隐含意义,这是英语典故性成语翻译中的常见情况。比如“Achilles'heel”,如果译成“阿克勒斯的脚后跟”,大多数中国人都无法理解。所以把它翻译成“唯一的致命弱点”更合适,因为这是它的隐含意。

9陈丽华,孙玥,张伟. 带边角裂纹悬臂矩形板的振动特性研究. 力学与工程应用, 2012,14:77~80 (Chen L H, Sun Y, Zhang W. Study of vibration characteristics of cantilever rectangular plate with side crack . Application of mechanics and Engineering, 2012,14:77~80 (in Chinese))

10 曹志远,杨昇田. 厚板动力学理论及其应用. 北京:国防工业出版社,1983:57~85 (Cao Z Y,Yang S T. Thick plate theory and its application. Beijng: Geological Publishing House, 1983:57~85 (in Chinese))

陈丽华,薛坚,张伟
《动力学与控制学报》2018年第01期文献

服务严谨可靠 7×14小时在线支持 支持宝特邀商家 不满意退款

本站非杂志社官网,上千家国家级期刊、省级期刊、北大核心、南大核心、专业的职称论文发表网站。
职称论文发表、杂志论文发表、期刊征稿、期刊投稿,论文发表指导正规机构。是您首选最可靠,最快速的期刊论文发表网站。
免责声明:本网站部分资源、信息来源于网络,完全免费共享,仅供学习和研究使用,版权和著作权归原作者所有
如有不愿意被转载的情况,请通知我们删除已转载的信息 粤ICP备2023046998号