粘弹性结构动力学分析的一种数值方法.docVIP

粘弹性结构动力学分析中的一种数值方法 彭凡 傅衣铭 （湖南大学工程力学系，长沙410082,中国） ：针对材料具积分型本构关系，及松弛模量为Prony级数形式的粘弹性结构动力学问题，本 文结合Newmark方法与Taylor方法，建立了计算该类问题的一种数值算法。且以简支梁为例， 应用该方法具体地分析了考虑线性与非线性粘弹性时梁的强迫振动响应。 关键词：粘弹性动力学数值方法响应 1引言 随着人们对结构材料物理与力学性质了解的不断深入，以及新型材料的广泛应用，粘弹性结构的动力学 研究受得了愈来愈多的重视，数值计算己成为一种主耍的分析手段。文［1］基于Newmark方法建立了粘弹性结 构动力学响应的有限元法，但只涉及到线性问题，而且在每一计算步，卷积积分的计算量较大。桂洪斌等⑵ 提出将粘弹性结构的动力学方程进行Laplace变换，然后在相域屮求解问题，显然这种处理方式同样只适应于 线性情况。当考虑几何，物理包括损伤等非线性因素时，粘弹性结构动力学的数值分析就变得十分复杂与困 难了。文［3,4］通过将微分一积分型非线性动力学方程化成高阶的微分方程，最终由Runge-Kutta法來获得数值 解，但只有当材料为标准线性固体或Prony级数取较少项数时，这种方法才比较容易实现。本文针对材料服 从积分型本构关系，且松弛模量为Prony级数形式的粘弹性结构动力学问题，建立了从时域内直接求解的数 值算法，它是基于Newmark方法与Taylor方法而得出的。其I1 Taylor方法为卷积积分的递规算法，能使计算 量显著降低［5］。文中通过对粘弹性梁的受迫振动分析来说明方法的应用。 2简支粘弹性梁受迫振动的动力学方程 (1)考虑一简支梁，其跨度为L,高为h,中点受横向周期激励H^Ot。设材料具非线性粘弹性，可由Leaderman161 本构关系描叙，则有 (1) 恥（￡（） + 式中Eo = E（0） , E（r）为松弛函数，g（w）为应变w的非线性函数： = fie2 + 冷 英屮0与7为常数。在小挠度情况下，梁的受迫振动方程为： 12 a%48() dx 12 a%4 8() dx2 I dx2 (3) dr = H 8 (x- L/2)sin dr = H 8 (x- L/2)sin 0t 3(r-r) 12 dx4 80 dx2 [ 式中o A分别为梁的质量密度及横截面面积，力为Dirac函数, 满足两个简支端条件，即w（0j） = %（0』） =w（厶『）=%（厶/） = 0的挠度w（x,/）取为 毗山）=工￡?⑴sin字 （4） 为说明问题起见，式⑷中只考虑R = 1的项，且令士⑴=/（/） 0将式⑷代入式⑶后，作Galerkin积分，并记 D(t) = E(r)/E0且取无量纲位移q(t) = f(t)/h ,经运算和整理后得到: q + gjtD *dq +a? = H() sin 价 兀\2pA\jL80丿Hq =2H 兀 \2pA\jL 80 丿 Hq = 2H LpAh D^dq- q(t) + 硕-力d(t~ T) D * 姑=q + j 響 _ ⑺ I d(t-T) 当材料为线粘弹性吋，简支梁的受迫振动方程为 q + co1 D * dq = Hq sin 0t (6) 式(5)与式(6)屮的￡)(/)可表为Prony级数的形式 K D ⑴= X()+》X/7 (7) k=\ 其中材料参量Xq、Xz 0,且X。+ ￡x严1。 *=1 3数值算法 TOC \o 1-5 \h \z 对式(5)进行数值求解，设时间增量步为,基于Newmark171方法，/ + 4时刻的运动方程为 ；讥》如皿+2 *窃)心=H() sin(r + A/) ⑻ 又 彳+少=--彳)— 0 -[亠-1彳 ⑼ oAt oAt \2a 丿 式中。和2是按积分精度及稳定性耍求而决定的参数。再市Taylor181的卷积积分数值递规算法有 (D*dg)g =““+[?+$ -4)+ 城?】 (10) (11)(D * dqy )心=亿+i I- 1 + 妙鹽 (11) K 1 -Q/ K K 式中 如= x° + ￡x』：m,其中 們)=匕1胡叭 心=险+￡器+”必*。(汀+；『 1 C. A/ A=1 1 其中 灯严二圍)+X用(4-必)]严“,勢)彳蹭+X卅((4)3_(皿)3)”2 而为时刻q的值。且当心1,即20时，= 0 o将式⑼~(12)代入式⑻后得到 A(G+zv)‘ + Bqg + F = 0 (12) 式屮 A =弘“+| ； B = y° + ； F 二—[h° sin(r + Ar) + 齐么 + y2q, + y? + e叽q( + 昭屮(q/ -阳—S、_ 在计算步/, 及F已知，通过Newton法解非线性方程(12)来求出