Implicit and explicit integration in the solution of the absolute nodal coordinate differential/algebraic equations

This investigation is concerned with the use of an implicit integration method with adjustable numerical damping properties in the simulation of flexible multibody systems. The flexible bodies in the system are modeled using the finite element absolute nodal coordinate formulation (ANCF), which can be used in the simulation of large deformations and rotations of flexible bodies. This formulation, when used with the general continuum mechanics theory, leads to displacement modes, such as Poisson modes, that couple the cross section deformations, and bending and extension of structural elements such as beams. While these modes can be significant in the case of large deformations, and they have no significant effect on the CPU time for very flexible bodies; in the case of thin and stiff structures, the ANCF coupled deformation modes can be associated with very high frequencies that can be a source of numerical problems when explicit integration methods are used. The implicit integration method used in this investigation is the Hilber–Hughes–Taylor method applied in the context of Index 3 differential-algebraic equations (HHT-I3). The results obtained using this integration method are compared with the results obtained using an explicit Adams-predictor-corrector method, which has no adjustable numerical damping. Numerical examples that include bodies with different degrees of flexibility are solved in order to examine the performance of the HHT-I3 implicit integration method when the finite element absolute nodal coordinate formulation is used. The results obtained in this study show that for very flexible structures there is no significant difference in accuracy and CPU time between the solutions obtained using the implicit and explicit integrators. As the stiffness increases, the effect of some ANCF coupled deformation modes becomes more significant, leading to a stiff system of equations. The resulting high frequencies are filtered out when the HHT-I3 integrator is used due to its numerical damping properties. The results of this study also show that the CPU time associated with the HHT-I3 integrator does not change significantly when the stiffness of the bodies increases, while in the case of the explicit Adams method the CPU time increases exponentially. The fundamental differences between the solution procedures used with the implicit and explicit integrations are also discussed in this paper.     

炸药冲击响应的二维细观离散元模拟

用离散元法在颗粒尺度模拟了塑料粘接炸药的冲击响应。炸药的细观结构基于Voronoi拼图建成。以含孔洞和不含孔洞的炸药模型为例,说明细观结构对热点的形成和分布有重要影响。在这些算例中,粘塑性变形是温度局域化的主要原因,而摩擦和体积粘性等机制可以忽略。     

弯曲波对含多圆孔薄板的散射与动应力集中

采用复变函数法和多极坐标方法,研究了弯曲波对含有多圆孔薄板的散射问题。通过板的弯曲波动方程和内力方程的推导,求出在入射弯曲波条件下该问题的一般解的函数逼近序列和边界条件的表达式。用展开正交函数的方法将待解的问题归结为对一组无穷代数方程组的求解。最后,给出了含3圆孔薄板的孔边动应力集中系数的结果,并分析了孔间距和波数对动应力分布的影响。     

一种磨损预测的优化算法研究

根据机械部件磨损机理复杂、磨损量预测难精确的特点,提出基于免疫粒子群参数优化的最小二乘支持向量机方法预测磨损量.该算法采用免疫粒子群优化最小二乘支持向量机建模参数,避免了算法陷入局部最优解,实现了精确度高、泛化能力强的磨损量预测模型.对轴承钢试件磨损进行了试验研究,试验数据分析结果表明,基于免疫粒子群的最小二乘支持向量机预测方法优于前向反馈神经网络算法、遗传算法及蚁群算法,预测误差较小,具有很好的预测能力.     

Mathematic modeling and characteristic analysis for dynamic system with asymmetrical hysteresis in vibratory compaction

We investigate dynamic characteristics of vibratory compaction system with asymmetrical hysteresis. An asymmetrical model derived from Bouc-Wen differential equation is employed to describe hysteretic behavior of vibration engineering. A practical polynomial expression for hysteretic restoring force is deduced to be substituted into standard equation of the system, assuming that the non-linearity of the restoring force is weak. An asymptotic method, which combines Krylov-Bogolyubov-Mitropolsky (KBM) method with harmonic balance (HB) method, is applied to analyze steady-state responses of the asymmetrical hysteretic system subjected to harmonic excitation. Dynamic responses, such as the restoring force time histories and frequency responses of the system for the first order approximate, are obtained. Furthermore, numerical solution obtained using Runge-Kutta method as well as results of experiments (asphalt compaction on the Beijing-Fuzhou highway) are compared with the asymptotic solution. These results investigated that asymmetrical hysteretic model and the methods applied in this paper are quite appropriate for engineering applications.     

Effect of the accommodation coefficient on transfer processes in a supersonic rarefied gas flow around a transversely aligned cylinder

A supersonic flow around a cylinder is studied by the direct statistical Monte Carlo method in a wide range of rarefaction: from regimes close to continuum to free-molecular flow. The effect of the accommodation coefficient on the flow near the cylinder and on heat transfer between the gas and the cylinder is examined. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 64–72, May–June, 2008.     

Stability of Poiseuille flow in the presence of a longitudinal magnetic field

The stability of the plane flow of an electrically conducting fluid with respect to small perturbations was studied at large Reynolds numbers in the presence of a longitudinal magnetic field. The dependence of the critical Reynolds number on the electrical conductivity is investigated. At large Reynolds numbers, a new branch of instability and a sudden change in the critical Reynolds numbers is found. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 45–53, May–June, 2008.     

Homoclinic orbits for some （2＋1）-dimensional nonlinear Schrodinger-like equations

Chaos is closely associated with homoclinic orbits in deterministic nonlinear dynamics. In this paper, analytic expressions of homoclinic orbits for some （2＋1）- dimensional nonlinear Schrodinger-like equations are constructed based on Hirota＇s bilinear method, including long wave-short wave resonance interaction equation, generalization of the Zakharov equation, Mel＇nikov equation, and g-Schrodinger equation are constructed based on Hirota＇s bilinear method.     

????й????????????????????????

在对室内液化石油气泄漏流动过程进行量级估算分析的基础上,提出了一种预测室 内泄漏液化石油气浓度分布的简化方法. 并结合一次居民室内液 化气爆炸事故,讨论了此方法在实际应用中的简易可行性.     

Application of nonsmooth transformations to analyze a vibroimpact duffing system

The forced vibrations of an impact Duffing oscillator with one degree of freedom are studied using Zhuravlev's nonsmooth transformation and the multiple-scales method. Detailed consideration is given to the forced vibrations in the domain of resonance family. The stability and bifurcations of periodic vibrations are analyzed Translated from Prikladnaya Mekhanika, Vol. 44, No. 10, pp. 114–121, October 2008.