The aim of the presented paper is to show the results of shape optimization of railway polynomial transition curves (TCs) of 5th, 7th, and 9th degrees through the use of the full vehicle model and new criteria of assessement concerning the jerk value. The search for the proper shape of TCs means that in this work, the evaluation of TC properties is based on select quantities and the generation of such a shape through the use of mathematically understood optimization methods. The studies presented have got a character of the numerical tests. For this work, advanced vehicle models describing dynamical track-vehicle and vehicle-passenger interactions as well as optimization methods were exploited. In the software vehicle model of a 2-axle freight car, the track discrete model, non-linear descriptions of wheel-rail contact are applied. This part of the software, the vehicle simulation software, is combined with a library optimization procedure into the final computer program.
This paper addresses the problem of road safety regarding barrier placement as relative to the curb. A short summary of existing regulations is presented. Numerical simulations using the explicit finite element system Ls-Dyna are shown. In the analysis, variable distance between the barrier and the curb is assumed. The obtained result reveals that the distance has little impact on the working width of the barrier.
The article presents an analysis of stresses in the current tool system of the die during the implementation of the third forging operation of the screw M12 class 10.9 with cylinder head and hexagonal socket. It was assumed that the level of negative cracking due to stress can be reduced by using a mounting interference between the die and the tube blank. Due to the design of the tool system value of the die, the interference value cannot be too large. Therefore, an analysis of the influence of the interference between the die and the tube blank in a die tool system on the value and distribution of stresses in the individual components. An analysis of the assembly stresses and the stresses occurring during the process of deformation of the shaped head of the screw was done. The calculations were performed using a commercial software package MARC / Mentat.
The object of the present study is to investigate the influence of damping uncertainty and statistical correlation on the dynamic response of structures with random damping parameters in the neighbourhood of a resonant frequency. A Non-Linear Statistical model (NLSM) is successfully demonstrated to predict the probabilistic response of an industrial building structure with correlated random damping. A practical computational technique to generate first and second-order sensitivity derivatives is presented and the validity of the predicted statistical moments is checked by traditional Monte Carlo simulation. Simulation results show the effectiveness of the NLSM to estimate uncertainty propagation in structural dynamics. In addition, it is demonstrated that the uncertainty in damping indeed influences the system response with the effects being more pronounced for lightly damped structures, higher variability and higher statistical correlation of damping parameters.
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.
A description of direct simulation of crosswind loads caused by critical vortex excitation and the response of the structure to these loads are presented in this paper. Tower-like structures of circular cross-sections are considered. A proposed mathematical model of vortex excitation has been numerically implemented and a selfserving computer program was created for the purpose. This software, cooperating with the FEM system, allows for a simulation of a crosswind load and lateral response in real time, meaning that at each time step of the calculations the load is generated using information regarding displacements seen beforehand. A detailed description of the mathematical model is neglected in this paper, which is focused on numerical simulations. WAWS and AR methods are used in simulations.
This study aims to design a novel air cleaning facility which conforms to the current situation in China, and moreover can satisfy our demand on air purification under the condition of poor air quality, as well as discuss the development means of a prototype product. Air conditions in the operating room of a hospital were measured as the research subject of this study. First, a suitable turbulence model and boundary conditions were selected and computational fluid dynamics (CFD) software was used to simulate indoor air distribution. The analysis and comparison of the simulation results suggested that increasing the area of air supply outlets and the number of return air inlets would not only increase the area of unidirectional flow region in main flow region, but also avoid an indoor vortex and turbulivity of the operating area. Based on the summary of heat and humidity management methods, the system operation mode and relevant parameter technologies as well as the characteristics of the thermal-humidity load of the operating room were analyzed and compiled. According to the load value and parameters of indoor design obtained after our calculations, the airflow distribution of purifying the air-conditioning system in a clean operating room was designed and checked. The research results suggested that the application of a secondary return air system in the summer could reduce energy consumption and be consistent with the concept of primary humidity control. This study analyzed the feasibility and energy conservation properties of cleaning air-conditioning technology in operating rooms, proposed some solutions to the problem, and performed a feasible simulation, which provides a reference for practical engineering.
The quantitative description of an airlift bioreactor, in which aerobic biodegradation limited by carbonaceous substrate and oxygen dissolved in a liquid takes place, is presented. This process is described by the double-substrate kinetics. Mathematical models based on the assumption of plug flow and dispersion flow of liquid through the riser and the downcomer in the reactor were proposed. Calculations were performed for two representative hydrodynamic regimes of reactor operation, i.e. with the presence of gas bubbles only within the riser and for complete gas circulation. The analysis aimed at how the choice of a mathematical model of the process would enable detecting the theoretical occurrence of oxygen deficiency in the airlift reactor. It was demonstrated that the simplification of numerical calculations by assuming the “plug flow” model instead of dispersion with high Péclet numbers posed a risk of improper evaluation of the presence of oxygen deficiency zones. Conclusions related to apparatusmodelling and process design were drawn on the basis of the results obtained. The paper is a continuation of an earlier publication (Grzywacz, 2012a) where an analysis of single-substrate models of the airlift reactor was presented.
CFD modelling of momentum and heat transfer using the Large Eddy Simulation (LES) approach has been presented for a Kenics static mixer. The simulations were performed with the commercial code ANSYS Fluent 15 for turbulent flow of three values of Reynolds number, Re = 5 000, 10 000 and 18 000. The numerical modelling began in the RANS model, where standard k−ε turbulence model and wall functions were used. Then the LES iterations started from the initial velocity and temperature fields obtained in RANS. In LES, the Smagorinsky–Lilly model was used for the sub-grid scale fluctuations along with wall functions for prediction of flow and heat transfer in the near-wall region. The performed numerical study in a Kenics static mixer resulted in highly fluctuating fields of both velocity and temperature. Simulation results were presented and analysed in the form of velocity and temperature contours. In addition, the surface-averaged heat transfer coefficient values for the whole insert length were computed and compared with the literature experimental data. Good compliance of the LES simulation results with the experimental correlation was obtained.
Rock excavation is a basic technological operation during tunnelling and drilling roadways in underground mines. Tunnels and roadways in underground mines are driven into a rock mass, which in the particular case of sedimentary rocks, often have a layered structure and complicated tectonics. For this reason, rock strata often have highly differentiated mechanical properties, diverse deposition patterns and varied thicknesses in the cross sections of such headings. In the field of roadheader technology applied to drilling headings, the structure of a rock mass is highly relevant when selecting the appropriate cutting method for the heading face. Decidedly differentiated values of the parameters which describe the mechanical properties of a particular rock layer deposited in the cross section of the drilled tunnel heading will influence the value and character of the load on the cutting system, generated by the cutting process, power demand, efficiency and energy consumption of the cutting process. The article presents a mathematical modelling process for cutting a layered structure rock mass with the transverse head of a boom-type roadheader. The assumption was made that the rock mass being cut consists of a certain number of rock layers with predefined mechanical properties, a specific thickness and deposition pattern. The mathematical model created was executed through a computer programme. It was used for analysing the impact deposition patterns of rock layers with varied mechanical properties, have on the amount of cutting power consumed and load placed on a roadheader cutting system. The article presents an example of the results attained from computer simulations. They indicate that variations in the properties of the rock cut – as cutting heads are moving along the surface of the heading face – may have, apart from multiple other factors, a significant impact on the value of the power consumed by the cutting process.
This article considers designing of a renewable electrical power generation system for self-contained homes away from conventional grids. A model based on a technique for the analysis and evaluation of two solar and wind energy sources, electrochemical storage and charging of a housing area is introduced into a simulation and calculation program that aims to decide, based on the optimized results, on electrical energy production system coupled or separated from the two sources mentioned above that must be able to ensure a continuous energy balance at any time of the day. Such system is the most cost-effective among the systems found. The wind system adopted in the study is of the low starting speed that meets the criteria of low winds in the selected region under study unlike the adequate solar resource, which will lead to an examination of its feasibility and profitability to compensate for the inactivity of photovoltaic panels in periods of no sunlight. That is a system with fewer photovoltaic panels and storage batteries whereby these should return a full day of autonomy. Two configurations are selected and discussed. The first is composed of photovoltaic panels and storage batteries and the other includes the addition of a wind system in combination with the photovoltaic system with storage but at a higher investment cost than the first. Consequently, this result proves that is preferable to opt for a purely photovoltaic system supported by the storage in this type of site and invalidates the interest of adding micro wind turbines adapted to sites with low wind resources.
The engine simulations have become an integral part of engine design and development. They are based on approximations and assumptions. The precision of the results depends on the accuracy of these hypotheses. The simplified models of frozen composition, chemical equilibrium and chemical kinetics provide the compositions of combustion products for engine cycle simulations. This paper evaluates the effects of different operating conditions and hypotheses on the exergetic analysis of a spark-ignition engine. The Brazilian automotive market has the highest number of flex-fuel vehicles. Therefore, a flex-fuel engine is considered for simulations in order to demonstrate the effects of these different hypotheses. The stroke length and bore diameter have the same value of 80 mm. The in-cylinder irreversibility is calculated for each case at the closed part of the engine cycle. A comparative analysis of these hypotheses provides a comprehensive evaluation of their effects on exergetic analysis. Higher values of accumulated irreversibility are observed for the oversimplified hypothesis.
The mathematical model and numerical simulations of the solidification of a cylindrical shaped casting, which take into account the process of filling the mould cavity by liquid metal and feeding the casting through the riser during its solidification, are presented in the paper. Mutual dependence of thermal and flow phenomena were taken into account because have an essential influence on solidification process. The effect of the riser shape on the effectiveness of feeding of the solidifying casting was determined. In order to obtain the casting without shrinkage defects, an appropriate selection of riser shape was made, which is important for foundry practice. Numerical calculations of the solidification process of system consisting of the casting and the conical or cylindrical riser were carried out. The velocity fields have been obtained from the solution of momentum equations and continuity equation, while temperature fields from solving the equation of heat conductivity containing the convection term. Changes in thermo-physical parameters as a function of temperature were considered. The finite element method (FEM) was used to solve the problem.
This work presents the methodology for analyzing the impact of ground vibrations induced during the drilling of gas/oil exploration wells on the surrounding constructions, as well as on humans and the natural environment. In the primary stage, this methodology is based on measurements of ground vibrations induced by a specific type of drilling system in the so-called reference site. In the next stage, ground vibrations are estimated in similar conditions to another design site, these conditions are assumed for a given drilling system, treated as a vibration source. In both sites, special seismic and geotechnical data are collected to construct numerical models for dynamic analyses. Finally, if it is required, a protection system is proposed with respect to the drilling technology and local conditions. The methodology presented has been tested on the terrain of an active natural gas mine used as the design site, and located in the southeastern part of Poland. The reference site was placed in the terrain of a working drilling system in similar conditions in the central part of Poland. Based on the results of numerical simulations, one may verify the different locations of the drilling rig in the design site with respect to the existing industrial structure. Due to the hazard from destructive ground vibrations, a certain vibroisolation system was proposed at the design site. Based on the results of numerical simulations one could rearrange the components of the drilling system in order to provide maximum security for the surrounding structures.
The article presents the use of computer graphics methods and computational geometry for the analysis on changes of geometrical parameters for a mixed zone in resistance-heated samples. To perform the physical simulation series of resistance heating process, the Gleeble 3800 physical simulator, located in the Institute for Ferrous Metallurgy in Gliwice, was used. The paper presents a description of the test stand and the method for performing the experiment. The numerical model is based on the Fourier-Kirchoff differential equation for unsteady heat flow with an internal volumetric heat source. In the case of direct heating of the sample, geometrical parameters of the remelting zone change rapidly. The described methodology of using shape descriptors to characterise the studied zone during the process allows to parametrise the heat influence zones. The shape descriptors were used for the chosen for characteristic timing steps of the simulation, which allowed the authors to describe the changes of the studied parameters as a function of temperature. Additionally, to determine the impact of external factors, the remelting zone parameters were estimated for two types of grips holding the sample, so-called hot grips of a shorter contact area with the sample, and so-called cold grips. Based on the collected data, conclusions were drawn on the impact of the process parameters on the localisation and shape of the mushy zone.
Scientific research discussed in the present article is focused on the determination of the vertical conveyance capacity in the process of mining minerals, while applying a mathematical calculation and verification of the calculation results by simulation. Input parameters for the capacity calculation include the transport cycle time. The article presents the results of measuring a transport cycle during the operation and a calculation of the transport cycle while using known formulas. On the basis of the observed findings, two methods of increasing the hoisting machine capacity were proposed. The first method is increasing the velocity from the original value of 6 m.s–1 to the velocity of 7 m.s–1. In this case, we achieved the daily capacity increase in 2-2.5%. The second method consisted in changing the hoisting machine acceleration and deceleration modes by which we achieved as much as 9% increase in the daily capacity. The article also describes a transport cycle simulation model, with its output being the number of work cycles parameter. The obtained parameter was used again in the capacity calculation. The simulation model was used in experiments for both, the current status as well as proposed solutions. The simulation model serves also for calculation verification.
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.
According to the European Environment Agency (EEA 2018), air quality in Poland is one of the worst in Europe. There are several sources of air pollution, but the condition of the air in Poland is primarily the result of the so-called low-stack emissions from the household sector. The main reason for the emission of pollutants is the combustion of low-quality fuels (mainly low-quality coal) and waste, and the use of obsolete heating boilers with low efficiency and without appropriate filters. The aim of the study was to evaluate the impact of measures aimed at reducing low-stack emissions from the household sector (boiler replacement, change of fuel type, and thermal insulation of buildings), resulting from environmental regulations, on the improvement of energy efficiency and the emission of pollutants from the household sector in Poland. Stochastic energy and mass balance models for a hypothetical household, which were used to assess the impact of remedial actions on the energy efficiency and emission of pollutants, have been developed. The annual energy consumption and emissions of pollutants were estimated for hypothetical households before and after the implementation of a given remedial action. The calculations, using the Monte Carlo simulation, were carried out for several thousand hypothetical households, for which the values of the technical parameters (type of residential building, residential building area, unitary energy demand for heating, type of heat source) were randomly drawn from probability distributions developed on the basis of the analysis of the domestic structure of households. The model takes the coefficients of correlation between the explanatory variables in the model into account. The obtained results were multiplied so that the number of hypothetical households was equal to 14.1 million, i.e. the real number of households in Poland. The obtained results allowed for identifying the potential for reducing the emission of pollutants such as carbon dioxide, carbon monoxide, dust, and nitrogen oxides, and improving the energy efficiency as a result of the proposed and implemented measures, aimed at reducing low-stack emission, resulting from the policy. The potential for emissions of gaseous pollutants is 94% for CO, 49% for NOx, 90% for dust, and 87% for SO2. The potential for improving the energy efficiency in households is around 42%.
Simplifications used in simulation program codes require the use of substitute parameters in the material databases (also called apparent or substitutive). On the one hand, they formally fit into the records used in the heat flow model, porosity, properties etc. and on the other hand they should be determined in conditions most similar to the real casting-mould system. The article presents results of a research on thermophysical parameters of gypsum mould used for precision casting moulds. Experiments were carried out on a cylindrical mould made of Plasticast gypsum, in which the heat source was a cylinder filled with liquid aluminium alloy of the temperature of 655°C. Energetic validation was carried out by using the NovaFlow&Solid ver. 6.3 simulation code. As a result of validation tests, substitute thermophysical parameters of gypsum were determined. For determined parameters, best-fit of solidification time from the experiment and simulation was obtained and the curves of gypsum mass heating were satisfactorily recreate.
The irregularity profiles of steel samples after vapour blasting were measured. A correlation analysis of profile parameters was then carried out. As the result, the following parameters were selected: Pq, Pt, PDq, Pp/Pt and Pku. Surface profiles after vapour blasting were modeled. The modeled surfaces were correctly matched to measured surfaces in 78% of all analyzed cases. The vapour blasting experiment was then carried out using an orthogonal selective research plan. The distance between the nozzle and sample d and the pressure of feed system p were input parameters; selected surface texture coefficients were output parameters. As the result of the experiment, regression equations connecting vapour blasting process parameters p and d with selected profile parameters were obtained. Finally, 2D profiles of steel samples were forecasted for various values of vapour blasting parameters. Proper matching accuracy of modeled to measured profiles was assured in 75% of analyzed cases.
This paper proposes a new approach for calculating the impulse response of room acoustics. Impulse response provides unique characterization of any discrete lineartime invariant (LTI) systems. We assume that the room is a linear time-invariant system and the impulse response is calculated simply by sending a Dirac Impulse into the system as input and getting the response from the output. Then, the output of the system is represented as a sum of time-shifted weighted impulse responses. Both mathematical justifications for the proposed method and results from simulation software developed to evaluate the proposed approach are presented in detail.
The paper presents results of the field tests on membrane biogas enrichment performed with the application of mobile membrane installation (MMI) with the feed stream up to 10 Nm3/h. The mobile installation equipped with four hollow fibre modules with polyimide type membranes was tested at four different biogas plants. Two of them were using agricultural substrates. The third one was constructed at a municipal wastewater plant and sludge was fermented in a digester and finally in the fourth case biogas was extracted from municipal waste landfill site. Differences in the concentration of bio-methane in feed in all cases were observed and trace compounds were detected as well. High selectivity polyimide membranes, in proper module arrangements, can provide a product of high methane content in all cases. The content of other trace compounds, such as hydrogen sulphide, water vapour and oxygen on the product did not exceed the values stated by standard for a biogas as a vehicle fuel. The traces of hydrogen sulphide and water vapour penetrated faster to the waste stream enriched in carbon dioxide, which could lead to further purification of the product – methane being hold in the retentate (H2O > H2S > CO2 > O2 > CH4 > N2). In the investigated cases, when concentration of N2 was low and concentration of CH4 higher than 50%, it was possible to upgrade methane to concentration above 90% in a two-stage cascade. To performsimulation ofCH4 andCO2 permeation through polyimide membrane,MATLABwas used. Simulation program has included permeation gaseous mixture with methane contents as observed at field tests in the range of 50 and 60% vol. The mass transport process was estimated for a concurrent hollow fibre membrane module for given pressure and temperature conditions and different values of stage cut. The obtained results show good agreement with the experimental data. The highest degree of methane recovery was obtained with gas concentrating in a cascade with recycling of the retentate.
Dynamic Mine disasters can be induced by the instability and failure of a composite structure of rock and coal layers during coal mining. Coal seam contains many native defects, severely affecting the instability and failure of the compound structure. In this study, the effects of coal persistent joint on the strength and failure characteristics of coal-rock composite samples were evaluated using PFC2D software. The results show that with the increase of included angle α between the loading direction and joint plane direction, the uniaxial compressive stress (UCS) and peak strain of composite samples first decrease and then gradually increase. The elastic moduli of composite samples do not change obviously with α. The peak strain at α of 45° is the lowest, and the UCS at α of 30° is the smallest. This is inconsistent with theoretical analysis of lowest UCS at α of 45°. This is because that the local stress concentration caused by the motion inconformity of composite samples may increase the average axial stress of upper wall in PFC2D software. Moreover, the coal persistent joint promotes the transformation from the unstable crack expansion to the macro-instability of composite samples, especially at α of 30° and 45°. The majority of failures for composite samples occur within the coal, and no obvious damage is observed in rock. Their failure modes are shear failure crossing or along the coal persistent joint. The failure of composite sample at α of 30° is a mixed failure, including the shear failure along the persistent joint in coal and tensile failure of rock induced by the propagation of coal persistent joint.
The cohesion and internal friction angle were characterized as quadratic functions of strain and were assumed to follow the Mohr-Coulomb criterion after the yield of peak strength. These mechanical parameters and their variations in post-peak softening stage can be exactly ascertained through the simultaneous solution based on the data points of stress-strain curves of triaxial compression tests. Taking the influence of the fault into account, the variation of strata pressure and roadway convergence with coal advancement, the temporal and spatial distribution of axial bolt load were numerically simulated by FLAC3D (Fast Lagrangian Analysis of Continua) using the ascertained post-peak mechanical parameters according to the cohesion weakening and friction strengthening model. The change mechanism of axial load of single rock bolt as abutment pressure changes was analyzed, through the comparison analysis with the results of axial bolt load by field measurements at a coal mine face. The research results show that the simulated results such as the period of main roof weighting, temporal and spatial distribution of axial bolt load are in accordance with field measurement results, so the validity of the numerical model is testified. In front of the working face, the front abutment pressure increases first and then decreases, finally tends to be stable. A corresponding correlation exists between the variation of axial bolt load and rock deformation along the bolt body. When encountered by a fault, the maximum abutment pressure, the influential range of mining disturbance and the roadway convergence between roof and floor before the working face are all increased. In the roadways along the gob, axial bolt loads on the side of the working face decrease, while the other side one increases after the collapse of the roof. As superficial surrounding rock mass is damaged, the anchoring force of rock bolts will transfer to inner rock mass for balancing the tensile load of the bolts.
The aim of this paper is to present the procedure test for calibration and validation of the numerical model for X22CrMoV12-1 steel multilayer welding. On the real multilayer weld was described how to arrange the whole experiment in order to obtain not only relevant input data but also verification data. Tests on a specially prepared specimen, welded with 8 beads in 4 layers, allows to determine the actual geometry of the single welded beads, registration of welding thermal cycles and the hardness distribution in successively deposited beads together with determining the heat influence of subsequent layers. The results of the real welding tests were compared with the results obtained from the numerical simulations and extended by the calculated stresses and distortions distributions of the tested specimen. A new, improved hardness prediction algorithm for high-alloy martensitic and bainitic steels was also proposed.