声散射问题的一种新的无限元方法

声波在无限大外域中的传播和衰减问题是研究声与结构相互作用的重要基础.本文提出了一种新的声无限单元耦合有限元方法于外域声散射问题的计算.应用加权残值法推导了声无限单元的刚度矩阵和质量矩阵，选取的形函数及权函数充分考虑了外域声传播的基本特点及无穷远辐射边界条件.计算实例表明这种方法对于求解无界域声场问题是非常有效的.     

港口内靠码头系泊船波浪力时域模型

为进行港口内靠码头系泊船波浪力的时域计算,提出了一种耦合数值方法。用包含船体的矩形柱面将整个港口区域划分为船体表面与柱面所围区域(内域)和柱面与整个港口边界面所围区域(外域)。内域流动由三维NavierStokes(NS)方程控制,通过VOF方法数值求解和描述其自由表面;外域流动由二维浅水波模型——改进的Boussinesq方程控制,应用有限差分三对角矩阵算法(TDMA)快速求解。内域解和外域解通过速度/压力和波面连续的匹配边界条件联系起来。应用所建立的NS/Boussinesq耦合模式同步和高效率地获得了港内波浪变形、系泊船附近三维流场、船表面波压力分布和波浪荷载。     

钢管混凝土拱极限承载力计算及相关参数分析

采用考虑材料非线性的钢管混凝土拱极限承载力计算方法,对两个模型拱进行了极限承载力计算及相关参数分析．在该方法中,对钢管混凝土拱结构采用纤维单元模型,该模型假定钢管与混凝土完全粘接,钢管对核心混凝土的套箍作用体现在以一维形式表达的核心混凝土的应力一应变关系曲线之中．针对材料非线性分析中单元内各点刚度参差不齐的特点,采用单元内设小元的方法(相当于子结构),编制了非线性有限元程序．在该程序中,计算模型完全是基于小元层次进行的,比如：单元刚度矩阵由小元刚度矩阵凝聚而成,单元节点的残余力由小元节点的残余力构成,故只需改变单元内小元个数这一个参数,就可实现对结构的重新划分,极大地降低了非线性方程组的阶数,且方便实用．在程序计算结果得到模型试验结果验证的基础上,还分析了套箍系数、截面含筋率及混凝土强度对结构极限承载力的影响。     

有限元与边界元耦合法在波浪力计算中的应用

利用有限元与边界元耦合法对三维无界区域中直立圆柱所受的波浪力进行计算，把整个求解区域分成内域和外域两部分，在内域采用有限元法，对外域采用边界元法计算，由加权余量法的理论知这种耦合在理论上是可行的，根据此耦合方法编制了完整的计算机程序，对海洋中直立圆柱的波浪力进行了分析．数值计算的结果与理论解吻合良好，表明该方法有效．     

基于耦合模型的椭圆型沉箱基础水动力分析

椭圆型沉箱是一种新型的重力墩式码头结构形式,其水动力特性一直是目前研究的热点问题。建立了适用于椭圆沉箱的波浪时域耦合模型以及模拟椭圆型沉箱的水动力特性,该模型有外域与内域之分,外域采用三维Boussinesq方程,该方程的求解基于有限差分法,内域采用二维Laplace方程,该方程的求解基于有限元法,将两者在交界面处进行匹配。耦合数值计算结果与模型实验结果对比可看出,两者之间吻合度良好,验证了该耦合模型的精度,也为下一步的工作奠定了基础。     

Sobolev方程的半离散混合有限元法

对Sobolev方程采用混合有限元法进行数值模拟，给出了相应的半离散格式及其误差估计，构造了几组简单的低阶元.与已有文献中的有限元方法相比，该方法所采用的变分形式较简单，计算量较小，精度较高.通过对单元刚度矩阵的分析，得出在一维和二维情形下通量函数选取某些不同模式得到的关于位移的单元刚度矩阵等同     

钢管混凝土拱极限承载力计算及相关参数分析

采用考虑材料非线性的钢管混凝土拱极限承载力计算方法,对两个模型拱进行了极限承载力计算及相关参数分析.在该方法中,对钢管混凝土拱结构采用纤维单元模型,该模型假定钢管与混凝土完全粘接,钢管对核心混凝土的套箍作用体现在以一维形式表达的核心混凝土的应力-应变关系曲线之中.针对材料非线性分析中单元内各点刚度参差不齐的特点,采用单元内设小元的方法(相当于子结构),编制了非线性有限元程序.在该程序中,计算模型完全是基于小元层次进行的,比如:单元刚度矩阵由小元刚度矩阵凝聚而成,单元节点的残余力由小元节点的残余力构成,故只需改变单元内小元个数这一个参数,就可实现对结构的重新划分,极大地降低了非线性方程组的阶数,且方便实用.在程序计算结果得到模型试验结果验证的基础上,还分析了套箍系数、截面含筋率及混凝土强度对结构极限承载力的影响.     

钢管混凝土拱极限承载力计算及相关参数分析

采用考虑材料非线性的钢管混凝土拱极限承载力计算方法,对两个模型拱进行了极限承载力计算及相关参数分析.在该方法中,对钢管混凝土拱结构采用纤维单元模型,该模型假定钢管与混凝土完全粘接,钢管对核心混凝土的套箍作用体现在以一维形式表达的核心混凝土的应力-应变关系曲线之中.针对材料非线性分析中单元内各点刚度参差不齐的特点,采用单元内设小元的方法(相当于子结构),编制了非线性有限元程序.在该程序中,计算模型完全是基于小元层次进行的,比如:单元刚度矩阵由小元刚度矩阵凝聚而成,单元节点的残余力由小元节点的残余力构成,故只需改变单元内小元个数这一个参数,就可实现对结构的重新划分,极大地降低了非线性方程组的阶数,且方便实用.在程序计算结果得到模型试验结果验证的基础上,还分析了套箍系数、截面含筋率及混凝土强度对结构极限承载力的影响.     

土与结构间一种新接触单元的数值解

推导了基于非线性弹性-理想塑性本构模型的接触带单元的单元刚度矩阵，并对接触带厚度的确定、非线性迭代方法等数值计算中一些重要的细节问题进行了论述．最后对某一典型算例进行了分析，分析结果表明，所提出的接触带单元能够合理地模拟接触面上的变形机理和受力状态，是一种比较适用的模型，在土与结构相互作用分析中具有一定的应用前景。     

波浪与斜坡堤护面块体相互作用的SPH-DEM数值模拟

基于SPH方法和离散单元法(DEM)建立了二维流固耦合数值模型,模拟了波浪作用下斜坡式防波堤护面块体的水动力变化特性.流体模型采用经黎曼接触算法和CSPM修正后的SPH方法,水槽固壁和造波边界以及结构物边界采用固壁粒子法处理.离散护面块体的运动和受力采用DEM方法来模拟.流体与固体的交界面处满足界面力平衡条件.应用所建立的二维SPH-DEM耦合数值模型模拟了波浪与斜坡堤护面块体的相互作用,分析了波浪爬坡过程中的流场变化以及波压力沿胸墙的分布规律.物理模型试验验证了所建立的数学模型.     

有限元与边界元耦合问题的探讨

本文阐述有限元与边界元耦合的基本方法,并对地下工程中线弹性问题的耦合法进行了尝试。     

有限元-边界元耦合法在二维电磁场数值计算中的应用及优化措施

在二维电磁场数值分析中,有限元-边界元耦合法将有限元法与边界元法结合,可以同时发挥两种方法的优势．笔者在介绍有限元-边界元耦合法原理与算法的基础上,分别研究子域逼近有限元法、Mortar元法和并行计算与有限元-边界元耦合法结合使用分别为关心区域相对整体区域尺寸较小、求解区域中包含运动导体和大规模高自由度电磁问题的数值计算提供优化途径．     

耦合分析轧件塑性变形与轧辊弹性变形的方法和应用

发展了一种耦合分析轧件塑性变形与轧辊弹性变形(压扁)的方法。方法考虑了轧件塑性变形与轧辊弹性变形的相互影响,通过在轧件与轧辊的接触面上建立力的平衡方程实现耦合;在耦合的基础上用刚塑性有限元法分析轧件的塑性变形,用边界元法分析轧辊的弹性变形。耦合分析所得单位压力,前滑,中性角等比刚体轧辊的要低。用耦合分析方法分析轧制过程能提高分析的精度。     

弹性力学平面问题的分域虚边界元法

把所研究的弹性域分成若干个子域，然后利用基本解根据边界条件和域与域之间的交界条件建立积分方程，最后采用虚边界元技术进行数值求解．文中除讨论两域耦合情况外，还讨论了任意多域耦合的情况，并给出了形成总体矩阵的一般规律．     

复杂结构应力分析的边界元重分析方法(英文)

本文提出了一种基于边界元技术的应力重分析方法。该方法的基本思路为：采用较粗糙的有限元网格确定局部重分析模型的边界条件,然后借助计算精度较高的边界元技术对局部应力进行重分析．本文推导了重分析边界元模型边界条件的计算公式及相关方程．数值算例表明,该方法对于大型复杂结构件的局部应力分析比较有效,且具有较高的数值精度．     

有限元与边界元耦合法在波浪力计算中的应用

利用有限元与边界元耦合法对三维无界区域中直立圆柱所受的波浪力进行计算,把整个求解区域分成内域和外域两部分,在内域采用有限元法,对外域采用边界元法计算,由加权余量法的理论知这种耦合在理论上是可行的,根据此耦合方法编制了完整的计算机程序,对海洋中直立圆柱的波浪力进行了分析.数值计算的结果与理论解吻合良好,表明该方法有效.     

平面理性八节点曲边四边形有限元RCQ8

由于理性有限元用弹性力学的解插值并且在物理域内直接列式,因而等参元相比,具有物理意义明确和精度高的优点,本文采用最高至４次的平面弹性力学的多项式解对单元位移插值,推出了八节点曲边平面理性元ＲＣＱ８,若干算例表明,本单元可以避免剪切自锁和体积自锁,具有很高的精度和效率以及数值稳定性。     

Sobolev方程的半有限元方法

对Sobolev方程采用半有限元法进行数值模拟.通过将空间变量和时间变量分离,得到Sobolev方程的离散格式.首先对空间变量应用有限元方法进行离散化,得到常微分方程组的初值问题;再对时间变量应用有限差分法进行离散化,得到一系列线性方程组,求解可得到Sobolev方程的数值解.本文从理论上推导出了本文所讨论的Sobolev方程半有限元算法的矩阵算法格式,分析了其可行性.在最后给出了数值例子,从数值例子中进一步验证了半有限元方法的可行性.     

飞机偏心受载加筋板结构的有限元简化建模方法

采用偏心梁、杆板组合两种有限元模型简化方法对偏心受载加筋板结构进行建模,得到的计算结果与实体三维有限元模型的计算结果进行对比。对比结果表明,杆板组合有限元模型能有效、合理地简化偏心受载的加筋板结构。     

几何边界变化的有限元,边界元耦合重分析理论

在一种新的变分原理的基础上,提出了利用边界元和有限元耦合进行几何边界变化的重分析理论。它不仅使分析人员的工作量大大减少,而且也节约了内存,提高了重分析的精度,使大型结构的重分析大为简化。