Perfectly matched layer boundary condition for two-dimensional Euler equations in generalized coordinate system

In the past, perfectly matched layer (PML) equations have been constructed in Cartesian and spherical coordinates. In this article, the focus is on the development of a PML absorbing technique for treating numerical boundaries, especially those with unbounded domains, in a generalized coordinate system for a flow in an arbitrary direction. The PML equations for two-dimensional Euler equations are developed in split form through a space–time transformation involving a complex variable transformation with the application of a pseudo-mean-flow in the PML domain. A numerical solver is developed using conventional numerical schemes without employing any form of filtering or artificial dissipation to solve the governing PML equations for two-dimensional Euler equations in a generalized coordinate system. Physical domains of arbitrary shapes are considered and numerical simulations are carried out to validate and demonstrate the effectiveness of the PML as an absorbing boundary condition in generalized coordinates.     

局部自然坐标与斜交笛卡儿坐标的对等性

本文通过分析法和几何法证明：在有限元法和边界元法中，平面三角形单元的面积坐标或四面体单元的体积坐标与由这些单元的二边或三边构成的斜交笛卡儿二维坐标或三维坐标，尽管定义不同，而其内容却是完全对等的。这样，完善地解决了单元中任意点在整体坐标与局部坐标的表示问题和基于局部自然坐标的面积分和体积分的数值积分的依据问题。     

考虑空气阻力和地球自转的远程抛射体的研究

针对地球表面坐标系计算远程抛射体落地时间和位置的不足,通过两种方法获得远程抛射体的落地信息. 首先通过在地球表面建立空间直角坐标系计算出考虑空气阻力和地球自转情况下任意抛射角的抛射体的解析解,然后获得了在地球表面所建立的空间直角坐标系与地心坐标系的转换关系式,以抛射体与地心的距离为地球半径作为落地点的条件来计算落地时间和落地位置. 另外,采用球面坐标系表示抛射体的位置、速度及加速度,通过数值计算方法获得抛射体的落地信息. 两种方法所得的结果比较接近,因此两种方法来获得远程抛射体的落地信息都是可行的.     

Visualization of three-dimensional incompressible flows by quasi-two-dimensional divergence-free projections in arbitrary flow regions

A visualization of three-dimensional incompressible flows by divergence-free quasi-two-dimensional projections of the velocity field onto three coordinate planes is revisited. An alternative and more general way to compute the projections is proposed. The approach is based on the Chorin projection combined with a SIMPLE-like iteration. Compared to the previous methodology based on divergence-free Galerkin–Chebyshev bases, this technique, formulated in general curvilinear coordinates, is applicable to any flow region and allows for faster computations. To illustrate this visualization method, examples in Cartesian and spherical coordinates, as well as post-processing of experimental 3D-PTV data, are presented.     

WGS84椭球下的UTM坐标与Clarke80椭球下的兰勃特坐标转换方法研究

介绍了两种椭球下的两种投影之间的坐标转换方法.在进行坐标转换时,WGS84椭球下的UTM坐标首先被转换成地理坐标,然后被转换成笛卡尔坐标;采用7参数法,可得到Clarke80椭球下的笛卡尔坐标,通过相应的坐标变换模型便可计算出该椭球下的大地坐标以及相应的兰勃特坐标.在给定了坐标系之间的7个转换参数以及单标准线的兰勃特投影(1SP)参数时,采用本方法WGS84椭球下的UTM和大地坐标可以很容易地被转换成Clarke80椭球下的大地坐标和兰勃特坐标,其转换精度优于0.1 m,适用于很多应用领域.     

两种坐标系下惯导传递对准效果比较

分别建立了发射点惯性坐标系和当地地理坐标系下的惯导误差传播模型,以及典型匹配模式下的传递对准基本方程。基于运载火箭主动段飞行,分析了姿态机动运动对两种坐标系下传递对准的不同影响。数学模型和仿真分析表明,两种坐标系下的传递对准在原理和系统特性上均有较大差异。典型匹配方案的对准效果对比表明,发射点惯性坐标系内的"角速度+加速度"匹配是一种能够实现子惯导快速、准确初始化的有效方法。     

Comparative analysis of the electroelastic thickness vibrations of layers with curved boundaries

The three-dimensional equations of electroelasticity in Cartesian, cylindrical, and spherical coordinates are represented in Hamiltonian form with respect to the thickness coordinate. The boundary-value problem with a harmonic potential difference and zero mechanical load given on the boundaries is solved numerically. The amplitude–frequency characteristics and natural frequencies are compared. The resonant and antiresonant frequencies of the current and the dynamic electromechanical coupling coefficient are determined     

Measurement of curved-surface deformation in cylindrical coordinates

Stereo vision is used to measure the strain field of a round tension test specimen in a cylindrical coordinate system. Initially, the displacement fields of the specimen are measured relative to a world coordinate system erected by the stereo vision. Through coordinate transformations, the measured displacement fields expressed in world coordinates are then converted to the displacement fields expressed in cylindrical coordinates. By differentiating the axial and circumferential displacements in the axial and circumferential directions, the axial, circumferential and shear strains are determined. Results indicate that the measured mean value of the axial strains is in good agreement with the measurements of the extensometer and the strain gage. The Poisson's ratio obtained by the circumferential and axial strains is close to .33 in the elastic state. The mean error of the computed shear strain is approximately .03 percent in the smaller elastic deformation and .08 percent in the larger plastic deformation.     

Analytical solution for the time-fractional heat conduction equation in spherical coordinate system by the method of variable separation

In this paper,using the fractional Fourier law,we obtain the fractional heat conduction equation with a time-fractional derivative in the spherical coordinate system.The method of variable separation is used to solve the timefractional heat conduction equation.The Caputo fractional derivative of the order 0 < α≤ 1 is used.The solution is presented in terms of the Mittag-Leffler functions.Numerical results are illustrated graphically for various values of fractional derivative.     

一种提高四边形四节点平面壳单元计算精度的新方法

平面壳单元是由平面应力单元和平板弯曲单元叠加组合而成,具有简单的理论表达,但是它在计算曲面壳体结构时误差较大。为了进一步提高平面壳单元的计算精度,本文提出了一种计算平面壳单元刚度矩阵的新方法。通过该方法在高斯积分点建立多个单元局部坐标系,并保证每个局部坐标系都位于单元在高斯点处的切平面上,从而可以有效适应曲面壳体形状,达到进一步提高平面壳单元计算精度的目的。为了在这种新坐标系下计算单元刚度矩阵,给出了求解形函数对局部坐标的导数、局部到自然坐标系积分转换的雅可比、以及局部到整体坐标系的转换矩阵的新型计算方法。通过将这些新坐标系以及新计算方法运用到平面壳单元DKQ24中,可以有效提高平面壳单元尤其是在计算曲面壳体时的精度。计算结果表明,本文方法和平面壳单元相结合,不仅具有平面壳单元简单的理论表达式,还能得到满意的精度。另外,本文方法还可以应用到其他类型的平面壳单元,为提高其他类型平面壳单元的计算精度提供了一种新的途径和思路,具有广阔的应用前景。     

On the formulation of the planar ANCF triangular finite elements

In this paper, new planar isoparametric triangular finite elements (FE) based on the absolute nodal coordinate formulation (ANCF) are developed. The proposed ANCF elements have six coordinates per node: two position coordinates that define the absolute position vector of the node and four gradient coordinates that define vectors tangent to coordinate lines (parameters) at the same node. To shed light on the importance of the element geometry and to facilitate the development of some of the new elements presented in this paper, two different parametric definitions of the gradient vectors are used. The first parametrization, called area parameterization, is based on coordinate lines along the sides of the element in the reference configuration, while the second parameterization, called Cartesian parameterization, employs coordinate lines defined along the axes of the structure (body) coordinate system. The fundamental differences between the ANCF parameterizations used in this investigation and the parametrizations used for conventional finite elements are highlighted. The Cartesian parameterization serves as a unique standard for the triangular FE assembly. To this end, a transformation matrix that defines the relationship between the area and the Cartesian parameterizations is introduced for each element in order to allow for the use of standard FE assembly procedure and define the structure (body) inertia and elastic forces. Using Bezier geometry and a linear mapping, cubic displacement fields of the new ANCF triangular elements are systematically developed. Specifically, two new ANCF triangular finite elements are developed in this investigation, namely four-node mixed-coordinate and three-node ANCF triangles. The performance of the proposed new ANCF elements is evaluated by comparison with the conventional linear and quadratic triangular elements as well as previously developed ANCF rectangular and triangular elements. The results obtained in this investigation show that in the case of small and large deformations as well as finite rotations, all the elements considered can produce correct results, which are in a good agreement if appropriate mesh sizes are used.     

A derivation of the Mehrabadi-Cowin equations

An alternative derivation to that given by Mehrabadi and Cowin (1978) is presented here for a pair of kinematic equations governing a certain class of flows in the plastic deformation of dilatant granular materials. This class has been described by Spencer (1981) as double shearing flows. In their derivation Mehrabadi and Cowin (1978), prior to presenting the equations relative to rectangular Cartesian coordinates, obtained an intermediate pair of equations relative to a non-orthogonal network of characteristic coordinates. The essential difference between the original and present derivation is that here, the flow rule, expressed relative to rotating, rectangular Cartesian coordinates, is transformed directly to obtain the kinematic equations relative to fixed rectangular Cartesian coordinate axes, without the need to obtain the characteristic equations.     

One-dimensional piston process in strong gravitational field

In this paper, the dynamic process driven by one-dimensional piston in strong gravitational field was studied on the Cartesian, cylinderical and spherical coordinates. The gasdynamic equations were numerically solved by the characteristic method. The solution which satisfies the velocity condition at piston and the boundary conditions connect the flow region and the quiet region is obtained. The present paper analyses especially the influence of coordinate systems on the field of compressible flow, uniform flow and rarefaction flow region, the shock velocity and the temperature distribution at the piston.     

平面问题极坐标下几何方程的一种解析法

经典弹性力学教材中都是通过单元的几何变形来推导平面问题极坐标下的几何方 程. 这里, 用解析法从直角坐标系下的几何方程推出极坐标下的几何方程.     

ANALYTICAL SOLUTIONS FOR THE LAYERED GEO-MATERIALS SUBJECTED TO AN ARBITRARY POINT LOAD IN THE CARTESIAN COORDINATE

This paper presents analytical solutions for the stress and displacement field in elastic layered geo-materials induced by an arbitrary point load in the Cartesian coordinate system. The point load solutions can be obtained by referring to the integral transform and the transfer matrix technique. However, former solutions usually exist in the cylindrical coordinate system subjected to axisymmetric loading. Based on the proposed solutions in the Cartesian coordinate, it is very easy to solve asymmetric problems and consider the condition with internal loads in multi-layered geo-materials. Moreover, point load solutions can be used to construct solutions for analytical examination of elastic problems and incorporated into numerical schemes such as boundary element methods. The results discussed in this paper indicate that there is no problem in the evaluation of the point load solutions with high accuracy and efficiency, and that the material non-homogeneity has a significant effect on the elastic field due to adjacent loading.     

Strain gradient theory in orthogonal curvilinear coordinates

In this short note, general formulations of the Toupin–Mindlin strain gradient theory in orthogonal curvilinear coordinate systems are derived, and are then specified for the cases of cylindrical coordinates and spherical coordinates. Expressions convenient for practical use are presented for the corresponding equilibrium equations, boundary conditions, and the physical components for strains and strain gradients in the two coordinate systems. The results obtained in this paper are general and complete, and can be useful for a wide range of applications, such as asymptotic crack tip field analysis, cylindrical and spherical cavity expansion problems, and the interpretation of micro/nano indentation tests and bending/twisting tests on small scales.     

New mixed plate-bending elements with shear strain interpolations

This paper is concerned with two mixed plate-bending elements with shear strain interpolations, a quadrilateral element and an 8-node serendipity-type element based on discussions on the element proposed in Ref.[1]. The shear strains and inner-forces in the natural coordinates are interpolated in an element and then transformed into Cartesian coordinates in accordance with covariant and contravariant tensor transformations, respectively. The quadrilateral element coincides with the element in Ref.[1] when it is rectangular. Numerical examples show that the two new elements are free from shear locking and spurious kinematic modes under regular and irregular meshes and have the advantages of being insensitive to element distortion and able to give fairly accurate results.The Project supported by National Natural Science Foundation of China.     

On Generating Isochoric Finite Deformations

New isochoric finite deformations may be generated from any such deformation described in rectangular Cartesian coordinates by changing coordinate systems.      

A velocity transformation method for the nonlinear dynamic simulation of railroad vehicle systems

The solution of the constrained multibody system equations of motion using the generalized coordinate partitioning method requires the identification of the dependent and independent coordinates. Using this approach, only the independent accelerations are integrated forward in time in order to determine the independent coordinates and velocities. Dependent coordinates are determined by solving the nonlinear constraint equations at the position level. If the constraint equations are highly nonlinear, numerical difficulties can be encountered or more Newton–Raphson iterations may be required in order to achieve convergence for the dependent variables. In this paper, a velocity transformation method is proposed for railroad vehicle systems in order to deal with the nonlinearity of the constraint equations when the vehicles negotiate curved tracks. In this formulation, two different sets of coordinates are simultaneously used. The first set is the absolute Cartesian coordinates which are widely used in general multibody system computer formulations. These coordinates lead to a simple form of the equations of motion which has a sparse matrix structure. The second set is the trajectory coordinates which are widely used in specialized railroad vehicle system formulations. The trajectory coordinates can be used to obtain simple formulations of the specified motion trajectory constraint equations in the case of railroad vehicle systems. While the equations of motion are formulated in terms of the absolute Cartesian coordinates, the trajectory accelerations are the ones which are integrated forward in time. The problems associated with the higher degree of differentiability required when the trajectory coordinates are used are discussed. Numerical examples are presented in order to examine the performance of the hybrid coordinate formulation proposed in this paper in the analysis of multibody railroad vehicle systems.     

Comparative study on multibody vehicle dynamics models based on subsystem synthesis method using Cartesian and joint coordinates

The subsystem synthesis method has been developed in order to improve computational efficiency for a multibody vehicle dynamics model. Using the subsystem synthesis method, equations of motion of the base body and each subsystem can be solved separately. In the subsystem synthesis method, various coordinate systems can be used and various integration methods can be applied in each subsystem, as long as the effective mass matrix and the effective force vector are properly produced. In this paper, comparative study has been carried out for the subsystem synthesis method with Cartesian coordinates and with joint relative coordinates. Two different integration methods such as an explicit integrator and an explicit implicit integrator are employed. In order to see the accuracy and computational efficiency from the different models based on the different coordinate systems and different integration methods, a rough terrain run simulations has been carried out with a 6 × 6 off-road multibody vehicle model.