弹性地基梁动力响应分析的新方法.pdfVIP

·4“· ‘工程力学》增刊 1999年 弹性地基梁动力响应分析的新方法车 黄伟江余勇潘小强罗恩 (中山大学力学系，广州510275， 提要 本文首先建立有阻尼弹性地基粱动力学的一类变量广义}hnin明型拟变分原理，它能反 映动力学初值一边值问题的全部特征，然后，以该拟变分原理为基础．提出时间域用五次H咖ite插值毒 项式插值的时间子域法与空间域用有限元法相结合的新方法，并用该法分析弹性地基粱的动力响应问题。 算例表明，这种新方法的精度和计算效率都明显高于目前国际上常用的wn啪．o法和Ne硝ark一日法 关键词 弹性地基粱，广义}IaIIIin咂型拟变分原理，有限元法，时间子域法，动力响f世 1 i I ton型拟变分原理 弹性地基梁一类变量Ham 按文【I】的思路．可以建立有阻尼弹性地基梁动力学的一类变量广义Hamilton掣拟变 分原理，对于边界条件x=o：w=万，尝=歹：x=，：M=厨，Q=虿和初始条件 凹 w(r，o)=茹。，丽掣=瓦的情形，其表达式为 d， 州=绯丽(铲警(爿一≯十晰詈～ 一孵w，哥瞎H豢一百剧。破+膀一厨觌…， +I{一咯劬+瓦西c删+pc蛐卜引莉笔竽卜=。 对于其他边界条件，上式中的划线项应作相应的改变。 当梁端支承分别为通常的简支、固支或自由、以及位移满足位移边界条件和位移初 羹伟虹．男，1975．6出生，研究生 +广东省自然科学基盒和田家自然科擘基金资助项目 《工程力学≥增刊 1999年 -48S· 始条件时．(11)式就变成为 mⅢ照而㈤一等㈢o：+加卜 一咖等～一多詈鼬=。 其中．丽、i和f分别是单位梁长的质量、粘性阻尼系数和弹性地基系数 2动力响应分析的新方法 2 1有限元法 对空问域用有限元法进行离散n按一般有限元的步骤．对任一单元P。．其节-_f立fj；勾 矿(，))=【w，(，)只(，)w∥)B(f)】7 12I I} l 位移模式为： w…(x，，)=【Ⅳ(x)】铮。’(，) (2 21 l (2 2)代入L12)．得： m阳I)=J{f[圭臼}『眇】臼)．圭白F医】臼)一圭矗)『医k}+泸】白)]竹} 瞳。：、 一』{c7}rFp白如一曲}rmp臼}=o 其中 沁】=∑￡而【Ⅳ(r)】7【Ⅳ(x)】出k】-∑￡田【Ⅳ’(x)r【Ⅳ”(r)】c 【(1】= k】-∑F【Ⅳ(x)】7[Ⅳ(x)】c打 【尸】=∑f，【Ⅳ(x)】出 白(f)}=∑驴(f) 以上各式均是按有限元的一般方法对号集总而成．没其最后的自由度数【即待定的 符点位移函数数目)为H。此外．在进行数值运算之前，还应对边界条件进行处理， 2．2时间子域法 ·486· 《工程力学≥增刊 I：J99年 时问子域法的思想是把所考察的整个时间响应历程划分成若干个时间f域．卉：任一 f域末 时问r域．f：．用适：§的插值函数逼近待定的函数．根据动力学的变分原理，求解÷B 端的状态值：然后．把前一个时间子域的术端状态值作为下‘小时间产域的韧婚：I《苍值． 重复t一步的计算；如此反复进行．直至最后一个时间于域， 对于任意时间子域[f，，f¨】，令f。=fH—f，，则局部时间坐标表示的时间于域为 ，∈[