混合元法及其在本质边界条件处理中的应用

为了解决无网格Galerkin(EFG)方法中本质边界条件的施加问题,提出了一种耦合径向基点插值无网格(RPIM)方法与EFG方法的混合元法.该方法将问题域分成两部分,在包含本质边界部分采用RPIM方法处理,其余部分采用EFG方法,然后利用两种无网格插值方法之间的关系进行耦合.数值算例结果表明:文中方法在简明有效地处理本质边界条件的同时,具有与EFG方法相当的计算精度,且易于编程计算.     

方形同轴电位数值计算的MATLAB实现

二维静态场对一些特殊对称的边界才能求它的解析解，而许多实际问题由于边界条件复杂需要借助数值法求场域的电位分布。文章以方形同轴线静态场域边值问题为例，论述了有限差分法和有限单元法数值计算的过程中，用MATLAB实现是一种行之有效的方法。     

SolvingPolyzoneProblemsofIndividualMediaUsingBoundaryElementMethod

文中在常规间接边界单元法的基础上，建立了多区域耦合边界单元法，以便求解矿山岩体中普遍存在的分区介质问题。用常规边界单元法解多域问题，要给定内外边界条件，但是，内部边界条件确定是困难的。本文试图用多域耦合边界单元法解决这一困难，即只需给定计算区域最外围边界条件，而各子域之间的边界条件借助于协调关系确定，并且各子域可用不同的材料参数     

径向基点插值无网格法与有限元耦合法

为了发挥无网格法和有限元法各自的优势,提出了径向基点插值无网格法与有限元直接耦合的计算方法.比较了点插值无网格方法与无网格Galerkin法(element-free Galerkin method, EFG)的优缺点:无网格法只需要结点信息,无需单元信息,克服了有限元计算中网格畸变和重新生成带来的困难,故其在分析裂纹扩展和局部大变形等问题方面具有优势.对悬臂梁受集中荷载问题进行了计算,两种计算结果一致,但是本算法效率较高.用本方法与FEM耦合方法对半无限平面受集中荷载问题进行了计算,结果表明:本文的耦合方法可行,计算效率较高.     

基于自组织规则的热力耦合解析

为解决平面微小结构设计问题中的可变边界条件、多物理场耦合效应等解析难题,利用平面微小结构连续性单元同桁架单元应变能相等的机理,提出了基于元胞自动机的热力耦合解析方法.以发光二极管(LED)基板热力耦合现象为研究对象,通过LED基板温度、位移及热应力耦合的数值结果与有限元解析结果的比较,证明了方法的有效性.适用于平面微小结构的热力耦合规则的提出,为复杂多物理场条件下的微结构工程设计提供了一个具有灵活性的性能解析途径.     

二维Biot固结方程的自然单元法求解

为了使自然单元法能够应用于土体等多孔介质的流固耦合计算,通过结合Biot固结理论及自然单元法自身特点,利用经典变分原理推导了固结微分方程的离散形式,并针对二维问题编制了相应的计算程序.算例结果表明,自然单元法的结果与解析解吻合良好,其精度高于有限单元法.从而验证了自然单元法在固结分析中的正确性,拓展了自然单元法的适用范围.     

无缝钢管纵连轧无网格法热力耦合分析

采用混合交换法修正无网格移动最小二乘近似函数,以便直接施加边界条件;采用伽辽金法构建速度场刚度矩阵;采用配点法构建温度场刚度矩阵;采用间接耦合法将速度场与温度场耦合求解;最终推导出了刚塑性无网格法热力耦合计算公式,实现了纵连轧过程的热力耦合模拟。通过数值模拟与实验研究对比表明,模拟结果与实验结果吻合度较高,充分验证了本文给出的无网格法的可靠性和正确性。     

杆的加模态式有限单元法

为了提高杆单元精度和降低自由度,将静态有限单元线性特征值问题和动态有限单元高精度的优点有机的结合起来形成一种高精度的动力有限单元法,将杆元固有振动方程分解为满足单元边界条件的平衡方程和固支边界条件的振动方程两部分,前者的解即静态元的位移函数,后者的解为为本方法添加部分,由此得出的有限单元法只对静态元刚度阵和质量阵做添加式修改,求解过程不变,理论和算例均表明其有效性。     

对二维温度场敞开边界计算问题的研究

对二维温度场敞开边界计算问题的有限单元处理方法进行了探讨。提出了扩大边界法。这种方法有效地解决了求解无限域温度场时边界的处理问题。解决了利用常规有限单元法无法克服的困难,具有很大的实用价值。     

拉格朗日无网格粒子法的声学计算及边界实现

针对现有拉格朗日方法处理声学问题时声学边界出现的匮乏现状,提出了耦合无网格-时域有限差分边界处理技术.边界粒子采用修正光滑粒子动力学方法(CSPM)求解,虚拟粒子采用时域有限差分方法(FDTD)求解;随后分别建立了声学软边界、硬边界和吸声边界的CSPM/FDTD实现方法,其中声传播控制方程采用拉格朗日形式的小振幅波动方程.通过管路声传播模型来对三种不同边界条件进行了验证,而后基于CSPM方法,着重对软边界弹性体在平面波中的时域声散射问题进行了模拟.计算结果表明:该方法能较好模拟出散射声场及其随时间的变化,同时不同初始粒子间距的计算结果也表明了其具有收敛性,验证了修正光滑粒子动力学方法和耦合声学边界处理技术的可行性.     

应用FDM与AWE技术快速求解导体目标宽带RCS

渐近波形估计(AWE)技术是分析目标宽带有效数值方法之一,但该方法需要多次存储高阶频率导数阻抗矩阵,内存消耗大.本文将快速偶极子法与AWE技术相结合,只需存储近区场阻抗矩阵及其高阶频率导数阻抗矩阵,并且大大加速了在迭代求解过程中的矩阵矢量乘积运算.与传统AWE技术相比,计算时间和内存消耗都得到了有效缩减,数值结果证明了本方法的有效性和精确性.     

基于最小二乘拟合和特征基函数法的目标宽带RCS快速计算

提出一种快速分析目标宽带雷达散射截面的方法,该方法将最小二乘拟合与特征基函数法相结合,通过计算选定的若干频率点的表面电流便可快速求解出整个频带内的表面电流.具体过程为利用特征基函数法求解选定频率点目标表面电流,进而利用最小二乘拟合实现表面电流和雷达散射截面的快速计算.数值计算结果表明:在不影响精度的前提下,该方法可大大提高计算效率、减少内存需求.     

铝型材挤压有限元和有限体积对比模拟

研究了有限元法(FEM)和有限体积法(FVM)模拟金属塑性成形的基本理论和求解方法．针对铝合金门窗型材制品截面形状十分复杂、壁薄而成形坯料厚的具体特点，分别采用FEM和FVM对其挤压成形过程进行了对比模拟．结果表明，FEM计算效率高，但在模拟终成形阶段时，由于频繁的网格再划分，造成模拟结果严重失真；FVM避免了网格再划分问题，模拟精度高，但计算时间长．提出了采用FEM模拟铝型材挤压预成形，FVM模拟终成形的复合模拟方法，解决了料厚而壁薄这类铝型材挤压成形的数值模拟问题。     

微分求积方法及其在力学应用中的若干新进展

简单地介绍了微分求积法（DQM）的基本原理及发展、影响微分求积方法数值解精度的因素、与有限差分方法和有限元方法比较而言DQM的优点,同时重点介绍了近年来DQM在固体力学和流体力学相关领域应用中的若干新进展和展望.     

Sramek公式及修正式计算每搏量与导纳法的比较

本文应用阻抗—导纳式心输出量机对40名正常成年人(男、女各半)进行每搏量(SV)的测定,并把Sramek公式及Sramek修正式和导纳法3种不同方法测得的SV进行比较。结果表明,应用Sramek公式和Sramek修正式测得的SV均大于标准对照值,而导纳法的SV值在标准对照值范围内。提示导纳法测得的SV值可能更为可靠。     

三维快速多极虚边界元配点法

以三维弹性力学问题为研究背景,提出了一种三维快速多极虚边界元配点法的求解思想,即将三维快速多极展开的基本思想和广义极小残值法运用于求解传统虚边界元配点法方程.文中将三维弹性问题的基本解推导为适合于虚边界元快速多极算法的展开格式,经数值计算格式的演变,使求解方程的计算量和储存量与所求问题的计算自由度数成线性比例,以达到数值模拟大规模自由度问题的目的.算例说明了该方法的可行性、计算效率和计算精度.此外,该方法的思想具有一般性,应用上具有扩展性.     

解Helmholtz方程的共轭梯度法三角阵预处理器

利用对应微分运算的差分运算矩阵,提出了差分求解Helmholtz方程的三角阵预处理器,阐述了在算子离散的过程中而非离散后构造高效率预处理器的基本思想.利用二阶频域Mur吸收边界条件下的二维导体柱散射的模型问题,结合法方程最小余量预处理共轭梯度法(PCGNR)验证了该预处理器的有效性.数值结果表明了该预处理器能够在网格精度提高和吸收边界趋远时,相对常规共轭梯度法具有降低计算复杂度的效果,而存储复杂度并没有提高.同时,也揭示了二维散射问题在网格精度与吸收边界距离一定的情况下,用不同共轭梯度法求解时,由于迭代次数变化较少,计算量几乎随未知量线性增长.     

六西格玛管理方法的应用

详细介绍了六西格玛管理方法的主要内容与一般推进过程,阐述了成功应用六西格玛管理方法的关键因素。     

时间序列非平稳检测方法的对比分析

对比分析了检测时间序列是否平稳的两种方法: 一种是基于简单的统计方法检测时间序列非平稳性的run-test (简称RT)方法(考虑了3种典型方法: 均值RT、方差RT和标准差RT), 另一种是利用动力非平稳检测的图示方法(space time-index, STI)。 通过分析人工生成的已知平稳性的时间序列和实际湍流时间序列, 对比了两种方法在时间序列非平稳检测上各自的优劣: RT方法简单、高效, 但精确度低; STI方法精确度高, 但复杂、低效。同时, 分析了3种不同的RT方法对不同类型时间序列的平稳性检测能力, 表明方差RT方法最适用于湍流时间序列非平稳性检测。     

电磁脉冲对数字电路的辐照法与注入法比较

为了研究电磁脉冲注入法在评价设备或系统中电子器件受电磁脉冲辐照作用的有效性,准确评价电子器件在工作时的抗电磁脉冲能力,选取了某导弹装备上信息处理模块的串行数据分配电路为研究对象,研究了2种雷电电磁脉冲对电路的辐照效应与注入效应,总结了耦合规律与作用机理,讨论了直接注入与辐射耦合的相关性,探讨了电磁脉冲抗扰度的实验方法.研究结果表明,电磁脉冲的注入法与辐照法是不相关的,注入法不能替代辐照法.