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.
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.
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.
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.
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 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.
The work reports on the development of random three-dimensional Laguerre-Voronoi computational models for open cell foams. The proposed method can accurately generate foam models having randomly distributed parameter values. A three-dimensional model of ceramic foams having pre-selected cell volumes distribution with stochastic coordinates and orientations was created in the software package ANSYSTM. Different groups of finite element models were then generated using the developed foam modeling procedure. The size sensitivity study shows that each of foam specimens at least contains 125 LV-cells. The developed foam models were used to simulate the macroscopic elastic properties of open cell foams under uni-axial and bi-axial loading and were compared with the existing open cell foam models in the literature. In the high porosity regime, it is found that the elastic properties predicted by random Laguerre-Voronoi foam models are almost the same as those predicted by the perfect Kelvin foam models. In the low porosity regime the results of the present work deviate significantly from those of other models in the literature. The results presented here are generally in better agreement with experimental data than other models. Thus, the Laguerre-Voronoi foam models generated in this work are quite close to real foam topology and yields more accurate results than other open cell foam models.
An alternative FEM algorithm of fi nding piston ring pressure distribution to a contact simulation is introduced. The method is basing on an analytical determining of required nodal displacement boundary conditions. Its several confi gurations are tested using APDL and compared to a no-separation contact simulation of a simple 2D fi nite element model of a two-stroke piston ring made of Titanium alloy. Each of the methods tested in the paper brings displacement result and Huber-Misses equivalent stresses close to each other. However, only one of those brings resulting contact pressure close to a no-separation contact simulation. Nonetheless, the obtained confi guration occurred to be less computationally effi cient than no- separation contact simulation performed in an ANSYS software.
The aim of this paper is to present methods of digitally synthesising the sound generated by vibroacoustic systems with distributed parameters. A general algorithm was developed to synthesise the sounds of selected musical instruments with an axisymmetrical shape and impact excitation, i.e., Tibetan bowls and bells. A coupled mechanical-acoustic field described by partial differential equations was discretized by using the Finite Element Method (FEM) implemented in the ANSYS package. The presented synthesis method is original due to the fact that the determination of the system response in the time domain to the pulse (impact) excitation is based on the numerical calculation of the convolution of the forcing function and impulse response of the system. This was calculated as an inverse Fourier transform of the system’s spectral transfer function. The synthesiser allows for obtaining a sound signal with the assumed, expected parameters by tuning the resonance frequencies which exist in the spectrum of the generated sound. This is accomplished, basing on the Design of Experiment (DOE) theory, by creating a meta-model which contains information on its response surfaces regarding the influence of the design parameters. The synthesis resulted in a sound pressure signal in selected points in space surrounding the instrument which is consistent with the signal generated by the actual instruments, and the results obtained can improve them.
Nominal strength reduction in cross ply laminates of [0/90]2s is observed in tensile tests of glass fiber composite laminates having central open hole of diameters varying from 2 to 10 mm. This is well known as the size effect. The extended finite element method (XFEM) is implemented to simulate the fracture process and size effect (scale effect) in the glass fiber reinforced polymer laminates weakened by holes or notches. The analysis shows that XFEM results are in good agreement with the experimental results specifying nominal strength and in good agreement with the analytical results based on the cohesive zone model specifying crack opening displacement and the fracture process zone length.
The paper deals with experimental investigations of a set of metal wave-ring gaskets of different thickness and different assembly interference. The gaskets were examined under assembly conditions, i.e. pressed in their seats with no operating pressure applied. The electric resistance wire strain gauges were used to measure the circumferential and axial strains at the inner surface of the gaskets. The traces of contact at the working surface of the gaskets were measured after disassembly the gaskets from their seats. The material tests were carried out to determine the real mechanical properties of materials applied for the gaskets and the seats. The results of experiment were verified by FEM calculations and compared with the analytical approach based on the simplified shell model proposed for the gasket.