Exploitation of hard coal seams by roadway system is applied by two coal mines in southern Poland in Upper Silesian Basin. It is a secondary mining exploitation carries out in safety pillars of urban areas and shafts within mining areas of closed coal mines. Roadway system is the excavation process of gateways which are made in parallel order leaving coal pillars between them. An optimal width of coal pillar makes roadway stable and reduces subsidence of terrain surface. The article presents results of subsidence simulation caused by partial extraction using empirical and numerical methods on the example of one exploitation field of “Siltech” coal mine. The asymptotic state of subsidence was considered after mining ceased in the study area. In order to simulate of subsidence, numerical model of rock mass and model of Knothe-Budryk theory were calibrated. Simulation of vertical displacements in numerical method was carried out using RS3 program by Rocscience based on finite element method. The assumption was made that model of rock mass is transversely isotropic medium, in which panels were designed according to order of extraction of coal seams. The results of empirical and numerical methods were compared with measured values of subsidence at benchmarks along drawn lines (subsidence profiles).
This paper presents new procedure modeling based on finite element method analysis of wood-framed timber structures. The fasteners linking boards of sheathing with the timber frame both modeled applying shell finite element, with individual material parameters, remain the main objective of this manuscript. Material parameters are obtained from experimental tests and numerical identification. The main objective of the paper is the elaboration of the numerical model with high precision of mapping, and, at the same time, diminishing the number of the unknown simplifying the process of the modeling of timber structures. The new presented method leads to a simplification of analysis of multistory wood-framed multifamily building structures.
The standard PN-EN_1993-1-5:_2008 (Eurocode 3) compared with the standard (PN-B-03200:_1990) used previously in Poland, introduces extended rules referring to the computations of the bearing capacity of the plated structural elements including the shear lag effect. The stress distribution in the width flanges is variable. Therefore in the case of the beam with the shear lag effect cannot be calculated by the classic beam theory. In this article a comparison of the results of the calculations of forces distribution, stresses and displacement according to the rule presented in PN-EN_1993 and results of the numerical computations for_3D model (using finite element method) is presented. The elastic shear lag effects, the elastic shear lag effects including effects of the plate buckling and the elastic-plastic shear lag effects including the local instabilities were analysed. The calculations were performed for beams with a small and a large span and an influence of stiffeners was analysed.
The paper presents an analysis of the influence of elevated temperature on thin-walled purlins restrained by sheeting. In the first part of the study the bearing capacity of purlins cooperating with sheeting is examined in normal and elevated temperature based on European Standards. Next, special attention is paid to creating a numerical FEM model of the restrained purlins in Abaqus program taking into account different materials properties with respect to temperature increase.
This paper presents a complex study of anhydrite interbeds influence on the cavern stability in the Mechelinki salt deposit. The impact of interbeds on the cavern shape and the stress concentrations were also considered. The stability analysis was based on the 3D numerical modelling. Numerical simulations were performed with use of the Finite Difference Method (FDM) and the FLAC3D v. 6.00 software. The numerical model in a cuboidal shape and the following dimensions: length 1400, width 1400, height 1400 m, comprised the part of the Mechelinki salt deposit. Three (K-6, K-8, K-9) caverns were projected inside this model. The mesh of the numerical model contained about 15 million tetrahedral elements. The occurrence of anhydrite interbeds within the rock salt beds had contributed to the reduction in a diameter and irregular shape of the analysed caverns. The results of the 3D numerical modelling had indicated that the contact area between the rock salt beds and the anhydrite interbeds is likely to the occurrence of displacements. Irregularities in a shape of the analysed caverns are prone to the stress concentration. However, the stability of the analysed caverns are not expected to be affected in the assumed operation conditions and time period (9.5 years).
The numerical algorithm of thermal phenomena is based on the solution of the heat conduction equations in Petrov-Galerkin’s formula using the finite element method. In the modeling of phase transformation in the solid state, the models based on the diagrams of continuous heating and continuous cooling (CHT and CCT). In the modeling of mechanical phenomena, equations of equilibrium and constitutive relationships were adopted in the rate form. It was assumed that the hardened material is elastic-plastic, and the plasticizing can be characterized by isotropic, kinematic or mixed strengthening. In the model of mechanical phenomena besides thermal, plastic and structural strains, the transformations plasticity was taken into account. Thermo-physical size occurring in the constitutive relationship, such as Young’s modulus and tangential modulus, while yield point depend on temperature and phase composition of the material. The modified Leblond model was used to determine transformation plasticity. This model was supplemented by an algorithm of modified plane strain state, advantageous in application to the modeling of mechanical phenomena in slender objects. The problem of thermoelasticity and plasticity was solved by the FEM. In order to evaluate the quality and usefulness of the presented numerical models, numerical analysis of temperature fields, phase fractions, stresses and strains was performed, i.e. the basic phenomena accompanying surface layer of progressive-hardening with a movable heat source of slender elements made of tool steel for cold work.
Deep excavation walls can be analyzed and calculated by using classical methods (currently rarely in use due to their many simplifications) or numerical methods. Among the numerical methods we can distinguish a simplified approach, in which the interaction between soil and a wall structure is modelled by a system of elasto-plastic supports, and the finite-element method (FEM) in which the soil is modelled with mesh of elements. It is a common view that if we want to analyze only wall constructions, the first, simplified method of calculation is sufficient. The second method, FEM, is required if we want to further analyze the stress and strain states in the soil and the influence of the excavation on the surrounding area. However, as it is demonstrated in the paper, important differences may appear in the calculation results of both methods. Thus, the safety design of a deep excavation structure depends very much on the choice of calculating method.
In the paper, the method of a numerical simulation concerning diagonal crack propagation in con-crete beams was presented. Two beams reinforced longitudinally but without shear reinforcement were considered during the Finite Element Method analysis. In particular, a nonlinear method was used to simulate the crack evaluation in the beams. The analysis was performed using the commercial program ANSYS. In the numerical simulation, the limit surface for concrete described by Willam and Warnke was applied to model the failure of concrete. To solve the FEM-system of equations, the Newton-Raphson method was used. As the results of FEM calculations, the trajectories of total stains and numerical images of smeared cracks were obtained for two analyzed beams: the slender beam S5 of leff = 1.8 m and the short beam S3k of leff = 1.1 m. The applied method allowed to generate both flexural vertical cracks and diagonal cracks in the shear regions. Some differences in the evaluation of crack patterns in the beams were observed. The greater number of flexural vertical cracks which penetrated deeper in the beam S5 caused the lower stiffness and the greater deformation in the beam S5 compared to the short beam S3k. Numerical results were compared with the experimental data from the early tests performed by Słowik . The numerical simulation yielded very similar results as the experiments and it confirmed that the character of failure process altered according to the effective length of the member. The proposed numerical procedure was successfully verified and it can be suitable for numerical analyses of diagonal crack propagation in concrete beams.
In general, uniform mixing of particles is desirable in the process of particle handling. However, during the charging of sinter feed and upper ore, size segregation must be induced to prevent heat imbalance, ensure bed permeability, and prevent the loss of fine ore. In this study, upper ore charging was simulated using a discrete element method (DEM) to find the optimal method for controlling particle size segregation, and the segregation characteristics in the upper ore bed were investigated when a deflector plate was applied to the charging machine. The degree of vertical segregation increased when a deflector plate was applied, and it was confirmed that the segregation direction in the upper ore bed can be controlled by adjusting the charging direction of the upper ore by using a deflector plate. In order to apply this method directly to the actual process, further study is needed to understand the influence of the characteristics of the deflector plate such as length and angle.
The objective of the paper is to analyse thermodynamical and operational parameters of the supercritical power plant with reference conditions as well as following the introduction of the hybrid system incorporating ORC. In ORC the upper heat source is a stream of hot water from the system of heat recovery having temperature of 90 °C, which is additionally aided by heat from the bleeds of the steam turbine. Thermodynamical analysis of the supercritical plant with and without incorporation of ORC was accomplished using computational flow mechanics numerical codes. Investigated were six working fluids such as propane, isobutane, pentane, ethanol, R236ea and R245fa. In the course of calculations determined were primarily the increase of the unit power and efficiency for the reference case and that with the ORC.
The main objective of this work was to present a successful stabilization action of a building structure in an active landslide. Firstly, history of the case and a FEM simulation explaining ensuing situation are presented. Then different structural measures to stabilize the whole system are discussed. The structural solution of the problem (pile system reaching solid rocky zone) is presented in more detailed way. The estimation of forces acting on the structure, caused by an unstable soil mass, being crucial for the design of stabilizing structure is described.
The paper presents thermodynamic analysis of the gas-steam unit of the 65 MWe combined heat and power station. Numerical analyses of the station was performed for the nominal operation conditions determining the Brayton and combined cycle. Furthermore, steam utilization for the gas turbine propulsion in the Cheng cycle was analysed. In the considered modernization, steam generated in the heat recovery steam generator unit is directed into the gas turbine combustion chamber, resulting in the Brayton cycle power increase. Computational flow mechanics codes were used in the analysis of the thermodynamic and operational parameters of the unit.
The paper presents a dynamic analysis of the damaged masonry building repaired with the Flexible Joint Method. Numerical analysis helped to determine the effect of the applied repairing method on natural frequencies as well as values of stresses and accelerations in the analyzed variants of numerical model. They confirmed efficiency of the proposed repair method.
The Dez dam was commissioned in 1963 and since sediments accumulated in the reservoir up to an elevation of approximately 15m below the intake of the power tunnel. One of the possible measures to improve operation of the reservoir is by heightening of the existing dam. This paper describes the conducted procedure for static and thermal calibration of this 203m dam in Iran based on micro geodesies measurements. Also the nonlinear response of existing dam is investigated under maximum credible earthquake ground motions considering joint behavior and mass concrete cracking and safety of dam is evaluated for possible heightening. For thermal calibration of provided numerical model, transient thermal analysis was conducted and results were compared with thermometers records installed in central block. In addition, for static calibration; thermal distribution within dam body, dam self weight, hydrostatic pressure and silt load applied on the 3D fi nite element model of dam-reservoir-foundation were considered. Results show that the distribution of stresses will be critical within dam for heightening case under seismic loads in MCL.
Because the heat release of plutonium material, the composite structure is heated and the stress and strain of the composite structure will increase, which will affect the thermodynamic properties of the structure. The thermodynamic analysis of complex structures, which are composed of concentric structures of plutonium, beryllium, tungsten, explosives, and steel, was carried out. The results showed that when the structure is spherical, the temperature is higher than that of the ellipsoid structure. Stress of the elliptical structure is greater than the spherical structure. This study showed that the more flat the shell is, the greater the stress concentration point occurs at the long axis, and the maximum stress occurs at the beryllium layer. These conclusions provide theoretical support for the plutonium composite component testing.
Thermo-chemical treatments are known to increase the fatigue life of industrial parts. Due to the imprecise consideration of residual stresses in predicting the durability of components subjected to cyclic loading and their effect on the fatigue life, the authors developed a numerical model combining the influence of residual stresses with stresses caused by bending. The authors performed the numerical simulation with the use of Finite Element Method to analyse material behaviour during cyclic loading. The residual stress state developed during nitriding was introduced onto cross-section of the numerical specimen. The goal of this work was better understanding of the real conditions of the nitride steel fatigue processes and improving the knowledge about numerical predicting of the fatigue life for parts with residual stresses. The results of simulation were compared with plane bending fatigue tests. The presented method indicates the possibility of increasing the accuracy of the fatigue analysis of elements after surface treatment, increasing its certainty and the ability to perform better optimization of service life.
The work deals with experimental and numerical thermodynamic analyses of cross-flow finned tube heat exchangers of the gas-liquid type. The aim of the work is to determine an impact of the gas non-uniform inlet on the heat exchangers performance. The measurements have been carried out on a special testing rig and own numerical code has been used for numerical simulations. Analysis of the experimental and numerical results has shown that the range of the non-uniform air inlet to the considered heat exchangers may be significant and it can significantly affect the heat exchanger efficiency.
The aim of the paper is to revalorizate of the Szewalski binary vapour cycle by analysing thermodynamical and operational parameters of this cycle. This was carried out by accessible numerical computational flow mechanics codes using the step-by-step modeling of separate elements. The binary vapour cycle is providing steam as the working fluid in the high temperature part of the cycle, while another fluid - a low boiling point fluid of low specific volume - as the working substance substituting conventional steam over the temperature range represented by the low pressure steam expansion. The steam cycle for reference conditions has been assumed. Four working fluids in the low temperature part of the binary cycle such as propane, isobutane, ethanol and ammonia have been investigated.
In this paper, thermodynamic analysis of a proposed innovative double Brayton cycle with the use of oxy combustion and capture of CO2, is presented. For that purpose, the computation flow mechanics (CFM) approach has been developed. The double Brayton cycle (DBC) consists of primary Brayton and secondary inverse Brayton cycle. Inversion means that the role of the compressor and the gas turbine is changed and firstly we have expansion before compression. Additionally, the workingfluid in the DBC with the use of oxy combustion and CO2 capture contains a great amount of H2O and CO2, and the condensation process of steam (H2O) overlaps in negative pressure conditions. The analysis has been done for variants values of the compression ratio, which determines the lowest pressure in the double Brayton cycle.
To solve the problem of large deformation soft rock roadway with complicated stress condition in Baluba copper mine, the characteristics of roadway deformation and failure modes are analyzed deeply on the basis of geological survey. Combined with the theoretical analysis and numerical simulation, the new reinforcement technology with floor mudsill and grouting anchor cable is proposed. Moreover, the three dimension numerical simulation model is established by the software FLAC-3D, the support parameter is optimized by it. The results show that the optical array pitch of the U-steel shelf arch is 0.8 m, and the optical array pitch of the grouting anchor cable is 2.4 m. At last, the field experiments are done all over the soft rock roadway. Engineering practice shows that the deformation of soft rock roadway in Baluba copper mine is effectively controlled by adopting the new reinforcement technology, which can provide certain references for similar engineering.
In order to optimize the stope structure parameters in broken rock conditions, a novel method for the optimization of stope structure parameters is described. The method is based on the field investigation, laboratory tests and numerical simulation. The grey relational analysis (GRA) is applied to the optimization of the stope structure parameters in broken rock conditions with multiple performance characteristics. The influencing factors include stope height, pillar diameter, pillar spacing and pillar array pitch, the performance characteristics include maximum tensile strength, maximum compressive strength and ore recovery rate. The setting of influencing factors is accomplished using the four factors four levels Taguchi experiment design method, and 16 experiments are done by numerical simulation. Analysis of the grey relational grade indicates the first effect value of 0.219 is the pillar array pitch. In addition, the optimal stope structure parameters are as follows: the height of the stope is 3.5 m, the pillar diameter is 3.5 m, the pillar spacing is 3 m and the pillar array pitch is 5 m. In-situ measurement shows that all of the pillars can basically remain stable, ore recovery rate can be ensured to be more than 82%. This study indicates that the GRA method can efficiently applied to the optimization of stope structure parameters.
A method of tensile testing of materials in dynamic conditions based on a slightly modified compressive split Hopkinson bar system using a shoulder is described in this paper. The main goal was to solve, with the use of numerical modelling, the problem of wave disturbance resulting from application of a shoulder, as well as the problem of selecting a specimen geometry that enables to study the phenomenon of high strain-rate failure in tension. It is shown that, in order to prevent any interference of disturbance with the required strain signals at a given recording moment, the positions of the strain gages on the bars have to be correctly chosen for a given experimental setup. Besides, it is demonstrated that - on the basis of simplified numerical analysis - an appropriate gage length and diameter of a material specimen for failure testing in tension can be estimated.