基于薄层单元与弹簧单元的滚动直线导轨副动力学建模

针对某大型光机装置中使用的LM型直线导轨-滑块的结合面,采用薄层单元与弹簧单元模拟其接触面,建立了大行程传输结构模拟件的有限元模型。通过模型修正方法,根据模态试验结果对薄层单元弹性模量以及弹簧刚度进行了识别。模型修正后,结构前三阶模态频率计算结果与试验结果最大差别为2.23%,地脉动载荷下各测点位移响应计算结果与试验结果最大差别为7.61%。计算结果与试验结果具有较好的一致性,验证了模型的有效性。     

强激光干扰下运动目标热损伤效应的有限元研究

为对空中飞行目标实施合理有效的强激光干扰,需要综合考虑采用的干扰方式、干扰效果以及考虑光斑偏移对干扰效果的影响。以热传导理论为基础,建立了激光辐照下移动目标温度场模型,采用伽辽金法建立了有限元分析模型,采用了等效比热容的方法处理材料的相变过程;并以Ansys软件为仿真平台,对目标在高斯分布的强激光作用下的温度场和热应力进行了综合仿真分析。结果表明,光斑与目标辐照区域的相对偏移对目标表面温度影响较大,对目标纵深破坏效应相对较小,应力变化主要发生在激光辐照区域,应力最大处也会随热源的移动而变化。     

Multi-focus design of underwater noise control linings based on finite element analysis

Varied, counter-demanding objectives in designing the underwater noise control linings are addressed using a finite element model based methodology. Four different kinds of designs are proposed to attend to diverse and conflicting requirements concerning echo reduction (ER) and transmission loss (TL) performance of these linings. In this regard a slightly modified hybrid type finite element based on the Pian and Tong (PT) formulation has been used to make the computational efforts less demanding as compared to the original one. The adequacy of this formulation has been shown by comparing its results with the analytical, finite element analysis based, and experimental results. Different unit cell representations for different types of distributions of air cavities on the linings are discussed with respect to their limitations and applicability. Effect of static pressure is studied by using a simplified technique which can be used to simulate deep sea testing environment. Performance variation of different designs is investigated under different water depths to study their applicability in such situations.     

A nodal discontinuous Galerkin finite element method for nonlinear elastic wave propagation

A nodal discontinuous Galerkin finite element method (DG-FEM) to solve the linear and nonlinear elastic wave equation in heterogeneous media with arbitrary high order accuracy in space on unstructured triangular or quadrilateral meshes is presented. This DG-FEM method combines the geometrical flexibility of the finite element method, and the high parallelization potentiality and strongly nonlinear wave phenomena simulation capability of the finite volume method, required for nonlinear elastodynamics simulations. In order to facilitate the implementation based on a numerical scheme developed for electromagnetic applications, the equations of nonlinear elastodynamics have been written in a conservative form. The adopted formalism allows the introduction of different kinds of elastic nonlinearities, such as the classical quadratic and cubic nonlinearities, or the quadratic hysteretic nonlinearities. Absorbing layers perfectly matched to the calculation domain of the nearly perfectly matched layers type have been introduced to simulate, when needed, semi-infinite or infinite media. The developed DG-FEM scheme has been verified by means of a comparison with analytical solutions and numerical results already published in the literature for simple geometrical configurations: Lamb's problem and plane wave nonlinear propagation.     

Laser cutting silicon-glass double layer wafer with laser induced thermal-crack propagation

This study was aimed at introducing the laser induced thermal-crack propagation (LITP) technology to solve the silicon-glass double layer wafer dicing problems in the packaging procedure of silicon-glass device packaged by WLCSP technology, investigating the feasibility of this idea, and studying the crack propagation process of LITP cutting double layer wafer. In this paper, the physical process of the 1064 nm laser beam interact with the double layer wafer during the cutting process was studied theoretically. A mathematical model consists the volumetric heating source and the surface heating source has been established. The temperature and stress distribution was simulated by using finite element method (FEM) analysis software ABAQUS. The extended finite element method (XFEM) was added to the simulation as the supplementary features to simulate the crack propagation process and the crack propagation profile. The silicon-glass double layer wafer cutting verification experiment under typical parameters was conducted by using the 1064 nm semiconductor laser. The crack propagation profile on the fracture surface was examined by optical microscope and explained from the stress distribution and XFEM status. It was concluded that the quality of the finished fracture surface has been greatly improved, and the experiment results were well supported by the numerical simulation results.     

A finite element model for laser-induced leaky waves at fluid–solid interfaces

The finite element method is firstly used to simulate the laser-induced leaky waves at fluid–solid interfaces. Corresponding models and arithmetic are developed, in that the fluid–solid interactions are described by a coupling matrix and the infinite boundary of fluid domain is modeled by acoustic absorption elements. Typical calculations are executed for air–aluminum plane and cylindrical interfaces. The results are in very good agreements with the experimental signals in available literatures, which verify the correctness of our finite element model for simulating laser-induced leaky wave at fluid–solid interfaces. And some elementary conclusions are obtained for the laser induced leaky waves.     

高速气流作用下激光加热金属平板数值模拟

数值模拟研究了高速气流作用下激光加热金属平板温度场。流体控制方程为三维雷诺平均Navier-Stokes方程,固体控制方程为能量方程,湍流粘性系数求解使用k-ε两方程模型。采用流固耦合计算方法,使用两相流模型模拟气流对烧蚀物的剥蚀,较完整地模拟了激光辐照金属材料的物理变化过程,计算得到了不同气流速度下金属平板的温度分布以及烧蚀形貌。分别使用两相流方法和动网格方法对高速气流作用下激光对金属板的烧蚀效应进行了计算,结果表明,两相流方法与动网格方法都能较好地模拟高速气流作用下激光加热金属平板的温度响应,由于两相流方法能够较全面地模拟对流换热、熔化与凝固过程以及金属液体在气流冲刷下的动力学过程,因此能获得比动网格方法更为合理的物理图像。     

Finite element modeling of heating phenomena of cracks excited by high-intensity ultrasonic pulses

A three-dimensional thermo-mechanical coupled finite element model is built up to simulate the phenomena of dynamical contact and frictional heating of crack faces when the plate containing the crack is excited by high-intensity ultrasonic pulses.In the finite element model,the high-power ultrasonic transducer is modeled by using a piezoelectric thermal-analogy method,and the dynamical interaction between both crack faces is modeled using a contact-impact theory.In the simulations,the frictional heating taking place at the crack faces is quantitatively calculated by using finite element thermal-structural coupling analysis,especially,the influences of acoustic chaos to plate vibration and crack heating are calculated and analysed in detail.Meanwhile,the related ultrasonic infrared images are also obtained experimentally,and the theoretical simulation results are in agreement with that of the experiments.The results show that,by using the theoretical method,a good simulation of dynamic interaction and friction heating process of the crack faces under non-chaotic or chaotic sound excitation can be obtained.     

Finite element analysis of ionic liquid gel soft actuator

A new type of soft actuator material-ionic liquid gel(ILG),which consists of HEMA,BMIMBF4,and TiO_2,can be transformed into gel state under the irradiation of ultraviolet(UV)light.In this paper,Mooney–Rivlin hyperelastic model of finite element method is proposed for the first time to study the properties of the ILG.It has been proved that the content of TiO_2 has a great influence on the properties of the gel,and Young’s modulus of the gel increases with the increase of its content,despite of reduced tensile deformation.The results in this work show that when the TiO_2 content is 1.0 wt%,a large tensile deformation and a strong Young’s modulus can be obtained to be 325%and 7.8 k Pa,respectively.The material parameters of ILG with TiO_2content values of 0.2 wt%,0.5 wt%,1.0 wt%,and 1.5 wt%are obtained,respectively,through uniaxial tensile tests,including C_(10),C_(01),C_(20),C_(11),C_(02),C_(30),C_(21),C_(12),and C_(03)elements.In this paper,the large-scaled general finite element software ANSYS is used to simulate and analyze the ILG,which is based on SOLID186 element and nonlinear hyperelastic Mooney–Rivlin model.The finite element simulation analysis based stress-strain curves are almost consistent with the experimental stress–strain curves,and hence the finite element analysis of ILG is feasible and credible.This work presents a new direction for studying the performance of soft actuator for the ILG,and also contributes to the design of soft robot actuator.     

激光支持等离子体流场的2维动态数值模拟

 采用高温气体（等离子体）电离度的一种近似计算方法，以及具有五阶精度的广义Godunov差分格式-加权本质无振荡格式，给出了高温气体状态方程的简便描述，并考虑激光与等离子体的相互作用和能量耦合模拟了强激光与固体靶相互作用时激光支持靶面等离子体流场的动态物理过程，并同理论模型结果进行了比较。结果表明：不考虑2维膨胀效应的理论模型结果比2维数值计算结果偏大，但两者只有接近一倍的差距，并无数量级的差别。考虑了膨胀效应的2维数值计算结果是可信的。     

Discrete element crowd model for pedestrian evacuation through an exit

A series of accidents caused by crowds within the last decades evoked a lot of scientific interest in modeling the movement of pedestrian crowds. Based on the discrete element method, a granular dynamic model, in which the human body is simplified as a self-driven sphere, is proposed to simulate the characteristics of crowd flow through an exit. In this model, the repulsive force among people is considered to have an anisotropic feature, and the physical contact force due to body deformation is quantified by the Hertz contact model. The movement of the human body is simulated by applying the second Newton's law. The crowd flow through an exit at different desired velocities is studied and simulation results indicated that crowd flow exhibits three distinct states, i.e., smooth state, transition state and phase separation state. In the simulation, the clogging phenomenon occurs more easily when the desired velocity is high and the exit may as a result be totally blocked at a desired velocity of 1.6 m/s or above, leading to faster-to-frozen effect.     

Simulating regular wave propagation over a submerged bar by boundary element method model and Boussinesq equation model

Numerical models based on the boundary element method and Boussinesq equation are used to simulate the wave transform over a submerged bar for regular waves.In the boundary-element-method model the linear element is used,and the integrals are computed by analytical formulas.The Boussinesq-equation model is the well-known FUNWAVE from the University of Delaware.We compare the numerical free surface displacements with the laboratory data on both gentle slope and steep slope,and find that both the models simulate the wave transform well.We further compute the agreement indexes between the numerical result and laboratory data,and the results support that the boundary-element-method model has a stable good performance,which is due to the fact that its governing equation has no restriction on nonlinearity and dispersion as compared with Boussinesq equation.     

激光烧结石墨烯铜复合材料温度场有限元模拟

在激光烧结石墨烯增强铜基复合材料的过程中,了解瞬时温度场分布对优化工艺参数、控制烧结质量有重要作用。建立了激光烧结预涂在42CrMo基板上的石墨烯铜的混合粉末的有限元模型。研究了激光烧结过程温度场分布,熔池的几何参数以及烧结层与基体的冶金结合宽度。为了验证模拟结果,使用与模拟相同的参数进行了单道激光烧结的实验。研究表明,热传导、热辐射和相变潜热在激光烧结过程的温度场分布中起重要作用。实验结果与模拟结果较为一致。所以可以依据模拟结果预测实验的温度场分布和熔池几何参数,同时也可以据此优化激光烧结参数。     

A general finite element model for numerical simulation of structure dynamics

A finite element model used to simulate the dynamics with continuum and discontinuum is presented. This new approach is conducted by constructing the general contact model. The conventional discrete element is treated as a standard finite element with one node in this new method. The one-node element has the same features as other finite elements, such as element stress and strain. Thus, a general finite element model that is consistent with the existed finite element model is set up. This new model is simple in mathematical concept and is straightforward to be combined into the existing standard finite element code. Numerical example demonstrates that this new approach is more effective to perform the dynamic process analysis in which the interactions among a large number of discrete bodies and continuum objects are included.     

贴片晶振在冲击环境下的损伤边界

贴片式石英晶体振荡器广泛应用于各类电子和通信设备系统中。针对晶振在冲击环境中容易出现结构破坏而导致系统工作异常的问题,通过分析单自由度系统在不同频率冲击载荷作用下的响应特点,建立了结构的应力响应水平与相关冲击响应谱谱值之间的联系,获得了较已有结论更合理的损伤边界形式。根据典型晶振结构易损组件的力学特性建立对应的简化分析模型,得到了贴片晶振在大频率范围内的结构损伤边界。利用有限元仿真软件,对晶振结构在0.5~30 kHz频率范围内冲击载荷下的响应进行仿真分析,验证了该损伤边界的有效性。这也为以贴片晶振为代表的微小元器件在冲击环境下的可靠性研究提供了一种可行的方法。     

重力作用下颗粒介质应力链的离散元模拟

用离散元的方法模拟了仅有重力作用的二维颗粒系统内部力的分布情况,并根据力的大小得到颗粒之间的应力链.模拟结果与颗粒介质研究中的两个著名模型q模型和α模型作了对比,并与光弹实验的结果作了比较.对比结果表明,模拟结果与实验相似,而与两个概率模型有一定的差异.另外计算结果还表明,颗粒介质中力大小的概率分布极为不均匀,较大的力概率呈指数衰减,应力链的分布具有分形特征. 关键词： 颗粒介质 离散元 应力链 光弹实验     

钛合金纳秒激光打孔数值模拟和实验研究

为了描述纳秒激光对钛合金打孔过程中孔的形貌及温度场的变化规律,建立激光打孔的物理模型,利用ANSYS中APDL语言进行编程,对温度场进行仿真分析,并利用单元生死技术模拟孔形貌的变化过程。从有限元数值模拟和实验两方面综合分析比较了激光工艺参数(脉冲能量和脉冲数量)对打孔质量(孔深和孔径)的影响,系统论述了钛合金纳秒激光打孔的一般规律,以达到工艺参数优化,提高打孔质量的目的。     

冲击载荷下混凝土本构模型构建研究

 在对混凝土动态力学性能和现有本构模型综合分析的基础上,构建了一个新的适用于冲击响应问题数值分析的混凝土本构模型。该本构模型全面考虑了压力、应力第三不变量、变形的硬化和软化、应变率强化以及拉伸损伤等各个影响因素。将其加入LTZ-2D程序,确定了本构模型参数,对混凝土靶板的穿透问题进行了数值验证分析。计算得到的弹体剩余速度同实验结果基本一致,同时得到了混凝土靶板破裂的计算图像。计算结果及其分析表明,所构建的本构模型能够较好地反映冲击载荷作用下混凝土动态响应的主要特性。     

光学元件微缺陷处电磁场分布特性的数值计算方法

 神光Ⅲ的光学系统具有高功率密度运行的特点，对光学元件的抗激光损伤能力等提出了很严格的要求。光学元件内部的微缺陷会引起局域场强增强，采用时域有限差分方法对亚波长量级的缺陷进行了电磁场的数值模拟，并对数值计算的参数选取给出了定量的判断。     

光机靶室钢混结构基于动刚度等效的简化建模

钢混结构分离式有限元建模方法由于考虑了钢筋的实际分布,其计算精度较高,但建模过程复杂繁琐。开展了基于动刚度等效的钢混结构简化建模研究,该方法将钢筋和混凝土考虑为一种等效材料。研究中首先以典型钢混组合结构为例,以等效模型与分离式模型的模态特性等效为目标,建立了分部件的动刚度等效模型,然后将各部件等效模型组合成整体结构等效模型,数值仿真结果显示在相同载荷作用下,整体参考模型与整体等效模型的响应分布非常接近,而等效模型的规模远小于参考模型,说明了该等效建模方式的有效性。进一步将该等效建模方法应用于神光Ⅲ主机靶室钢混结构的动力学分析中,并计算了该巨型光机结构等效模型在实测基础加速度激励下的随机振动响应,靶球上各方向的响应实测结果与计算结果较为吻合。所采用的动刚度等效建模方法可应用在巨型光机结构动力学稳定性的进一步研究中。