基于逆向工程的三维复杂地质体精细建模及ADINA前处理在FLAC3D 建模中的应用

逆向工程技术作为一门新兴学科在复杂模型建模方面具有独特优势,并在许多领域得到了广泛的应用,在此基础上本文首次提出了基于逆向工程的三维复杂地质体精细建模技术,使得建立模型达到既“可视”又“可算”的目的。利用现行大型通用有限元仿真平台(ABAQUS,AD INA,ANSYS等)良好的网格划分功能且具有与CAD /CAM软件有着良好的接口这一优点,对模型进行网格划分,然后转换为FLAC^3D能够接受的格式,从而克服了FLAC^3D在前处理建模方面的不足。以AD INA为例对相应的接口程序开发过程进行了阐述,将上述技术方法应用于工程实例,取得了很好的效果。     

一类复杂建筑外围装饰结构的风荷载参数研究

在建筑外围布置复杂装饰结构能够进一步提高建筑的视觉效果,但复杂装饰结构对风荷载的敏感性往往要超过主体建筑,因此对建筑结构的抗风设计安全性提出了更高的要求.建筑外围装饰结构由于其体量小数量大而导致直接进行物理风洞试验变得几乎不可能,因此本文研究利用数值风洞技术来解决外围装饰结构的风荷载问题.结合高性能的硬件计算平台,通过设置合理的网格剖分形式,解决了流苏结构多个移动位置下的三维空间高质量网格的生成难题.数值模拟了各层流苏在不同风向角下的风荷载体型系数分布规律,并统计分析了所有风向角下流苏结构体型系数极值分布,给出了设计所需的各层流苏极值风参数取值范围.数值模拟结果显示出了流苏结构风荷载的极值区域位于建筑外形的拐角区域,当结构面临可能的大风天气时,流苏结构的停止位置宜避免将长条的流苏杆件置于上述区域,以提高流苏结构的抗风安全性.     

大变形问题分析的局部Petrov-Galerkin法

在微机电系统(MEMS)的建模和模拟研究中,大变形或大移动要充分予以考虑.用有限元法分析这类问题,由于难以避免的网格畸变,使模拟效率精度降低甚至失效,无网格方法(Meshless Method)则能在分析这类问题时显示出明显的优势,无网格局部Petrov-Galerkin(MLPG)法被誉为是一种有发展前景的真正无网格法.本文进一步发展了MLPG法,通过对任意的离散分布节点采用局部径向基函数构造插值形函数和Heaviside权函数,分析方程采用局部加权弱形式离散,建立了变量仅依赖于初始构型的完全Lagrange分析格式,最后用Newton-Raphson法迭代求解.文中分析了悬臂梁典型算例和微机电开关非线性大变形问题,通过与有限元结果的比较,表明本文提出的大变形问题无网格局部Petrov-Galerkin法具有稳定性好及收敛性快等优点.     

重叠网格法应用于模拟高速列车隧道气动效应

高速列车通过隧道时所诱发的压力波不仅会使乘客耳感舒适性变差,而且会造成车辆零部件疲劳。基于计算流体力学软件中有限体积方法,本文使用重叠网格法模拟了复杂外形高速列车完全通过隧道时引起的隧道内的压力波动。描述了重叠网格法及使用该方法进行列车过隧道建模的过程。利用文献提供的数据,分别采用滑移网格法和重叠网格法建模计算,并进行对比。通过与京沪线高速列车实车试验数据的对比,结果表明:对于模拟高速列车隧道气动效应,重叠网格法不仅网格设计简单,对车体、隧道复杂外形适应性强,且计算结果具有较高的计算精度及可靠性,具有较高的工程价值,可为今后高速铁路隧道气动效应的相关研究提供参考。     

不可压粘性流体浸入边界直接力法及软件实现

浸入边界法通过在N-S方程中施加体积力模拟不可滑移固壁边界及动边界，避免生成复杂贴体网格及动网格，极大地节省了网格建模时间及动网格计算消耗。本文提出一种新型附加体积力简化计算方法，将简化附加体积力以源项形式嵌入动量方程迭代中，通过用户自定义函数对CFD软件FLUENT二次开发，实现了浸入边界法和通用流体力学求解器的耦合计算。通过静止圆柱和动圆柱绕流数值模拟进行了验证，并探讨了插值函数对计算精度的影响。研究表明，通过引入浸入边界模型，能够提高计算效率，并实现结构网格背景下复杂边界和动边界的高效建模。     

不可压粘性流体浸入边界直接力法及软件实现

浸入边界法通过在N-S方程中施加体积力模拟不可滑移固壁边界及动边界,避免生成复杂贴体网格及动网格,极大地节省了网格建模时间及动网格计算消耗。本文提出一种新型附加体积力简化计算方法,将简化附加体积力以源项形式嵌入动量方程迭代中,通过用户自定义函数对CFD软件FLUENT二次开发,实现了浸入边界法和通用流体力学求解器的耦合计算。通过静止圆柱和动圆柱绕流数值模拟进行了验证,并探讨了插值函数对计算精度的影响。研究表明,通过引入浸入边界模型,能够提高计算效率,并实现结构网格背景下复杂边界和动边界的高效建模。     

HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics,CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度. 其次,通过NASATurbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶). 最后,通过NASA CommonResearchModel标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.      

HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics, CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度.其次,通过NASA Turbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶).最后,通过NASA Common Research Model标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.     

两套节点格林元嵌入式离散裂缝模型数值模拟方法

对于原始嵌入式离散裂缝模型(EDFM), 在计算包含裂缝单元的基质网格内的压力分布时采用了线性分布假设, 这导致了油藏开发早期对非稳态窜流量的计算精度不足. 因此, 本文提出了一种两套节点格林元法的EDFM数值模拟方法. 两套节点格林元法的核心思想是将压力节点与流量节点区分开, 一套压力节点设置在单元顶点, 另一套流量节点设置在网格边的中点, 满足局部物质守恒、具有二阶精度的同时, 可适用于任意网格类型. 本文将两套节点格林元法与EDFM耦合, 采用了非稳态渗流控制方程的边界积分形式推导了基质网格与裂缝网格之间传质量的新格式, 代替了线性分布假设以提高模拟精度; 此外, 修正后的EDFM能适应任意形态的基质网格剖分, 拓展了原始EDFM仅适用于矩形基质网格、难以考虑复杂油藏边界的局限性. 研究表明: 通过对比商业模拟软件tNavigator? LGR模块与原始EDFM, 验证了本文模型具有较高的早期计算精度; 以复杂油藏边界?缝网?SRV分区模型为例, 通过对比SFEM-COMSOL商业模拟软件, 验证了本文模型处理复杂问题的适应性. 本文研究可用于裂缝性油藏开发动态的精确模拟.      

扰动风作用下微型旋翼流场的数值模拟及流动控制研究

对于处在低雷诺数下的微型旋翼,扰动风对其流场影响很大。本文基于动态嵌套网格技术的双时间法,采用网格速度法模拟扰动风,研究了不同频率和幅值的正弦式扰动风下三维微型旋翼上的拉力系数响应情况,并采用变形网格方法模拟了旋翼桨叶的周期性边距运动对扰动风的流动控制,为今后旋翼机的主动控制研究提供了线索。研究结果表明,低雷诺数下,旋翼在周期性扰动风作用下的拉力响应曲线拟合函数与扰动风函数形式相同,响应曲线函数同样具有周期性变化趋势,波动幅值与扰动风幅值成正比例关系;对不同飞行状态下的旋翼,采用周期等角变距方式,能够找到一个适合的控制参数使得旋翼在正弦式扰动风环境下旋翼的拉力波动在一定范围内得到有效抑制。     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).     

Dynamics and Breakup of Viscoelastic Liquids (A Review)

The phenomena of hydrodynamic breakup of liquid jets, drops, films, bridges, and filaments are reviewed for liquids with viscoelastic properties. The reasons for breakup are capillary instabilities, collisions with rigid obstacles, and other forms of dynamic action. The relationship between the properties of the liquids and the features of the breakup process is discussed.     

Flow stability analysis and excitation using pulsating jets

Classical flow stability applied to transition from laminar to turbulent flow may also describe the behavior of vorticity fluctuations created by a pulsating jet placed along a solid boundary. A numerical laminar flow experiment involving a pulsating jet placed along the surface of a duct with flow separation downstream, resulted in eliminating most part of the separated flow region. Applying the same approach to a turbulent flow, it was possible to develop a turbulent stability flow formulation and apply successfully turbulent pulsating jet flow separation control. To cite this article: D. Skamnakis, K. Papailiou, C. R. Mecanique 333 (2005).