Search

Micro-Electro-Mechanical System (MEMS)-based pressure sensors operating on the principle of piezoresistivity have found profound application in various fields like automobile, aerospace, aviation, biomedical and consumer electronics. Various research studies have been conducted to optimize the design of MEMS-based pressure sensors to meet specific requirements of different fields. Modification in the structure of the piezoresistors placed on these sensors has shown great effect in this regard. However, most of these improvements have been validated through fabrication and measurement, but there has been a lack of significant studies developing analytical models to explain these improvements. This paper studies the performance of a single-turn piezoresistor design on a square silicon diaphragm. The analytical model relates the dimensions of the single-turn piezoresistor on a square diaphragm to the output voltage, and hence, sensor sensitivity is laid out. The correctness of the relation is also validated through Finite Element Analysis (FEA) performed using COMSOL Multiphysics software. Hence, an optimized single-turn design is presented which achieves a sensitivity of 203.57 mV/V/MPa over a pressure range of 0–1 MPa. These results are then compared to work from existing literature. The comparison shows an improved performance which was achieved by optimizing the design through its derived analytical model. The proposed sensor can be utilized in disposable blood pressure measurement system where high sensor sensitivity is required. [ABSTRACT FROM AUTHOR]

The paper presents the results of work leading to the construction of a spatial hybrid model based on finite element (FE) and Monte Carlo (MC) methods allowing the computer simulation of physical phenomena accompanying the steel sample testing at temperatures that are characteristic for soft-reduction process. The proposed solution includes local density variations at the level of mechanical solution (the incompressibility condition was replaced with the condition of mass conservation), and at the same time simulates the grain growth in a comprehensive resistance heating process combined with a local remelting followed by free/controlled cooling of the sample tested. Simulation of grain growth in the entire computing domain would not be possible without the support of GPU processors. There was a 59-fold increase in the computing speed on the GPU compared to single-threaded computing on the CPU. The study was complemented by examples of experimental and computer simulation results, showing the correctness of the adopted model assumptions. [ABSTRACT FROM AUTHOR]

This paper presents a quasistatic problem of an elastic body in frictional contact with a moving foundation. The model takes into account wear of the contact surface of the body caused by the friction. We recall existence and uniqueness results obtained in Sofonea et al. (2017). The main aim of this paper is to present a fully discrete scheme for numerical approximation together with an error estimation of a solution to this problem. Finally, computational simulations are performed to illustrate the mathematical model. [ABSTRACT FROM AUTHOR]

Simulating multi-body dynamics with both rigid and flexible parts and with frictional contacts is a hard problem. We solve this by expressing the couplings between the bodies as position level constraints. The implicit treatment of the constraint directions gives us improved stability over velocity based methods. Then by employing regularization of nonlinear constraints and a convex minimization formulation, we bridge constraint-based methods to traditional force-based methods. In fact, the former are just a dual variables formulation of the latter. We solve this dual problem using position based dynamics (PBD). We show how PBD is a completely valid modeling technique and we extend it with an accurate contact and Coulomb friction model. We further show for the first time how the same solver can be used to simulate both rigid and deformable solids with two way coupling. For the soft bodies we introduce a novel form of linear finite elements expressed as constraints, that is more accurate than PBD mass-spring systems. More of our results include the energy conserving Newmark integrator and the accelerated Jacobi solver suitable for parallel architectures. Note that this paper is an extended and revised version of the conference paper published in [1]. [ABSTRACT FROM AUTHOR]

Metallic double-ring metamaterial without splits in both of the rings was recently studied to show its negative refraction with magnetic resonance [Z.-G. Dong et al., Appl. Phys. Lett. 92, 064101 (2008)]. As a consequent work, systematic numerical results are presented in this paper to investigate the dependence of the magnetically resonant transmission on the structural parameters. The origin of the transmission band with magnetic resonance is shown rather explicitly with these numerical results. In addition, the optimized structure of the double-ring metamaterial is used to explore its terahertz response with all the sizes scaled down simultaneously. [ABSTRACT FROM AUTHOR]

Abstract In this work, we present and analyze a new weak Galerkin (WG) finite element method for the Brinkman equations. The finite element spaces which are made up of piecewise polynomials are easy to be constructed. Especially, the variational form considered in this work is based on two gradient operators. The stability, priori error estimates and L 2 error estimates for velocity are proved in this paper. In addition, we prove that the new method also yields globally divergence-free velocity approximations. The convergence rates are independent of the Reynolds number, thus the new WG finite element method is efficient for both the Stokes and Darcy equations dominated. Finally, numerical results illustrate the performance of the method and support the theoretical properties of the estimator. [ABSTRACT FROM AUTHOR]

Numerical models for vibro-acoustic analyses of complex mechanical systems are becoming more and more popular, in particular in the field of virtual product development. Therefore, reliable, comprehensive, and validated modeling methodologies remain crucial. However, system characteristics such as elasticities of the drive train components and nonlinear characteristics can lead to complex, and costly numerical models with a huge number of degrees of freedom. This may raise not only the need for novel and reasonable modeling strategies, but also exacerbates validation process, due to the wide scope in terms of operating conditions. In practice, structure-borne noise signals, for example, from accelerometers, are often used for the validation of mechanical systems. By choice of a sufficient number of measurement points, the interpretation becomes more complex. A lot of vibration response curves then need to be compared and interpreted over a wide operating range. In general, the interpretation focuses on deviations in quality and quantity. In this paper, to overcome these mentioned challenges, a validation methodology is proposed allowing a fast and transparent check of a number of captured signals. Therefore, it is shown how the original information can be reduced in a meaningful manner, making it possible to run a fast and accurate validation. The method is demonstrated on a real application with high mechanical complexity and it is shown that the chosen parameters are reliable. [ABSTRACT FROM AUTHOR]

Abstract Downburst produces strong winds near the ground and causes severe damage to buildings and people. The non-stationarity of downburst is characterized by moving downdraft and pulsed jet to the ground. In practice, both characteristics have often been observed simultaneously. These characteristics significantly affect the wind loads on buildings, particularly on low-rise buildings. However, the previous experimental studies of downbursts focused on each of these characteristics separately. No study has been made of their coupled effect. Therefore, we have developed a downburst simulator, which can generate a pulsed jet and a moving downdraft either separately or simultaneously. In the present paper, we have experimentally investigated the effects of non-stationary characteristics of downbursts on the wind field and wind loads on a flat-roofed low-rise building using this downburst simulator. The results may become a useful database for validating numerical simulation models of downbursts. Highlights • Effects of movement and impinging jet on the flow field of downburst is investigated. • Non-stationary flow fields of moving and pulsed jet are experimentally reproduced. • Downburst-induced pressure on a building is measured with a downburst simulator. • Downburst-induced wind forces are different from those by turbulent boundary layers. [ABSTRACT FROM AUTHOR]

Due to recent advances in genetic manipulation, transgenic mosquitoes may be a viable alternative to reduce some diseases. Feasibility conditions are obtained by simulating and analyzing mathematical models that describe the behavior of wild and transgenic populations living in the same geographic area. In this paper, we present a reaction–diffusion model in which the reaction term is a nonlinear function that describes the interaction between wild and transgenic mosquitoes, considering the zygosity, and the diffusive term that represents a nonuniform spatial spreading characterized by a random diffusion parameter. The resulting nonlinear system of partial differential equations is numerically solved using the sequential operator splitting technique, combining the finite element method and Runge-Kutta method. This scheme is numerically implemented considering uncertainty in the diffusion parameters of the model. Several scenarios simulating spatial release strategies of transgenic mosquitoes are analyzed, demonstrating an intrinsic association between the transgene frequency in the total population and the strategy adopted. [ABSTRACT FROM AUTHOR]

Numerical schemes provably preserving the positivity of density and pressure are highly desirable for ideal magnetohydrodynamics (MHD), but the rigorous positivity-preserving (PP) analysis remains challenging. The difficulties mainly arise from the intrinsic complexity of the MHD equations as well as the indeterminate relation between the PP property and the divergence-free condition on the magnetic field. This paper presents the first rigorous PP analysis of conservative schemes with the Lax--Friedrichs (LF) flux for 1D and multidimensional ideal MHD. The significant innovation is the discovery of the theoretical connection between the PP property and a discrete divergence-free (DDF) condition. This connection is established through the generalized LF splitting properties, which are alternatives to the usually expected LF splitting property that does not hold for ideal MHD. The generalized LF splitting properties involve a number of admissible states strongly coupled by the DDF condition, making their derivation very difficult. We derive these properties via a novel equivalent form of the admissible state set and an important inequality, which is skillfully constructed by technical estimates. Rigorous PP analysis is then presented for finite volume and discontinuous Galerkin schemes with the LF flux on uniform Cartesian meshes. In the 1D case, the PP property is proved for the first-order scheme with proper numerical viscosity, and also for arbitrarily high-order schemes under conditions accessible by a PP limiter. In the 2D case, we show that the DDF condition is necessary and crucial for achieving the PP property. It is observed that even slightly violating the proposed DDF condition may cause failure to preserve the positivity of pressure. We prove that the 2D LF type scheme with proper numerical viscosity preserves both the positivity and the DDF condition. Sufficient conditions are derived for 2D PP high-order schemes, and extension to 3D is discussed. Numerical examples provided in the supplementary material further confirm the theoretical findings. [ABSTRACT FROM AUTHOR]

This paper aims to gain fundamental understanding of the microscopic mechanisms that control the transition between secondary and tertiary creep around salt caverns in typical geological storage conditions. We use a self-consistent inclusion-matrix model to homogenize the viscoplastic deformation of halite polycrystals and predict the number of broken grains in a Representative Elementary Volume of salt. We use this micro-macro modeling framework to simulate creep tests under various axial stresses, which gives us the critical viscoplastic strain at which grain breakage (i.e., tertiary creep) is expected to occur. The comparison of simulation results for short-term and long-term creep indicates that the initiation of tertiary creep depends on the stress and the viscoplastic strain. We use the critical viscoplastic deformation as a yield criterion to control the transition between secondary and tertiary creep in a phenomenological viscoplastic model, which we implement into the Finite Element Method program POROFIS. We model a 850-m-deep salt cavern of irregular shape, in axis-symmetric conditions. Simulations of cavern depressurization indicate that a strain-dependent damage evolution law is more suitable than a stress-dependent damage evolution law, because it avoids high damage concentrations and allows capturing the formation of a damaged zone around the cavity. The modeling framework explained in this paper is expected to provide new insights to link grain breakage to phenomenological damage variables used in Continuum Damage Mechanics. [ABSTRACT FROM AUTHOR]

Due to the nature of salient poles in both the stator and rotor, the new class of flux-controllable stator-permanent-magnet brushless (FC-SPMBL) motors still suffers from severe torque ripples. In this paper, a new torque ripple minimization approach, namely, the harmonic current injection method, is proposed and implemented in the FC-SPMBL motor. By injecting proper harmonic current components into the fundamental current, the variation of instantaneous power is kept to minimum. Thus, the torque ripple is minimized. Both simulation and experimental results show that the proposed approach can effectively suppress the torque ripple. [ABSTRACT FROM AUTHOR]

Despite the tremendous effort of researchers and manufacturing engineers in improving the predictability of the air bending process, there is still a strong need for comprehensive and dependable prediction models. Currently, available modeling approaches all present some relevant limitations in practical applications. In this paper, we propose a new method, which represents an improvement over all existing modeling and prediction techniques. The proposed method can be used for accurate prediction of the main response variables of the air bending process: the angle a after springback and the bend deduction BD. The metamodeling method is based on the hierarchical fusion of different kinds of data: the deterministic low-fidelity response of numerical finite element method (FEM) simulations and the stochastic high fidelity response of experimental tests. The metamodel has been built over a very large database, unprecedented in the scientific literature on air bending, made of more than 500 numerical simulations and nearly 300 experimental tests. The fusion is achieved first by interpolating the FEM simulations with a kriging predictor; then, the hierarchical metamodel is built as a linear regression model of the experimental data, using the kriging predictor among the regressors. The accuracy of the method has been proved using a variant of the leave-one-out cross validation technique. The quality of the prediction yielded by the proposed method significantly over-performs the current prediction of the press brake on-line numerical control. [ABSTRACT FROM AUTHOR]

In this paper, the authors determine the shear modulus using a pure tensile test without the application of pure shear external load. In order to validate the computer-aided theoretical model, a measurement setup was designed and built by the authors. A theoretical model was realized by finite element model in order to simulate the torsion effect exerted by a pure shear mode. The authors have validated the theoretical and experimental results using a new and innovative photostress measurement. [ABSTRACT FROM AUTHOR]

The paper has considered the fractal flow in a dual media with fractal properties, where the media could be elastic, heterogeneous and visco-elastic. We argued that, the fractal flow within a geological formation with elastic property cannot be accurately described with the concept of differentiation with local operator, as this operator is unable to include into mathematical formulation the effect of elasticity. Thus to include into mathematical formula the observed facts, we have modified the model by replacing the local derivative with the non-local operator with power. A more complex problem was considered where the geological formation is considered to have visco-elastic and heterogeneity properties. We argued that, the flow within a matrix rock with these two properties cannot either be described with local derivative nor a non-local derivative with power law. In this case two non-local operators were considered, an operator with Mittag-Leffler kernel and Mittag-Leffler-Power law kernel [F. Ali et al., J. Magn. Magn. Mater. 423, 327 (2017); F. Ali et al., Eur. Phys. J. Plus 131, 310 (2016); F. Ali et al., Eur. Phys. J. Plus 131, 377 (2016); F. Ali et al., Nonlinear Sci. Lett. A 8, 101 (2017); N.A. Sheikh et al., Neural Comput. Appl. (2016) https://doi.org/10.1007/s00521-016-2815-5]. For each model, a detailed study of existence and uniqueness of the system solutions was presented using the fixed point theorem. We solved numerically each model using a more acculturate numerical scheme known as Upwind. Some numerical simulations are presented to underpin the effect of the suggested fractional differentiation. [ABSTRACT FROM AUTHOR]

This paper presents a mathematical description of the electromechanical energy converters with permanent magnets. Unlike well-known mathematical models the mathematical description is obtained by the joint solution of Maxwell and Fourier equations. This allows to determine the geometric dimensions of the active part of electromechanical energy converters taking into account thermal and electromagnetic processes on the stage of the design calculations. Finite element method is used for the final calculations of the full geometry of electromechanical energy converters. Computer simulation and experimental study of electric machines with high-coercivity permanent magnets were investigated to evaluate the effectiveness of the obtained results. The results of experimental studies have been compared with mathematical model and its high accuracy was obtained. [ABSTRACT FROM AUTHOR]

This paper presents a multipurpose numerical model that can be used to identify limiting values of the input process parameters for various laser material processing applications. The simulations layout is designed based on a central composite design (CCD) with four factors and five levels. Laser power, scanning speed, spot diameter and plate thickness are considered as input variables; while maximum temperature of the plate is considered as output variable of the model. Based on the simulation output a regression equation is developed to estimate the maximum temperature as a function of the laser process parameters (power, scanning speed, spot diameter) and the sheet thickness. The results obtained from the developed mathematical model and numerical simulation shows good agreement, which gives an assurance to use the developed model for laser processing of Ti with satisfactory precision. [ABSTRACT FROM AUTHOR]

Let M be a compact connected 3-submanifold of the 3-sphere S with one boundary component F such that there exists a collection of n pairwise disjoint connected orientable surfaces S = { S , · · ·, S } properly embedded in M, ∂ S = {∂ S , · · ·, ∂ S } is a complete curve system on F. We call S a complete surface system for M, and ∂ S a complete spanning curve system for M. In the present paper, the authors show that the equivalent classes of complete spanning curve systems for M are unique, that is, any complete spanning curve system for M is equivalent to ∂ S. As an application of the result, it is shown that the image of the natural homomorphism from the mapping class group M( M) to M( F) is a subgroup of the handlebody subgroup H . [ABSTRACT FROM AUTHOR]

Free bulging process is an experimental technique which can be used to characterize a sheet material under conditions of biaxial tension during hot forming. Analytical and semi-analytical models of this process are usually based on the hypothesis offering certain relations between the geometrical characteristics of a bulge during forming. The paper presents an original relation between a specimen thickness at the dome pole and the dome height which is used by the semi-analytical method for simulation of free bulging process. The finite-element computer simulation results are generalized to obtain this relation. The influence of the material constants on the geometrical parameters of the bulge is studied. It is shown that the sheet thickness corresponding to a specific dome height is dictated by the strain rate sensitivity index of the material. The equation describing the influence of the strain rate sensitivity index on the dome apex thickness is presented. [ABSTRACT FROM AUTHOR]

In this paper a novel modal converter based on an asymmetric dual-core photonic crystal fiber is proposed and numerically analyzed by using the full-vector finite element method. This converter allows mode conversion between the LP01 and LP11 modes, and the LP01 and LP21 modes. These modes are obtained at 1.55 μm. In addition, it was found that the operating wavelength of this device has a great dependence on the different geometric parameters of the structure, such as the diameter of the holes and the pitch as well. Finally, a compact device that can be used in the O + S + C + L + U bands with efficiency greater than 80% and a total length of 2 mm was obtained. This is a very interesting alternative to fabricate new all-fiber optic devices that could be implemented in mode division multiplexing systems. [ABSTRACT FROM AUTHOR]

Abstract: FE model updating techniques are used to update dynamic FE models of structures in the light of modal test data. Iterative methods of model updating that update a set of chosen parameters of the model, so as to reduce the difference between the natural frequencies and the mode shapes of the FE model and the corresponding quantities obtained through a modal test on the structure, are probably the most widely used methods. Once experimental modal data has been identified, a necessary prior step, before updating can be carried out, is that of establishing the correspondence between the FE model modes and the experimentally identified modes. It is however experienced that, many a times a situation is encountered where not all of the modes identified through an experiment can be correlated with certainty with those predicted by the FE model and some experimental modes may be left uncorrelated. There could be several reasons for this lack of correlation as identified in the paper. But the consequence is that such uncorrelated modes cannot be used in FE model updating using existing iterative methods based on modal data even when they form valid known pieces of information about the structure. This is a disadvantage since it reduces the quantity of experimental data available for model updating and hence makes the updating process less effective in yielding an updated model that is a closer representation of the structure. This paper identifies this as a limitation of the existing iterative methods of model updating based on modal data and puts forward a notion of FE model updating using uncorrelated modes. The paper proposes a solution to overcome this limitation in the form of a new method of FE model updating that accepts both correlated as well as uncorrelated modes for updating. This is in contrast to all the current iterative modal data based methods of model updating that are based on the assumption of availability of correlated mode pairs and hence cannot use uncorrelated mode shapes and corresponding natural frequencies in the updating process. Formulation of the proposed method is described followed by a couple of numerical examples based on a beam structure to validate the method. The robustness of the method in the presence of simulated noise is also studied. Another numerical example of a more complex F-shape structure is also presented. The method is then validated though an experimental study. The proposed method is found to successfully update an FE model yielding correct estimates of the updating parameters in the presence of uncorrelated modes. [Copyright &y& Elsevier]

A new class of history-dependent quasivariational inequalities was recently studied in [M. Sofonea and A. Matei, History-dependent quasivariational inequalities arising in contact mechanics. Eur. J. Appl. Math. 22 (2011) 471–491]. Existence, uniqueness and regularity results were proved and used in the study of several mathematical models which describe the contact between a deformable body and an obstacle. The aim of this paper is to provide numerical analysis of the quasivariational inequalities introduced in the aforementioned paper. To this end we introduce temporally semi-discrete and fully discrete schemes for the numerical approximation of the inequalities, show their unique solvability, and derive error estimates. We then apply these results to a quasistatic frictional contact problem in which the material’s behavior is modeled with a viscoelastic constitutive law, the contact is bilateral, and friction is described with a slip-rate version of Coulomb’s law. We discuss implementation of the corresponding fully-discrete scheme and present numerical simulation results on a two-dimensional example. [ABSTRACT FROM AUTHOR]

This paper deals with the mathematical modelling of confinement of paralic ecosystems. It is based on the recent paper (Frénod and Goubert in Ecol. Model. 200(1-2):139-148, ) that presents a modelling procedure in order to compute the confinement field of a lagoon. Here, we improve the existing model in order to account for tide oscillations in any kind of geometry such as non-rectangular lagoons with a non-flat bottom. The new model, that relies on PDEs rather than ODEs, is then implemented thanks to the finite element method. Numerical results confirm the feasibility of confinement studies thanks to the introduced model. [ABSTRACT FROM AUTHOR]

This paper proposes an effective method to predict the thermal conductivity of plain woven blended fabric to optimize woven fabric structure, and to evaluate thermal comfort. The unit cell model of fabric is established for numerical simulation of heat transfer through thickness. The thermal conductivity of blended yarns is calculated by a series model. The temperature and heat flux distributions are verified experimentally. [ABSTRACT FROM AUTHOR]

Purpose This paper aims to introduce a multiscale computational method for structural failure analysis with inheriting simulation of moving trans-scale boundary (MTB). This method is motivated from the error in domain bridging caused by cross-scale damage evolution, which is common in structural failure induced by damage accumulation.Design/methodology/approach Within the method, vulnerable regions with high stress level are described by continuum damage mechanics, while elastic structural theory is sufficient for the rest, dividing the structural model into two scale domains. The two domains are bridged to generate mixed dimensional finite element equation of the whole system. Inheriting simulation is developed to make the computation of MTB sustainable.Findings Numerical tests of a notched three-point bending beam and a steel frame show that this MTB method can improve efficiency and ensure accuracy while capturing the effect of material damage on deterioration of components and structure.Originality/value The proposed MTB method with inheriting simulation is an extension of multiscale simulation to structural failure analysis. Most importantly, it can deal with cross-scale damage evolution and improve computation efficiency significantly. [ABSTRACT FROM AUTHOR]

The capacity to quickly and accurately simulate extracellular stimulation of neurons is essential to the design of next-generation neural prostheses. Existing platforms for simulating neurons are largely based on finite-difference techniques; due to the complex geometries involved, the more powerful spectral or differential quadrature techniques cannot be applied directly. This paper presents a mathematical basis for the application of a spectral element method to the problem of simulating the extracellular stimulation of retinal neurons, which is readily extensible to neural fibers of any kind. The activating function formalism is extended to arbitrary neuron geometries, and a segmentation method to guarantee an appropriate choice of collocation points is presented. Differential quadrature may then be applied to efficiently solve the resulting cable equations. The capacity for this model to simulate action potentials propagating through branching structures and to predict minimum extracellular stimulation thresholds for individual neurons is demonstrated. The presented model is validated against published values for extracellular stimulation threshold and conduction velocity for realistic physiological parameter values. This model suggests that convoluted axon geometries are more readily activated by extracellular stimulation than linear axon geometries, which may have ramifications for the design of neural prostheses. [ABSTRACT FROM AUTHOR]

Model-based damage detection entails the calibration of damage-indicative parameters in a physics-based computer model of an undamaged structural system against measurements collected from its damaged counterpart. The approach relies on the premise that changes identified in the damage-indicative parameters during calibration reveal the structural damage in the system. In model-based damage detection, model calibration has traditionally been treated as a process, solely operating on the model output without incorporating available knowledge regarding the underlying mechanistic behavior of the structural system. In this paper, the authors propose a novel approach for model-based damage detection by implementing the Extended Constitutive Relation Error (ECRE), a method developed for error localization in finite element models. The ECRE method was originally conceived to identify discrepancies between experimental measurements and model predictions for a structure in a given healthy state. Implementing ECRE for damage detection leads to the evaluation of a structure in varying healthy states and determination of discrepancy between model predictions and experiments due to damage. The authors developed an ECRE-based damage detection procedure in which the model error and structural damage are identified in two distinct steps and demonstrate feasibility of the procedure in identifying the presence, location and relative severity of damage on a scaled two-story steel frame for damage scenarios of varying type and severity. [ABSTRACT FROM AUTHOR]

In this paper we extensively analyze the magnetization of superconducting bulk cylinders in parallel and transverse magnetic field. We use analytic formulations with different J c ( B ) dependencies, and compare the magnetization to 2D and 3D numerical computations for finite geometries with different aspect ratios. This subject is a basic study for understanding the hysteretic losses in a cylindrical geometry, which is the typical geometry occurring in superconducting wires and cables. The shift of the peak of the magnetization cycle is always reproduced by FEM 2D and 3D, as a consequence of considering the local field in the J c ( B ). Using different J c ( B ) dependencies and different cylinder aspect ratio, we have reproduced the behavior of the magnetization in parallel fields, checking simulation results versus analytical calculations. In this framework we also developed the analytic behavior of the magnetization in parallel fields using the generalized Kim model. Finally, for a cylinder in a perpendicular magnetic configuration, we found that 2D numerical simulations, once suitably scaled, are in agreement to 3D simulations. [ABSTRACT FROM AUTHOR]

In this work, an effective field-modulated permanent-magnet (FMPM) machine was investigated, in which the spoke-magnet outer rotor and open-slot stator were employed. The objective of this paper is to provide the mathematical modeling analysis that was performed for the purpose of control research on this type of FMPM machine. The simulation results by means of finite element analysis (FEA) are given to verify the theoretical analysis and the validity of mathematical model. A prototype machine was also fabricated for experimentation. Both the analytical model and the FEA results are validated by experimental tests on the prototype machine. [ABSTRACT FROM AUTHOR]

Shape-memory polymers (SMPs) represent a class of smart materials able to store a temporary shape and to recover the original shape upon an external stimulus, such as temperature. The present paper proposes a three-dimensional finite-strain phenomenological model for thermo-responsive SMPs, which distinguishes between two material phases presenting different properties and is based on a rule of mixtures. The proposed model is motivated by the earlier work of Reese et al. (2010) and it considers several significant material features that had not been addressed in previous phenomenological approaches. Specifically, the model reproduces both heating-stretching-cooling and cold drawing shape-fixing procedures and it takes into account the non-ideal behavior of realistic SMPs (i.e. imperfect shape-fixing and incomplete shape-recovery). Several numerical tests are reported to assess model performances, from simple uniaxial and biaxial tests to complex simulations of biomedical devices. Comparisons with experimental data taken from the literature are also provided to validate the model. The proposed improvements increase the model applicability over a wide range of polymer types and operating conditions. [ABSTRACT FROM AUTHOR]

Abstract: This paper concerns the mathematical and numerical modeling of solidification process of pure metals. Such process takes place without development of “mushy zone”, which very often occurs during solidification of alloys. The model is based on the finite element method (FEM) and takes into account the existence of a sharp solidification front which moves according to time. The paper discusses how to include the conditions of continuity on the moving interface as well as the technique of front tracking based on the level set method (LSM). The paper also presents the results of computer simulations of the solidification process of pure copper. Two in-home solvers for one- and two-dimensional problems were built. Both of them are based on FEM. The correctness of the front tracking technique was checked by comparing the results obtained from 2D simulation to 1D simulation. Furthermore, the final calculation process was done for the two-dimensional area at different temperatures on the boundaries to demonstrate the effectiveness of the solver for a complex shape of the solidification front. [Copyright &y& Elsevier]

Abstract: Stress peening forming is widely used in the aeronautics industry to induce curvatures in wing skins. Most of the investigations of stress peen forming are empirical and experimental. In this paper, a three step numerical model that can simulate this process was developed. First, an implicit Finite Element Analysis (FEA) with ANSYS where a prebending moment along the spanwise direction of the component was performed. Then, an explicit FEA with LS-DYNA simulating shot impacts on the pre-stressed component was executed in order to obtain the resulting stresses inside the component. Finally, an implicit FEA with ANSYS was performed for calculating the arc heights and the curvature radii of the component in chordwise and spanwise directions. Numerical analysis of the process shows that the prebending moments have an influence not only on the residual stress profiles but also on the curvatures of the deformed component in chordwise and spanwise directions. This model was used to establish a relationship between the prebending moment and the resulting arc heights and residual stress profiles. The numerical strategies developed in this paper supply a useful tool for studying and optimizing the stress peening process. [Copyright &y& Elsevier]

Abstract: In the comparison with the theoretical analyses of wire strands reported in the literature where obviously single-layered strands with a construction of the 1+6 wires were modelled and analysed, this paper is focused on a multi-layered strand with a construction of the 1+6+12+18 wires. The geometric parametric equations developed in the first part of this paper are implemented in CATIA V5 software code for geometric modelling of the multi-layered strand. The methodology of their implementation and the approach for the generation of the strand geometric model are demonstrated. To predict the behaviour of the multi-layered strand under tensile loads, the mathematical geometric model is further implemented in a finite element program. For this purpose ABAQUS/Explicit software is used. The derived 3D geometric models of the multi-layered strands and the results of the finite element elastic behaviour analyses of the strand under tension loads are validated through comparisons with experimental and theoretical data available. The results obtained confirm the correctness of the derived parametric equations and mathematical and physical importance of the finite element model developed. [Copyright &y& Elsevier]

Abstract: In this paper the mathematical geometric models of the single-lay wire strands and double-lay wire ropes with defined initial parameters are presented. The present geometric models fully consider the single-helix configuration of individual wires in the strand and the double-helix configuration of individual wires within the wound strands of the ropes. The mathematical representation of the single and double helixes is in form of parametric equations with variable input parameters which determine the centreline of an arbitrary circular wire of the right hand lay and left hand lay strands and ropes of the Lang lay and regular lay construction. The concrete forms of the parametric equations are derived and presented. The application of the derived geometric analytical model is illustrated by numerical examples. Techniques for the implementation of the derived mathematical models in CATIA V5 software and procedures for the generation of the rope model are briefly presented. Correctness of the derived parametric equations and a performance of the generated rope model are controlled by visualizations. The application of the derived mathematical model and the development of a finite element model for the numerical simulation of the multi-layered strand under tension tests are treated in the second part of the paper . [Copyright &y& Elsevier]

Abstract: The paper reports a preliminary study of the behavior of a high performance controlled-clearance piston gauge (CCPG) in the pressure range up to 1 GPa through finite elemental analysis (FEA). The details of the experimental characterization of this CCPG has already been published (Yadav et al., 2007 ). We have already pointed out that the use of Heydemann–Welch (HW) model for the characterization of any CCPG, has some limitation due to the fact that the linear extrapolation of the cube root of the fall rate versus jacket pressure (v 1 / 3–pj ) curve is assumed to be independent of the rheological properties of the pressure transmitting fluids. The FEA technique addresses this problem through simulation and optimization with a standard ANSYS program where the material properties of the piston and cylinder, pressure dependent density and viscosity of the pressure transmitting fluid etc. are to be used as the input parameters. Thus it provides characterization of a pressure balance in terms of effective area and distortion coefficient of the piston and cylinder. The present paper describes the results obtained on systematic studies carried out on the effect of gap profile between piston and cylinder of this controlled-clearance piston gauge, under the influence of applied pressure (p) from 100MPa to 1000MPa, on the pressure distortion coefficient (λ) of the assembly. The gap profile is also studied at different applied jacket pressure (pj ) such that pj/p varied from 0.3, 0.4 and 0.5. [ABSTRACT FROM AUTHOR]

Purpose - The purpose of this paper is to investigate the transient three-dimensional temperature distribution for a laser sintered duraform fine polyamide part by a moving Gaussian laser beam. The primary objective of the present paper is to develop computationally efficient numerical simulation technique with the commercially available finite element software domain for the accurate prediction of the temperature history and heat-affected zones of the laser sintered parts so as to finally obtain the density of the sintered sample. Design/methodology/approach - The paper proposes a mathematical model of scanning by moving laser beam and sintering sub-model. Based on the mathematical models, a simulation model was developed by using author written subroutines in ANSYS® 11.0, a general purpose finite element software. The simulation model was then run at experimental designed points using two-level factorial design of experiments (DOE) approach. The data thus generated were used to predict the equation for the density of sintered part in terms of process parameters using Design Expert software in order to analyse the designed experiments. Findings - Laser power and scan spacing were found to be significant parameters affecting the part density. Amongst the interaction terms, significant effect of laser power was found on the part density at the lower settings of the scan velocity. Temperature-time plots were generated to study the transient temperature distribution for the sintering process and with further applicability to study the thermal stresses. Research limitations/implications - The simulation model hence developed can be used for only simple part geometries and cannot be generalised for any complex geometry. Originality/value - The paper presents a simulation model which is integrated with a DOE approach so as to develop a robust as well as simple and fast approach for the optimization of quality objective. [ABSTRACT FROM AUTHOR]

This paper provides a comprehensive methodology for accurate analysis and design of miniaturized radio frequency (RF) surface acoustic wave (SAW) duplexer package. Full-wave analysis based on the three dimensional (3-D) finite element method (FEM) is applied to get the package model. The die model is obtained by combining the parasitics and acoustics models. The modeling of bonding wire is also discussed. The models of package, die, and bonding wires are assembled together to get the total response. Based on this methodology, several novel ideas are proposed to significantly improve the isolation. Simulation and measurement results are compared, and excellent agreement is found. The technique developed in this paper reduces the design cycle time greatly and can be applied to various RF SAW device packages. [ABSTRACT FROM AUTHOR]

Purpose – Numerical solution of Plateau’s problem of minimal surface using non-variational finite element method. The paper aims to discuss this issue. Design/methodology/approach – An efficient algorithm is proposed for the computation of minimal surfaces and numerical results are presented. Findings – The solutions obtained here are examined for different cases of non-linearity and are found sufficiently accurate. Originality/value – The manuscript provide the non-variational solution for Plateau’s problem. Thus it has a good value in engineering application. [ABSTRACT FROM AUTHOR]

Stiffened plates are the most common structural units used in different structural applications ranging from buildings and bridges to ships and offshore structures. The present paper aiming at investigating the influence of using Y-stiffeners profile instead of conventional stiffeners on the load carrying capacity of ship deck panels under vertical hull girder bending moment. ANSYS (Release 15.0) nonlinear finite element software was used to verify the ultimate strength results obtained from applying the International Association Classification Societies-Common Structural Rules (IACS-CSR) on both conventional and novel deck panels. The developed numerical model with conventional stiffeners profile was validated in a good agreement with a benchmark study results included in the International Ship Structures Committee (ISSC). Furthermore, both analytical and numerical simulation results were in good agreement that demonstrates the capability of using IACS-CSR for designing the new stiffeners profile in practical applications. In comparison with the conventional stiffeners profile, the proposed Y-stiffeners resulted in a considerable increase in the panel safety margin to weight (area) ratio and higher ultimate strength. Moreover, the numerical simulations provided that implementing the Y-stiffeners profile can also rise the first yielding stress point of the structural response with nearly the same strain absorption in comparison with T-stiffeners profile. [ABSTRACT FROM AUTHOR]

This paper presents an application of the finite element method and boundary element method to determine the distribution of the elongation. Computer simulations were performed using the computation of numerical algorithms according to a mathematical structure of the model and taking into account the values of all other elements of the fiber Bragg grating (FBG) sensor. Experimental studies were confirmed by elongation measurement system using one uniform FBG. [ABSTRACT FROM AUTHOR]

Finite element simulations are being more and more applied when studying the crash-worthiness of vehicles during impact. This paper deals with setting up such a simulation and discusses several ways to simplify and verify a simulated crash. For this purpose, a notch impact-testing machine will be released from a certain angle and crash into a model constructed with three different wall thicknesses. The plastic and elastic deformation is measured in the front of the model and is then used for validation of the simulation. In the end, the simulation was found to be in good agreement with the real crash data. [ABSTRACT FROM AUTHOR]

The modeling of the friction heat generation has become increasingly important in product design process including areas such as electronics, automotive, aerospace, railway (e. g. wheel and rail rolling contact, braking systems, and so on), medical industries, etc. Determination of generated friction heat in the contact of wheel and rail is important for understanding the damage mechanisms on these two bodies such as wear. This paper presents a method to determine the friction generated heat in contact of wheel and rail during normal operation using transient structural-thermal analysis in ANSYS software. [ABSTRACT FROM AUTHOR]

Purpose – The purpose of this paper is to extend the bridge scale method (BSM) developed for granular materials with only the solid phase to that taking into account the effects of wetting process in porous continuum. The granular material is modeled as partially saturated porous Cosserat continuum and discrete particle assembly in the coarse and fine scales, respectively. Design/methodology/approach – Based on the mass and momentum conservation laws for the three phases, i.e. the solid skeleton, the pore water and the pore air, the governing equations for the unsaturated porous Biot-Cosserat continuum model in the coarse scale are derived. In light of the passive air pressure assumption, a reduced finite element model for the model is proposed. According to the decoupling of the fine and coarse scale calculations in the BSM, the unsaturated porous Cosserat continuum model using the finite element method and the discrete element model using the discrete element method for granular media are combined. Findings – The numerical results for a 2D example problem of slope stability subjected to increasing rainfall along with mechanical loading demonstrate the applicability and performance of the present BSM. The microscopic mechanisms of macroscopic shear band developed in the slope are demonstrated. Research limitations/implications – Do not account for yet the effects of unsaturated pore water in the fine scale. Originality/value – The novel BSM that couples the Biot-Cosserat porous continuum modeling and the discrete particle assembly modeling in both coarse and fine scales, respectively, is proposed to provide a micro-macro discrete-continuum two-scale modeling approach for numerical simulations of the hydro-mechanical coupling problems in unsaturated granular materials. [ABSTRACT FROM AUTHOR]

This paper presents the results of calculations to determine the influence of the exciting force application point on the value of the vibration diagnostic parameter to identify the presence of a breathing crack at sub- and superharmonic resonances of the first flexural vibration mode in a cantilever beam of constant rectangular cross section using a three-dimensional finite element model. The investigation results are compared with the known data obtained by the analytical and numerical methods. [ABSTRACT FROM AUTHOR]

The goal of this paper is to prove the well-posedness for the governing equations which are used for cylindrical cloaking simulation. A new time-domain finite element scheme is developed to solve the governing equations. Numerical results demonstrating the cloaking phenomenon with the cylindrical cloak are presented. We finally extend the analysis and simulation to an elliptical cloak model. [ABSTRACT FROM AUTHOR]

In this paper the numerical simulation of aeroelastic interactions of flexibly supported two-degrees of freedom (2-DOF) airfoil in two-dimensional (2D) incompressible viscous turbulent flow subjected to a gust (sudden change of flow conditions) is considered. The flow is modeled by Reynolds averaged Navier-Stokes equations (RANS), and by k - ω turbulence model. The considered flow problem is discretized in space using the fully stabilized finite element (FE) method implemented in the developed in-house program, which allows to solve interaction problems. In order to treat the time dependent inlet boundary condition the standard stabilization procedure was modified. Further, the under relaxation procedure was introduced in order to overcome the artificial instability of the coupling algorithm. The aeroelastic response to a sudden gust is numerically analyzed with the aid of the developed FE code. [ABSTRACT FROM AUTHOR]

This paper deals with the computational modeling and numerical simulation of contact problems at finite deformations using the finite element method. Quasi-static and dynamic problems are considered and two particular frictional conditions, full stick friction and frictionless cases, are addressed. Lagrange multipliers and regularized formulations of the contact problem, such as penalty or augmented Lagrangian methods, are avoided and a new direct elimination method is proposed. Conserving algorithms are also introduced for the proposed formulation for dynamic contact problems. An assessment of the performance of the resulting formulation is shown in a number of selected benchmark tests and numerical examples, including both quasi-static and dynamic contact problems under full stick friction and frictionless contact conditions. Conservation of key discrete properties exhibited by the time stepping algorithm used for dynamic contact problems is also shown in an example. [ABSTRACT FROM AUTHOR]

SUMMARY This paper presents an embedded beam formulation for discretization independent finite element (FE) analyses of interactions between pile foundations or rock anchors and the surrounding soil in geotechnical and tunneling engineering. Piles are represented by means of finite beam elements embedded within FEs for the soil represented by 3D solid elements. The proposed formulation allows consideration of piles and pile groups with arbitrary orientation independently from the FE discretization of the surrounding soil. The interface behavior between piles and the surrounding soil is represented numerically by means of a contact formulation considering skin friction as well as pile tip resistance. The pile-soil interaction along the pile skin is considered by means of a 3D frictional point-to-point contact formulation using the integration points of the beam elements and reference points arbitrarily located within the solid elements as control points. The ability of the proposed embedded pile model to represent groups of piles objected to combined axial and shear loading and their interactions with the surrounding soil is demonstrated by selected benchmark examples. The pile model is applied to the numerical simulation of shield driven tunnel construction in the vicinity of an existing building resting upon pile foundation to demonstrate the performance of the proposed model in complex simulation environments. Copyright © 2014 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

This paper investigates the robustness of a parameter estimation procedure for nonlinear Finite Element (FE) model updating. Through this procedure, polynomial Response Surface (RS) models are constructed to approximate the response of a nonlinear FE model at every time step of the analysis. Subsequently, the optimization problem of model updating is solved iteratively in time which results in histograms of the updating parameters. With the assumption of White Gaussian measurement noise, it is shown that this parameter estimation technique has low sensitivity to the standard deviation of the measurement noise. In order to validate this, a parametric sensitivity study is performed through numerical simulations of nonlinear systems with single and multiple degrees of freedom. The results show the least sensitivity to measurement noise level, selected time window for model updating, and location of the true model parameters in RS regression domain, when vibration frequency of the system is outside the frequency bandwidth of the load. Further application of this method is also presented through a case study of a steel frame with bilinear material model under seismic loading. The results indicate the robustness of this parameter estimation technique for different cases of input excitation, measurement noise level, and true model parameters. [ABSTRACT FROM AUTHOR]

Raked piles are believed to behave better than vertical piles in a laterally flowing liquefied ground. This paper aims at numerically simulating the response of raked pile foundations in liquefying ground through nonlinear finite element analysis. For this purpose, the OpenSees computer package was used. A range of sources have been adopted in the definition of model components whose validity is assessed against case studies presented in literature. Experimental and analytical data confirmed that the backbone force density-displacement (p-y) curve simulating lateral pile response is of acceptable credibility for both vertical and raked piles. A parametric investigation on fixed-head piles subject to lateral spreading concluded that piles exhibiting positive inclination impart lower moment demands at the head while those inclined negatively perform better at liquefaction boundaries (relative to vertical piles). Further studies reveal substantial axial demand imposed upon negatively inclined members due to the transfer of gravity and ground-induced lateral forces axially down the pile. Extra care must be taken in the design of such members in soils susceptible to lateral spreading such that compressive failure (i.e. pile buckling) is avoided. [ABSTRACT FROM AUTHOR]