Feature Extraction and Recognition for Rolling Element Bearing Fault Utilizing Short-Time Fourier Transform and Non-negative Matrix Factorization

Due to the non-stationary characteristics of vibration signals acquired from rolling element bearing fault, thc time-frequency analysis is often applied to describe the local information of these unstable signals smartly. However, it is difficult to classitythe high dimensional feature matrix directly because of too large dimensions for many classifiers. This paper combines the concepts of time-frequency distribution（TFD） with non-negative matrix factorization（NMF）, and proposes a novel TFD matrix factorization method to enhance representation and identification of bearing fault. Throughout this method, the TFD of a vibration signal is firstly accomplished to describe the localized faults with short-time Fourier transform（STFT）. Then, the supervised NMF mapping is adopted to extract the fault features from TFD. Meanwhile, the fault samples can be clustered and recognized automatically by using the clustering property of NMF. The proposed method takes advantages of the NMF in the parts-based representation and the adaptive clustering. The localized fault features of interest can be extracted as well. To evaluate the performance of the proposed method, the 9 kinds of the bearing fault on a test bench is performed. The proposed method can effectively identify the fault severity and different fault types. Moreover, in comparison with the artificial neural network（ANN）, NMF yields 99.3% mean accuracy which is much superior to ANN. This research presents a simple and practical resolution for the fault diagnosis problem of rolling element bearing in high dimensional feature space.

Effects of rare earth element lanthanum on the microstructure of copper matrix diamond tool materials

Effects of rare earth element La on the microstructure of Cumatrix diamond tools were researched under the conditions of variousmaterials components and the process parameters in order to improvematerials properties. SEM, XPS and X-ray were used to investigate thefracture section, microstructure and the element valence inmaterials. The Results shown that the combination of rare earthelement La and transition element Ti is advantageous to the bondingstate Between diamond particles and matrix, so it can improve thematerials properties. Suitable sintering temperature is 790 deg. C.

Damage prediction for magnesium matrix composites formed by liquid-solid extrusion process based on finite element simulation

A damage prediction method based on FE simulation was proposed to predict the occurrence of hot shortness crocks and surface cracks in liquid-solid extrusion process. This method integrated the critical temperature criterion and Cockcroft ＆ Latham ductile damage model, which were used to predict the initiation of hot shortness cracks and surface cracks of products, respectively. A coupling simulation of deformation with heat transfer as well as ductile damage was carried out to investigate the effect of extrusion temperature and extrusion speed on the damage behavior of Csf/AZ91D composites. It is concluded that the semisolid zone moves gradually toward deformation zone with the punch descending. The amplitude of the temperature rise at the exit of die from the initial billet temperature increases with the increase of extrusion speed during steady-state extrusion at a given punch displacement. In order to prevent the surface temperature of products beyond the incipient melting temperature of composites, the critical extrusion speed is decreased with the increase of extrusion temperature, otherwise the hot shortness cracks will occur. The maximum damage values increase with increasing extrusion speed or extrusion temperature. Theoretical results obtained by the Deform^TM-2D simulation agree well with the experiments.

THE STRESS SUBSPACE OF HYBRID STRESS ELEMENT AND THE DIAGONALIZATION METHOD FOR FLEXIBILITY MATRIX H

The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix; 2) The equivalent assumed stress modes lead to the identical hybrid element. The Hilbert stress subspace of the assumed stress modes is established. So, it is easy to derive the equivalent orthogonal normal stress modes by Schmidt's method. Because of the resulting diagonal flexibility matrix, the identical hybrid element is free from the complex matrix inversion so that the hybrid efficiency, is improved greatly. The numerical examples show that the method is effective.

MatDEM-fast matrix computing of the discrete element method

Discrete element method can effectively simulate the discontinuity,inhomogeneity and large deformation and failure of rock and soil.Based on the innovative matrix computing of the discrete element method,the highperformance discrete element software MatDEM may handle millions of elements in one computer,and enables the discrete element simulation at the engineering scale.It supports heat calculation,multi-field and fluidsolid coupling numerical simulations.Furthermore,the software integrates pre-processing,solver,postprocessing,and powerful secondary development,allowing recompiling new discrete element software.The basic principles of the DEM,the implement and development of the MatDEM software,and its applications are introduced in this paper.The software and sample source code are available online(http://matdem.com).

Wigner Function：from Ensemble Average of Density Operator to Its Matrix Element in Entangled Pure States

We show that the Wigner function (an ensemble average of the density operator ρ, Δ is the Wigner operator) can be expressed as a matrix element of ρ in the entangled pure states. In doing so, converting from quantum master equations to time-evolution equation of the Wigner functions seems direct and concise. The entangled states are defined in the enlarged Fock space with a fictitious freedom.

THE GENERATION OF NON－LINEAR STIFFNESS MATRIX OF TRIANGLE ELEMENT WHENCONSIDERING BOTH THE BENDING AND IN－PLANEME MBRANE FORCES

Using Stricklin Melhod￣[5],we have this paper has derived the formulas for the ge-neration of non-linear element stiffness matrix of a triangle element when considering both the bending and the in-plane membrane forces. A computer programme for the calculation of large deflection and inner forces of shallow shells is designed on theseformulas. The central deflection curve computed by this programme is compared with other pertaining results.

Finite Element Analysis of Push-Out Test of SiC Fiber Reinforced Titanium Matrix Composites

Based on the interphase layer model and the spring layer model, an improved interface model was developed to evaluate the interfacial shear strength of Titanium matrix composites(TMCs) and to analyze the effects of various parameters on the interfacial properties. The results showed that the improved interface model is more suitable for calculating the interfacial properties of SiC fiber reinforced titanium matrix composites. The interfacial shear strength of SiC/Timetal-834 predicted is 500 MPa. In addition, in order to better understand the interfacial properties of composites, some push out phenomenon were analyzed.

Wigner Function:from Ensemble Average of Density Operator to Its One Matrix Element in Entangled Pure States

We show that the Wigner function W = Tr(△ρ) (an ensemble average of the density operator ρ, △ is theWigner operator) can be expressed as a matrix element of ρ in the entangled pure states. In doing so, converting fromquantum master equations to time-evolution equation of the Wigner functions seems direct and concise. The entangledstates are defined in the enlarged Fock space with a fictitious freedom.

Analytical Expressions of Matrix Elements of Physical Quantities for Dirac Oscillator

The analytical expressions of the matrix elements for physical quantities are obtained for the Dirac oscillator in two and three spatial dimensions. Their behaviour for the case, of operator's square is discussed in details. The two-dimensional Dirac oscillator has similar behavior to that for three-dimensional one.

Finite element model for arch bridge vibration dynamics considering effect of suspender length adjustment on geometry stiffness matrix

In this paper, we established a finite element (FEM) model to analyze the dynamic characteristics of arch bridges. In this model, the effects of adjustment to the length of a suspender on its geometry stiffness matrix are stressed. The FEM equations of mechanics characteristics, natural frequency and main mode are set up based on the first order matrix perturbation theory. Applicantion of the proposed model to analyze a real arch bridge proved the improvement in the simulation precision of dynamical characteristics of the arch bridge by considering the effects of suspender length variation.

THE DYNAMIC STIFFNESS MATRIX OF THE FINITE ANNULAR PLATE ELEMENT

The dynamic deformation of harmonic vibration is used as the shape functions of the finite annular plate element, and sonic integration difficulties related to the Bessel's functions are solved in this paper. Then the dynamic stiffness matrix of the finite annular plate element is established in closed form and checked by the direct stiffness method. The paper has given wide convcrage for decomposing the dynamic matrix into the power series of frequency square. By utilizing the axial symmetry of annular elements, the modes with different numbers of nodal diameters at s separately treated. Thus some terse and complete results are obtained as the foundation of structural characteristic analysis and dynamic response compulation.

Analytic Expression of Arbitrary Matrix Elements for Boson Exponential Quadratic Polynomial Operators

Making use of the transformation relation among usual, normal, and antinormal ordering for the multimode boson exponential quadratic polynomial operators （BEQPO＇s）I we present the analytic expression of arbitrary matrix elements for BEQPO＇s. As a preliminary application, we obtain the exact expressions of partition function about the boson quadratic polynomial system, matrix elements in particle-number, coordinate, and momentum representation, and P representation for the BEQPO＇s.

A Pragmatic Approach to Coping with Matrix Effects during ICP-MS Analysis of Trace Elements in Silicate Rocks and Calibration of REE Interferences

Operating an Agilent 7700X ICP-MS spectrometer under robust plasma conditions (1550 W) with a He-filled octopole collision cell and analysing solutions (?1 total dissolved solids) still suffered analyte peak suppression due to matrix effects. International reference rocks BCR-1, BHVO-1, AGV-1, G-2 and BCR-2 all showed count rate reductions for 36 elements (mass range 7Li to 238U) averaging ~10% but with no dependence on isotope mass. Use of an internal standard (103Rh) and/or using a ten-fold dilution of sample solutions reduced these effects but problems with reduced count rates combined with larger errors for some elements introduced other problems. The best approach was to normalise the count rates for each element in the other samples against those for BCR-1 as an external standard;thus the count suppression due to the matrix effect is corrected for each individual element. This approach provides standardization “traceability” in line with the ERM ISO/IEC requirement. Experiments are also reported on quantifying the proportions of Ba and selected REE oxide/hydroxide components versus parent isotopes (XO/X and XOH/X). This information is essential for correcting peak interferences on higher mass number REE for the rock samples, and equations are developed to use measured CeO/Ce and CeOH/Ce ratios to predict such values for any other member of the REE suite. Concentrations obtained show excellent agreement with recommended values for the international reference materials especially for the REE. Robust data are also provided for two other standard rocks: nepheline syenite STM-1 and quartz syenite CAAS-1;the latter shows exceptional enrichments of Zr, REE, Th, and U.

General calculation formulas and recurrence relations of radial matrix elements for relativistic n-dimensional hydrogen atom of spin S=0

In this paper, the general calculation formulas of radial matrix elements for relativistic n-dimensional hydrogen atom of spin S=0 are obtained, and the recurrence relation of different power order radial matrix elements are also derived.

考虑剪切变形的欧拉梁单元一致质量矩阵

为明确考虑剪切变形的欧拉梁单元一致质量矩阵推导方法,以形函数为基础,基于虚功原理,区分伸缩、扭转以及是否考虑剪切变形的弯曲受力状态,给出了欧拉梁单元一致质量矩阵表达式。研究表明:不考虑剪切变形时,欧拉梁受弯状态的一致质量分析中一般不计单元上微元体水平方向的惯性力,此时仅需弯曲变形引发的竖向位移形函数即可;考虑剪切变形时,则需计入水平方向的惯性力,且同时需要完整的竖向位移形函数和纯弯曲转角位移形函数;对于变截面欧拉梁单元,其一致质量矩阵表达式过于复杂,可在等截面梁的基础上,根据元素位置对矩阵元素匹配左右端或平均截面面积及截面极惯性矩近似给出实用的质量矩阵表达式;考虑剪切变形时,欧拉梁的刚度矩阵亦可经竖向位移和纯弯曲转角位移形函数计算得到,这一过程与使用纯弯曲竖向位移和纯剪切竖向位移形函数的过程在本质上是一致的。

NORMALIZED ELEMENT NODE TOPOLOGICAL MATRIX OF FINITE ELEMENT MESHES

On the basis of concept of element node topological analysis, the normalized element node topological matrices for finite element meshes are presented in the paper, including 3-node and 6-node triangular element, 4-node and 8-node quadrilateral element, 8-node and 20-node hexahedral element. It is beneficial to further analyzing topological characteristics of finite element models and automatic generation of meshes

超声振动辅助研磨SiCp/Al材料去除特性的研究

SiCp/Al因具有高强度,低密度,高耐磨性等优点被应用在许多领域,如航空航天,汽车和电子领域等。但SiCp/Al是典型的难加工材料,在研磨过程中极易造成表面损伤,从而影响工件的表面精度和材料去除率。基于金属基复合材料的研磨去除机理,建立有限元仿真模型。仿真模拟超声去除SiCp/Al的表面损伤形式、颗粒的破碎方式与材料去除过程。通过超声振动辅助研磨SiCp/Al实验可知,超声振动能够有效降低表面粗糙度。材料去除率随主轴转速和研磨深度的增加而增大,随进给速度的增加而减小。

埃尔米特梁单元的分块对角与高阶质量矩阵

埃尔米特梁单元常用的集中质量矩阵,是由挠度自由度对应的一致质量矩阵元素通过行求和或节点积分构造.然而,数值结果表明该集中质量矩阵在求解包含自由端的梁振动问题时,会出现频率精度掉阶现象.本文首先从保障质量矩阵最优收敛性的数值积分精度出发,分别针对三次和五次梁单元,发展了质量矩阵的梯度增强节点积分方案.利用梯度增强节点积分方案,可以得到具有分块对角形式的单元质量矩阵,而其组装的整体质量矩阵除边界节点外仍然呈现对角形式.对于两种单元,其分块对角质量矩阵分别具有4阶最优精度和6阶次优精度.再者,将标准一致质量矩阵和具有同阶精度的梯度增强节点积分质量矩阵进行优化组合,建立了具有超收敛特性的高阶质量矩阵.最后,通过数值算例系统验证了三次和五次单元的分块对角与高阶质量矩阵的频率计算精度.

判定广义严格对角占优矩阵的新条件

仅利用矩阵自身的元素,得到了广义严格对角占优矩阵的几个新的判定条件,扩大了判别范围.给出了数值算例.