三角波加载下金属铝动态破坏现象的微观模拟

采用嵌入原子势模型和分子动力学方法, 模拟研究了三角波加载下金属铝动态破坏的微观过程和动力学性质. 根据原子中心对称参数变化给出了样品微结构演化过程, 解读了熔化前后破坏过程的形态差异; 基于Virial定理统计了样品中压力和温度等力学量波形, 分析了熔化前后材料的强度变化. 通过不同碰撞速度的模拟, 讨论了破碎区内物质形态和密度分布的变化, 给出了材料破坏深度的变化规律. 研究还发现, 熔化后材料的动态拉伸强度已显著降低, 而此时由声学近似推算的材料拉伸强度已明显高于内部应力直接计算结果. 关键词： 破坏 分子动力学 冲击     

舰船水下接触爆炸数值研究

计及船体外板的材料强度,基于SPH算法建立水下接触爆炸数值模型,提出混合密度、广义光滑长度等解决数值模拟中计算崩溃的问题,对舰船的水下接触爆炸进行数值模拟,计算结果和冲击波基本理论及接触爆炸试验吻合良好.结果表明,接触爆炸时钢结构会产生圆形冲塞破口,其破口大小和炸药的厚度接近,当板与炸药的正则化厚度达1.8时,板仅会发生较大的塑性变形而不会冲塞破坏.并得到钢的反射系数约为1.2,及接触爆炸钢板的近似变形公式.     

激光驱动X光背光照相技术在金属靶微层裂研究中的应用探索

利用神光Ⅱ装置发射强激光一方面对铝靶进行烧蚀加载,另一方面驱动铜背光靶发射X光对铝靶微层裂产物进行针孔背光照相,结果获得了可忽略动态模糊的高清晰照片,标准丝阵检验表明分辨率优于40 m。铝靶微层裂产物照片清晰显示,靶破碎存在不同分区, 这直观揭示了靶破碎的不同机制。实验成功验证了激光驱动X光背光照相技术在金属靶微层裂研究中存在重要的应用前景。     

Enhancement of stability and accuracy of the moving particle semi-implicit method

As a Lagrangian gridless particle method, the MPS (Moving Particle Semi-implicit) method has been proven useful in a wide range of engineering applications. Up to now, most of MPS applications have been limited to problems with compressive stress states. This paper investigates the performance and stability of MPS method in simulation of more general hydrodynamic problems, including those characterized by tensile stresses or by changes in the stress states. It is shown that MPS-based simulations are prone to become destabilized in presence of attractive interparticle forces. This instability appears to be similar to the so-called tensile instability in SPH method. Two new modifications, namely, a modified Poisson Pressure Equation and a corrective matrix for pressure gradient model, are proposed to resolve this shortcoming. These two new modifications together with two previously proposed improvements are shown to stabilize and enhance the performance of MPS method.     

Potential of a piezoelectric energy harvester from sea waves

A sea wave energy harvester from the longitudinal wave motion of water particles is developed. The harvester consisting of a cantilever substrate attached by piezoelectric patches and a proof mass is used to collect electrical energy owing to the electromechanical coupling effect of the piezoelectric patches from the longitudinal wave motion. To describe the energy harvesting process, a mathematical model is developed to calculate the output charge and voltage from the piezoelectric patches according to the Airy linear wave theory and classical elastic beam model. Results show that the mean value of the generated power increases with the increase in the ratio of the width to the thickness of the cantilever, the wave height, the sea depth (which equals to the cantilever height in this study), the ratio of the proof mass to the cantilever mass, and the ratio of the sea depth to the wave length. A value of the power up to 55 W can be realized for a practical sea wave with the values of the sea depth, wave height and wave length to be 3 m, 2 m, and 15 m, respectively. The collected power harvesting with respect to different categories of the sea waves are provided. Our simulations also show the generated electric power can be further increased by an increase in dimensions of the harvester considering the scale effect. This research develops a new technique for energy harvesting from sea waves by piezoelectric energy harvesters.     

基于周向SH0导波的管道轴向缺陷定量表征*

随着管道服役时间的增加,其损伤逐渐累积,最终导致泄漏或爆炸事故。引入导波在线监测管道损伤是确保其安全运行的重要保障。周向零阶水平剪切(CSH0)波具有不易频散、对管道轴向缺陷敏感等优秀特性,非常适用于管道轴向缺陷的定量表征。该文 基于周向SH0导波在管道的传播特性,通过有限元模拟和试验对管道周向SH0模态导波与轴向缺陷的定量关系展开研究。结果表明：管道曲率半径越大,周向SH0模态的传播特性越接近板中SH0模态的传播特性;导波信号的反射系数和透射系数对轴向缺陷的长度和深度均呈现单调变化。依此建立了缺陷定量表征的关系云图,可用于评估缺陷尺寸。     

SPH方法对金属表面微射流的数值模拟

 光滑粒子流体动力学（Smoothed Particle Hydrodynamics）方法是一种无网格的拉氏计算方法,具有计算格式简单、易于计算大变形问题等优势。用二维SPH程序模拟了冲击加载下金属铝表面的沟槽状缺陷产生微射流的过程,得到了喷射物的总质量、最大速度和质量-速度曲线。还计算了相同深度、不同夹角的沟槽微射流,并与实验结果进行了比较。结果表明：计算得到的喷射物总质量和最大速度与实验结果符合较好,SPH方法对于估算金属表面的微射流提供了一种有效的数值模拟手段。     

非等温非牛顿黏性流体流动问题的修正光滑粒子动力学方法模拟

为准确、有效地模拟非等温非牛顿黏性流体的流动问题,本文基于一种不含核导数计算的核梯度修正格式和不可压缩条件给出了一种改进光滑粒子动力学(SPH)离散格式,它较传统SPH离散格式具有较高精度和较好稳定性.同时,为准确地描述温度场的演化过程,建立了非牛顿黏性的SPH温度离散模型.通过对等温Poiseuille流、喷射流和非等温Couette流、4:1收缩流进行模拟,并与其他数值结果作对比,分别验证了改进SPH方法模拟非牛顿黏性流动问题的可靠性和提出的SPH温度离散模型求解非等温流动问题的有效性和准确性.随后,运用改进SPH方法结合SPH温度离散模型对环形腔和C形腔内非等温非牛顿黏性流体的充模过程进行了试探性模拟研究,分析了数值模拟的收敛性,讨论了不同位置处热流参数对温度和流动的影响.     

Experimental and numerical study of the propagation of focused wave groups in the nearshore zone

The propagation of focused wave groups in intermediate water depth and the shoaling zone is experimentally and numerically considered in this paper. The experiments are carried out in a two-dimensional wave flume and wave trains derived from Pierson-Moskowitz and JONSWAP spectrum are generated. The peak frequency does not change during the wave train propagation for Pierson-Moskowitz waves; however, a downshift of this peak is observed for JONSWAP waves. An energy partitioning is performed in order to track the spatial evolution of energy. Four energy regions are defined for each spectrum type. A nonlinear energy transfer between different spectral regions as the wave train propagates is demonstrated and quantified. Numerical simulations are conducted using a modified Boussinesq model for long waves in shallow waters of varying depth. Experimental results are in satisfactory agreement with numerical predictions, especially in the case of wave trains derived from JONSWAP spectrum.     

爆炸焊接斜碰撞过程的数值模拟

 借助动力学分析软件ANSYS/LS-DYNA 11.0,建立了考虑材料强度的流体弹塑性力学模型,运用光滑粒子流体动力学（SPH）方法,对爆炸焊接斜碰撞过程进行了数值模拟。界面波和射流的模拟结果与爆炸焊接生产实际中的现象较为吻合,说明采用SPH方法对爆炸焊接斜碰撞过程进行模拟研究是有科学价值的。     

循环爆炸作用下地下洞室的动态响应及损伤累积

为研究循环爆炸对地下洞室的影响,基于相似模型试验,采用通用有限元软件ABAQUS对比研究了洞室拱顶高水平单次爆炸和低水平10次循环爆炸作用下地下洞室围岩的应力波衰减规律、损伤累积规律及洞壁位移和环向应变分布特征。结果表明:循环爆炸中,洞室围岩的应力波衰减速度随着爆炸次数的增加先减小后增大。单次爆炸中,洞壁环向峰值应变从拱顶至直墙脚由拉应变转为压应变;循环爆炸中,随着爆炸次数的增加,拱顶环向峰值应变由压应变转为拉应变。爆炸荷载总水平相同时,低水平循环爆炸中洞室围岩的损伤面积和程度比高水平单次爆炸大。循环爆炸中,围岩的损伤累积呈现不可逆的逐级增加趋势,且累积损伤和爆炸次数之间呈明显的非线性关系。     

Study on the structural characteristics of HgCdTe photodiodes using laser beam-induced current

The structural characteristics of typical n⁺-on-p HgCdTe photodiodes have been studied by laser beam-induced current (LBIC). The dependence of LBIC on laser wavelength, junction depth and localized leakage has been presented. The diffusion length of minority carrier of p-type region is extracted by the exponential decay fitting of the curve of LBIC. It is found that the peak magnitude of LBIC and junction depth approximate to a linear relationship for practical values of device fabrication. The diffusion length monotonously increases with the junction depth. A notable shift of LBIC profile is observed when localized leakage exists. This provides a useful explain for LBIC applying to characterize the structure and process uniformity of HgCdTe infrared detectors.     

Search for two-scale localization in disordered wires in a magnetic field

A recent paper [A. V. Kolesnikov and K. B. Efetov, Phys. Rev. Lett. 83, 3689 (1999)] predicts a two-scale behavior of wave function decay in disordered wires in the crossover regime from preserved to broken time-reversal symmetry. We have tested this prediction by a transmission approach, relying on the Borland conjecture that relates the decay length of the transmittance to the decay length of the wave functions. Our numerical simulations show no indication of two-scale behavior.     

完全变光滑长度SPH法及其实现

提出完全变光滑长度SPH法及其算法实现.方程组基于对称形式核函数近似,SPH密度演化方程与变光滑长度方程隐式关联;在Springel提出的全守恒SPH方程组基础上,通过将分散核近似形式改进为对称核近似形式得到SPH动量方程和能量方程.采用迭代求解密度演化方程和变光滑长度方程,显式求解SPH动量方程和能量方程,增加的计算量相对很少.给出三个1D激波管算例和2D Sedov算例验证方法的有效性.数值结果表明,算法保持动量和能量的守恒律,解决了传统SPH法中由于变光滑长度影响带来的计算误差,且在模拟2DSedov问题时能得到比Springel方法更准确的压强峰值位置和中心压强值.特别适合于模拟爆炸与冲击、大变形大扭曲等密度梯度和光滑长度梯度剧烈变化的问题.     

碳纳米管-聚乙烯复合材料界面力学特性分析

本文采用分子动力学模拟办法对碳纳米管-聚乙烯复合材料的界面力学特性进行了模拟和分析. 通过对单壁碳纳米管从无定形聚乙烯中抽出过程进行模拟, 研究了界面剪切应力随碳管滑移速度、聚乙烯分子链长和碳纳米管管径之间的变化关系, 并对界面的滑移机理进行了讨论. 模拟结果发现, 随着聚合物分子链长的增加, 界面临界剪切应力有显著增大, 而滑移剪切应力略显增加; 界面临界剪切应力和滑移剪切应力随着碳纳米管管径的增大而明显增加. 本文同时对界面应力的变化机理进行了模拟和讨论.     

Mechanism of deformation-twin formation in nanocrystalline metals

A theoretical model is proposed to describe the nucleation of deformation twins at grain boundaries in nanocrystalline materials under the action of an applied stress and the stress field of a dipole of junction or grain-boundary wedge disclinations. The model is used to consider pure nanocrystalline aluminum and copper with an average grain size of about 30 nm. The conditions of barrier-free twinning-dislocation nucleation are studied. These conditions are shown to be realistic for the metals under study. As the twin-plate thickness increases, one observes two stages of local hardening and an intermediate stage of local flow of a nanocrystalline metal on the scale of one nanograin. In all stages, the critical stress increases with decreasing disclination-dipole strength. The equilibrium thickness and shape of the twin plate are analyzed and found to agree well with the well-known results of experimental observations.     

Simulating coupled dynamics of a rigid-flexible multibody system and compressible fluid

As a subsequent work of previous studies of authors, a new parallel computation approach is proposed to simulate the coupled dynamics of a rigid-flexible multibody system and compressible fluid. In this approach, the smoothed particle hydrodynamics (SPH) method is used to model the compressible fluid, the natural coordinate formulation (NCF) and absolute nodal coordinate formulation (ANCF) are used to model the rigid and flexible bodies, respectively. In order to model the compressible fluid properly and efficiently via SPH method, three measures are taken as follows. The first is to use the Riemann solver to cope with the fluid compressibility, the second is to define virtual particles of SPH to model the dynamic interaction between the fluid and the multibody system, and the third is to impose the boundary conditions of periodical inflow and outflow to reduce the number of SPH particles involved in the computation process. Afterwards, a parallel computation strategy is proposed based on the graphics processing unit (GPU) to detect the neighboring SPH particles and to solve the dynamic equations of SPH particles in order to improve the computation efficiency. Meanwhile, the generalized-alpha algorithm is used to solve the dynamic equations of the multibody system. Finally, four case studies are given to validate the proposed parallel computation approach.     

Absolute peak slope time based thickness measurement using pulsed thermography

Several methods have been reported in the literature using pulsed thermography for quantitative measurement of defect depth or sample thickness. In this paper, based on the analysis of a theoretical one-dimensional solution of pulsed thermography, we proposed to use the absolute peak slope time (APST) for quantitative measurement of defect depth. APST is the peak slope time of the curve which is obtained by multiplying the original temperature decay with the square root of the corresponding time. The theoretical model shows that APST has linear relation with square of defect depth, which was verified with the experimental results of an aluminum and a steel specimen with six flat-bottom wedges and holes as simulated defects respectively.     

Oscillatory disturbance in force calibration of optical tweezers

In the calibration of the optical trap stiffness, it is found that there appears an attenuating oscillation as an oscillatory disturbance added to the trapped bead movement, when the scanner is driven by a triangular wave input.An equivalent oscillator model is put forward to explain the mechanism of the oscillatory disturbance. Both the measurements and calculations show that the attenuating oscillation comes from the oscillation of the scanner and the triangular wave drive causes this additional oscillation of the scanner. Furthermore, the analysis indicates that the oscillatory disturbance will become stronger, when the stiffness of the trap increases or the natural frequency of the scanner decreases. We adopt another driving way, i.e. a sinusoidal wave input is used instead of the triangular wave input. Our experiment has verified that in this case the oscillatory disturbance is eliminated completely.