双频自动网格法研究

传统的自动网格法常采用均匀布点法，这对应变集中区域的变形测量极为不利。节点过多，检测精度会下降，节点少了，会降低应变集中区域的空间分辨率。本文提出了双频自动网格法，在应变集中区域增加节点。文中研究了双频网格点的制作技术及对双频网格图的自动处理与分析。并应用双频自动网格法测量了中心带孔铝试件在拉伸加载过程中的孔边位移场。     

网格法及其在大变形测量中的应用

网格法具有操作简单、测量范围大、易于自动处理等特点。本文首先介绍了网格法测量的基本原理,其次用仿真的实验对网格法进行了检验,最后用皮革材料进行了拉伸实验。由于对镜头误差得到了修正、变形前后的网格点得到了自动匹配,从而实现了大变形的自动高精度测量。     

应用自动网格法测量带孔洞大样品的变形

为了能在现场检测带孔洞大样品的变形，本文提出了运用环形编码法的改进自动网格法。文中还研究了适用于大样品测量的网格制作技术，编制了网格自动处理及分析软件，进行了计算机模拟和实验标定，并分析了测试精度。     

自动网格法在轮胎橡胶力学行为测试中的应用

本文利用自动网格法光学测量技术的非接触、高精度以及可以获得全场变形信息等优点，建立了一套包含图像采集功能的载荷－变形同步测试,实验轮胎橡胶单向拉伸较大变形力学行为的测试。试验结果表明自动网格法适用于轮胎橡胶等软试件较大变形力学行为的测试，同时也验证了加和卸载过程中，胎面碳黑填充橡胶（以下简称胎面胶）具有明显的迟滞效应、Mullins效应以及残余变形；而纵横向平均应变的数据则有力的证明了在单向拉伸的较大变形范围内，胎面胶是可压缩的，且体积比随轴向拉伸应变增大而逐渐增大。     

应用自动网格法测量带孔洞大样品的变形

为了能在现场检测带孔洞大样品的变形,本文提出了运用环形编码法的改进自动网格法。文中还研究了适用于大样品测量的网格制作技术,编制了网格自动处理及分析软件,进行了计算机模拟和实验标定,并分析了测试精度。     

有限宽板孔边弹塑性变形测试

测量应变集中区域的弹塑性变形，对于研究材料损伤或微裂纹的产生，以预防宏观裂纹的产生及扩展具有重要的意义。为给航空发动机轮盘的损伤容限设计研究提供相关材料的力学实验参数，做了两个有限宽中心带孔试件的拉伸试验。试验中，采用单调逐级加、卸栽的循环加栽方式，并应用自动网格法和数字图像相关技术测量了孔用的全场位移分布。再应用最小二乘法将离散的位移分量拟合成二元多项式函数，求解出拉伸方向的应变分布，并给出孔边应力σ与孔边应变εy的关系曲线。     

超细网格法的频域分析技术

本文把通讯理论中的"系统"的观点和数学上的傅立叶分析方法引入网格法的变形分析中,形成了一种全场自动处理的网格图像分析新方法,它具有直观性好和灵敏度可由网格栅线密度控制的特点,为观察细观的实验现象和测量微区域内的位移场和应变场提供了一种有效的实验手段.     

无网格法研究进展及其应用

从加权残量法的角度出发，系统地总结了现有各种无网格法的基本格式，阐明了无网格法的特点，论述了无网格法的研究进展，给出了无网格法在碰撞、动态裂纹扩展、金属加工成型、流体力学以及其它领域中的应用。     

用记忆材料模拟流变大变形问题

本文采用具有记忆特性的有机玻璃材料结合网格法测量技术，以软岩巷道围岩流变形为例展开了有限变形模拟试验的研究工作，并以有限变形理论为基础结合现代计算机技术进行了实验数据处理。实例分析表明，所用方法和技术具有广泛的灵活性、可靠性的实用性。     

无网格法及其最新进展

无网格法具有许多独特的优点,因此有人认为无网格法将成为继有限元法之后新一代的数值方法. 到目前为止,已提出了不下十几种无网格法, 这些无网格法各有不同的优缺点.本文评述几种主要无网格法(着重应用于固体力学中的基于弱式的无网格法),将它们适当地分类, 论述典型的无网格形状函数的构造方法,介绍无网格法的发展现状, 评价已提出无网格法的优缺点,并比较典型无网格法的相同和不同点, 以及讨论无网格法发展所面临的问题等.最后就无网格法的发展趋势进行了展望.      

THE PROGRESS OF THE OVERLAPPING GRID TECHNIQUES

从网格装配和插值计算两个主要方面对现有的重叠网格方法进行了综述。首先,从挖洞方法和建立嵌入网格关系环节的寻点技术出发归纳和介绍了网格装配方法;其次,介绍了数值迭代过程中的插值计算方法,并特别讨论了插值守恒性以及插值计算精度等问题;另外,对重叠网格方法的并行计算和应用成果也作了介绍;最后,通过总结认为重叠网格方法在改进网格装配方法、改善插值和并行计算效率等方面仍需进一步研究。     

7阶WENO-S格式的计算效率研究

WENO-S格式是一类适合于含间断问题数值模拟的加权本质无振荡格式.这类格式的光滑因子满足对单频波为常数,这使得其近似色散关系与线性基底格式一致,并且具有良好的小尺度波动模拟能力．计算效率是数值方法性能指标的一个重要方面.由于WENO-S格式的光滑因子在各子模板上的计算公式除下标不同外形式一致,在计算线性对流方程相邻数值通量时,部分光滑因子完全相同．为此提出一种消除WENO-S格式冗余光滑因子计算的方法.该方法要求一条网格线上用于重构或插值的量可以表示为一个序列.基于此要求分析其对于几种不同物理问题的可行性和使用方法.以7阶WENO-S格式为例介绍了格式性质和去除冗余光滑因子计算的方法.该方法中预先计算和存储一条网格线上的所有光滑因子,在网格点较多的情况下,光滑因子计算次数约为原7阶WENO-S格式的1/4.对一维对流问题、球面波传播问题、二维旋转问题、二维小扰动传播问题及一维和二维无黏流动问题进行了数值模拟.结果表明该格式对多种流动结构具有良好的捕捉能力,并且同时具有良好的计算效率,去除冗余计算后又降低了约20%的计算时间.     

OVERLAPPING GRIDS FOR THE SIMULATION OF PARTICLE INTERACTION AT THE MICRO SCALE

This paper describes the use of overlapping grids for the calculation of flow around single and multipleparticle configurations at the micro scale. The basic equations for calculation are those for conservation of mass and momentum which are solved using a common Finite-Volume formulation. The hydrodynamic particle-particle and particle-wall interaction can be calculated by using an overlapping or Chimera grid scheme. With the grid structuring procedure it is possible to use simple and structured grids around the particles and the overall main grid geometry. The particle grids are lapped over the main grid such that they can move independently after each time step without remeshing the whole geometry. The paper gives results for the validation of the code developed for general test cases, for a rotating ellipsoid in simple shear flow, the flow around particles attached to a wall, the motion of a particle in the vicinity of a wall and some results for the flow through a packed bed configuration.     

纳米级变形的扫描隧道显微镜测量研究

提出了应用扫描隧道显微镜及利用物质原子结构测量纳米级变形的新方法，对高定向裂解石墨、单晶硅材料的纳米级变形进行了测量研究，得到了强激光辐照石墨的残余应变场，对拨原子后单晶硅形成的纳米裂纹附近区域的残余变形进行了分析     

动网格生成技术及非定常计算方法进展综述

对应用于飞行器非定常运动的数值计算方法(包括动态网格技术和相应的数值离散格式)进行了综述.根据网格拓扑结构的不同,重点论述了基于结构网格的非定常计算方法和基于非结构/混合网格的非定常计算方法,比较了各种方法的优缺点.在基于结构网格的非定常计算方法中,重点介绍了刚性运动网格技术、超限插值动态网格技术、重叠动网格技术、滑移动网格技术等动态结构网格生成方法,同时介绍了惯性系和非惯性系下的控制方程,讨论了非定常时间离散方法、动网格计算的几何守恒律等问题.在基于非结构/混合网格的非定常计算方法中,重点介绍了重叠非结构动网格技术、重构非结构动网格技术、变形非结构动网格技术以及变形/重构耦合动态混合网格技术等方法,以及相应的计算格式,包括非定常时间离散、几何守恒律计算方法、可压缩和不可压缩非定常流动的计算方法、各种加速收敛技术等.在介绍国内外进展的同时,介绍了作者在动态混合网格生成技术和相应的非定常方法方面的研究与应用工作.     

线化欧拉方程的高阶间断有限元数值解法研究

采用高阶间断有限元法于非结构网格上针对复杂外形数值求解声学控制方程------线化欧拉方程. 背景流场采用有限体积法于结构网格求得, 一种高精度数据传递方法将基于有限体积法的背景流场数据传递到声场计算所采用的较为稀疏的非结构网格上, 保证了背景流场信息的完整和精确. 为提高计算效率, 采用了一种更为直接的Quadrature-FreeImplementation技术以及网格分区并行技术. 数值结果表明采用高阶的情况下即使在稀疏的网格上也可以捕捉到细微的声场结构.      

THE INVERSE ANALYSIS OF INTERNAL DEFECTS (INCLUSIONS) IN STRUCTURES USING A COMBINED XFEM AND GA METHOD

建立扩展有限元法与遗传算法相结合的结构缺陷反演分析模型.扩展有限元法通过引入不连续位移模式使得网格剖分无需依赖结构内部的不连续界面.通过改变水平集函数表征结构缺陷（夹杂）的位置和大小，遗传算法在每次迭代过程中具有全局和局部搜索能力，通过评估响应测点的响应量适应度值决定是否进一步迭代.对带有单个圆形缺陷（夹杂）和多个缺陷（夹杂）的结构进行了反演分析，并就响应测点的布置进行了讨论.结果表明，建立的反演分析模型能准确地探测结构存在的单个甚至多个缺陷（夹杂）.     

水下发射航行体尾涡不稳定性分析

在水下连续发射过程中前一发航行体尾流会对后一发航行体运动姿态稳定性产生流动干扰现象. 因此, 研究尾流中涡旋结构演变机理对解决多弹体水下连续发射流动干扰难题具有重要的意义. 本文采用改进型分离涡模型与能量方程, VOF多相流模型与重叠网格技术相结合方法, 对航行体水下发射尾流演变过程开展精细化模拟研究, 其中模拟结果和实验吻合度较好, 验证了本文数值方法的有效性. 以航行体尾流区域为重点研究对象, 分析了尾流区瞬态流场分布, 讨论了横流强度和雷诺数对尾涡结构演变以及脉动压力分布特性的影响. 结果表明: 由于尾流区高速流体核心区与低速自由流相互作用导致Kelvin-Helmholtz不稳定现象出现, 可以清晰地发现涡旋结构在剪切力的作用下发生脱落. 在横流条件下, 航行体尾端脱落的涡环与涡腿形成发卡涡, 而多个发卡涡沿轴向间隔排列组成发卡涡包存在于尾流中. 随着横流强度增大, 形成多级发卡涡包结构, 而导致脉动压力二次峰值均出现的主要原因是尾流涡旋流场演变引起的. 随着雷诺数的增大, 尾流中由圆柱形涡和U型涡组成的二次涡结构逐渐明显, 不稳定性加强.      

THE NUMERICAL STUDY OF LATTICE BOLTZMANN METHOD BASED ON DIFFERENT GRID STRUCTURE

传统的格子波尔兹曼方法(lattice-Boltzmann method, LBM)通常基于标准均匀网格, 这主要取决于速度的空 间离散格式.均匀网格结构的特点, 使LBM在处理具有复杂边界的问题时遇到较大的困难, 从而限制了它的应用.另外, 对于较为复杂的流动, 其流场存在流动变化剧烈和平缓的区域, 在流动变化剧烈的区域, 往往需要足够的网格点才能更好地捕捉到流场信息, 而均匀网格会使得网格数量过多, 这会增加计算量, 但网格数量过少又无法获得必要的流场信息, 使LBM的计算效率降低.为了解决上述问题, 用不同的网格结构, 以顶盖驱动的腔体内流、柱体绕流和翼型绕流为例, 探讨了提高LBM算法的计算效率和适用性问题.     

裂缝在非均匀岩层内扩展机理研究

薄差储层有效开发的关键问题是明确水力裂缝在砂/泥岩界面上的扩展特性. 为了分析裂缝在砂/泥岩界面上的扩展机理, 以ABAQUS 为平台, 利用扩展有限元框架下的分离裂缝模型模拟了裂缝在砂/泥岩界面上的扩展过程, 该模型提出预设虚节点法, 通过在位移的插值函数中引入附加增强项, 实现裂缝扩展独立于网格边界, 获得了裂缝通过砂/泥岩界面的扩展规律. 并采用白光散斑实验对低渗透储层砂/泥岩界面的裂缝扩展进行了实时跟踪, 裂缝扩展过程与数值模拟过程吻合较好, 说明扩展有限元法是定量分析裂缝扩展的有效手段, 同时分析了裂缝穿过界面扩展的影响因素.