锥柱形膜片后屈曲的材料特性预示研究

使用正交表设计数值试验, 分析弹性模量、屈服强度、延伸率等材料特性参数对锥柱形膜片结构后屈曲状态的影响, 预示材料特性对膜片翻转的影响. 正交数值试验在不影响结果科学性的基础上, 极大地降低了试验次数. 数值试验结果表明: 在一定结构情况下, 弹性模量与屈服强度对膜片正向翻转影响非常显著, 而延伸率的影响与前两个因素比较而言, 几乎可以忽略. 并且, 结构顶点位移随材料的弹性模量增大而增大; 随材料的屈服强度减小而增大.因此, 在结构固定的情况下, 为了提高材料的翻转性能, 可以选择弹性模量较大而屈服强度较小的材料, 且延伸率为20%左右最好.     

柱形汇聚几何中内爆驱动金属界面不稳定性

采用自研的高保真度爆轰与冲击动力学程序，对柱形汇聚几何中内爆驱动金属材料界面不稳定性的动力学行为，进行了数值模拟研究。结果表明，首次冲击后至约12 μs，界面发展以RM（Richtmyer-Meshkov）不稳定性为主；12 μs后至冲击波聚心反弹加载前，界面聚心运动处于加速减速状态，界面发展由RT (Rayleigh-Taylor)不稳定性主导；冲击波聚心反弹加载后，界面发展又由RM不稳定性主导。另外，还研究了初始条件（初始振幅、初始波长、钢壳初始厚度和几何构型）对柱形内爆驱动金属材料界面不稳定性的影响。结果显示:初始振幅较大时振幅增长也较大；初始波长较小（模数较大）时振幅增长较小，而且存在一个截止波长；钢壳厚度会抑制扰动增长，也存在一个截止厚度；几何汇聚效应会使扰动增长速度更快。     

内爆压缩多层密绕螺线管的数值模拟

利用AUTODYN二次开发接口建立三维多层密绕螺线管数值模型,并实现周期性边界条件,应用光滑粒子动力学方法对内爆压缩多层密绕螺线管过程及其界面不稳定性开展数值模拟。计算结果表明,内爆压缩螺线管结构过程存在扰动快速增长至后期的界面失稳现象,与对应的实验结果较为相符。同时,计算显示螺线管结构参数对界面不稳定性发展具有显著影响,螺旋角度减小,结构压缩后期的界面不稳定性趋于严重;铜线直径减小,结构压缩后期的界面不稳定性趋于减弱。     

黏性对内空泡诱导柱状液滴界面不稳定性的影响

Rayleigh-Taylor不稳定性存在于爆炸、液滴形成和液体喷雾等工程应用过程中,是流体力学关注的经典问题之一.内空泡振荡诱导液滴界面演化问题是其研究中基本模型之一,空泡振荡作用下液滴界面发生扰动并发展,其特征形态主要表现为破碎、通气和稳定.液体黏性是影响界面不稳定性发展的重要因素,文章通过建立高精度的数值模拟方法,开展液体黏性对内空泡诱导柱状液滴界面不稳定性的影响研究.在数值模拟中,基于开源OpenFOAM框架的多相可压缩求解器直接求解Navier-Stokes方程,采用isoAdvector的几何流体体积法捕捉界面演化特征.结果表明,液体黏性的增加会减缓空泡的收缩,进而减缓液滴界面扰动的发展,该影响下通气工况液滴通气发生时间增加,而稳定工况最大扰动幅值减小.最大扰动幅值的减小直接影响了液滴的特征形态,基于一系列数值模拟结果归纳得到液滴不稳定性相图.在文章讨论的参数范围内,随着黏性增加,小液滴(R_d0 <2 mm)的形态从破碎转变为通气进而变成稳定;中液滴(2 mm d0 <3 mm)的形态从通气转变为稳定,不出现...     

沟槽表面填充材料对界面摩擦振动噪声的影响

采用机械加工方式,在基体(蠕墨铸铁)表面加工出平行且等间距的沟槽,并在沟槽中分别填满铝材料和锻钢材料,对该填充材料试样和基体试样进行界面摩擦噪声试验,并采用数值模拟方法加以分析,得出界面摩擦振动噪声行为在不同沟槽填充材料下的变化特性,并将不同界面摩擦磨损行为与振动噪声动态关联起来.试验及数值模拟结果表明,在沟槽中填充铝材料加剧界面磨损,增大摩擦系数,并加速系统不稳定振动的产生,增大噪声强度,而在沟槽中填充钢材料能有效延缓界面不稳定振动和摩擦尖叫噪声信号的产生.此外,计算结果表明,界面磨擦磨损行为导致的界面接触非线性特性是产生多频率不稳定振动和噪声的重要因素.     

岩体—界面系统剪切不稳定性分析

基于界面的刚塑性应变软化假设，分析了岩体－界面系统在端部剪力和岩体中分布剪切载荷共同作用下的变形、应力和损伤演化。利用位能原理和稳定性的能量准则，得到了岩体界面系统的不稳定性条件。分析结果表明，损伤区达到边界之前，系统可能是稳定的或者是不稳定的。依赖于分布剪切载荷和界面摩擦力的比较，若均布剪切载荷大于界面摩擦力，则系统不稳定，否则系统稳定；当损伤区达到边界之后，系统的不稳定性决定于载荷及界面材料性质，界面软化刚度系数和界面强度对于不稳定性有明显影响。     

不同初始磁场对激波冲击R22重气柱过程影响的数值模拟

基于磁流体动力(magneto-hydrodynamic,MHD)方程,采用CTU+CT方法,对在不同初始磁场作用下的平面入射激波与磁化R22重质气柱作用过程进行了数值研究。数值结果清晰地描述了不同初始磁场条件时激波诱导R22气柱界面不稳定性的过程,揭示了磁场控制界面不稳定性的机理。另外,还分析了磁感应强度对界面不稳定性的影响,发现在磁场较小时,涡层附着于界面,但随着磁感应强度的增大,平均涡量随之增大,涡层与界面逐渐分离,最终更好地抑制了界面不稳定性。同时,还发现平均涡度拟能随着磁感应强度的增大而减小,而垂直磁场比平行磁场更能降低平均涡度拟能,因而平均涡度拟能可较好地反映磁场对不稳定性的影响效果。     

等热流密度加热垂直下降液膜的研究

利用等热流密度加热条件下降膜流动的三维模型方程进行线性稳定性分析和数值模拟。线性稳定性分析表明,模型方程在小到中等Reynolds数下都适用,并且流向不稳定性增长率随着Reynolds数和Marangoni数增加而增加,展向不稳定性增长率则随着Marangoni数增加而增加,随着Reynolds数增加而减小,流向和展向对扰动波数都存在一个不稳定区间。三维数值模拟表明,在等热流密度加热条件下,液膜在随机扰动的情况下最终会形成带孤立波的三维溪流状结构,液膜与气体的换热也因溪流状结构的出现而加强;在随机扰动的基础上引入占优势地位的展向最不稳定扰动会使得换热增强,液膜会提前破裂;在随机扰动的基础上引入占优势地位的流向最不稳定扰动时,液膜的换热会增强,但不会提前破裂;在随机扰动的基础上同时引入占优势地位的流向和展向最不稳定扰动时,换热会加强且液膜会提前破裂。     

圆柱形界面相裂纹在扭转载荷作用下的动态断裂

本文研究了位于界面相中的圆柱形界面裂纹的扭转冲击问题.采用Laplace、Fourier变换和位错密度函数将混合边值问题转化为求解Cauchy核奇异积分方程,利用Laplace数值反演技术计算了动态应力强度因子.讨论了材料特性和结构的几何尺寸对动态应力强度因子的影响.结果表明,随着界面相厚度的增加,无量纲化的动态应力强度因子减小.当裂纹靠近剪切弹性模量大的材料时,无量纲化的动态应力强度因子增大,反之减小.界面相两侧不同的材料组合对裂尖动态应力强度因子的影响是随着剪切弹性模量和质量密度的比值的增加而减小.界面相中裂纹长度对裂尖动态应力强度因子的影响比其他因素的影响大.     

多次激波诱导正弦扰动预混火焰界面失稳的数值研究

激波诱导火焰失稳是实际中常见的现象,为深入研究火焰失稳特性,采用三维单步化学反应的Navier-Stokes方程和9阶weighted essentially non-oscillatory (WENO)的高精度格式,对不同马赫数的入射激波及其反射激波多次诱导正弦型预混火焰界面失稳的现象进行了三维数值模拟,并对计算结果的可靠性进行了验证。研究结果显示,在激波的多次作用下,火焰界面的演变过程主要受Richtmyer-Meshkov (RM)不稳定特性和化学反应特性的双重影响,且随着入射激波强度的增强,上述2种特性均得到进一步强化。为构造体现反应性RM不稳定特性的参数,根据火焰界面混合区平均涡量和化学反应速率,提出了表征界面受不稳定性和化学反应影响的量纲一参数。通过分析发现,在同一入射激波强度下,该参数的对数形式随入射激波和反射激波的多次作用呈基本相同的线性增长趋势;对不同马赫数的入射激波,该参数对数形式的线性增长率也基本一致。这样的变化表明该量纲一参数能够反映反应性RM不稳定过程中火焰界面发展的内在规律。     

Elastic-plastic crack growth on plane strain bimaterial interface

Finite element computation are carried out to simulate plane strain crack growth on a bimaterial interface under the assumption of small scale yielding. The modified Gurson constitutive equation and the element vanish technique introduced by Tvergaard et al. are used to model the final formation of an open crack. It is found from the calculation that the critical fracture toughness for crack growth is much lower in bimaterials than that in homogeneous material. The critical fracture toughness is strongly dependent on material properties of the bimaterial pair and the mixed mode of remote loads. The interface crack grows in the more compliant (lower hardening) material or in the weaker (lower yield strength) material. In Mode-I loading, the crack grows zigzag along the interface. Project supported by Fok Ying-Tung Education Foundation and National Natural Science Foundation of China.     

拉压异性线性等向强化材料厚壁球壳极限分析

运用Mohr屈服准则对承受内压的拉压屈服强度不同的线性等向强化材料的厚壁球壳进行了极限载荷分析，得到了依赖于拉压比和强化模量的厚壁球壳极限载荷解析式和依赖于拉压比和强化模量的厚壁球壳极限载荷分析式。结果表明，材料拉压屈服强度的不同和强化特性对厚壁球壳极限载荷均有一定的影响。     

An analytical study of the non-linear vibrations of cylindrical shells

The primary resonance response of simply supported circular cylindrical shells is investigated using the perturbation method. Donnell's non-linear shallow-shell theory is used to derive the governing partial differential equations of motion. The Galerkin technique is then employed to transform the equations of motion into a set of temporal ordinary differential equations. Considering only the primary resonance case, the method of multiple scales is used to study the periodic solutions and their stability. The necessary and sufficient conditions for appearance of the so-called companion mode are also discussed. To this end, a range of the possible multi-mode solution is obtained for response and excitation amplitudes and also excitation frequency as a function of damping, geometry and material properties of the shell. Parametric studies are performed to illustrate the effect of different values of thickness, length and material composition on the possibility of the companion mode participation in primary resonance response.     

Method of successive approximations for solving boundary-value problems of plasticity with allowance for the stress mode

A method for solving boundary-value problems of plasticity with allowance for the stress mode is developed. To describe the elastoplastic deformation of an isotropic material, use is made of constitutive equations that include two nonlinear functions dependent on the stress mode and determined experimentally. The elastoplastic state of a thin cylindrical shell under axisymmetric loading is calculated as a numerical example. The numerical results demonstrate good convergence of the method. The effect of the stress mode on the strain distribution in a cylindrical shell is assessed     

Linear Rayleigh-Taylor instability analysis of double-shell Kidder＇s self-similar implosion solution

This paper generalizes the single-shell Kidder＇s self-similar solution to the double-shell one with a discontinuity in density across the interface. An isentropic implosion model is constructed to study the Rayleigh-Taylor instability for the implosion compression. A Godunov-type method in the Lagrangian coordinates is used to compute the one-dimensional Euler equation with the initial and boundary conditions for the double-shell Kidder＇s self-similar solution in spherical geometry. Numerical results are obtained to validate the double-shell implosion model. By programming and using the linear perturbation codes, a linear stability analysis on the Rayleigh-Taylor instability for the double-shell isentropic implosion model is performed. It is found that, when the initial perturbation is concentrated much closer to the interface of the two shells, or when the spherical wave number becomes much smaller, the modal radius of the interface grows much faster, i.e., more unstable. In addition, from the spatial point of view for the compressibility effect on the perturbation evolution, the compressibility of the outer shell has a destabilization effect on the Rayleigh-Taylor instability, while the compressibility of the inner shell has a stabilization effect.     

Bifurcation phenomena of a biphasic compressible hyperelastic spherical continuum

A basic linear bifurcation problem is solved for a representative biphasic composite spherical cell. Setting is that of an internal inclusion that expends with an external shell made of a different material. Exact rate equations are derived in the framework of hyperelastic continuum mechanics. Sample examples are solved numerically revealing sensitivity of critical strain levels and bifurcated mode pattern to strength differential between both phases. Mathematical framework can be adapted to other types of constitutive responses.     

Decay of confined,two-dimensional,spatially periodic arrays of vortices: A numerical investigation

The disarrangement of a perturbed lattice of vortices was studied numerically. The basic state is an exponentially decaying, exact solution of the Navier-Stokes equations. Square arrays of vortices with even numbers of vortex cells along each side were perturbed and their evolution was investigated. Whether the energy in the perturbation grows somewhat before it decays or decays monotonically depends on the initial strength of the vortices of the basic state, the extent of lateral confinement and the structure of the perturbation. The critical condition for temporally local instability, i.e. the critical amplitude of the basic state that must be exceeded to allow energy transfer from the basic state to the perturbation, is discussed. In the strongly confined case of a square lattice of four vortices the appearance of enchancement of global rotation is the result of energy transfer from the basic state to a temporally local unstable mode. Energy is transferred from the basic state to larger-scaled structures (inverse cascade) only if the scales of the larger structures are inherently contained in the initial structure of the perturbation. The initial structure of the double array of vortices is not maintained except for a very special form of perturbation. The facts that large scales decay more slowly than small scales and that, when non-linearities are sufficiently strong, energy is transferred from one scale to another explain the differences in the disarrangement process for different initial strengths of the vortices of the basic state. The stronger vortices, i.e. the vortices perturbed in a manner that increases their strength, tend to dominate the weaker vortices. The pairing and subsequent merging (or capture) of vortices of like sense into larger-scale vortices are described in terms of peaks in the evolution of the square root of the palinstrophy divided by the enstrophy.     

Analysis of dynamic growth of a tensile crack in an elastic-plastic material

The influence of inertia on the stress and deformation fields near the tip of a crack growing in an elastic-plastic material is studied. The material is characterized by the von Mises yield criterion and J₂ flow theory of plasticity. The crack grows steadily under plane strain conditions in the tensile opening mode. Features of the stress and deformation state at points near the moving crack tip are described for elastic-perfectly plastic response and for several crack propagation speeds. It is found that inertia has a significant effect on the elastic-plastic response of material particles near the crack tip, and that elastic unloading may occur behind the crack tip for higher speeds. The relationship between the applied crack driving force, represented by a remote stress intensity factor, and the crack tip speed is examined on the basis of a critical crack tip opening angle growth criterion. The calculated result is compared with dynamic fracture toughness versus crack speed data for a 4340 steel.     

Distribution of dislocation source length and the size dependent yield strength in freestanding thin films

A method is proposed to estimate the size-dependent yield strength of columnar-grained freestanding thin films. The estimate relies on assuming a distribution of the size of Frank-Read sources, which is then translated into a log-normal distribution of the source strength, depending on film thickness, grain size and theoretical strength of the material, augmented with a single fit parameter. Two-dimensional discrete dislocation plasticity (DDP) simulations are carried out for two sets of Cu films and the fit parameter is determined from independent experiments. Subsequent DDP predictions of the stress-strain curves in comparison with the corresponding experimental data show excellent agreement of initial yield strength and hardening rate for films of varying film thickness and grain size. Having thus demonstrated the power of the proposed strength distribution, it is shown that the mode of this distribution governs the most effective source strength. This is then used to suggest a method to estimate the yield strength of thin films as a function of film thickness and grain size. Simple maps are presented that are in very good agreement with recent experimental results for Cu thin films.