大学物理《力学》中的自由度分析

对《力学》中的物体自由度进行多方面分析,以深化教学、提高学生正 确分析物理问题的能力.使用实际教学分析的研究方法,在《力学》范围内讨论自由度与坐标、 自由与约束的关系并得以下结论: (1) 同一物体的自由度随其所在的``空间'不同而不同, 不因坐标系的选取不同而 异, 在同类参考系中不因参考系的动静而有别;(2)自由度遵循叠加原理. 讨论了质点系的总自由度及相关计算问题,并指出研究《力学》中自由度的意义.     

Direct modeling of flow of FENE fluids

Direct simulations of macromolecular fluids are carried out for flows between parallel plates and in expanding and contracting channels. The macromolecules are modeled as FENE dumbbells with soft disks or Lennard-Jones dumbbell-dumbbell interactions. The results are presented in terms of profiles and contour plots of velocity, pressure, temperature, density, and flow fields. In addition the data for potential energy, shear stress, and the normal components of the stress tensor are collected. In general, an excellent agreement is found between the simulated profiles and the well-known flow structures, such as flow separation and formation of viscous eddies, indicating that micro-hydrodynamics is a viable tool in linking macroscopic phenomena with the underlying physical mechanisms. The simulations are performed in the Newtonian regime, for medium-size systems comprising up to 3888 dumbbells. This number is sufficiently large to control boundary and particle number effects. The flow is induced by gravity. The traditional stochastic (thermal) and periodic boundary conditions are employed. Also, diffusive boundary conditions, which could include a stagnant fluid layer and repulsive potential walls, are developed. The scaling problems, which are related to the application of a large external force in a microscopic system (of the size of the order 100 Å), result in extreme pressure and temperature gradients. In addition, the viscosity and thermal conductivity coefficients obtained from velocity and temperature profiles of the channel flow are presented. These results are confirmed independently from modeling of Couette flow by the SLLOD equations of motion and from the Evans algorithm for thermal conductivity.     

圆柱壳的非线性自由振动分析

本文分析了各向同性封闭圆柱壳的非线性自由振动。文中采用经典的非线性弹性力学方法推导了圆柱壳的大振幅运动方程,这些方程的静态形式与冯·卡门的板理论方程具有同样的精度。文中讨论了四种基本振动模态,并且还以数学公式的形式给出了一般的最终结果,一些例子以曲线给出结果,并进行了比较。结果还表明线性振动可以作为非线性振动的一种特例。     

长杆弹对混凝土的侵爆效应

对适用于长杆弹在混凝土介质中侵彻和爆炸全过程的三维数值模拟方法和技术进行了研究。描述了进行相应数值模拟的有关方法和关键技术。确定了靶体C30混凝土材料所使用的本构模型及其相应的参数。对卵形头长杆弹在C30混凝土中侵彻到一定深度再爆炸的全物理过程进行了三维数值模拟,分别给出了C30混凝土靶体在侵彻和爆炸作用下的破坏效应图像。将侵彻计算图像与实验结果进行了比较,两者定性符合。     

Identification of regions of fastest mixing in a system of point vortices

This paper describes a numerical method for efficiently identifying the regions of fastest mixing of a passive dye in a flow due to a system of point vortices. Results obtained from computations are presented for systems of three and four point vortices, both in the unbounded domain and inside a circular cylinder. The flow is two‐dimensional and the fluid is incompressible. The regions where mixing is possible are found by studying the largest Lagrangian Lyapunov exponent distribution with respect to various initial positions of tracer particles. The regions of fastest mixing are then identified from the Lyapunov exponent distribution at small times. The results of the method are verified by quantifying the mixing by using a traditional box counting technique. The technique is then applied to several different initial configurations of vortices and some interesting results are obtained. Some numerical findings about the nature of the exponents computed are also discussed. Copyright © 2002 John Wiley Sons, Ltd.     

Natural harmonic oscillations of a heavy fluid in basins of complex shape

The natural harmonic oscillations of a heavy fluid in uniform-depth basins of complex shape, including those with three or more symmetry axes, are investigated within the shallow-water approximation. The waves are assumed to be gently sloping. The modes found are compared with the similar modes in an elliptic basin representing the class of basins with two symmetry axes. The mode characteristics associated with the number of basin symmetry axes and the basin shape are explored. Basins with two symmetry axes whose shape differs considerably from the elliptic, in particular, nonconvex basins, are considered. Both rotating and non-rotating basins are studied. The possibility of approximating the amplitudes of certain rotating-basin mode classes by Bessel functions is discussed.     

Nonlinear study of the dynamic behavior of a string-beam coupled system under combined excitations

In this paper,the nonlinear dynamic behavior of a string-beam coupled system subjected to external,parametric and tuned excitations is presented.The governing equations of motion are obtained for the nonlinear transverse vibrations of the string-beam coupled system which are described by a set of ordinary differential equations with two degrees of freedom.The case of 1:1 internal resonance between the modes of the beam and string,and the primary and combined resonance for the beam is considered.The method of multiple scales is utilized to analyze the nonlinear responses of the string-beam coupled system and obtain approximate solutions up to and including the second-order approximations.All resonance cases are extracted and investigated.Stability of the system is studied using frequency response equations and the phase-plane method.Numerical solutions are carried out and the results are presented graphically and discussed.The effects of the different parameters on both response and stability of the system are investigated.The reported results are compared to the available published work.     

Experimental validation of the constitutive equations of thermoelastoplasticity including the third invariant of the stress deviator

The equations describing thermoelastoplastic deformation along nonstraight paths and taking into account the third invariant of the stress deviator are experimentally validated. The equations contain two scalar functions that are specified in base tests on tubular specimens. The test data are tabulated. The values of the scalar functions for strains, temperature, and stress mode are found by using nonlinear interpolation of the data and the temperature similarity of the functions. The stresses in elements of the body are calculated from given strains by the method of successive approximations     

Multiscale modeling of fracture of MgO: Sensitivity of interatomic potentials

Three different interatomic potentials, namely, B-G I Model, B-G II Model and L-C Model, are used in multiscale modeling and simulation of a center-cracked specimen made of magnesia subjected to monotonically increasing loading. The specimen is decomposed into a far field, a near field and a crack-tip region. The analytical solution in the far field from linear elastic fracture mechanics (LEFM) is utilized. The solution of the near field is based on a multiscale field theory. In the crack-tip region, molecular dynamics (MD) simulation is employed. These methodologies are integrated to simulate mixed mode fracture of magnesia (MgO). Three different interatomic potentials are examined and the interatomic potential and interatomic force between Mg-Mg, Mg-O and O-O are shown. The numerical results of crack propagation demonstrate that (1) crack closure is witnessed in B-G I Model but not in B-G II Model and L-C Model, (2) B-G II Model and L-C Model diverge in the early stage. The cause of instability and the remedy are also discussed.     

张拉膜结构力密度法混合找形分析

论述了张拉膜结构力密度法混合找形的基本理论;所谓力密度法混合找形,即部分单元力密度控制,部分单元弹性控制;力密度控制采用线性求解,弹性控制采用非线性求解,通过迭代计算混合找形求出各结点的坐标值。据此编制了相应的计算软件;对工程实例进行了验算,结果表明,本文给出的计算结果与德国著名软件easy的计算结果相吻合。     

矩形板非线性蠕变损伤屈曲特性分析

基于Von Karman非线性板理论和Kachanov-Rabotnov损伤理论,建立了在横向和面内载荷共同作用下考虑蠕变损伤效应的矩形板的非线性控制平衡方程,采用有限差分法和时间增量算法对未知变量进行离散,对整个问题进行迭代求解,分析了几何非线性、荷载等因素对板非线性蠕变损伤特性的影响。     

Numerical analysis of nonlinear modes of piecewise linear systems torsional vibrations

The nonlinear modes of essentially nonlinear piecewise-linear finite degrees of freedom systems are calculated by the numerical methods, which are suggested in this paper. The basis of these methods is numerical solutions of the equations of the systems motions in configuration space. The numerical method for the nonlinear modes of essentially nonlinear piecewise-linear systems forced vibrations is suggested. The basis of this approach is the combination of the Rauscher method and the calculations of the autonomous system nonlinear modes. The nonlinear modes of the diesel engine transmission torsional vibrations are analyzed numerically. The vibrations are described by essentially nonlinear piecewise-linear system with fifteen degrees of freedom. The NNMs of this system forced vibrations are observed in the resonance regions. Both NNMs and the motions, which are essentially differ from NNMs, are observed in the distance from the resonances. NNMs of the forced vibrations of the systems with dissipation are close to NNMs of the system without dissipation.     

A review of autonomous agricultural vehicles (The experience of Hokkaido University)

Robotic farming will play an undeniably significant role in future sustainable agriculture. Autonomous agricultural vehicles for arable crops and their components are reviewed herein, and their differing possible components, advantages and disadvantages are discussed. The autonomous agricultural vehicles are qualified from technical points of view, including each vehicle's hardware unit (platform development, platform, transporter system, operation functions, communication, sensors [positional sensor, attitude sensor, and safety sensor], and control unit), the physical environment, and the control algorithm. The development process for different vehicles is described. As an operational case study, all investigations conducted at Hokkaido University between 1990 and 2018 are discussed and evaluated based on the described classification. The development procedure, the process of component selection, developmental challenges, and the performance indicators of the vehicles are discussed in detail. The development process and applicability of each component are then presented, and recommendations for future studies are noted. Finally, the most important experiences and lessons, and some recommendations are outlined.     

Non-linear temporal dynamics of two-mode interactions of magnetized flow

This paper considers, in the frame work of the model of two superposed layers of viscous-potential incompressible magnetic fluids, the problem on formation of resonant waves of two modes on the interface between fluids that arisen as a result of second-harmonic resonance. The fluids moving with uniform velocities parallel to their interface are stressed by a tangential magnetic field. The analysis includes the linear, as well as the non-linear effects where the analytical solutions are constructed using the method of multiple scales, in both space and time, and hence the solvability conditions correspond to the uniform (convergent) solutions are obtained. The solvability conditions are then exploited to derive a more general system of non-linear partial differential equations with complex coefficients governing the amplitudes of the resonant waves. These equations are examined for solutions corresponding to sinusoidal wavetrains consequently different kinds of instabilities are demonstrated. The stability criterion in each case is derived and discussed both analytically and graphically.     

Symmetries of boundary layer equations of power-law fluids of second grade

A modified power-law fluid of second grade is considered. The model is a combination of power-law and second grade fluid in which the fluid may exhibit normal stresses, shear thinning or shear thickening behaviors. The equations of motion are derived for two dimensional incompressible flows, and from which the boundary layer equations are derived. Symmetries of the boundary layer equations are found by using Lie group theory, and then group classification with respect to power-law index is performed. By using one of the symmetries, namely the scaling symmetry, the partial differential system is transformed into an ordinary differential system, which is numerically integrated under the classical boundary layer conditions. Effects of power-law index and second grade coefficient on the boundary layers are shown and solutions are contrasted with the usual second grade fluid solutions.     

Thermal bending of thick rectangular plates of bimodulus composite materials

Wiighted residual solutions are presented for thermal bending of laminated composite plates. The material of each layer is assumed to be clastically and thermoelastically orthotropic and bimodular. The formulations are based on the thermoclastic version of the theory of Whitney-Pagano laminated plate, which includes thickness shear detormations. The results are obtained for deflections and neutral-surface positions and are found to be in good agreement with the closed-form solution.     

Skew-gradient representations of constrained mechanical systems

The characteristics of stationary and non-stationary skew-gradient systems are studied. The skew-gradient representations of holonomic systems, Birkhoffian systems, generalized Birkhoffian systems, and generalized Hamiltonian systems are given. The characteristics of skew-gradient systems are used to study integration and stability of the solution of constrained mechanical systems. Examples are given to illustrate applications of the result.     

A study of the perforation of stiffened plates by rigid projectiles

In the present paper, a four-stage perforation model that accurately predicts the residual velocity is developed by adopting an energy method. The four stages are plug formation, dishing formation, petal formation and projectile exit. In addition, some important experimental results are presented and analyzed to validate the present perforation model. In the experiments, high speed camera system is used to record the perforation process. Observations on target damage and measurements of initial velocities and residual velocities with the aid of the system are presented. Numerical simulations are carried out for projectiles against single and layered plates adopted in the experiments. The perforation process is studied and the deformation and failure modes are obtained. The predictions of numerical simulations and analytical model are found in reasonably good agreement with those of experiments, and can be used to predict the ballistic limit and residual velocity of stiffened plates perforated by rigid projectiles. The project supported by the National Natural Science Foundation of China (90305018) The English text was polished by Yunming Chen     

预扭转钨合金杆弹侵彻钢靶的数值模拟

在实验的基础上 ,对预扭转钨合金长杆弹提出了一个唯象的本构模型和破坏准则 ,并针对预扭转钨合金长杆弹侵彻厚钢靶进行二维有限元计算 ,得到了侵彻过程中的主要物理图象和曲线 ,进而对计算结果进行了分析讨论。计算结果表明 ,侵彻过程中的主要物理数据与实验测量结果基本一致。     

预处理方法在低速弹箭流场计算中的应用

基于Weiss-Smith预处理矩阵和全局截断预处理参数,采用有限体积方法对雷诺平均Navier-Stokes方程进行离散。对流项离散采用二阶线性重构和AUSM ＋-up格式,时间推进方法采用多重网格下的LU-SGS方法。结合M PI消息传递方法,建立了一套计算低速流动的并行数值方法。计算了低速椭球体的流场和气动力,压力系数和切应力系数计算结果与文献实验结果对比吻合度较好。生成了末敏弹的流场计算网格,对绕末敏弹流场进行了数值模拟。对多重网格下多进程的加速比和并行效率进行了测试,显示了程序良好的并行效率。计算的气动力结果与实验结果吻合。综合结果表明：本文的数值方法能够用于低速弹箭流场和气动力计算,为新型弹箭的设计和定型提供保证。