【摘 要】
：
验证与确认是评估数值模拟精确性和可信度的主要技术手段。本文着眼于验证环节中的解验证，探讨了如何对数值模拟中离散误差造成的性能参数不确定度进行后验估计。基于Richardson 外推公式，本文研究了几种收敛解的估算方法，尝试使用收敛解的估计值量化性能参数的不确定度，为可靠性评估提供输入信息。
【机 构】
：
北京应用物理与计算数学研究所,北京100094
【出 处】
：
中国力学大会2011暨钱学森诞辰100周年纪念大会
论文部分内容阅读
验证与确认是评估数值模拟精确性和可信度的主要技术手段。本文着眼于验证环节中的解验证，探讨了如何对数值模拟中离散误差造成的性能参数不确定度进行后验估计。基于Richardson 外推公式，本文研究了几种收敛解的估算方法，尝试使用收敛解的估计值量化性能参数的不确定度，为可靠性评估提供输入信息。
其他文献
研究用pFFT快速边界元方法模拟声散射问题的关键技术。采用Burton-Miller方程消除外问题解的不唯一现象。为此，文中研究了采用常量元时该模型中超奇异积分的计算方法。通过平面声波的刚性圆球声散射的数值模拟验证了所建立的声学pFFT 快速边界元方法。
光滑粒子流体动力学方法(SPH )是一种基于核估计的无网格拉格朗日数值方法.DSPH (非连续的SPH )是相对于传统的SPH 方法提出来的一种可用于模拟非连续物理现象的方法，并且能够解决边界缺陷问题，有很好的应用前景.本文研究了一维情况下界面处函数阶跃大小，函数变化率及粒子分布的均匀性对DSPH 模拟结果的影响.对函数及其导数估计式的分析及ID情况下不连续函数的数值计算均显示了DSPH方法的优越
One of major difficulties in the implementation of meshfree methods using the moving least square (MLS)approximation is the imposition of essential boundary conditions as the approximations do not pas
等几何分析方法(IGA-Isogeometric Analysis )采用NURBS 技术实现了计算机辅助几何设计(CAGD )、有限元分析(FEA )和结构优化的无缝结合。在此框架下，本文采用等几何分析方法对带孔结构进行了结构分析与孔形优化设计研究。在研究等几何分析的基本原理和方法的基础上，给出了多片NURBS 模型的等几何分析和含NURBS 修剪技术的等几何分析的一般步骤，并采用这两种分析方法
自然单元法主要是基于给定结点的Voronoi 图，利用自然相邻插值进行形函数的构造，其形函数满足Kronecker delta 性质，便于施加本质边界条件，这使得自然单元法同时兼有有限单元法和无网格法的优点。在小应变假设的条件下，利用几何非线性的应变-位移关系，推导了考虑材料几何非线性问题的自然单元法模型。算例表明，该模型在处理几何非线性问题时，具有一定的合理性和可行性，是一种有效的数值方法。
本文从高玉臣院士提出的基面力概念出发,针对一种基于基面力概念的新型有限元法―基面力元法(BaseForce Element Method, BFEM )的研究及应用进行综述。分析结果表明：基面力元法的计算所得结果与理论解吻合较好；基于余能原理的基面力元法较传统平面4节点等参元(Q4 单元)具有较高的计算性能；在高度几何非线性分析中,基于余能原理的基面力元法较平面4 节点减缩积分单元(Q4R 单元)
自适应无网格法是针对有限元法无法求解或者不易求解的复杂问题，利用无网格法结点排布灵活、易于增删结点、便于自适应分析等优点发展起来的。在对自适应方法理论基础和发展进行总结基础上，采用基于应变能梯度的自适应无网格-有限元耦合方法，对HA 涂层等离子喷涂制备的生物涂层材料的无摩擦接触问题进行分析、计算精度和效率较好。
本文研究了针对频域有限元直接动态分析的面向目标误差估计以及误差范围估计计算方法。面向目标的误差估计方法就是专门针对如何准确和经济地估算特定值误差的一种方法，利用原问题的共轭偶问题进行计算。频域有限元的直接动态分析是模拟频域扫描实验的一种计算方法，专门针对谐振激励的线性动态响应问题，利用将原自由度分解为实部和虚部描述频率的变化，从而计算变形体的动态响应。本文利用扩展针对有限元的面向目标误差估计的自由
以典型螺栓联接结构为实例，讨论了有限元模型确认的技术思路。选取两块钢板组成的螺栓联接结构作为应用实例，初步提出了有限元模型确认流程。针对有限元模型确认过程中存在的若干问题，如目标结构中不确定性的来源，分层确认，有限元模型参数不确定性的非参数核密度估计，模型确认准则和代理模型的建立等，进行了讨论。
在升阶谱有限元法中，低阶升阶谱元的自由度是高阶升阶谱元的自由度的一个子集，这为节省计算量和编制自适应分析程序提供了极为有利的条件。但升阶谱有限元要用到高阶甚至很高阶多项式作为附加自由度的基底函数，从而将出现数值稳定性问题。诸德超已解决了一维问题和正规域内的二维或三维问题，但对于不规则域的问题目前还没有解决。本文在单元边界上采用满足C0或C1连续的形函数，在单元内部采用正交升阶谱形函数，用微分求积原