Optimal design of butterfly-shaped linear ultrasonic motor using finite element method and response surface methodology

A new method for optimizing a butterfly-shaped linear ultrasonic motor was proposed to maximize its mechanical output. The finite element analysis technology and response surface methodology were combined together to realize the optimal design of the butterfly-shaped linear ultrasonic motor. First, the operation principle of the motor was introduced. Second, the finite element parameterized model of the stator of the motor was built using ANSYS parametric design language and some structure parameters of the stator were selected as design variables. Third, the sample points were selected in design variable space using latin hypercube Design. Through modal analysis and harmonic response analysis of the stator based on these sample points, the target responses were obtained. These sample points and response values were combined together to build a response surface model. Finally, the simplex method was used to find the optimal solution. The experimental results showed that many aspects of the design requirements of the butterfly-shaped linear ultrasonic motor have been fulfilled. The prototype motor fabricated based on the optimal design result exhibited considerably high dynamic performance, such as no-load speed of 873 ram/s, maximal thrust of 27.5 N, maximal efficiency of 43%, and thrust-weight ratio of 45.8.

Prolonged simulation of near-free surface underwater explosion based on Eulerian finite element method

n the area of naval architecture and ocean engineering,the research about the underwater xplosion problem is of great significance.To achieve prolonged simulation of near-free surface underwater explosion,the underwater explosion transient numerical model is established in this paper based on compressible Eulerian finite element method(EFEM).Compared with Geers Hunter formula,EFEM is availably validated by simulating the free-field underwater xplosion case.Then,the bubble pulsation and flow field dynamic characteristics of the cases with different underwater explosive depth are compared in this work.Lastly,the height of the water hump and the pressure of flow flied are analyzed quantitatively through the simulation results.

Modeling of contact surface morphology and dust particles by using finite element method

The effect of dust particles on electric contacts and a hazardous size range of hard dust particles using a rigid model were discussed before. As further research, elastic-plastic model of finite element analysis was established in this work, which is closer to real condition. In this work, the behavior of large size and small size particles, and the influence of particles hardness were investigated. The calculating result of small-size particles presents a general hazardous size coefficient for different contact surface morphology; for large-size particles, it presents a hazardous size coefficient for complicated composition of the dust. And the effect of the dust shape is also discussed.

Finite element simulation for mechanical response of surface mounted solder joints under different temperature cycling

Nonlinear finite element simulation for mechanical response of surface mounted solder joint under different temperature cycling was carried out. Seven sets of parameters were used in order to evaluate the influence of temperature cycling profile parameters. The results show that temperature cycling history has significant effect on the stress response of the solder joint. Based on the concept of relative damage stress proposed by the authors, it is found that enough high temperature holding time is necessary for designing the temperature cycling profile in accelerated thermal fatigue test.

Finite Element Simulation of Processes Involving Moving Heat Sources. Application to Welding and Surface Treatment

A wide range of welding and surface treatment processes involve the use of a heat source which is moving at a constant speed over the component. The numerical simulation of such processes implies a transient analysis using a very refined mesh in order to follow properly the path of the heat source. The 3D-mesh size can be very large if one consider the welds length or the heat-treated surface size in industrial components. To reduce the computational time to acceptable values, several techniques have been investigated. The first type is to use analytical methods such as Rosenthal equations. The second type of solutions consists in performing a transient analysis using adaptive meshing. But, for a large proportion of the involved processes, practical experience demonstrates the existence of quasi steady state conditions over the major part of the heat source path. Numerical algorithms have therefore been developed to directly compute the steady temperature, metallurgical phase proportion and stress distributions. This paper gives a general overview of the different numerical methods used to simulate welding and surface treatment processes with a special emphasis on the steady state calculation. The benefits and limitations of each of them are discussed and applications are presented.

Single Track Laser Surface Hardening Model for AISI 4340 Steel Using the Finite Element Method

Laser surface hardening becomes one of the most effective techniques used to enhance wear and fatigue resistance of mechanical parts. The characteristics of the hardened surface depend on the physicochemical properties of the material as well as the heating system parameters. To adequately exploit the benefits presented by the laser heating method, it is necessary to develop a comprehensive strategy to control the process parameters in order to produce desired hardened surface attributes without being forced to use the traditional and fastidious trial and error procedures. This study presents a comprehensive approach used to build a simplified model for predicting the hardness profile. A finite element method based prediction model for AISI 4340 steel is investigated. A circular shape with a Gaussian distribution is used for modeling the laser heat source. COMSOL MULTIPHYSICS software is used to solve the heat transfer equations, estimate the temperature distribution in the part and consequently predict the hardness profile. A commercial 3 kW Nd:Yag laser system is combined to a structured experimental design and confirmed statistical analysis tools for conducting the experimental calibration and validation of the model. The results reveal that the model can effectively lead to a consistent and accurate prediction of the hardness profile characteristics under variable hardening parameters and conditions. The results show great concordance between predicted and measured values for the dimensions of hardened and melted zones.

SEMI-ELLIPTIC SURFACE CRACK IN AN ELASTIC SOLID WITH FINITE SIZE UNDER IMPACT LOADING

In this paper a semi-elliptic surface crack problem in an elastic solid of finite size under impact loading is investigated. An analysis is performed by means of fracture dynamics and the finite element method, and a three-dimensional finite element program is developed to compute the dynamic stress intensity factor. The results reveal that the effects of the solid＇s boundary surface, crack surface, material inertia and stress wave interactions play significant roles in dynamic fracture.

Dynamic finite element model updating using meta-model and genetic algorithm

Current dynamic finite element model updating methods are not efficient or restricted to the problem of local optima. To circumvent these, a novel updating method which integrates the meta-model and the genetic algorithm is proposed. Experimental design technique is used to determine the best sampling points for the estimation of polynomial coefficients given the order and the number of independent variables. Finite element analyses are performed to generate the sampling data. Regression analysis is then used to estimate the response surface model to approximate the functional relationship between response features and design parameters on the entire design space. In the fitness evaluation of the genetic algorithm, the response surface model is used to substitute the finite element model to output features with given design parameters for the computation of fitness for the individual. Finally, the global optima that corresponds to the updated design parameter is acquired after several generations of evolution. In the application example, finite element analysis and modal testing are performed on a real chassis model. The finite element model is updated using the proposed method. After updating, root-mean-square error of modal frequencies is smaller than 2%. Furthermore, prediction ability of the updated model is validated using the testing results of the modified structure. The root-mean-square error of the prediction errors is smaller than 2%.

Finite Element Analysis for Magneto-Convection Heat Transfer Performance in Vertical Wavy Surface Enclosure:Fin Size Impact

The goal of this paper is to represent a numerical study of magnetohydrodynamic mixed convection heat transfer in a lid-driven vertical wavy enclosure with a fin attached to the bottomwall.We use a finite elementmethod based on Galerkin weighted residual(GWR)techniques to set up the appropriate governing equations for the present flow model.We have conducted a parametric investigation to examine the impact of Hartmann and Richardson numbers on the flow pattern and heat transmission features inside a wavy cavity.We graphically represent the numerical results,such as isotherms,streamlines,velocity profiles,local and mean Nusselt numbers,and average surface temperature.Comparisons between the results of this work and previously published work in a literature review have been produced to examine the reliability and consistency of the data.The different sizes of the fin surface significantly impact flow creation and temperature fields.Additionally,the long fin size is necessary to enhance the heat transfer rate on the right surface at large Richardson numbers and low Hartmann numbers.Fin surfaces can significantly increase the mixing of fluid inside the enclosure,which can mean reductions in reaction times and operating costs,along with increases in heat transfer and efficiency.

Effect of Interconnect Linewidth on Evolution of Intragranular Microcracks Due to Electromigration Analyzed by Finite Element Method

The effect of interconnect linewidth on the evolution of intragranular microcracks due to surface diffusion induced by electromigration is analyzed by finite element method.The numerical results indicate that there exists critical values of the linewidth hc,the electric fieldχc and the aspect ratioβc.When h>hc,χ<χc orβ<βc,the microcrack will evolve into a stable shape as it migrates along the interconnect line.When h≤hc,χ≥χc orβ≥βc,the microcrack will split into two smaller microcracks.The critical electric field,the critical aspect ratio and the splitting time have a stronger dependence on the linewidth when h≤6.In addition,the decrease of the linewidth,the increase of the electric field or the aspect ratio is beneficial to accelerate microcrack splitting,which may delay the open failure of the interconnect line.

Investigation of FE Model Size Definition for Surface Coating Application

An efficient prediction mechanical performance of coating structures has been a constant concern since the dawn of surface engineering. However, predictive models presented by initial research are normally based on traditional solid mechanics, and thus cannot predict coating performance accurately. Also, the high computational costs that originate from the exclusive structure of surface coating systems （a big difference in the order of coating and substrate） are not well addressed by these models. To fill the needs for accurate prediction and low computational costs, a multi-axial continuum damage mechanics （CDM）-based constitutive model is introduced for the investigation of the load bearing capacity and fracture properties of coatings. Material parameters within the proposed constitutive model are determined for a typical coating （TIN） and substrate （Cu） system. An efficient numerical subroutine is developed to implement the determined constitutive model into the commercial FE solver, ABAQUS, through the user-defined subroutine, VUMAT. By changing the geometrical sizes of FE models, a series of computations are carried out to investigate （1） loading features, （2） stress distributions, and （3） failure features of the coating system. The results show that there is a critical displacement corresponding to each FE model size, and only if the applied normal loading displacement is smaller than the critical displacement, a reasonable prediction can be achieved. Finally, a 3D map of the critical displacement is generated to provide guidance for users to determine an FE model with suitable geometrical size for surface coating simulations. This paper presents an effective modelling approach for the prediction of mechanical performance of surface coatings.

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.

Propagations of Rayleigh and Love waves in ZnO films/glass substrates analyzed by three-dimensional finite element method

Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass substrates, the simulation results confirm that the Rayleigh waves along the [0001] direction and Love waves along the [1ˉ100] direction are successfully excited in the multilayered structures. Next, the crystal orientations of the ZnO films are rotated, and the influences of ZnO films with different crystal orientations on SAW characterizations, including the phase velocity, electromechanical coupling coefficient, and temperature coefficient of frequency, are investigated. The results show that at appropriate h/λ, Rayleigh wave has a maximum k^2 of 2.4% in（90°, 56.5°, 0°） ZnO film/glass substrate structure; Love wave has a maximum k^2 of 3.81% in（56°, 90°, 0°） ZnO film/glass substrate structure. Meantime, for Rayleigh wave and Love wave devices, zero temperature coefficient of frequency（TCF） can be achieved at appropriate ratio of film thickness to SAW wavelength. These results show that SAW devices with higher k^2 or lower TCF can be fabricated by flexibly selecting the crystal orientations of ZnO films on glass substrates.

Analysis and Comparison of Safety and Fatigue Comfort of Postal Bicycles Based on Finite Element Method and sEMG

The safety and the fatigue comfort were compared between a domestic and a Japanese postal bicycle. Firstly, the fatigue comfort of these two kinds of bicycles was evaluated by surface electromyographic signal (sEMG) experiment, in which human lower limb muscle groups were research objects, and the average EMG (AEMG) index and median frequency (MF) were chosen as the evaluation indexes. Secondly, the safety of these two kinds of bicycle frames was analyzed and compared by using the finite element analysis. The results show that the riding fatigue comfort of the Japanese postal bicycle is better, and the Japanese postal bicycle frame is more safe and reasonable although both the postal bicycles meet the requirement for strength. Finally, based on the above analysis, the frame structure and related parameters of the domestic postal bicycle were improved with reference to the Japanese postal bicycle and biomechanics theory.

Nonlinear finite-element-based structural system failure probability analysis methodology for gravity dams considering correlated failure modes

The structural system failure probability(SFP) is a valuable tool for evaluating the global safety level of concrete gravity dams.Traditional methods for estimating the failure probabilities are based on defined mathematical descriptions,namely,limit state functions of failure modes.Several problems are to be solved in the use of traditional methods for gravity dams.One is how to define the limit state function really reflecting the mechanical mechanism of the failure mode;another is how to understand the relationship among failure modes and enable the probability of the whole structure to be determined.Performing SFP analysis for a gravity dam system is a challenging task.This work proposes a novel nonlinear finite-element-based SFP analysis method for gravity dams.Firstly,reasonable nonlinear constitutive modes for dam concrete,concrete/rock interface and rock foundation are respectively introduced according to corresponding mechanical mechanisms.Meanwhile the response surface(RS) method is used to model limit state functions of main failure modes through the Monte Carlo(MC) simulation results of the dam-interface-foundation interaction finite element(FE) analysis.Secondly,a numerical SFP method is studied to compute the probabilities of several failure modes efficiently by simple matrix integration operations.Then,the nonlinear FE-based SFP analysis methodology for gravity dams considering correlated failure modes with the additional sensitivity analysis is proposed.Finally,a comprehensive computational platform for interfacing the proposed method with the open source FE code Code Aster is developed via a freely available MATLAB software tool(FERUM).This methodology is demonstrated by a case study of an existing gravity dam analysis,in which the dominant failure modes are identified,and the corresponding performance functions are established.Then,the dam failure probability of the structural system is obtained by the proposed method considering the correlation relationship of main failure modes on the basis of the mechanical mechanism analysis with the MC-FE simulations.

Modelling of Wideband Concentric Ring Frequency Selective Surface for 5G Devices

Frequency selective surfaces(FSSs)play an important role in wireless systems as these can be used as filters,in isolating the unwanted radiation,in microstrip patch antennas for improving the performance of these antennas and in other 5G applications.The analysis and design of the double concentric ring frequency selective surface(DCRFSS)is presented in this research.In the sub-6 GHz 5G FR1 spectrum,a computational synthesis technique for creating DCRFSS based spatial filters is proposed.The analytical tools presented in this study can be used to gain a better understanding of filtering processes and for constructing the spatial filters.Variation of the loop sizes,angles of incidence,and polarization of the concentric rings are the factors which influence the transmission coefficient as per the thorough investigation performed in this paper.A novel synthesis approach based on mathematical equations that may be used to determine the physical parameters ofDCRFSSbased spatial filters is presented.The proposed synthesis technique is validated by comparing results from high frequency structure simulator(HFSS),Ansys electronic desktop circuit editor,and an experimental setup.Furthermore,the findings acquired from a unit cell are expanded to a 2×2 array,which shows identical performance and therefore proves its stability.

Performance Comparison of Two Meta-Model for the Application to Finite Element Model Updating of Structures

To investigate the application of meta-model for finite element( FE) model updating of structures,the performance of two popular meta-model,i. e.,Kriging model and response surface model( RSM),were compared in detail. Firstly,above two kinds of meta-model were introduced briefly. Secondly,some key issues of the application of meta-model to FE model updating of structures were proposed and discussed,and then some advices were presented in order to select a reasonable meta-model for the purpose of updating the FE model of structures. Finally,the procedure of FE model updating based on meta-model was implemented by updating the FE model of a truss bridge model with the measured modal parameters. The results showed that the Kriging model was more proper for FE model updating of complex structures.

PRINCIPAL COMPONENT DECOMPOSITION BASED FINITE ELEMENT MODEL UPDATING FOR STRAIN-RATE-DEPENDENCE NONLINEAR DYNAMIC PROBLEMS

Thin wail component is utilized to absorb impact energy of a structure. However, the dynamic behavior of such thin-walled structure is highly non-linear with material, geometry and boundary non-linearity. A model updating and validation procedure is proposed to build accurate finite element model of a frame structure with a non-linear thin-walled component for dynamic analysis. Design of experiments （DOE） and principal component decomposition （PCD） approach are applied to extract dynamic feature from nonlinear impact response for correlation of impact test result and FE model of the non-linear structure. A strain-rate-dependent non-linear model updating method is then developed to build accurate FE model of the structure. Computer simulation and a real frame structure with a highly non-linear thin-walled component are employed to demonstrate the feasibility and effectiveness of the proposed approach.

FINITE ELEMENT ANALYSIS ON EVOLUTION PROCESS FOR DAMAGE MICROCRACK HEALING

Based on the thermal kinetic and mass conservation, a series of controlling equations for the finite element are derived and related programs are developed to simulate the damage microcrack healing process controlled by surface diffusion. Two kinds of typical models for microcrack splitting are proposed, i.e., the grain boundary energy existing on the crack surface and residual stresses applying on the crack surface. And the conditions of microcrack splitting in the two models are given as a function of the microcrack aspect ratio. The microcrack with traction-free surfaces will directly evolve into a spheroid.