The present paper is dedicated to presentation and energy verification of the methods of stabilization the strain energy by penalty coefficients. Verification of the methods is based on the consistency and ellipticity conditions to be satisfied by the finite elements. Three methods of stabilization are discussed. The first does not satisfy the above requirements. The second is consistent but cannot eliminate parasitic energy terms. The third method, proposed by the author, is based on the decomposition of the element stiffness matrix. The method can help to eliminate locking of the finite elements. For two-noded beam element with linear shape functions and exact integration a stabilized free of locking (and elliptical) element is received (equivalent to reduced integration element). Two plate finite elements are analyzed: four-noded rectangular element and DSG triangle. A new method of stabilization with the use of four independent parameters is proposed. The finite elements with this kind of stabilization satisfy the consistency condition. In the rectangular element it was not possible to eliminate one parasitic term of energy which appears during the procedure. For DSG triangle all parasitic terms of energy are eliminated. The penalty coefficients depends on the geometry of the triangle.
A numerical analysis of the initially clamped bolt joint subject to the working pressure is presented in the paper. Special, hexahedral 21- and 28-node isoparametric finite elements have been employed to model the contact zone. In this model, one takes into account loading due to the working pressure in the gap between the gasket and the flange arising as an effect of the progressing joint opening, what has not been considered in recent papers. Nonlinear stiffness characteristics of the bolt and the flange with the gasket are developed. Working pressure corresponding to the critical bolt force resulting in the joint leakage (complete opening between the gasket and the flange) is determined. FE computational results are compared with the available experimental results. The numerical results are presented using the authors' own graphical postprocessor.
This paper aims at providing a framework for comprehensive steady-state time-domain analysis of rotating machines considering motion. The steady-state waveforms of electromagnetic and circuit quantities are computed via iterative solution of the nonlinear field-circuit-and-motion problem with constraints of time periodicity. The cases with forced speed and forced load torque are considered. A comparison of execution times with a conventional time-stepping transient model is carried out for two different machines. The numerical stability of a time-periodic model with forced speed is shown to be worse than that of traditional transient time-stepping one, although the model converges within a reasonable number of iterations. This is not the case if forced load via equation of mechanical balance is accounted for. To ensure convergence of the iterative process the physical equation of motion is replaced by the fixed-point equation. In this way the model delivers time-periodic solutions regarding not only the electromagnetic quantities but also the rotational speed.
The paper presents an analysis of the effect of shape of primary silicon crystals on the sizes of stresses and deformations in a surface layer
of A390.0 alloy by Finite Elements Method (FEM). Analysis of stereological characteristics of the studied alloy, performed based on a
quantitative metallographic analysis in combination with a statistical analysis, was used for this purpose. The presented simulation tests
showed not only the deposition depth of maximum stresses and strains, but also allowed for determining the aforementioned values
depending on the shape of the silicon crystals. The studied material is intended for pistons of internal combustion engines, therefore the
analysis of the surface layer corresponded to conditions during friction in a piston-cylinder system of an internal combustion engine having
power of up to 100 kW. The obtained results showed important differences in the values of stresses and strains up to 15% between various
shape of the silicon crystals. Crystals with sharp edges caused higher stresses and deformation locally than those with rounded shapes.
The paper presents the methodology that makes it possible to evaluate computational model and introduce current corrections to it. The methodology ensures proper interpretation of nonlinear results of numerical analyses of thin-walled structures. The suggested methodology is based on carrying out, in parallel to nonlinear numerical analysis, experimental research on some selected crucial zones of loadcarrying structures. Attention is drawn to the determinants concerning the performance of an adequate experiment. The author points out on indicating the role of model tests as a fast and economically justified research instruments practicable when designing thin-walled load-carrying structures.
The presented considerations are illustrated by an example of a structure whose geometrical complexity and ranges of deformation are characteristic for modern solutions applied in the load-carrying structures of airframes. As the representative example, one selected the area of the load-carrying structure that contains an extensive cut-out, in which the highest levels and stress gradients occur in the conditions of torsion evoking the post-buckling states within the permissible loads. The stress distributions within these ranges of deformations were used as the basis for determining the fatigue life of the structure.
This paper presents the possibility to apply numerical simulation in static analysis of reinforcedconcrete structure strengthened with carbon fibre reinforced polymer composite strips (CFRP).Reinforced concrete beams, with strengthening in form values CFRP made of carbon fibres andepoxy resin, featuring various width, as well as non-strengthened bent beams, were analysed. Thesimply supported beams arranged in a free support scheme were subjected to two concentratedforces within full range of loading (until collapse). The numerical analysis was performed throughapplication of the Finite Elements Method (FEM), and the calculation model applied took intoaccount the geometric and physical nonlinearity. The problem was solved by application of thequasi-staticstrategy method of calculations using ABAQUS software. While analysing the results,we focused on the run of changes in structure displacement and development of material damage,up to the point of destruction of the beam.
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 multi-phase permanent-magnet machines with a fractional-slot concentrated-winding (FSCW) are a suitable choice for certain purposes like aircraft, marine, and electric vehicles, because of the fault tolerance and high power density capability. The paper aims to design, optimize and prototype a five-phase fractional-slot concentrated-winding surface-mounted permanent-magnet motor. To optimize the designed multi-phase motor a multi-objective optimization technique based on the genetic algorithm method is applied. The machine design objectives are to maximize torque density of the motor and maximize efficiency then to determine the best choice of the designed machine parameters. Then, the two-dimensional Finite Element Method (2D-FEM) is employed to verify the performance of the optimized machine. Finally, the optimized machine is prototyped. The paper found that the results of the prototyped machine validate the results of theatrical analyses of the machine and accurate consideration of the parameters improved the acting of the machine.
By simulating the actual working conditions of a cable, the temperature variation rule of different measuring points under different load currents was analyzed. On this basis, a three-dimensional finite element model (FEM) was established, and the difference and influence factors between the simulation temperature and the experimental measured value were discussed, then the influence of thermal conductivity on the operating temperature of the conductor layer was studied. Finally, combined with the steady-state thermal conductivity model and the experimental measured data, the relation between thermal conductivity and load current was obtained.
This paper deals with the modelling of traction linear induction motors (LIMs) for public transportation. The magnetic end effect inherent to these motors causes an asymmetry of their phase impedances. Thus, if the LIM is supplied from the three-phase symmetrical voltage, its phase currents become asymmetric. This effect must be taken into consideration when simulating the LIMs’ performance. Otherwise, when the motor phase currents are assumed to be symmetric in the simulation, the simulation results are in error. This paper investigates the LIM performance, considering the end-effect induced asymmetry of the phase currents, and presents a comparative study of the LIM performance characteristics in both the voltage and the current mode.
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.
In this paper, finite element modelling is employed for simulating and analysing seepage and slope stability of earthfill dam via GeoStudio software. Two products are employed, which are SLOPE/W for slope stability and SEEP/W for seepage analysis. The behaviour of earthfill dam with four different types of sandy soils having different values of hydraulic conductivity (K) has been studied. Different upstream (US) slopes of 1:2, 1:2.5, 1:3 and 1:3.5 for the earthfill dam are simulated. The downstream (DS) slope is constant at 1:2. The results showed for all the four types of soils that when the US slope is increased, the amount of seepage from the dam increases and the factor of safety (F) decreases. For each US slope, when K (type of soil) increases, both seepage and F increase. Fine sand soil is associated with less seepage and less F. Sixteen equations are obtained to predict both seepage and F with respect to US slope for each type of soil and K of the soil for US slope. An experimental model for earthfill dam is constructed in the laboratory of hydraulics, Benha University to investigate the seepage of water through earthfill dams. It is concluded that seepage decreased when K decreased, and when the US slope for each type of soil decreased. The seepage increased when K increased for each US slope. Seven equations are obtained to predict seepage with respect to US slope for each type of soil, and K for each US slope.
The paper aims was assessing risks of mandible fractures consequent to impacts or sport accidents. The role of the structural stiffness of mandible, related to disocclusion state, was evaluated using the finite element method. It has been assumed, that the quasi-static stress field, due to distributed forces developed during accidents, could explain the common types of mandibular fractures. Mandibular condyles were supposed jammed in the maxillary fossae. The force of 700 N, simulating an impact on mandible, has been sequentially applied in three distinct areas: centrally, at canine zone and at the mandibular angle. Clinically most frequent fractures of mandible were recognized through the analysis of maximal principal stress/strain fields. It has been shown that mandibular fracture during accidents can be analyzed at satisfactory level using linear quasi-static models for designing protections.
This article presents a sequential model of the heating-remelting-cooling of steel samples based on the finite element method (FEM) and the smoothed particle hydrodynamics (SPH). The numerical implementation of the developed solution was completed as part of the original DEFFEM 3D package, being developed for over ten years, and is a dedicated tool to aid physical simulations performed with modern Gleeble thermo-mechanical simulators. Using the developed DEFFEM 3D software to aid physical simulations allows the number of costly tests to be minimized, and additional process information to be obtained, e.g. achieved local cooling rates at any point in the sample tested volume, or characteristics of temperature changes. The study was complemented by examples of simulation and experimental test results, indicating that the adopted model assumptions were correct. The developed solution is the basis for the development of DEFFEM 3D software aimed at developing a comprehensive numerical model allows the simulation of deformation of steel in semi solid state.