Reduced-order finite element method based on POD for fractional Tricomi-type equation

The reduced-order finite element method （FEM） based on a proper orthogo- nal decomposition （POD） theory is applied to the time fractional Tricomi-type equation. The present method is an improvement on the general FEM. It can significantly save mem- ory space and effectively relieve the computing load due to its reconstruction of POD basis functions. Furthermore, the reduced-order finite element （FE） scheme is shown to be un- conditionally stable, and error estimation is derived in detail. Two numerical examples are presented to show the feasibility and effectiveness of the method for time fractional differential equations （FDEs）.

Reduced-order proper orthogonal decomposition extrapolating finite volume element format for two-dimensional hyperbolic equations

This paper is concerned with establishing a reduced-order extrapolating fi- nite volume element （FVE） format based on proper orthogonal decomposition （POD） for two-dimensional （2D） hyperbolic equations. For this purpose, a semi discrete variational format relative time and a fully discrete FVE format for the 2D hyperbolic equations are built, and a set of snapshots from the very few FVE solutions are extracted on the first very short time interval. Then, the POD basis from the snapshots is formulated, and the reduced-order POD extrapolating FVE format containing very few degrees of freedom but holding sufficiently high accuracy is built. Next, the error estimates of the reduced-order solutions and the algorithm procedure for solving the reduced-order for- mat are furnished. Finally, a numerical example is shown to confirm the correctness of theoretical conclusions. This means that the format is efficient and feasible to solve the 2D hyperbolic equations.

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

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.

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.

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.

Finite element analysis of stress and strain distributions in mortise and loose tenon furniture joints

We studied the effect of loose tenon dimensions on stress and strain distributions in T-shaped mortise and loose tenon （M＆amp;LT） furni-ture joints under uniaxial bending loads, and determined the effects of loose tenon length （30, 45, 60, and 90 mm） and loose tenon thickness （6 and 8 mm） on bending moment capacity of M＆amp;LT joints constructed with polyvinyl acetate （PVAc） adhesive. Stress and strain distributions in joint elements were then estimated for each joint using ANSYS finite element （FE） software. The bending moment capacity of joints increased significantly with thickness and length of the tenon. Based on the FE analysis results, under uniaxial bending, the highest shear stress values were obtained in the middle parts of the tenon, while the highest shear elastic strain values were estimated in glue lines between the tenon sur-faces and walls of the mortise. Shear stress and shear elastic strain values in joint elements generally increased with tenon dimensions and corre-sponding bending moment capacities. There was consistency between predicted maximum shear stress values and failure modes of the joints.

Modal analysis of prestressed draft pad of freight wagons using finite element method

The present work intends to investigate dynamic behaviour of draft gear using finite element method. The longitudinal force that the draft gear absorbs usually leads to the failure of its components, especially, the load bearing draft pads. Dynamic behaviour of an individual draft pad and a draft gear is determined and characterized with exciting frequen- cies and corresponding mode shapes. The effect of compressive prestress load on the dynamic behaviour of an individual draft pad is also determined as the draft pads in assembled state are under constant axial compressive force in the draft gear. The vibration characteristics of individual draft pad are compared with draft pads that are part of draft gear. The modal analysis gives us a basis for subjecting a draft pad to higher frequency loading for determining its fatigue behaviour.

Study of Axisymmetric Infinite Guide Lined with Locally Reacting Material without Flow Using DtN Operators

The present work proposed a new method for the modeling by the finite element method of the acoustic propagation problems in infinite axisymmetric cylindrical guides lined with locally reacting absorbent materials without flow. The method deals with the development of an efficient transparent boundary condition based on DtN operators. The method developed in this study is successfully applied to a straight axisymmetric lined guide by imposing a mode on one of the artificial boundaries of the truncated guide. The results are in good agreement with analytical solutions. Applying the method for a non-uniform axisymmetric lined guide which is a complex case, proved its effectiveness and the results compared to those of PML layers are in very good agreement.

A NEW HYBRID QUADRILATERAL FINITE ELEMENT FOR MINDLIN PLATE

In this paper a new quadrilateral plate element concerning the effect of transverse shear strain has been presented. It is derived from the hybrid finite element model based on the principles of virtual work. The outstanding advantage of this element is to use more rational trial functions of the displacements. For this reason, every variety of plate deformation can be simulated really whilc the least degrees of freedom is employed.A wide range of numerical tests was conducted and the results illustrate that this element has a very wide application scope to the thickness of plates and satisfactory accuracy can be obtained by coarse mesh for all kinds of examples.

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors （SIFs） directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.

A reduced-order extrapolation algorithm based on CNLSMFE formulation and POD technique for two-dimensional Sobolev equations

A reduced-order extrapolation algorithm based on Crank-Nicolson least-squares mixed finite element （CNLSMFE） formulation and proper orthogonal decomposition （POD） technique for two-dimensional （2D） Sobolev equations is established. The error estimates of the reduced-order CNLSMFE solutions and the implementation for the reduced-order extrapolation algorithm are provided. A numerical example is used to show that the results of numerical computations are consistent with theoretical conclusions. Moreover, it is shown that the reduced-order extrapolation algorithm is feasible and efficient for seeking numerical solutions to 2D Sobolev equations.

Reduced-order modeling and vibration transfer analysis of a fluid-delivering branch pipeline that consider fluid–solid interactions

The efficient dynamic modeling and vibration transfer analysis of a fluid-delivering branch pipeline(FDBP)are essential for analyzing vibration coupling effects and implementing vibration reduction optimization.Therefore,this study proposes a reduced-order dynamic modeling method suitable for FDBPs and then analyzes the vibration transfer characteristics.For the modeling method,the finite element method and absorbing transfer matrix method(ATMM)are integrated,considering the fluid–structure coupling effect and fluid disturbances.The dual-domain dynamic substructure method is developed to perform the reduced-order modeling of FDBP,and ATMM is adopted to reduce the matrix order when solving fluid disturbances.Furthermore,the modeling method is validated by experiments on an H-shaped branch pipeline.Finally,transient and steady-state vibration transfer analyses of FDBP are performed,and the effects of branch locations on natural characteristics and vibration transfer behavior are analyzed.Results show that transient vibration transfer represents the transfer and conversion of the kinematic,strain,and damping energies,while steady-state vibration transfer characteristics are related to the vibration mode.In addition,multiple-order mode exchanges are triggered when branch locations vary in frequency-shift regions,and the mode-exchange regions are also the transformation ones for vibration transfer patterns.

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.

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 ModeⅡfracture behavior of rocks under hydrostatic pressure using the short core in compression(SCC)method

The shear failure of rocks under both a static triaxial stress and a dynamic disturbance is common in deep underground engineering and it is therefore essential for the design of underground engineering to quantitively estimate the dynamic ModeⅡfracture toughness KⅡCof rocks under a triaxial stress state.However,the method for determining the dynamic KⅡCof rocks under a triaxial stress has not been developed yet.With an optimal sample preparation,the short core in compression(SCC)method was designed and verified in this study to measure the dynamic KⅡCof Fangshan marble(FM)subjected to different hydrostatic pressures through a triaxial dynamic testing system.The formula for calculating the dynamic KⅡCof the rock SCC specimen under hydrostatic pressures was obtained by using the finite element method in combination with secondary cracks.The experimental results indicate that the failure mode of the rock SCC specimen under a hydrostatic pressure is the shear fracture and the KⅡCof FM increases as the loading rate.In addition,at a given loading rate the dynamic rock KⅡCis barely affected by hydrostatic pressures.Another important observation is that the dynamic fracture energy of FM enhances with loading rates and hydrostatic pressures.

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.

Texture evolution and slip mode of a Ti-5.5Mo-7.2Al-4.5Zr-2.6Sn-2.1Cr dual-phase alloy during cold rolling based on multiscale crystal plasticity finite element model

The complex micromechanical response among grains remains a persistent challenge to understand the deformation mechanism of titanium alloys during cold rolling.Therefore,in this work,a multiscale crystal plasticity finite element method of dual-phase alloy was proposed and secondarily developed based on LS-DYNA software.Afterward,the texture evolution and slip mode of a Ti-5.5Mo-7.2Al-4.5Zr-2.6Sn-2.1Cr alloy,based on the realistic 3D microstructure,during cold rolling(20%thickness reduction)were systematically investigated.The relative activity of the■slip system in theαphase gradually increased,and then served as the main slip mode at lower Schmid factor(<0.2).In contrast,the contribution of the■slip system to the overall plastic deformation was relatively limited.For theβphase,the relative activity of the<111>{110}slip system showed an upward tendency,indicating the important role of the critical resolved shear stress relationship in the relative activity evolutions.Furthermore,the abnormally high strain of very fewβgrains was found,which was attributed to their severe rotations compelled by the neighboring pre-deformedαgrains.The calculated pole figures,rotation axes,and compelled rotation behavior exhibited good agreement to the experimental results.

MODES COMPUTATION AND ANALYSIS FOR LONG MULTI-SPAN BUNDLE CONDUCTORS OF POWER TRANSMISSION LINES WITH GALLOPING CONTROL DEVICES

A new type of element which is suitable for solving the modes of thegalloping long multi-span bundle conductor structures is presented. The element is composed of allsub-conductor segments between two spacers. Based on the linearized governing differential equationsof the conductors, the mass matrix and stiffness matrix of the element in consideration of theconstrained relations imposed on the conductors by spacers are derived. The dynamic characteristicsof the galloping control devices can be directly added to the element. The modes for an actual powerline structure are computed by using the element formula and FEM procedures, where seven cases ofdifferent galloping control device allocations are considered. Compared with the measured data, themethod is shown to be reliable and effective. Analysis and discussions of the computational resultsare given. Some hints that are helpful to further investigation of galloping are also obtained .

Rock fracture under anti-plane shear (Mode Ⅲ) loading

Anti-plane punch-through shear test and anti-planefour-point bending test are used to study the crack initiation and propagation under anti-plane shear (Mode Ⅲ) loading. The tensile and shear stresses at the crack tip are calculated by finite element method. The results show that under Mode Ⅲ loading the maximum principal stress σ1 at crack tip is smaller or a little larger than the maximum shear stress τmax. Since the tensile strength of brittle rock is much lower than its shear strength, σ1 is easy to reach its critical value before τmax reaches its critical value and thus results in Mode I fracture. The fracture trajectory is helicoid and the normal direction of tangential plane with the fractured helicoid is along the predicted direction of the maximum principal stress at the notch tip. It is further proved that Mode Ⅰ instead of Mode Ⅲ fracture occurs in brittle rock under Mode Ⅲ loading. The fracture mode depending on the fracture mechanism must be distinguished from the loading form.