单集中阻尼张紧弦动力时程分析方法对比研究.pdf
《单集中阻尼张紧弦动力时程分析方法对比研究.pdf》由会员分享,可在线阅读,更多相关《单集中阻尼张紧弦动力时程分析方法对比研究.pdf(6页珍藏版)》请在咨信网上搜索。
1、Jun.20232023年6 月JOURNALOFDONGCUANERSTTYOFTECHNOLOGYVol.30No.3第30 卷第3期学院学报东莞理单集中阻尼张紧弦动力时程分析方法对比研究王宪东1王弘阳?(1.重庆交通大学土木工程学院,重庆400074;2.兰州理工大学理学院,甘肃兰州730050)摘要:集中阻尼张紧弦动力时程分析方法对比研究,采用精细时程积分法与纽马克(Newmark)法两种方法进行算例分析,主要探讨精细时程积分法与纽马克法对积分步长的敏感度。研究发现,精细时程积分法和纽马克法,可以无条件稳定,当积分步长逐渐增大时,精细时程积分法依旧可以精确计算,得到准确结果,纽马克法计
2、算精度受限,导致出现误差,故在进行单集中阻尼张紧弦的动力时程分析时,精细时程积分法优于纽马克法。关键词:张紧弦;阻尼;振动;状态空间;积分步长;精细时程积分法;Newmark法中图分类号:0 30 2文献标志码:A文章编号:1 0 0 9-0 31 2(2 0 2 3)0 3-0 0 9 5-0 6在桥梁工程中,斜拉桥上斜拉索的长期性能问题一直是备受关注的热点1-3。斜拉索自身的细长使得它们容易振动,为了研究斜拉索振动特性,可以将其简化为一根两端固定的弦。外置阻尼器在实际工程中应用广泛,通常用来对斜拉索进行减振,因此与动力学有关的工程索结构的减振问题必然涉及到阻尼混合弦的动力时程分析问题。从1
3、 9 8 1 年开始,关于集中阻尼弦的问题外国学者 Carne4已经有研究。关于阻尼混合弦的本征问题备受研究界关注5-6 ,Main7研究了单阻尼混合弦本征问题,作者利用复特征值问题的解析公式,详细研究了带有线性粘滞阻尼器的张拉索的自由振动。BM Pacheco、YFu j i n o、A Su l e-kh8为了简化桥梁工程中拉索的粘性阻尼器的设计程序,提出了一种通用的估计曲线,并通过几个数值例子说明了这种设计辅助工具的方便性。郑罡9】利用Dirac函数,在全域里建立并求解了集中阻尼弦的动力学方程,讨论并且得出了混合动力学系统的特性,文献9 所讨论的这些问题,对本研究有重要的指导作用。该领域
4、关于集中阻尼弦的动力时程分析问题研究较少,学者陈建兵、曾小树、彭勇波1 0 首次分析了粘滞阻尼器系统的刚性特征,在此基础上引人向后差分格式,数值分析结果表明此方法具有足够的精度和效率。李寿英、顾明、陈政清川研究了阻尼器对拉索风雨激振的控制效果,详细分析了阻尼器各种参数对控制效果的影响。Chen L 等1 2 提提出了一种弦系统自由振动和强迫振动分析的通用数值模型,作者使用该模型对比较典型的弦系统进行了自由振动分析,表明了该方法的有效性和准确性。针对动力时程分析方法,学者王元丰、储德文1 3】扶提出一种将精细积分法与纽马克-(Ne w m a r k-)法耦合起来的结构动力学时程积分方法,文献【
5、1 3中通过理论推导和算例验证,表明了该方法在结构动力分析中的有效性。李靖、李彬、唐小微1 4分析了在进行动力时程分析时优化时间步长与计算精度的相关性,并且基于收敛性、稳定性、计算精度和计算效率四个方面对时间步优化方法的优劣进行了评价。精细积分方法是钟万1 51】于1 9 9 3年提出的一种高效精确的逐步积分方法,它是显式且无条件稳定的,并适用于非正交阻尼体系。结构动力方程的精细时程积分法1 5】对于线性的定常结构动力系统的分析提供了精细时程积分方法,对于本文研究动力时程问题提供了方法思路。综上可知,该领域运用两种及两种以上动力时程分析方法,针对集中阻尼弦系统的动力时程问题进行的对比研究较少,
6、同时此类研究多集中于单一方法的探究。动力时程分析的部分研究虽然运用多种方法进行对比分析,但并不针对阻尼弦系统,也未能给出影响计算精度的因素。本文针对单集中阻尼张紧弦进行动力时程分析,运用两种方法进行有限元算例分析,其中精细时程积收稿日期:2 0 2 2-0 8-1 8作者简介:王宪东(1 9 9 7 一),男,甘肃天水人,硕士生,主要从事工程力学研究,Email:1 2 8 7 548 9 0 7 q q.c o m。962023年东莞理工学院学报分法在状态空间中求解。通过与理论解进行对比分析,证明了精细时程积分法精度更高,更贴近真实解,同时探讨了影响两种方法计算精度的因素,以及阻尼位置改变时
7、动力时程的变化情况。1数学模型1.1研究模型的建立本研究对象为一根两端固结的带集中阻尼器的张紧弦,弦两端承受张力T,除此之外不受任何外荷载作用。此研究对象作为单阻尼系统,阻尼所在位置可任意改动,研究对象模型示意图如图1 所示。LTmT图1带集中阻尼器的张紧弦其中:L为弦长,m为弦的线质量密度,T为张力常数,c为阻尼系数,以弦的左端点为坐标原点,以该点指向弦右端点的方向为正方向建立坐标,阻尼在弦上的位置为x。1.2集中阻尼弦初值问题的有限元动力学控制方程根据有限元方法,阻尼弦系统的振动控制方程为Mu+C u+Ku=O,(1)其中【M,【C ,【K】分别为结构的总质量,阻尼,总刚矩阵,【M,【C】
8、,【K】均为对称矩阵,节点位移列阵为 u。u ,【u ,u 分别为加速度、速度、位移向量。1.3有限元模型本文有限元研究模型为单位阻尼弦9 ,单位阻尼弦为具有单位长度、单位张力和单位线密度的阻尼弦,本文所研究的单位阻尼弦带有单项阻尼c。将张紧弦在进行有限元分析时划分为1 0 0 0个桁架单元,从左至右共有1 0 0 1 个节点,令弹性模量E=2.110ll。1)质量矩阵:单元质量矩阵为集中质量矩m;0阵,M.=之,将单元质量矩阵进行20m;组装,得到最终的总质量矩阵。2)刚度矩阵:本文中的刚度矩阵是通过考虑几何非线性由张力产生的矩阵,因此单元刚度-1,将单元刚度1矩阵为K.=二.阵进行组装,得
9、到最终的总刚度矩阵。3)阻尼矩阵:本文中阻尼矩阵是由集中阻尼器产生的非经典阻尼矩阵,因此阻尼矩阵为n行n列的方阵Cnxn(n 为系统总的自由度数),矩阵中元素为Ci。假设阻尼位于第N个节点处fc,i=j=NCo,其它若此时阻尼器安装位置未设置节点,阻尼等效到阻尼器相邻的两个节点处。2分析方法2.1结构动力方程的精细时程积分法精细积分方法是钟万魏1 5】于1 9 9 3年提出的一种高效精确的逐步积分方法,它是显式且无条件稳定的,并适用于非正交阻尼体系。精细时程积分法求解此问题,需要通过降阶的办法转化形式,在状态空间中求解。引人状态变量(2)V结合恒等式=,(3)方程(1)可降阶为下列一阶微分方程
10、(iu)=Hul,(4)其中0=H,B=-M-K,G=-M-CLBG(5)方程的齐次解为V(t)=T(T)(C),(6)其中【T()=exp(H),在积分步长内tEt k,t k+内,T=t-Tk。则在本步长的初始时刻t=tk,即T=0时,得到方程(4)通解(V(t)=T(T)/V(t),(7)(V(t)根据所假定初始方程得出,此时问题归结为精细地计算T(),具体求解方法见结构动力方程的精细时程积分法1 52.2Newmark法Newmark法所采用的数值积分格式,其实质就是线性加速度法的推广,Newmark法的递推公97第3期王宪东,等:单集中阻尼张紧弦动力时程分析方法对比研究式为式为1M+
11、-C+KWi+AttAt(B-2)c)FF+42(B-1)cli,+BM+M+LAtAtAt(8)在本文的研究中,参数取0.5,参数取0.5,且已知当0.5,0.2 5(0.5+)时,此算法无条件稳定,具体求解方法见有限单元法原理与应用1 6 O3算例分析本算例有限元研究模型主要参数如表1 所示。表1模型参数弦长单位密度阻尼C张力T弹性模量11112.110ll3.1算例1:中点阻尼弦此时系统固有振型存在理论解,理论解有两种情况,针对这两种情况,分为工况1 与工况2展开研究。为了分析精细时程积分法与Newmark法对积分步长的敏感度,该算例考虑了两个积分步长:0.0 1 s,0.1 s,与理论
12、解来进行对比分析,数值解按理论解给定初始条件。坐标轴标签的物理含义:t(s)为时间坐标,时间单位“秒”;x为阻尼在弦上的位置。3.1.11工况1设理论解为u=sin2Txcos2t,有限元分析设定阻尼位于中点,此时阻尼不起作用,模型两端固接。以图2 中位移均为无量纲位移。1)积分步长取0.0 1 s,总时长1 0 s,共1 0 0 0步,图2 给出了0.1 s时系统位移时程曲线(积分步长第1 0 步)与理论解曲线的对比图像,图3给出了1 0 号自由节点位移时程曲线与理论解曲线的对比图像。通过以上图像对比可知,当积分步长取0.0 1s时,精细积时程分法与Newmark法的系统位移时程曲线、1 0
13、 号自由节点位移时程曲线与理论解曲线重合,两种方法均可得到与理论解高度吻合的结果。2)积分步长取0.1 s,总时长1 0 s,共1 0 0步,图4给出了1 s时系统位移时程曲线(积分步长第1 0 步)与理论解曲线的对比图像,图5给1.0理论解0.5一一精细时程积分法一Newmark法0.04-0.5-1.00.00.20.40.60.81.0X图2第0.1 s时系统的位移时程图理论解一精细时程积分法一Newmark法0.03-0.02-0.070246810t/s图310号自由节点的位移时程图出了1 0 号自由节点位移时程曲线与理论解曲线的对比图像。1.0理论解一会一精细时程积分法0.5一Ne
14、wmark法0.0-0.5-1.00.00.20.40.60.81.0X图4第1 s时系统的位移时程图鑫A鑫鑫八鑫0.03-0.02-0.070246810t/s图510号自由节点的位移时程图通过以上图像对比可知,当积分步长取0.1 s时,精细时程积分法的系统位移时程曲线、1 0 号自由节点位移时程曲线与理论解曲线重合,与理论解高度吻合;Newmark法计算精度受限,导致出现误差,时程曲线与理论解出现偏移,本算例所研究的初值问题,不存在荷载模拟导致的误差,由此可见是算法本身导致的误差。以系统的位移982023年东莞理工学院学报时程图4为例,Newmark法的最大相对误差为11.6%,此误差使系
15、统的动力时程分析结果比实际情况要低。用显式算法求解该问题是导致Newmark法出现误差的主要原因,显式算法基于动力学方程,直接根据当前时间步的状态推导出下一时间步的状态,但此算法的计算结果对时间步长的敏感性较高。3.1.2工况2设理论解为u=(A+B)e(a+i)+(A-B)e(a-ji),12其中A=sinh(x)cos(wx),B=jcosh(ox)sin(wx)1u=(A+B)e(a+in)+(A-B)e(a-im),2其中A=sinh(1-x)cos(w(1-x)B=jcosh(1-x)sin(w(1-x)有限元分析设定阻尼位于中点,此时阻尼起作用,模型两端固接。根据文献9 本征函数形
16、式,此时工况2 理论解中参数=-1.09861229,取1 阶本征函数时,理论解中参数w=。以下图像中位移均为无量纲位移。1)积分步长取0.0 1 s,总时长1 0 s,共1000步,图6 给出了0.1 s时系统位移时程曲线(积分步长第1 0 步)与理论解曲线的对比图像,图7 给出了1 0 号自由节点位移时程曲线与理论解曲线的对比图像。0.03-0.02-0.070246810t/s图6第0.1 s时系统的位移时程图通过以上图像对比可知,精细积时程分法与Newmark法的系统位移时程曲线、1 0 号自由节点位移时程曲线与理论解曲线重合,两种方法均可得到与理论解高度吻合的结果。2)积分步长取0.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 集中 阻尼 张紧弦 动力 分析 方法 对比 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。