非结构网格中子输运方程的离散纵标节块数值解法

基于离散节块方法（DNTM）的思想，研究了一种新的适用于二维三角形节块的离散坐标数值解法．使用坐标变换将任意三角形变换为正三角形，源项和通量的空间分布使用二次正交多项式展开，横向泄漏项的空间二次分布通过使用二元二次多项式近似正三角形节块的通量分布来获得、针对小尺寸易出现不收敛的情况，提出了一种改进的节块平衡有限差分方法．数值计算结果表明，该方法能满足非结构几何区域问题的求解需求，与细网SN方法相比，该方法在较粗的网格下具有高的精度，同时还适用于较细网格下的中子输运计算．     

非线性节块展开法

非线性节块展开法是一种快速计算反应堆物理参数的节块方法。与有限差分法相比，它具有速度快、精度高的特点；与另一种节块法—Ｇｒｅｅｎ函数展开法相比，它又保留了有限差分法的优点：不但可以得到有效倍增因子、功率分布等参数，还可以得到差分方程的系数矩阵，因而十分适合与广义微扰理论（ＧＰＴ）相结合以获得更快的计算速度，应用于堆芯燃料管理优化中。通过两个二维基准问题的检验计算，比较了非线性节块展开法和Ｇｒｅｅｎ函数展开法程序的速度与精度，结果表明两个程序基本相当，非线性节块展开法可以满足实际应用的需要。     

基于有限元法的二维裂纹应力强度因子研究

基于有限元分析方法，对有限大平板中存在的中心穿透裂纹，分别用不同的方法分析其裂纹尖端应力、应变场分布，计算出裂纹尖端的应力强度因子．通过对求得的应力强度因子值与解析解的比较，表明用有限元方法计算应力强度因子具有相当高的精度，并且操作简便。     

Cauchy微分方程的一种新的数值解法

众所周知，求解微分方程（组）常用的数值方法有有限差分法，有限元素法等，这些方法都是将微分方程（组）分离散化后求解．若将网格划分得粗了，则求解精度不高，不能满足工程实际需要，若将网格划分得细了，则所需计算机内存量和计算量都太大．为解决上述问题，本文给出微分方程（组）的解的概率表达式的一种新的数值解法──概率数值解法．     

基于序列转换的p型有限元后误差估计

提出了基于序列转换的ｐ型有限元后误差估计外推算法，当问题的解是光滑的时候，可采用ε算法对能量模进行估计；当问题的解是非光滑的时候采用多项式外推算法，这种情况下，如果采用均匀或近似均匀网格，用ｈ２或（ｈ２，ｈ４）外推，如奇异点附近采用强梯度网格，则ｈ３或（ｈ３，ｈ６）外推更加合适．计算结果表明所提方法有好的计算精度和较高的可靠性     

渤海湾10年来化学耗氧量的变化趋势与污染风险

基于1995年到2005年渤海湾化学耗氧量（COD）的监测数据序列，研究了渤海湾10年来COD的季节（5月、8月及10月）及年际间的变化趋势，运用极值理论研究了COD的极值分布规律，计算了COD分别超过不同污染水平的风险概率．研究结果表明：渤海湾COD季节之间变化不明显，5月、8月及10月的COD具有指数分布特征；COD年际间变化呈下降趋势，与同期COD排放变化趋势吻合；COD的污染风险概率符合Pareto型超阈值分布函数．     

通量展开节块法求解六角形几何三维多群中子扩散方程

提出了一种在三维六角形几何节块内数值求解多群中子扩散方程的节块法，该方法把节块内各群中子通量分布用解析基函数近似展开．为了改善节块耦合关系，采用了一种新的节块边界条件：面平均偏流0次矩和1次矩同时保持连续．将响应矩阵技术应用于迭代求解过程，使得该方法具有较高的计算效率．通过对基准问题的校验计算表明，该方法能准确地给出有效增值系数以及节块功率；对于二维多群问题，所有基准题的组件最大功率偏差均小于1％．     

改进的混合界面动态子结构法及其在车身结构动力分析中的应用

本文针对现行的混合界面动态子结构法进行改进，提高了模态综合效率和和精度．将改进的方法运用到承载式车身结构（ＳＶＷ７１８１）动力分析中，采用空间梁单元和板壳单元进行车身结构有限元建模，对车身结构进行了子结构划分计算．计算结果表明，该法计算量小，综合效率高，避免了以往方法计算量大、精度不够高、对计算机内存要求高的缺陷．而且为该型车的动态性能评估和进一步研究定型提供了理论计算依据，表明该方法具有实用推广价值和应用前景．     

计算中子扩散方程的全堆芯格林函数快速节块方法

导出了全堆芯格林函数快速节块方法理论模型，提出了对各个节块逐次扰动进行正向扩散计算建立格林函数库的思想，研制相应软件ＣＧＦＭ，对一系列基准问题和我国泰山核电厂的问题进行了校核，计算结果表明，本文方法在保证一定的计算精度下，计算速度比先进节块方法ＮＧＦＭ快１０倍左右。     

结构系统中疲劳失效途径失效概率计算的二阶可靠性方法

针对结构系统疲劳可靠性分析中失效途径极限状态函数的非线性性质,提出用二阶可靠性方法（ＳＯＲＭ）计算疲劳失效途径的失效概率．文中介绍了疲劳失效途径的可靠性模型及结构系统二阶可靠性方法的基本原理,并通过例题说明了计算过程．计算结果表明,采用二阶可靠性方法可使失效概率的计算精度得到显著的改善．     

一种三角形表面网格质量综合度量方法

在计算机图形学中,有限元等科学计算与分析方法得到了越来越广泛的应用．而有限元单元网格质量的好坏将直接影响有限元计算过程的收敛性以及有限元解的误差．通过对计算机图形学中常用的三个表面网格模型进行分析,利用模糊数学的方法,给出了一种三角形网格质量综合度量方法．利用该方法得到的评价指标可以综合、快速地评价一个有限元网格模型的质量．     

一种消除线性独立性问题的单位分解三角形单元

为建立高效的消除线性独立性问题的单位分解有限单元公式,采用多项式基函数用来作局部近似。普通的线性三角形形函数作单位分解函数,提出有限元无网格耦合三角形单元。该三角形单元形函数具有Kronecker delta性质,能够直接施加位移边界条件。数值算例表明,该三角形单元能够消除普通单位分解有限元的线性相关问题,并且具有较高的计算精度,结果优于普通线性三角形单元和线性四边形等参元。     

基于BB型对偶剖分的抛物方程有限体积元法

讨论了抛物方程基于三角形剖分和BB型对偶剖分的有限体积元法, 给出半离散及全离散有限体积元格式的最佳阶L²和H¹误 差估计.     

有限元逆算法及其在板料成形工艺优化中的应用

采用三角形单元开发了精度较高的有限元逆算法,算法中考虑了摩擦、压边力、拉深筋、背压力、曲压料面等实际工艺条件．方程组求解引入快速求解算法,比一般的LDLT解法快几倍,尤其对于大型复杂冲压件问题来说求解效率提高更多．用此算法对典型的L型件进行了模拟,并与增量算法和实验结果进行了比较．可以看出该算法已经达到很高的精度,而且模拟过程迭代收敛速度快,计算时间很短．在此基础上以工件的厚向应变作为目标函数,采用一种实用性很强的工艺参数快速优化算法,即黄金分割法对压边力进行优化．仍以L型件为例,优化过程的计算时间说明快速优化算法相对传统增量算法具有明显的优势．     

一种新的矩形网格生成等值线算法

提出了一种利用数据关联表生成矩形网格等值线的算法。该算法的优点是计算效率高，避免了以往等值线追踪算法起始点选取困难、网格出口边判断复杂的问题。算法的追踪结果精度取决于等值点的计算方法，其精度是可控的。这种算法原理可以扩展应用到三角网格的等值线方法中。     

一种基于代数图论的有限元模型节点排序方法

提出了一种基于有限元模型中节点自由度构造赋权单元团图的方法,根据代数图论的理论,应用赋权单元团图的拉普拉斯矩阵的Fiedler向量,对有限元模型的节点进行排序,以达到减少结构刚度矩阵的半带宽和外形的目的,该方法不但能适用于一般有限元模型,而且适用于包含不同类型单元、具有不同自由度节点的混合节点模型,对于混合节点模型,该方法比基于单元团图的拉普拉斯矩阵的代数图论方法能够取得更加满意的结果,据此编制的前处理程序,可以对任意编号的模型进行优化处理,数值算例结果表明本方法是有效的.     

弹性波动方程数值解的有限元并行算法

在求解弹性波动方程中，有限元法的高内存量和巨大运算量的需求在基于单CPU串行算法中一直难于满足，制约其优势的发挥。根据有限元法的“化整为零、集零为整”的基本思想与并行处理技术的“分而治之”的原则基本一致，采用基于多CPU的并行算法，从有限元参数矩阵计算和线性方程组求解两个方面人手，把求解区域分到多个CPU上并行计算参数矩阵，对线性方程组采用循环块三对角线方程组进行并行求解。对比了不同大小空间和不同CPU个数下的加速比，证实了多CPU的并行算法能够克服基于单CPU串行算法的物理限制，满足了有限元法的巨大空间量和运算量的需求。此算法具有理论上的正确性和实践上的可行性。     

液化气船液罐鞍座温度场的有限元计算

本文采有四边形９节点等参单元有限元法编制了一通用程序,它适用于任何二维“场”问题分析．对任意形状、不同边界条件,它都能较精确地进行处理,具有网格自动生成和编号的功能,能打印出每个节点的温度值,并画出各种鞍座的等温线簇．它的精度远高于常规的三角形单元,在３８６微机上的运行时间在５ｍｉｎ左右,一次最多可计算５３０个点．     

抛物问题的一类变网格矩形有限元方法

传统有限元方法要求区域剖分满足正则性条件.对一类抛物问题,首先可利用相应的矩形元,绕开区域剖分中正则性条件的限制,再结合变网格思想,导出所讨论问题的全离散有限元格式;其次,利用Riesz投影算子,通过一些新的技巧和方法,得到与传统有限元相同的最优误差估计.     

非稳态方程的一类Crouzeix-Raviart型有限元方法

传统有限元方法要求区域剖分满足正则性条件.对一类非稳态方程,首先利用相应的Crouzeix-Raviart型元,绕开区域剖分中正则性条件的限制,导出所讨论问题的全离散格式;其次,利用Riesz投影算子,通过一些新的技巧和方法,得到最优误差估计.