Magnetic properties of silicon iron electrical steel are determined by using standardized measurement setups and distinct excitation parameters. Characteristic values for magnetic loss and magnetization are used to select the most appropriate material for its application. This approach is not sufficient, because of the complex material behavior inside electrical machines, which can result in possible discrepancies between estimated and actual machine behavior. The materials’ anisotropy can be one of the problems why simulation and measurement are not in good accordance.With the help of a rotational single sheet tester, the magnetic material can be tested under application relevant field distribution. Thereby, additional effects of hysteresis and anisotropy can be characterized for detailed modelling and simulation.
Underground mining extraction causes the displacement and changes of stress fields in the surrounding rock mass. The determination of the changes is extremely important when the mining activity takes place in the proximity of post-flotation tailing ponds, which may affect the stability of the tailing dams. The deterministic modeling based on principles of continuum mechanics with the use of numerical methods, e.g. finite element method (FEM) should be used in all problems of predicting rock mass displacements and changes of stress field, particularly in cases of complex geology and complex mining methods. The accuracy of FEM solutions depends mainly on the quality of geomechanical parameters of the geological strata. The parameters, e.g. young modulus of elasticity, may require verification through a comparison with measured surface deformations using geodetic methods. This paper presents application of FEM in predicting effects of underground mining on the surface displacements in the area of the KGHM safety pillar of the tailing pond of the OUOW Żelazny Most. The area has been affected by room and pillar mining with roof bending in the years 2008-2016 and will be further exposed to room-and-pillar extraction with hydraulic filling in the years 2017–2019.
The accurate prediction of iron losses has become a prominent problem in electromagnetic machine design. The basis of all iron loss models is found in the spatial field-locus of the magnetic flux density (B) and magnetic field (H). In this paper the behavior of the measured BH-field-loci is considered in FEM simulation. For this purpose, a vector hysteresis model is parameterized based on the global measurements, which then can be used to reproduce the measurement system and obtain more detailed insights on the device and its local field distribution. The IEM has designed a rotary loss tester for electrical steel, which can apply arbitrary BH-field-loci occurring during electrical machine operation. Despite its simplicity, the proposed pragmatic analytical model for vector hysteresis provides very promising results.
The paper presents a new geotechnical solution indicating a possibility of effective building structures protection. The presented solutions enable minimization of negative effects of underground mining operations. Results of numerical modelling have been presented for an example of design of preventive ditches reducing the influence of mining operations on the ground surface. To minimize the mining damage or to reduce its reach it is reasonable to look for technical solutions, which would enable effective protection of building structures. So far authors concentrated primarily on the development of building structure protection methods to minimize the damage caused by the underground mining. The application of geotechnical methods, which could protect building structures against the mining damage, was not considered so far in scientific papers. It should be noticed that relatively few publications are directly related to those issues and there are no practical examples of effective geotechnical protection. This paper presents a geotechnical solution indicating a possibility of effective protection of building structures. The presented solutions enable minimization of negative effects of underground mining operations. Results of numerical modelling have been presented for an example of design of preventive ditches reducing the influence of mining operations on the ground surface. The calculations were carried out in the Abaqus software, based on the finite element method.
The paper presents a numerical model of the novel design of the axial magnetic bearing with six cylindrical poles. The motivation behind this idea was to eliminate vibrations in rotating machinery due to the axial load. Common conception of such a bearing provides a single component of the electromagnetic force, which is not enough to reduce transverse and lateral vibrations of the armature. The proposed design allows for avoiding wobbling of the disc with the use of a few axial force components that are able to actively compensate the axial load and stabilise the disc in a balanced position. Before a real device is manufactured, a virtual prototype should be prepared. The accurate numerical model will provide essential knowledge about the performance of the axial magnetic bearing.
Experimental design and computational model for predicting debonding initiation and propagation are of interest of scientists and engineers. The design and model are expected to explain the phenomenon for a wide range of loading rates. In this work, a method to measure and quantify debonding strength at various loading rates is proposed. The method is experimentally verified using data obtained from a static test and a pulse-type dynamic test. The proposed method involves the cohesive zone model, which can uniquely be characterized with a few parameters. Since those parameters are difficult to be measured directly, indirect inference is deployed where the parameters are inferred by minimizing discrepancy of mechanical response of a numerical model and that of the experiments. The main finding suggests that the design is easy to be used for the debonding characterization and the numerical model can accurately predict the debonding for the both loading cases. The cohesive strength of the stress-wave case is significantly higher than that of the static case; meanwhile, the cohesive energy is twice larger.
The present paper is devoted to the discussion and review of the non-destructive testing methods mainly based on vibration and wave propagation. In the first part, the experimental methods of actuating and analyzing the signal (vibration) are discussed. The piezoelectric elements, fiber optic sensors and Laser Scanning Doppler Vibrometer (SLDV) method are described. Effective detecting of the flaws needs very accurate theoretical models. Thus, the numerical methods, e.g. finite element, spectral element method and numerical models of the flaws in isotropic and composite materials are presented. Moreover, the detection of the damage in structures, which are subjected to cyclic or static loads, is based on the analyzing of the change in natural frequency of the whole structure, the change of internal impedance of the material and the change in guided waves propagating through the investigated structure. All these cases are characterized in detail. At the end of this paper, several applications of the structural health monitoring systems in machine design and operation are presented.
The paper presents the dynamic behaviour of three-layer annular plates with damaged facings. The plate is composed of thin laminated, fibre-reinforced composite facings and thicker, foam core. Failure of the plate facings is modelled as fibre or matrix cracks. The plate loaded in the plane of facings with quickly increasing radially compressed forces loses its dynamic stability. Evaluation of the critical state of the plate with failures was carried out using both analytical and numerical solutions. The comparison of results between plates with material properties treated as isotropic, quasi-isotropic and composite has been conducted. Presented tables and figures create the image of dynamic responses of examined composite plates with structure failures.
A two-scale numerical homogenization approach was used for granular materials. At small-scale level, granular micro-structure was simulated using the discrete element method. At macroscopic level, the finite element method was applied. An up-scaling technique took into account a discrete model at each Gauss integration point of the FEM mesh to derive numerically an overall constitutive response of the material. In this process, a tangent operator was generated with the stress increment corresponding to the given strain increment at the Gauss point. In order to detect a loss of the solution uniqueness, a determinant of the acoustic tensor associated with the tangent operator was calculated. Some elementary geotechnical tests were numerically calculated using a combined DEM-FEM technique.
Based on the example of the pilot area in Kiev the influence of the increased static load on the superstructure of the stress-strain state of the slope was studied. The efficiency of the proposed methodology when considering the work of "home-slope-retaining structure" depending on natural and anthropogenic factors was demonstrated.
Some materials-related microstructural problems calculated using the phase-field method are presented. It is well known that the phase field method requires mesh resolution of a diffuse interface. This makes the use of mesh adaptivity essential especially for fast evolving interfaces and other transient problems. Complex problems in 3D are also computationally challenging so that parallel computations are considered necessary. In this paper, a parallel adaptive finite element scheme is proposed. The scheme keeps the level of node and edge for 2D and level of node and face for 3D instead of the complete history of refinements to facilitate derefinement. The information is local and exchange of information is minimized and also less memory is used. The parallel adaptive algorithms that run on distributed memory machines are implemented in the numerical simulation of dendritic growth and capillary-driven flows.
Numerical methods are mostly used to predict the acoustic pressure inside duct systems. In this paper, the development of a numerical method based on the convected Helmholtz equation to compute the acoustic pressure inside an axisymmetric duct is presented. A validation of the proposed method was done by a comparison with the analytical formulation for simple cases of hard wall and lined ducts. The effect of the flow on the acoustic pressure inside these ducts was then evaluated by computing this field with different Mach numbers.
The present article investigates the dynamic behavior of a fully assembled turbogenerator system influenced by misalignment. In the past, most of the researchers have neglected the foundation flexibility in the turbogenerator systems in their study, to overcome this modelling error a more realistic model of a turbogenerator system has been attempted by considering flexible shafts, flexible coupling, flexible bearings and flexible foundation. Equations of motion for fully assembled turbogenerator system including flexible foundations have been derived by using finite element method. The methodology developed based on least squares technique requires forced response information to quantify the bearing–coupling–foundation dynamic parameters of the system associated with different faults along with residual unbalances. The proposed methodology is tested for the various level of measurement noise and modelling error in the system parameters, i.e., 5% deviation in E (modulus of elasticity) and ρ (density), respectively, for robustness of the algorithm. In a practical sense, the condition analyzed in the present article relates to the identification of misalignment and other dynamic parameters viz. bearing and residual unbalance in a rotor integrated with flexible foundation.
The rigid finite element method (RFEM) has been used mainly for modelling systems with beam-like links. This paper deals with modelling of a single set of electrodes consisting of an upper beam with electrodes, which are shells with complicated shapes, and an anvil beam. Discretisation of the whole system, both the beams and the electrodes, is carried out by means of the rigid finite element method. The results of calculations concerned with free vibrations of the plates are compared with those obtained from a commercial package of the finite element method (FEM), while forced vibrations of the set of electrodes are compared with those obtained by means of the hybrid finite element method (HFEM) and experimental measurements obtained on a special test stand.
Th This paper presents a simple method of evaluating the load resistance and stiffness of corrugatedsheets locally loaded with suspended technical fixtures. As a part of this research, parametricnumerical analyses of corrugated sheets of different span, and with differently located concentratedforces, were carried out. Stress distributions in the individual folds in the elastic range and inthe ultimate limit state were identified. On their basis, equivalent concentrated load factors forthe individual folds in the sheet were determined. The load factors enable analyses of the loadresistance and stiffness of corrugated sheets loaded with concentrated forces, which can be helpfulin design practice.
A three Dimensional finite element model (FEM) incorporating the anisotropic properties and temperature profile of hot mix asphalt (HMA) pavement was developed to predict the structural responses of HMA pavement subject to heavy loads typically encountered in the field. In this study, ABAQUS was adopted to model the stress and strain relationships within the pavement structure. The results of the model were verified using data collected from the Korean Highway Corporation Test Road (KHCTR). The results demonstrated that both the base course and surface course layers follow the anisotropic behavior and the incorporation of the temperature profile throughout the pavement has a substantial effect on the pavement response predictions that impact pavement design. The results also showed that the anisotropy level of HMA and base material can be reduced to as low as 80% and 15% as a result of repeated loading, respectively.
The main idea of this work is to demonstrate an application of the generalized perturbation-based Stochastic Finite Element Method for a determination of the reliability indicators concerning elastic stability for a certain spectrum of the civil engineering structures. The reliability indicator is provided after the Eurocode according to the First Order Reliability Method, and computed using the higher order Taylor expansions with random coefficients. Computational implementation provided by the hybrid usage of the FEM system ROBOT and the computer algebra system MAPLE enables for reliability analysis of the critical forces in the most popular civil engineering structures like simple Euler beam, 2 and 3D single and multi-span steel frames, as well as polyethylene underground cylindrical shell. A contrast of the perturbation-based numerical approach with the Monte-Carlo simulation technique for the entire variability of the input random dispersion included into the Euler problem demonstrates the probabilistic efficiency of the perturbation method proposed.
The paper presents a tool for accurate evaluation of high field concentrations near singular lines, such as contours of cracks, notches and grains intersections, in 3D problems solved the BEM. Two types of boundary elements, accounting for singularities, are considered: (i) edge elements, which adjoin a singular line, and (ii) intermediate elements, which while not adjoining the line, are still under strong influence of the singularity. An efficient method to evaluate the influence coefficients and the field intensity factors is suggested for the both types of the elements. The method avoids time expensive numerical evaluation of singular and hypersingular integrals over the element surface by reduction to 1D integrals. The method being general, its details are explained by considering a representative examples for elasticity problems for a piece-wise homogeneous medium with cracks, inclusions and pores. Numerical examples for plane elements illustrate the exposition. The method can be extended for curvilinear elements.
The aim of the study is to identify the relevant aspects of numerical analysis of impact of projectiles with soft cores into a package composed of thin flexible plies located on the plastic backing. In order to illustrate the problem, normal impact of 7.62 mm TT projectile into an unclamped package comprising 36 plies of Dyneema SB71 supported on the plastic backing was selected. The problem was solved with the use of the finite element method (FEM) with the explicit integration scheme (central difference method) of motion equations in the matrix form. Based on the conducted numerical computations, it was revealed that obtaining the extreme deformations of a projectile soft core and the backing material in Lagrangian description requires employment of adaptive methods. The proposed R-adaptive method performs its role but must be used carefully due to the mass loss which may appear during calculations.
The paper presents a multi-scale mathematical model dedicated to a comprehensive simulation of resistance heating combined with the melting and controlled cooling of steel samples. Experiments in order to verify the formulated numerical model were performed using a Gleeble 3800 thermo-mechanical simulator. The model for the macro scale was based upon the solution of Fourier-Kirchhoff equation as regards predicting the distribution of temperature fields within the volume of the sample. The macro scale solution is complemented by a functional model generating voluminal heat sources, resulting from the electric current flowing through the sample. The model for the micro-scale, concerning the grain growth simulation, is based upon the probabilistic Monte Carlo algorithm, and on the minimization of the system energy. The model takes into account the forming mushy zone, where grains degrade at the melting stage – it is a unique feature of the micro-solution. The solution domains are coupled by the interpolation of node temperatures of the finite element mesh (the macro model) onto the Monte Carlo cells (micro model). The paper is complemented with examples of resistance heating results and macro- and micro-structural tests, along with test computations concerning the estimation of the range of zones with diverse dynamics of grain growth.
The microscale deformation behaviour of the Al-4.5Cu-2Mg alloy has been studied to understand the influence of various processing routes and conditions, i.e. the gravity casting with and without grain refiner, the rheocast process and the strain induced melt activation (SIMA) process. The micromechanics based simulations have been carried out on the optical microstructures of the alloy by 2D representative volume elements (RVEs) employing two different boundary conditions. Microstructural morphology, such as the grain size, the shape and the volume fraction of α-Al and binary eutectic phases have a significant effect on the stress and strain distribution and the plastic strain localization of the alloy. It is found that the stress and strain distribution became more uniform with increasing the globularity of the α-Al grain and the α-Al phase volume fraction. The simulated RVEs also reveals that the eutectic phase carries more load, but least ductility with respect to the α-Al phase. The SIMA processed alloy contains more uniform stress distribution with less stress localization which ensures better mechanical property than the gravity cast, grain refined and rheocast alloy.
The field of mechanical manufacturing is becoming more and more demanding on machining accuracy. It is essential to monitor and compensate the deformation of structural parts of a heavy-duty machine tool. The deformation of the base of a heavy-duty machine tool is an important factor that affects machining accuracy. The base is statically indeterminate and complex in load. It is difficult to reconstruct deformation by traditional methods. A reconstruction algorithm for determining bending deformation of the base of a heavy-duty machine tool using inverse Finite Element Method (iFEM) is presented. The base is equivalent to a multi-span beam which is divided into beam elements with support points as nodes. The deflection polynomial order of each element is analysed. According to the boundary conditions, the deformation compatibility conditions and the strain data measured by Fiber Bragg Grating (FBG), the deflection polynomial coefficients of a beam element are determined. Using the coordinate transformation, the deflection equation of the base is obtained. Both numerical verification and experiment were carried out. The deflection obtained by the reconstruction algorithm using iFEM and the actual deflection measured by laser displacement sensors were compared. The accuracy of the reconstruction algorithm is verified.
In order to enhance the acoustical performance of a traditional straight-path automobile muffler, a multi-chamber muffler having reverse paths is presented. Here, the muffler is composed of two internally parallel/extended tubes and one internally extended outlet. In addition, to prevent noise transmission from the muffler’s casing, the muffler’s shell is also lined with sound absorbing material. Because the geometry of an automotive muffler is complicated, using an analytic method to predict a muffler’s acoustical performance is difficult; therefore, COMSOL, a finite element analysis software, is adopted to estimate the automotive muffler’s sound transmission loss. However, optimizing the shape of a complicated muffler using an optimizer linked to the Finite Element Method (FEM) is time-consuming. Therefore, in order to facilitate the muffler’s optimization, a simplified mathematical model used as an objective function (or fitness function) during the optimization process is presented. Here, the objective function can be established by using Artificial Neural Networks (ANNs) in conjunction with the muffler’s design parameters and related TLs (simulated by FEM). With this, the muffler’s optimization can proceed by linking the objective function to an optimizer, a Genetic Algorithm (GA). Consequently, the discharged muffler which is optimally shaped will improve the automotive exhaust noise.
The paper presents an approach of numerical modelling of alloy solidification in permanent mold and transient heat transport between the casting and the mold in two-dimensional space. The gap of time-dependent width called "air gap", filled with heat conducting gaseous medium is included in the model. The coefficient of thermal conductivity of the gas filling the space between the casting and the mold is small enough to introduce significant thermal resistance into the heat transport process. The mathematical model of heat transport is based on the partial differential equation of heat conduction written independently for the solidifying region and the mold. Appropriate solidification model based on the latent heat of solidification is also included in the mathematical description. These equations are supplemented by appropriate initial and boundary conditions. The formation process of air gap depends on the thermal deformations of the mold and the casting. The numerical model is based on the finite element method (FEM) with independent spatial discretization of interacting regions. It results in multi-mesh problem because the considered regions are disconnected.