This study was carried out on the background of Sutong Bridge project based on fracture mechanics, aiming at analyzing the growth mechanism of fatigue cracks of a bridge under the load of vehicles. Stress intensity factor (SIF) can be calculated by various methods. Three steel plates with different kinds of cracks were taken as the samples in this study. With the combination of finite element analysis software ABAQUS and the J integral method, SIF values of the samples were calculated. After that, the extended finite element method in the simulation of fatigue crack growth was introduced, and the simulation of crack growth paths under different external loads was analyzed. At last, we took a partial model from the Sutong Bridge and supposed its two dangerous parts already had fine cracks; then simulative vehicle load was added onto the U-rib to predict crack growth paths using the extended finite element method.
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 present paper reports the results of theoretical and experimental studies of the process of die forging a bimetallic door handle intended for the production of a helicopter. The aim of the studies was to develop and implement a technology for die forging of a product with a specific mass similar to that of magnesium alloys which will have, however higher corrosion resistance. Numerical modelling and industrial tests were carried out based on the previously forging processes for an AZ31 alloy door handle. The material for the tests was a bimetallic bar produced by the explosive welding method, in which the core was of alloy AZ31, and the cladding layer was made of 1050A grade aluminium. The studies were conducted for two variants: Variant I – the forging process was mapped by numerical modelling and industrial tests for the die shape and parameters used in the forging of the AZ31 alloy door handle, Variant II – the tool shape was optimized and process parameters were selected so as to obtain a finished product characterized by a continuous Al layer. From the theoretical studies and experimental tests carried out it has been found that the application of the Variant I does not assure that a finished door handle characterized by a continuous cladding layer will be produced. Within this study, a novel method of bimetallic door handle die forging (Variant II) has been developed, which limits the amount of the flash formed and assures the integrity of the cladding layer.
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.
The paper presents the results of experimental-numerical tests of firing at aluminum composite materials. The test materials were manufactured by pressure infiltration of porous ceramic preforms made of -Al2O3 particles in the amount of 30% and 40% by volume. The EN AW-7075 alloy was chosen as the material matrix, and the steel 7.62×39 mm (M 43) FMJ (Full Metal Jacket) intermediate ammunition was selected for firing. In the result of the experiment, the samples were perforated with a clear difference in the muzzle diameter. The projectile with fragments caused damage to up to three reference plates placed behind the samples (witness plates) in composites with 40% of particles by volume. The mechanics of crack propagation during ballistic impacts of the projectile was characterized based on microstructure studies. Then, using numerical analysis of impact load, the examination of composite materials puncture in the ABAQUS environment was carried out. The Finite Element Method (FEM) was employed for the discretization of geometric models using Hex elements. The Johnson-Cook constitutive model describing the relationship between stress and strain in metal-ceramic composites was applied for the analyses. Numerical models were then subjected to numerical verification using smoothed particle hydrodynamics (SPH). Based on the obtained results, it was found that the hybrid FEM/SPH method correlates significantly with the experimental results.
Percutaneous RF ablation is one of alternative treatment for non-surgical liver tumors. Ablative changes in hepatic tissue can be successfully estimated using the finite element method. The authors created a 3D model of a multi-tine applicator immersed in liver tissue, and then determined the optimal values of voltage applied to such an RF electrode, which do not exceed the therapeutic temperature range valid during thermal ablation procedure. Importantly, the simulations were carried out for the RF electric probes with 2 to 5 evenly spaced arms. Additionally, the thermal damage of hepatic tissue for multi-armed applicators working at pre-defined limit values of voltages was established based on the Arrhenius model.
In this paper a scaling approach for the solution of 2D FE models of electric machines is proposed. This allows a geometrical and stator and rotor resistance scaling as well as a rewinding of a squirrel cage induction machine enabling an efficient numerical optimization. The 2D FEM solutions of a reference machine are calculated by a model based hybrid numeric induction machine simulation approach. In contrast to already known scaling procedures for synchronous machines the FEM solutions of the induction machine are scaled in the stator-current-rotor-frequency-plane and then transformed to the torque- speed-map. This gives the possibility to use a new time scaling factor that is necessary to keep a constant field distribution. The scaling procedure is validated by the finite element method and used in a numerical optimization process for the sizing of an electric vehicle traction drive considering the gear ratio. The results show that the scaling procedure is very accurate, computational very efficient and suitable for the use in machine design optimization.
This research presents a method for the simulation of the magneto-mechanical system dynamics taking motion and eddy currents into account. The major contribution of this work leans on the coupling the field-motion problem considering windings as the current forced massive conductors, modelling of the rotor motion composed of two conductive materials and the torque calculation employing the special optimal predictor combined with the modified Maxwell stress tensor method. The 3D model of the device is analysed by the time stepping finite element method. Mechanical motion of the rotor is determined by solving the second order motion equation. Both magnetic and mechanical equations are coupled in the iterative solving process. Presented method is verified by solving the TEAM Workshop Problem 30.
Dry electrostatic precipitators (ESPs) are widely used for purification of exhaust gases in industrial applications. Maintenance of their high efficiency depends primarily on periodical cleaning of the collecting electrodes (CEs). Dust removal (regeneration of CEs) is realized by inducing periodical vibrations of the electrodes. The paper presents results of vibration modelling of a system of CEs; the results were obtained by means of the finite element method, the hybrid finite element method, the finite strip method and a model formulated using Abaqus. Numerical results are compared with those obtained from experimental measurements. Conclusions concerning numerical effectiveness and exactness of the methods are formulated and reasons for differences are discussed.
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 finite element method (FEM) is one of the most frequently used numerical methods for finding the approximate discrete point solution of partial differential equations (PDE). In this method, linear or nonlinear systems of equations, comprised after numerical discretization, are solved to obtain the numerical solution of PDE. The conjugate gradient algorithms are efficient iterative solvers for the large sparse linear systems. In this paper the performance of different conjugate gradient algorithms: conjugate gradient algorithm (CG), biconjugate gradient algorithm (BICG), biconjugate gradient stabilized algorithm (BICGSTAB), conjugate gradient squared algorithm (CGS) and biconjugate gradient stabilized algorithm with l GMRES restarts (BICGSTAB(l)) is compared when solving the steady-state axisymmetric heat conduction problem. Different values of l parameter are studied. The engineering problem for which this comparison is made is the two-dimensional, axisymmetric heat conduction in a finned circular tube.
The paper presents a model of a rapping system of an electrostatic precipitator. The rapping system consists of a set of collecting electrodes hanging on a suspension bar and braced together in a brushing bar. The suspension and brushing bars are modeled using the rigid finite element method, while the collecting plates are modeled using the hybrid method. The method combines the rigid finite element method with the classical finite element method. As a result, the mass matrix is diagonal. Some results of numerical simulations concerning free vibrations of the collecting plates and the influence of the number of elements, into which the plate is divided, on the vibrations of the rapping system are presented.
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.
Various methods for steady-state and transient analysis of temperature distribution and efficiency of continuous-plate fins are presented. For a constant heat transfer coefficient over the fin surface, the plate fin can be divided into imaginary rectangular or hexangular fins. At first approximate methods for determining the steady-state fin efficiency like the method of equivalent circular fin and the sector method are discussed. When the fin geometry is complex, thus transient temperature distribution and fin efficiency can be determined using numerical methods. A numerical method for transient analysis of fins with complex geometry is developed. Transient temperature distributions in continuous fins attached to oval tubes is computed using the finite volume - finite element methods. The developed method can be used in the transient analysis of compact heat exchangers to calculate correctly the heat flow rate transferred from the finned tubes to the fluid.
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.
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 article presents a numerical model of the concrete heat accumulator for solar heating systems. Model uses control volume finite element method with an explicit solution method for time integration. The use of an explicit method is an essential advantage in the simulation of time-dependent changes in temperature of the air at the accumulator inlet. The study compares the results of numerical model calculations of the accumulator heating with experimental measurements and with computational fluid dynamics modeling. The comparison shows a good correlation between the results of calculation using the model and the results of measurements.
Miniature heat exchangers are used to provide higher cooling capacity for new technologies. This means a reduction in their size and cost but the identical power. The paper presents the method for determination of boiling heat transfer coefficient for a rectangular minichannel of 0.1 mm depth, 40 mm width and 360 mm length with asymmetric heating. Experimental research has focused on the transition from single phase forced convection to nucleate boiling, i.e., the zone of boiling incipience. The ‘boiling front’ location has been determined from the temperature distribution of the heated wall obtained from liquid crystal thermography. The experiment has been carried out with R-123, mass flux 220 kg/(m2s), pressure at the channel inlet 340 kPa. Local values of heat transfer coefficient were calculated on the basis of empirical data from the experiment following the solution of the two-dimensional inverse heat transfer problem. This problem has been solved with the use of the finite element method in combination with Trefftz functions. Temperature approximates (linear combinations of Trefftz functions) strictly fulfill the governing equations. In presented method the inverse problem is solved in the same way as the direct problem. The results confirmed that considerable heat transfer enhancement takes place at boiling incipience in the minichannel flow boiling. Moreover, under subcooling boiling, local heat coefficients exhibit relatively low values.