【摘 要】
：
多次透射公式(multi-transmitting formula, MTF)是在近场波动数值模拟中一种广泛应用的人工边界条件, 具有形式简单、精度可控和通用性好的优点, 但高阶MTF与有限元方法相结合有时会出现飘移问题. 现有的几种MTF消飘方法往往会显著地影响边界精度, 为此本文提出一种新的消飘因子修正MTF格式, 能够在较高精度水平下实现对飘移问题的有效控制. 该方法保持MTF的一次透射项不变, 仅对各高次透射误差项进行修正, 从而大幅降低了因消飘因子造成的精度损失. 消飘因子设置格式确保在零频和零
【机 构】
：
南京工业大学工程力学研究所,南京 211816;中国地震局地球物理研究所,北京 100081;南京工业大学交通运输工程学院,南京 211816
论文部分内容阅读
多次透射公式(multi-transmitting formula, MTF)是在近场波动数值模拟中一种广泛应用的人工边界条件, 具有形式简单、精度可控和通用性好的优点, 但高阶MTF与有限元方法相结合有时会出现飘移问题. 现有的几种MTF消飘方法往往会显著地影响边界精度, 为此本文提出一种新的消飘因子修正MTF格式, 能够在较高精度水平下实现对飘移问题的有效控制. 该方法保持MTF的一次透射项不变, 仅对各高次透射误差项进行修正, 从而大幅降低了因消飘因子造成的精度损失. 消飘因子设置格式确保在零频和零波数情形下能够满足GKS准则, 从理论上保证了消飘目标的实现. 进一步给出该方法的高阶统一形式, 并将传统的消飘因子修正MTF的方法归结为该统一形式的一个特例. 通过反射系数分析, 证明本文方法不仅具有精度优势, 而且消飘因子选取的适应性更强、取值范围更广. 数值算例表明, 本文的消飘因子对波动能量比较集中的法向和小角度透射波动的模拟精度影响很小, 在控制高阶MTF飘移问题和保持模拟精度方面, 均能够取得明显效果.
其他文献
活性流体在用于开发新材料方面具有巨大潜力, 满足这一需求就要定量掌握活性流体所表现的特殊力学行为, 特别是流变行为. 扩展布朗运动方程, 建立自驱动活性粒子的运动模型, 基于反向非平衡法确定活性流体的黏度, 考察活性粒子体积分数、直行速度和转向扩散系数对活性流体流变行为的影响规律, 确定活性流体特殊流变行为的形成机理. 结果表明, 活性流体的流变曲线可被划分为黏度下降区、过渡区和牛顿区; 活性粒子体积分数越高, 活性流体的非牛顿特性越显著, 活性粒子的直行运动引起活性流体在低剪切速率区域黏度下降, 直行运
建立含裂隙均匀基质中的稳态渗流模型,通过数值模拟得到渗流速度和压力分布,分析不同分布型式的平行和交叉裂隙对基质渗流的影响,并结合油气开采和驱替评价不同裂隙可能带来的后果.研究发现,在裂隙端部压力发生突变,相应的渗流速度也出现急剧变化;对于基质而言,裂隙入口端等同于汇,而开口端等同于源.平行对称分布的裂隙导致周围基质中流速降低,裂隙流量占截面流量的主要部分,而不对称分布的平行裂隙有利于基质中的流动.
通过分析实际设计过程及高温管道中各应力分量的影响因素,发现管系的热应力大小主要由管系的走向和各管段的实际长度决定.由此,提出了管系管段长度比的概念.通过力法对管段长度和热应力进行分析发现:(1)管系相邻管段长度比越大,各种排布方式中最大的热应力越大;(2)以目前常用的P91管道为例,给定管系的实际尺寸和运行参数,得出了管系所允许的最小管段长度;(3)在管系最小管段长度确定的前提下,进一步得出管系所允许的最大管段长度比;(4)随着管段数目的增加,管系允许的最短管长趋近于某个确定值;(5)在管系最小管段长度确
提出了一种从离散分子动力学模拟(MD)到连续介质弹性有限元计算分析(FEA)的过渡方法, 简称MD-FEA方法. 首先通过MD计算获得晶体材料原子的移动位置, 然后根据晶体结构的周期性特征构造连续介质假设下的有限单元变形模型, 进一步结合材料的力学行为本构关系获得应变和应力场. 为了检验MD-FEA方法的有效性, 将该方法应用于详细分析Al-Ni软硬组合两相材料纳米柱体的拉伸变形问题和基底材料为Al球形压头材料为金刚石的纳米压痕问题. 采用MD-FEA方法获得了上述两种问题的应力?应变场, 并将计算结果分
混沌和分岔使得多稳态俘能系统的非线性动力学响应对系统结构参数非常敏感, 导致了系统的非线性特性正向设计比较困难. 为了定量地表征非线性恢复力与结构参数的关系, 提出了一种多稳态俘能系统的准确磁力建模方法. 推导了多稳态俘能系统端部磁铁和外部磁铁的相对距离和转角位置, 并采用磁荷理论建立了多稳态系统的非线性磁力模型. 通过搭建实验平台测量了不同结构参数条件下多稳态系统的非线性磁力, 并对比了本方法与传统方法和实验测量的结果. 结果表明: 本方法的磁力计算结果与实验测量值吻合较好, 双稳态系统和三稳态系统的磁
本文运用Green函数法求解了多裂纹Euler–Bernoulli曲梁(Euler–Bernoulli curved beam,ECB)强迫振动下的解析解,并且考虑了阻尼效应.采用分离变量法、Laplace变换法和矩阵传递法得到了两端简支的多裂纹Euler–Bernoulli曲梁的Green函数解.通过研究表明,将半径R设置为无穷大,可以简化为Euler–Bernoulli直梁(Euler–Bernoulli beam,EB)振动模型.数值计算中,通过与已有文献中的解析解做对比,验证了解的有效性,而且进一
参变气动弹性建模是可变后掠翼气动弹性研究中的难点之一. 当地建模技术是一种构造线性参变 (LPV)模型的实用方法, 但一直缺乏有效途径来解决当地气动弹性模型的不一致问题. 气动弹性模型的不一致性体现在当地结构动力学模型和非定常气动力模型随参数变化的不连续性. 本文提出了一种自下而上的方法, 对变后掠翼不一致的当地气动弹性模型进行了一致性处理. 首先, 采用匈牙利算法跟踪结构模态并按模态分支进行排序, 使得匹配后的模态能够保证结构动力学模型的一致性; 其次, 对有理函数拟合表达式中的系数矩阵进行缩放处理,
波浪在大尺寸结构表面产生不可忽略的散射波, 该散射波在多柱体体系中继续传播, 并在同体系中的其他柱体上产生高次散射波. 本文基于椭圆坐标系和绕射波理论首先推导了波浪作用下椭圆单柱体产生的散射波压力公式, 随后考虑该散射波在多柱体系中的传播, 将其视为第二次入射波, 推导出柱体上第二次散射波压力公式, 同理可以推导出高次散射波压力公式, 最后得到椭圆多柱体波浪力解析解, 并用数值解验证了本文解析方法的正确性. 本文以双柱体和四柱体体系为例, 分析了不同参数(波数、净距、波浪入射角度等)下, 高次散射波对柱体
问题的简化分析往往基于假设条件,说明假设的合理性不仅是结论适用范围的基础,也是培养学生解决复杂系统工程问题的必要环节,本文采用量级估计的方式,说明材料力学基本变形量之间的关系以及常见简化假设的基础,为教学活动中解释相关问题提供参考.
以往针对结构物垂直贯穿冰层破裂的研究多不考虑水的作用, 与实际应用场景不符. 本文应用 LS-DYNA 有限元软件建立了基于结构化-任意拉格朗日欧拉(S-ALE)流固耦合方法及罚函数接触算法的冰?水?结构物耦合作用数值模拟方法. 采用欧拉算法描述空气域和水域, 采用拉格朗日算法描述圆柱体结构和冰层结构, 使用弹塑性应变率模型表征冰材料力学性质. 自主搭建了圆柱体垂直贯穿冰层试验台架, 验证了有限元方法计算结构物?冰层相互作用问题的可行性. 通过模拟圆柱体垂直出水破冰过程, 并与无水环境下圆柱体垂直贯穿冰层