显式波动有限元分析的高精度连分式人工边界方法

来源 :北京工业大学 | 被引量 : 2次 | 上传用户:eric_vl
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
动力荷载作用下重大土木基础设施工程需要考虑结构与无限域地基的动力相互作用,形成无限域介质中波动传播问题,可以采用有限元法结合人工边界条件的时域数值方法求解。显式有限元法适合求解波传播问题,动力有限元方程的显式时间积分算法通常应为二阶精度单步法,能够考虑非对角阻尼矩阵,并且具有良好的稳定性。高精度时域人工边界条件能够精确模拟被截去无限域的辐射阻尼特性,是工程数值计算领域的一个热点和难点课题,目前的主要困难在于,基于频域解析解的精确人工边界条件直接转化到时域面临计算量和存储量巨大的时间卷积运算。一个现有的解决方案是,在频域内采用有理函数逼近动力刚度,利用最小二乘优化算法确定有理函数的系数;然后将有理函数描述的动力刚度关系等价的转化为时间二阶结构动力学方程系统,与有限元方程耦合采用显式时间积分方法求解。上述方案的不足在于,当有理函数系数较多时,优化算法效率低甚至难于收敛。最近的有效方法是动力刚度的连分式逼近,通过代数计算由动力刚度所满足的方程确定连分式系数,避免优化计算,适合系数较多以及矩阵值问题。本文针对半空间出平面SH波(二维标量波)问题,提出一种新的显式时间积分算法以及一种新型的动力刚度连分式,发展了显式有限元结合高精度人工边界条件的时域波动数值模拟方法。具体研究工作如下:1、提出一种新的求解动力有限元方程的显式时间积分方法。该方法的位移公式为前一时刻反应的二阶泰勒展开;速度公式为预估-修正形式,预估速度为前一时刻反应的一阶泰勒展开,预估加速度为预估速度在当前时刻的动力方程,修正速度为采用预估加速度的梯形求积公式;加速度公式为当前时刻的动力方程。当质量矩阵对角时该方法为二阶精度显式方法,能够考虑非对角阻尼矩阵和复杂的非线性内力项;该方法为适合变时步分析和简单起步的单步法,稳定性优于现有几类方法。通过线性和非线性单和多自由度结构动力学问题算例,与现有几类方法对比,验证了提出方法的有效性。编写了基于该方法的FORTRAN有限元程序,与ABAQUS有限元程序对比,验证了有限元程序的正确性。2、提出一种新型的动力刚度连分式,建立了相应的高精度时域人工边界条件。证明了现有几类高精度时域人工边界条件等价于在频域内对动力刚度进行高频连分式逼近,指出高频连分式具有零频(静力)刚度不准确的问题;指出Song等提出的双渐近连分式具有待定参数难于确定的不足;提出一种具有精确零频刚度的连分式。通过将连分式转化为有理函数,证明了连分式的稳定性。基于已经发展的将有理函数描述的动力刚度关系等价的转化为时间二阶结构动力学方程的处理方法,建立了基于动力刚度连分式的高精度时域人工边界条件。该边界条件与有限元方程耦合后可以采用提出的显式时间积分方法求解。通过波传播问题数值算例,与大区域有限元模型结果进行对比,验证了建立的高精度连分式人工边界条件的有效性。
其他文献
天冬氨酸激酶(Aspartate kinase,AK)是一种控制生物合成工业中重要天冬氨酸族氨基酸的关键酶,受支路代谢终产物调控。本研究通过改善AK活力,达到解除其在天冬氨酸族氨基酸代
当今世界充满了各种各样的关系,这些关系用拓扑学的原理可以简化为一张一张网络图。因此复杂的世界可以用一张张网络图表示。关系紧密的节点可以划分在同一个社区中,因此可以
额外维理论为理解物理世界提供了一种新的思路。在解决诸如层次问题等方面取得了巨大的进展。另一方面为了解决广义相对论面临的诸如奇点问题、不可重整化问题、暗能量问题、
作为一门古老而基础的数学分支,早期的概率论主要研究的是生活中的可能性问题。随着其理论的不断深入和研究方法的革新,逐渐在金融、医学、军事、甚至社会科学中有了广泛的应
本实验依据精子顶体表面存在大量促性腺激素释放激素受体(GnRHR)的特点和绿脓杆菌外毒素A转膜区(TLD)有将外源蛋白跨膜主动转运至细胞内的能力以及超嗜热古菌组蛋白(HphA)高
在工程领域中,常发生起源于断裂或终止于裂纹扩展的灾难性破坏事故。从离散裂缝断裂力学角度加以数值模拟研究结构的断裂破坏过程,对于提高工程设计水平有着重要的意义。扩展
《魅力汉语》综合系列教材是由广西师范大学出版的一套精品教材。分为初级、中级和高级三个阶段,在我西南地区使用较多。本文参照《高等学校外国留学生汉语教学大纲(长期进修)》从课文、生词、语言点和练习四个部分对《魅力汉语》(高级)综合教材进行分析。经考察发现,此套教材从整体上看:针对性较强,循序渐进、教学难度逐步增加,版式设计精美、图文并茂等特点。从内容上看:课文题材丰富,内容生动形象。课文之间生词数量差
本文主要运用整函数的相关理论和亚纯函数的Nevanlinna基本理论,来研究几类整系数线性微分方程解的复振荡性质,全文分为以下四章.第一章,作为全文的预备知识,简要介绍了整函
我国推行公司注册资本认缴制改革以来,公司的财务自治得到扩展,筹资效率得以提升。但同时,制度变革催生了股东出资义务加速到期问题。然而,我国现行法律仅仅规定了破产和清算情境下的加速到期规则;现有的学理讨论也大多从维护交易安全出发研究公司债权人诉请股东出资义务加速到期问题,鲜有从公司自治角度出发分析公司财务自治和股东期限利益的冲突。实际上,拒绝或者不当限制股东出资义务加速到期规则的适用,不但可能损害公司
本文创新性地将工业上广泛应用的质量控制图方法应用到对我国各省市宏观经济数据质量的研究上来,并对传统的质量控制图进行改进,提出了两种全新的用于检测省市宏观经济数据的