无网格MPM法三维爆炸焊接数值模拟

物质点法MPM(Material Point Method)是无网格方法之一.它是在质点网格法(PIC)基础上发展而来的一种新数值方法,它利用了欧拉法和拉格朗日法两者的优点,计算物质点在冲击载荷下的应力和应交;通过物质点来跟踪材料体的变形和破损,而在整个计算过程中背景网格始终固定不变,避免了重新划分网格.本文应用MPM法计算三维爆炸焊接问题,在爆轰载荷作用下的飞板和基板的金属动态变形过程进行了三维数值模拟,并且对飞板的碰撞点速度和爆轰压力变化进行了计算分析.     

滑移爆轰问题无网格MPM法数值模拟

利用MPM(material point method)计算方法,以爆炸焊接作为实例,对炸药滑移爆轰过程和金属飞板与基板之间的碰撞变形进行了数值模拟.并将无网格MPM法的数值计算结果与近似解析公式的计算结果进行了比较,二者基本吻合,有力地证明了无网格MPM法在求解爆炸冲击问题中的有效性和健壮性.     

爆轰波碰撞聚能无网格MPM法数值模拟

爆轰聚能经常被应用于工程爆破切割,为了探索新的聚能爆破形式,本文应用无网格MPM法对圆柱筒装药和导爆索多点起爆形式的爆轰波碰撞聚能问题进行三维数值模拟,采用显式积分算法完成爆轰波形成和汇聚过程以及爆轰压力场分布的数值计算,并将数值模拟结果与实验结果做了对比分析。由此验证了无网格MPM法的准确性,也为炸药聚能爆破提供了新的思路和高效的操作方法。     

裂隙岩体局部化破坏分析的多重势面分岔模型

采用多重势面弹塑性分岔理论对裂隙岩体的局部化破坏进行计算模拟．基于多重势面弹塑性理论分析局部化问题，构造了适用于裂隙岩体破坏的多重势面不连续分岔模型，并使用数值方法求解局部化方阶在有限元方法的基础上，使用该模型计算裂隙岩体的局部化破坏条带．算例分析表明这一模型用于分析裂隙岩体的局部化破坏是有效的．     

裂隙岩体损伤局部化破坏分岔模型及其应用

采用概率统计方法分析节理裂隙岩体的几何特征，定义了反映裂隙岩体几何特征的组构张量．基于不可逆热力学理论，通过裂纹扩展的细观分析，得出了损伤的发展机理和演化方程，把损伤演化和裂隙的几何特征的变化联系起来，建立了弹塑性损伤本构关系．为分析含有节理裂隙岩体在发生局部化破坏时的特征，通过对发生局部化时的裂隙岩体的分析。构造了适用于节理裂隙岩体局部化分析的不连续分岔模型．利用非线性规划数值解法，可以得出局部化破坏的方向特征．在有限元方法中，根据该模型给出了节理裂隙岩体相关的算例，分析表明该模型用于分析裂隙岩体的局部化破坏是有效的．     

页岩气藏压裂缝网扩展数值模拟

为探究页岩气藏水力压裂复杂裂缝网络的形成机理,开展了缝网扩展的数值模拟研究.考虑应力阴影和天然裂缝作用,建立了井筒和裂缝中流体流动模型,利用位移不连续方法求解应力与位移不连续量,然后构建了压力与裂缝宽度的迭代方程,并采用牛顿迭代法求解.通过比较数值解经典模型解析解,验证了模型和迭代解法的正确性.多簇裂缝同步扩展时裂缝间距越小,压裂液分配到各条裂缝越不均匀,靠近井筒跟部的裂缝的分流量越大,从而裂缝宽度越大;考虑天然裂缝作用时,逼近角越小或者应力各向异性越弱,水力裂缝越容易发生转向扩展,裂缝网络越复杂.      

巴西圆盘试样中心斜裂纹疲劳扩展轨迹的边界元模拟

提出了一个平面弹性体多裂纹疲劳扩展模型. 它主要涉及到复合型加载情况下多裂纹尖端疲劳扩展的数学模型及杂交位移不连续法(一种边界元法). 在数值模拟中, 对每一裂纹扩展增量分析时,在其先前的边界上增添裂纹扩展增量, 且只对新增添的裂纹扩展增量划分单元, 同时, 按照这种边界元法的实施方法对一些单元特征进行调整, 就可以方便地模拟裂纹扩展. 用这种数值方法模拟了巴西圆盘试样中心斜裂纹疲劳扩展轨迹,数值结果说明了预报模型的有效性, 揭示了裂纹体几何对疲劳扩展的影响.      

无网格MPM法三维前处理系统设计

无网格法在解决爆炸冲击、裂纹扩展、高速碰撞以及大变形等问题时,由于与网格无关,能够避免有限元法因网格畸变计算精度降低或计算终止的矛盾,是目前求解固体力学问题的有效方法之一,也是近年来研究的热点。前处理对于无网格法数值计算是非常重要的环节,基于面向对象程序设计思想,应用C 和OpenGL开发了三维图形交互方式的MPM法前处理系统,提出了应用递归四叉树法、扫描表示法和CSG法划分任意形状非等间距的MPM法质点单元,为建立通用的无网格前处理平台提出了有效的解决方案。     

混凝土三点弯曲试件破坏过程的数值模拟

根据混凝土试件三点弯曲试验的物理模型,用MFPA数值模拟软件分别对砂浆和混凝土的三点弯曲试件的破坏全过程进行模拟．给出两个试件的载荷与加载点位移关系曲线,并进行了对比分析．     

含有序孔洞的多孔氧化铝薄膜的力学性能测量和数值模拟

高纯铝箔在特定的溶液下经过电化学阳极氧化腐蚀,可在其表面生成一层多孔的非晶氧化铝层,其孔径分布非常均匀,孔大致呈六方密排布。由于此类薄膜具有规则的纳米级孔径,大的比表面积,良好的自组织排列性,所以其日益受到人们的关注。然而,到目前为止,对于此类薄膜力学性能的研究还很少,所以在一定程度上限制其功能的开发和应用。为了获得此类多孔膜的弹性常数,本文首先由实验出发,通过光力学检测(双光束散斑干涉)的方法得到薄膜拉伸时的整体表观弹性模量。然后对薄膜建立二维有限元模型,运用均匀化理论反推出其基体(无孔结构)的弹性模量,同时考察了不同的基体泊松比对模型整体表观模量的影响,并且用一般有限元方法验证了沿特定方向拉伸时均匀化模型计算的有效性。     

Numerical studies of the fingering phenomena for the thin film equation

We present a new interpretation of the fingering phenomena of the thin liquid film layer through numerical investigations. The governing partial differential equation is h_t + (h²?h³)_x = ??·(h³?Δh), which arises in the context of thin liquid films driven by a thermal gradient with a counteracting gravitational force, where h = h(x, y, t) is the liquid film height. A robust and accurate finite difference method is developed for the thin liquid film equation. For the advection part (h²?h³)_x, we use an implicit essentially non‐oscillatory (ENO)‐type scheme and get a good stability property. For the diffusion part ??·(h³?Δh), we use an implicit Euler's method. The resulting nonlinear discrete system is solved by an efficient nonlinear multigrid method. Numerical experiments indicate that higher the film thickness, the faster the film front evolves. The concave front has higher film thickness than the convex front. Therefore, the concave front has higher speed than the convex front and this leads to the fingering phenomena. Copyright © 2010 John Wiley & Sons, Ltd.     

薄膜润滑中的微极流体效应

利用考虑微结构与微旋效应的微极流体理论来研究薄膜润滑的特性。微极效应将引起等效黏度的增加,从而影响油膜厚度。它对压力分布和膜厚形状的影响很小。数值模拟结果表明,等效黏度的增加效应与实验值有很好的类似性。     

Controlled buckling of thin film on elastomeric substrate in large deformation

Electronic systems with large stretchability have many applications. A precisely controlled buckling strategy to increase the stretchability has been demonstrated by combining lithographically patterned surface bonding chemistry and a buckling process. The buckled geometry was assumed to have a sinusoidal form, which may result in errors to determine the strains in the film. A theoretical model is presented in this letter to study the mechanics of this type of thin film/substrate system by discarding the assumption of sinusoidal buckling geometry. It is shown that the previous model overestimates the deflection and curvature in the thin film. The results from the model agree well with finite element simulations and therefore provide design guidelines in many applications ranging from stretchable electronics to micro/nano scale surface patterning and precision metrology.     

Beating capillarity in thin film flows

The combination of substrate unevenness and capillarity is known to induce far‐reaching perturbations at the free surface of thin liquid films. These might be undesired and this paper explores the possibility to control the free surface of thin liquid films to give it a prescribed profile by a suitable design of the underlying substrate. This corresponds to the inverse of the widely studied forward problem, which considers the effect of substrate unevenness on a free surface. Assuming that the steady free surface profile can be described by the lubrication approximation, this optimal control problem is shown to be governed by a first‐order partial differential equation, which is solved numerically using the method of characteristics. The proposed method is successfully tested for a range of desired free surface profiles and the domain of existence of a solution to the inverse problem is probed. Expectedly, it is shown that, owing to surface tension, not all free surface profiles can be achieved but in some cases capillarity can be beaten and a prescribed free surface profile obtained. Copyright © 2009 John Wiley & Sons, Ltd.     

A parametric study of laser induced thin film spallation

We report parametric studies of elastic wave generation by a pulsed laser and associated spalling of thin surface films by the corresponding high stresses. Two different substrate materials, single crystal Si (100) and fused silica, are considered. Spallation behavior of Al thin films is investigated as a function of substrate thickness, film thickness, laser energy, and various parameters governing the source. Surface displacement due to the stress wave is measured by Michaelson interferometry and used to infer the stresses on the film interface. Consistent with previous studies, the maximum stress in the substrate and at the film/substrate interface increases with increasing laser fluence. For many of the conditions tested, the substrate stress is large enough to damage the Si. Moreover, the maximum interface stress is found to increase with increasing film thickness, but decrease with increasing substrate thickness due to geometric attenuation. Of particular significance is the development of a decompression shock in the fused sillica substrates, which results in very high tensile stresses at the interface. This shock enhances the failure of thin film interfaces, especially in thicker samples.     

铁电薄膜在红外辐射作用下的响应

研究了红外探测器中红外焦平面列阵的象元———铁电薄膜微桥在红外辐射作用下的输出信号。用层合板壳作为微桥结构的力学模型,中间一层为压电材料,上下两层为金属材料电极。采用了力-热-电耦合的控制方程和变分原理,推导出了基于Mindlin假设的压电材料层合板有限元公式。以红外探测器在夜间从飞机上探测地面的坦克为例,用有限单元法模拟了铁电薄膜微桥在红外辐射作用下的力、热、电输出信号,并对结果作了分析比较。     

压电薄膜表面贯穿裂纹的有限元分析

基于有限元软件ANSYS数值模拟,计算了激光作用下的压电薄膜表面贯穿裂纹外场应力强度因子和电位移强度因子,并且研究了90°畴变所诱致的畴变增韧行为。首先,求解无裂纹压电薄膜在激光作用下的热-力-电响应,将求得的应力和电位移场反向作用于裂纹面,求解裂纹尖端处的外场应力和电位移强度因子,然后基于小范围畴变理论求解了90°畴变所致的屏蔽应力强度因子。讨论了薄膜表面裂纹的外场应力强度因子、电位移强度因子及屏蔽应力强度因子随激光作用时间和裂纹位置的变化关系,从而预测压电薄膜体系在加热工作状况下的裂纹扩展和断裂行为。     

嵌入薄膜传感器悬臂梁三向力测量技术研究

以薄膜传感器悬臂梁作为等效模型,通过传感器的应变效应对三向力测量技术进行了研究。为提高薄膜传感器的应变输出响应,对悬臂梁上布放薄膜传感器的位置加设弹性结构,研究了三向力测力模型输出电压与传感器所在位置应变的关系;分析了受力位置对测力模型输出响应的影响关系,结合实验验证了其工作原理、测力模型应变输出响应与可控尺寸参数的关系。研究表明:该测力模型可实现三向力测量,各个方向最大测量误差均在9%以内,悬臂梁宽度方向x和高度方向z的交叉干涉误差分别为2.84%和3.37%;当悬臂梁自由端受力位置发生变化时,测力模型输出响应只在梁长度方向y上发生变化。     

Numerical simulation of thin airfoil stall by using a modified DES approach

The detached‐Eddy simulation (DES) method was applied to calculate pre‐ and post‐stall aerodynamic characteristics of airfoil stall. A discrepancy between numerical and experimental data was observed near the stall regime for the airfoil NACA64_A‐ 006 which is a thin airfoil stall type. The reason of this discrepancy and one possible way for improvement of the numerical model are discussed here. It is shown that the use of the Baldwin–Lomax model in the RANS region improves the DES results in this case. If the relevant factors (grid density, time step, turbulence model, etc.) are appropriately taken into account, the DES approach could reliably predict stall aerodynamical characteristics. Copyright © 2006 John Wiley & Sons, Ltd.