A New Classification Method of Stress Modes in Hybrid Stress Finite Element

The paper presents a new method for classifying the stress modes in hybrid stress finite element in terms of natural stress modes in finite element and the rank analysis of matrix G in forming element It reveals the relation among the different assumed stress field, and gives the general method in forming stress field Comparing with the method of eigenvalue analysis, the new method is more efficient

Modes of Occurrence and Cleaning Potential of Trace Elements in Coals from the Northern Ordos Basin and Shanxi Province, China

Based on the analyses of 43 elements in 16 samples of the raw coal and feed coal collected from the northern Ordos basin and Shanxi Province, the modes of occurrence of these elements were studied using the method of cluster analysis and factor analysis, and the cleaning potential of the hazardous elements relatively enriched in the coals was discussed by analyzing six samples of the cleaned coal from the coal-washing plants and coal cleaning simulation experiments. The results shows that the elements Br and Ba show a strong affinity to the organic matter, Cs, Cd, Pb, Zn and Hg partly to the organic matter, and the other trace elements are mainly associated with the mineral matter. Cs, Mo, P, Pb, Zn and S have positive correlations with the two principal factors, reflecting the complexity of their modes of occurrence. Some elements that were thought to show a faint relationship (Be with S and Sb with carbonates) in other rocks are found to have a strong interrelation in the coals. Clay minerals (mainly kaolinite) dominate in the coals, and Ta, Th, Ti, Sc, REE, Hf, U, Se, W, V, Nb, Mo, Al, P, Cr, Pb and Zn are distributed mostly in kaolinite, while K, Rb, Cs, and Na have much to do with illite. Conventional cleaning can reduce the concentrations of most hazardous elements in various degrees. The hazardous elements S, As, Sb, Se, Mo, Pb, Cd and Hg relatively enriched in some coals from the area studied have a relatively high potential of environmental risks. However, by physical coal cleaning processes, more than 60% of As and Hg were removed, showing a high degree of removal, more than 30% of Sb, as well as S, Pb and Cd partly associated with the inorganic matter were removed. Se and Mo showing a relatively low degree of removal could be further removed by deep crushing of the coal during physical cleaning processes, and the concentrations of S, Pb, Cd and Hg with a partial association with the organic matter could be decreased in such ways as the coal blending. Cluster analysis together with factor analysis is a rapid and effective way to deduce the mode of occurrence of an element from bulk samples, and the removability data of most hazardous elements are basically consistent with their modes of occurrence suggested, which indicates that the statistical analysis could predict the cleaning potential of hazardous elements during the physical coal cleaning.

Geochemistry and Modes of Occurrence of Hazardous Trace Elements in the No.11 Coal Seam,Antaibao Surface Mine,Shanxi Province

Concentrations of seventeen hazardous trace elements including As, Pb, Hg, Se, Cd, Cr, Co, Mo, Mn, Ni, U, V, Th, Be, Sb, Br and Zn in the No.ll coal seam, Antaibao surface mine, Shanxi Province were determined using Instrumental Neutron Activation Analysis （INAA）, Inductively Coupled Plasma Atomic Emission Spectrometry （ICP-AES）, Cold-Vapor Atomic Absorption Spectrometry （CV-AAS） and Graphite Furnace Atomic Absorption Spectrometry （GF-AAS）. Comparisons with average concentrations of trace elements in Chinese coal show that the concentrations of Hg and Cd in the No. 11 coal seam, Antaibao surface mine are much higher. They may be harmful to the environment in the process of utilization. The variations of the trace elements contents and pyritic suffur in vertical section indicated that： （a） the concentrations of As, Pb, Mn, and pyritic sulfur decrease from roof to floor; （b） the concentrations of Cr, Zn and Mo are higher in roof, floor and lower in coal seam; （c） the concentration of Br, Sb, and Hg are higher in coal seam and lower in roof and floor; （d） the concentrations of Mo, V, Th and AI vary consistently with the ash yield. Cluster analysis of trace elements, pyritic sulfur, ash yield and major elements, such as AI, Fe, P, Ca shows that： （a） pyritic sulfur, Fe, As, Mn, Ni, Be are closely associated and reflect the influence of pyrite; （b） Mo, Se, Pb, Cr, Th, Co, Ca and A! are related to clay mineral, which is the main source of ash; （c） U, Zn, V, Na, P maybe controlled by phosphate or halite; （d） Hg, Br, Sb and Cd may be mainly organic-associated elements which fall outside the three main groups. The concentration distribution characteristics of trace elements in coal seam and the cluster analysis of major and trace elements showed that the contents of trace elements in the No. 11 coal seam, Antaibao surface mine, are mainly controlled by detrital input and migration from roof and floor.

Finite element analysis of keyhole plasma arc welding based on an adaptive heat source mode

An adaptive heat source mode is proposed to account for the keyhole effect and the characteristics of volumetric distribution along the direction of the workpiece thickness. Finite element analysis of the temperature field in keyhole plasma arc welding is conducted and the weld geometry is obtained. The predicted results are in agreement with the measured ones.

Effect of occurrence mode of heavy metal elements in a low rank coal on volatility during pyrolysis

The harmful trace elements will be released during coal utilization, which can cause environment pollution and further endangering human health, especially for heavy metal elements. Compared to combustion, the release of heavy metal elements during coal pyrolysis process, as a critical initial reaction stage of combustion, has not received sufficient attention. In the present paper, a low rank coal, from Xinjiang province in China, was pyrolyzed in a fixed bed reactor from room temperature, at atmospheric pressure, with the heating rate of 10 °C/min, and the final pyrolysis temperature was from 400 to 800℃ with the interval of 100℃. The volatility of heavy metal elements (including As, Hg, Cd and Pb) during pyrolysis process was investigated. The results showed the volatility of all heavy metal elements increased obviously with increasing temperature, and followed the sequence as Hg > Cd > As > Pb, which was mainly caused by their thermodynamic property and occurrence modes in coal. The occurrence modes of heavy metals were studied by sink-andfloat test and sequential chemical extraction procedure, and it can be found that the heavy metal elements were mainly in the organic and residual states (clay minerals) in the raw coal. And most of the organic heavy metals escaped during the pyrolysis process, the remaining elements were mainly in the residual state, and the elements in Fe-Mn state also tended to remain in the char.

Boundary Element Analysis forModeⅢCrack Problems of Thin-Walled Structures from Micro-to Nano-Scales

This paper develops a new numerical framework for modeⅢcrack problems of thin-walled structures by integrating multiple advanced techniques in the boundary element literature.The details of special crack-tip elements for displacement and stress are derived.An exponential transformation technique is introduced to accurately calculate the nearly singular integral,which is the key task of the boundary element simulation of thin-walled structures.Three numerical experiments with different types of cracks are provided to verify the performance of the present numerical framework.Numerical results demonstrate that the present scheme is valid for modeⅢcrack problems of thin-walled structures with the thickness-to-length ratio in the microscale,even nanoscale,regime.

A B-spline active contour model based on finite element method

A B-spline active contour model based on finite element method is presented, into which the advantages of a B-spline active contour attributing to its fewer parameters and its smoothness is built accompanied with reduced computational complexity and better numerical stability resulted from the finite element method. In this model, a cubic B-spline segment is taken as an element, and the finite element method is adopted to solve the energy minimization problem of the B-spline active contour, thus to implement image segmentation. Experiment results verify that this method is efficient for B-spline active contour, which attains stable, accurate and faster convergence.

DEFORMATION RIGIDITY OF ASSUMED STRESS MODES IN HYBRID ELEMENTS

The new methods to determine the zero-energy deformation modes in the hybrid elements and the zero-energy stress modes in their assumed stress fields are presented by the natural deformation modes of the elements. And the formula of the additional element deformation rigidity due to additional mode into the assumed stress field is derived. Based on, it is concluded in theory that the zero-energy stress mode cannot suppress the zero-energy deformation modes but increase the extra rigidity to the nonzero-energy deformation modes of the element instead. So they should not be employed to assume the stress field. In addition, the parasitic stress modes will produce the spurious parasitic energy and result the element behaving over rigidity. Thus, they should not be used into the assumed stress field even though they can suppress the zero-energy deformation modes of the element. The numerical examples show the performance of the elements including the zero-energy stress modes or the parasitic stress modes.

Orthogonal basic deformation mode method for zero-energy mode suppression of hybrid stress elements

A set of basic deformation modes for hybrid stress finite elements are directly derived from the element displacement field. Subsequently, by employing the so-called united orthogonal conditions, a new orthogonalization method is proposed. The result- ing orthogonal basic deformation modes exhibit simple and clear physical meanings. In addition, they do not involve any material parameters, and thus can be efficiently used to examine the element performance and serve as a unified tool to assess different hybrid elements. Thereafter, a convenient approach for the identification of spurious zero-energy modes is presented using the positive definiteness property of a flexibility matrix. More- over, based on the orthogonality relationship between the given initial stress modes and the orthogonal basic deformation modes, an alternative method of assumed stress modes to formulate a hybrid element free of spurious modes is discussed. It is found that the orthogonality of the basic deformation modes is the sufficient and necessary condition for the suppression of spurious zero-energy modes. Numerical examples of 2D 4-node quadrilateral elements and 3D 8-node hexahedral elements are illustrated in detail to demonstrate the efficiency of the proposed orthogonal basic deformation mode method.

Calculation of the Coupling Coefficient of Twin-Core Fiber Based on the Supermode Theory with Finite Element Method

Currently, coupled mode theory (CMT) is widely used for calculating the coupling coefficient of twin-core fibers (TCFs) that are used in a broad range of important applications. This approach is highly accurate for scenarios with weak coupling between the cores but shows significant errors in the strong coupling scenarios, necessitating the use of a more accurate method for coupling coefficient calculations. Therefore, in this work, we calculate the coupling coefficients of TCFs using the supermode theory with finite element method (FEM) that has higher accuracy than CMT, particularly for the strong coupling TCF. To investigate the origin of the differences between the results obtained by these two methods, the modal field distributions of the supermodes of TCF are simulated and analyzed in detail.

A Comparison of Active and Passive Metamaterials from Equivalent Lumped Elements Modes

With ever-increasing operating frequencies and complicated artificial structures, loss effects become more and more important in applications of metamaterials. Based on circuit theory and transmission line principle, the design equations for effective electromagnetic (EM) parameters (attenuation constant α, phase constant β, characteristic impedance Z0) of general active and passive metamaterial are compared and derived from the equivalent lumped circuit parameters (R, G, LL, CL, LR, CR). To verify the design equations, theα, βand Z0 indifferent cases, including balanced, unbalanced, lossless, passive and active, are shown by numerical simulations. The results show that using the active method can diminish the loss effects. Meantime, it also has influence on phase constant and real part of characteristic impedance.

REFINED HYBRID MINDLIN PLATE ELEMENT FOR GEOMETRICALLY NON LINEAR ANALYSIS

Based upon a generalized variational principle, which relaxed the inter element continuity requirements, a novel refined hybrid Mindlin plate element is developed, its non linear element stiffness matrices are decomposed into a series of matrices with respect to the assumed strain modes. The formulation presented in this paper is different from any other non linear mixed/hybrid element formulation all successful experience of linear hybrid formulation is absorbed into the formulation(adding non conforming modes and realizing orthogonalization) Numerical results show that the present approach is more effective than any other non linear hybrid element formulation over the accuracy and computational efficiency. In addition, non conforming modes can also overcome the shear locking effect.

Investigation of Microstructure, Natural Frequencies and Vibration Modes of Dragonfly Wing

In the present work, a thorough investigation on the microstructural and morphological aspects of dragonfly wings was carried out using scanning electron microscope. Then, based on this study and the previous reports, a precise three-dimensional numerical model was developed and natural frequencies and vibration modes of dragonfly forewing were determined by finite element method. The results shown that dragonfly wings are made of a series of adaptive materials, which form a very complex composite structure. This bio-composite fabrication has some unique features and potential benefits. Furthermore, the numerical results show that the first natural frequency of dragonfly wings is about 168 Hz and bending is the predominant deformation mode in this stage. The accuracy of the present analysis is verified by comparison of calculated results with experimental data. This paper may be helpful for micro aerial vehicle design concerning dynamic response.

Assessment of structural stability in Bohai Sea area based on AHP-GDM model

The AHP-GDM model is used for the assessment of structural stability, with the Bohai Sea area as an exam- ple. In this model, the credit degree of each expert is calculated through the assessment matrix based on the similarity and diversity of vector. The comprehensive opinions of expert panel are quantitatively obtained by considering the effect ofcredit degree. According to the geological structural setting, the Bohai Sea is di- vided into twelve assessment zones of structural stability by non-uniform element method. The structural stability grade of each zone is obtained on the basis of the latest geophysical data, earthquake statistical data, and the information of fault activities, current stress field and crustal deformation. The results show that there are one relatively stable area, three relatively sub-stable areas, six relatively sub-unstable areas and two relatively unstable areas. The assessment results of non-uniform element method are very close with those of uniform grid method with size of 0.25 in longitude direction and 0.14 in latitude direction. However the workload of non-uniform element method is only 1 / 16 of the latter. Compared with traditional assessment methods of structural stability, a more objective and reliable assessment result can be obtained by combining non-uniform element method and AHP-GDM model.

Generalized mixed finite element method for 3D elasticity problems

Without applying any stable element techniques in the mixed methods, two simple generalized mixed element(GME) formulations were derived by combining the minimum potential energy principle and Hellinger–Reissner(H–R) variational principle. The main features of the GME formulations are that the common C0-continuous polynomial shape functions for displacement methods are used to express both displacement and stress variables, and the coefficient matrix of these formulations is not only automatically symmetric but also invertible. Hence, the numerical results of the generalized mixed methods based on the GME formulations are stable. Displacement as well as stress results can be obtained directly from the algebraic system for finite element analysis after introducing stress and displacement boundary conditions simultaneously. Numerical examples show that displacement and stress results retain the same accuracy. The results of the noncompatible generalized mixed method proposed herein are more accurate than those of the standard noncompatible displacement method. The noncompatible generalized mixed element is less sensitive to element geometric distortions.

Elemental characteristics and paleoenvironment reconstruction： a case study of the Triassic lacustrine Zhangjiatan oil shale, southern Ordos Basin, China

Using trace elements to reconstruct paleoenvironment is a current hot topic in geochemistry. Through analytical tests of oil yield, ash yield, calorific value, total sulfur, major elements, trace elements, and X-ray diffraction, the quality, mineral content, occurrence mode of elements, and paleoenvironment of the Zhangjiatan oil shale of the Triassic Yanchang Formation in the southern Ordos Basin were studied. The analyses revealed relatively high oil yield(average 6.63%) and medium quality. The mineral content in the oil shale was mainly clay minerals,quartz, feldspar, and pyrite; an illite–smectite mixed layer comprised the major proportion of clay minerals. Compared with marine oil shale in China, the Zhangjiatan oil shale had higher contents of quartz, feldspar, and clay minerals, and lower calcite content. Silica was mainly in quartz and Fe was associated with organic matter, which is different from marine oil shale. The form of calcium varied. Cluster analyses indicated that Fe, Cu, U, V, Zn, As,Cs, Cd, Mo, Ga, Pb, Co, Ni, Cr, Sc, P, and Mn are associated with organic matter while Ca, Na, Sr, Ba, Si, Zr, K,Al, B, Mg, and Ti are mostly terrigenous. Sr/Cu, Ba/Al, V/(V+ Ni), U/Th, AU, and δU of oil shale samples suggest the paleoclimate was warm and humid, paleoproductivity of the lake was relatively high during deposition of the shale—which mainly occurred in fresh water—and the paleo-redox condition was dominated by reducing conditions. Fe/Ti ratios of the oil shale samples suggest clear hydrothermal influence in the eastern portion of the study area and less conspicuous hydrothermal influence in the western portion.

Identification of the anomaly component using BEMD combined with PCA from element concentrations in the Tengchong tin belt, SW China

Concentration of elements or element groups in a geological body is the result of multiple stages of rockforming and ore-forming geological processes.An ore-forming element group can be identified by PCA(principal component analysis)and be separated into two components using BEMD(bi-dimensional empirical mode decomposition):(1)a high background component which represents the ore-forming background developed in rocks through various geological processes favorable for mineralization(i.e.magmatism,sedimentation and/or metamorphism);(2)the anomaly component which reflects the oreforming anomaly that is overprinted on the high background component developed during mineralization.Anomaly components are used to identify ore-finding targets more effectively than ore-forming element groups.Three steps of data analytical procedures are described in this paper;firstly,the application of PCA to establish the ore-forming element group;secondly,using BEMD on the o re-forming element group to identify the anomaly components created by different types of mineralization processes;and finally,identifying ore-finding targets based on the anomaly components.This method is applied to the Tengchong tin-polymetallic belt to delineate ore-finding targets,where four targets for Sn(W)and three targets for Pb-Zn-Ag-Fe polymetallic mineralization are identified and defined as new areas for further prospecting.It is shown that BEMD combined with PCA can be applied not only in extracting the anomaly component for delineating the ore-finding target,but also in extracting the residual component for identifying its high background zone favorable for mineralization from its oreforming element group.

Optimization of Bearing Locations for Maximizing First Mode Natural Frequency of Motorized Spindle-Bearing Systems Using a Genetic Algorithm

This paper has developed a genetic algorithm (GA) optimization approach to search for the optimal locations to install bearings on the motorized spindle shaft to maximize its first-mode natural frequency (FMNF). First, a finite element method (FEM) dynamic model of the spindle-bearing system is formulated, and by solving the eigenvalue problem derived from the equations of motion, the natural frequencies of the spindle system can be acquired. Next, the mathematical model is built, which includes the objective function to maximize FMNF and the constraints to limit the locations of the bearings with respect to the geometrical boundaries of the segments they located and the spacings between adjacent bearings. Then, the Sequential Decoding Process (SDP) GA is designed to accommodate the dependent characteristics of the constraints in the mathematical model. To verify the proposed SDP-GA optimization approach, a four-bearing installation optimazation problem of an illustrative spindle system is investigated. The results show that the SDP-GA provides well convergence for the optimization searching process. By applying design of experiments and analysis of variance, the optimal values of GA parameters are determined under a certain number restriction in executing the eigenvalue calculation subroutine. A linear regression equation is derived also to estimate necessary calculation efforts with respect to the specific quality of the optimization solution. From the results of this illustrative example, we can conclude that the proposed SDP-GA optimization approach is effective and efficient.

3D modeling and deformation analysis for electromagnetic sheet forming process

Electromagnetic forming （EMF） is a high-speed forming method which can be quite effective in increasing the forming limits of metal sheet. However, the EMF process is complicated due to magnetic-structure coupling analysis. Numerical simulation offers an opportunity to overcome the problem. Nevertheless, most present models for EMF process are limited to 2D axisymmetric model. So, a three-dimensional （3D） finite element model was established to analyze the electromagnetic sheet bulging. The contact between the sheet and the die and the effect of sheet deformation on the magnetic field analysis were both taken into consideration during the forming process. The simulation results of deflection at the sheet center and 20 mm away from the center were in agreement with the experimental ones. The plastic strain energy and plastic strain were analyzed.

Dynamic analysis of beam-cable coupled systems using Chebyshev spectral element method

The dynamic characteristics of a beam-cable coupled system are investigated using an improved Chebyshev spectral element method in order to observe the effects of adding cables on the beam. The system is modeled as a double Timoshenko beam system interconnected by discrete springs. Utilizing Chebyshev series expansion and meshing the system according to the locations of its connections, numerical results of the natural frequencies and mode shapes are obtained using only a few elements, and the results are validated by comparing them with the results of a finite-element method. Then the effects of the cable parameters and layout of connections on the natural frequencies and mode shapes of a fixed-pinned beam are studied. The results show that the modes of a beam-cable coupled system can be classified into two types, beam mode and cable mode, according to the dominant deformation. To avoid undesirable vibrations of the cable, its parameters should be controlled in a reasonable range, or the layout of the connections should be optimized.