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 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.
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.
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.
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.
Technological development offers a wide range of new possibilities for implementation of production processes. Continual production development is the main key to success and competitiveness improvement, labour productivity and image-building for all manufacturing companies. The article deals with designing of new workplace with implementation and utilization of automated robot for faster and safer handling of cast stock. The new layout of workplace is created in software Process Simulate.
The main aim of the presented research was to check mechanical response of human body model under loads that can occur during airplane accidents and compare results of analysis with some results of experimental tests described in literature. In simulations, new multi-purpose human body model, the VIRTHUMAN, was used. The whole model, as well as its particular segments, was earlier validated based on experimental data, which proved its accuracy to simulate human body dynamic response under condition typical for car crashes, but it was not validated for loads with predominant vertical component (loads acting along spinal column), typical for airplane crashes. Due to limitation of available experimental data, the authors focused on conducting calculations for the case introduced in 14 CFR: Parts 23.562 and 25.562, paragraph (b)(1), knowing as the 60 pitch test. The analysis consists in comparison of compression load measured in lumbar section of spine of the FAA HIII Dummy (experimental model) and in the Virthuman (numerical model). The performed analyses show numerical stability of the model and satisfactory agreement between experimental data and simulated Virthuman responses. In that sense, the Virthuman model, although originally developed for automotive analyses, shows also great potential to become valuable tool for applications in aviation crashworthiness and safety analyses, as well.
Eutectoid growth, as the important reaction mechanism of the carbon steel heat treatment, is the basis to control the microstructure and performance. At present, most studies have focused on lamellar growth, and did not consider the nucleation process. Mainly due to the nucleation theory is inconclusive, a lot of research can support their own theory in a certain range. Based on the existing nucleation theory, this paper proposes a cooperative nucleation model to simulate the nucleation process of eutectoid growth. In order to ensure that the nucleation process is more suitable to the theoretical results, different correction methods were used to amend the model respectively. The results of numerical simulation show that when the model is unmodified, the lateral growth of single phase is faster than that of longitudinal growth, so the morphology is oval. Then, the effects of diffusion correction, mobility correction and ledges nucleation mechanism correction on the morphology of nucleation and the nucleation rate were studied respectively. It was found that the introduction of boundary diffusion and the nucleation mechanism of the ledges could lead to a more realistic pearlite.
An on-line optimising control strategy involving a two level extended Kalman filter (EKF) for dynamic model identification and a functional conjugate gradient method for determining optimal operating condition is proposed and applied to a biochemical reactor. The optimiser incorporates the identified model and determines the optimal operating condition while maximising the process performance. This strategy is computationally advantageous as it involves separate estimation of states and process parameters in reduced dimensions. In addition to assisting on-line dynamic optimisation, the estimated time varying uncertain process parameter information can also be useful for continuous monitoring of the process. This strategy ensures that the biochemical reactor is operated at the optimal operation while taking care of the disturbances that are encountered during operation. The simulation results demonstrate the usefulness of the two level EKF assisted dynamic optimizer for on-line optimising control of uncertain nonlinear biochemical systems.
There are certain well-known methods of diminishing concentrations of nitrogen compounds, but they are ineffective in case of nitrogen-rich wastewater with a low content of biodegradable carbon. Partial nitritation followed by anaerobic ammonium oxidation (Anammox) process appear to be an excellent alternative for traditional nitrification and denitrification. This paper presents the feasibility of successful start-up of Anammox process in a laboratory-scale membrane bioreactor (MBR). It was shown that the combination of membrane technology and Anammox process allowed to create a new highly efficient and compact system for nitrogen removal. It was possible to achieve average nitrogen removal efficiency equal to 76.7 ± 8.3%. It was shown that the start-up period of 6 months was needed to obtain high nitrogen removal efficiency. The applied biochemical model of the Anammox process was based on the state-of-the-art Activated Sludge Model No.1 (ASM 1) which was modified for accounting activity of autotrophs (nitrite-oxidising bacteria and nitrateoxidising bacteria) and anammox bacteria. In order to increase the predictive power of the simulation selected parameters of the model were adjusted during model calibration. Readjustment of the model parameters based on the critically evaluated data of the reactor resulted in a satisfactory match between the model predictions and the actual observations.
This study, describing computer simulation of a glider crash against a non-deformable ground barrier, is a part of a larger glider crash modeling project. The studies were intended to develop a numerical model of the pilot - glider - environment system, whereby the dynamics of the human body and the composite cockpit structure during a crash would make it possible to analyze flight accidents with focus on the pilot's safety. Notwithstanding that accidents involving glider crash against a rigid barrier (a wall, for example) are not common, establishing a simulation model for such event may prove quite useful considering subsequent research projects. First, it is much easier to observe the process of composite cockpit structure destruction if the crash is against a rigid barrier. Furthermore, the use of a non-deformable barrier allows one to avoid the errors that are associated with the modeling of a deformable substrate, which in most cases is quite problematic. Crash test simulation, carried out using a MAYMO package, involved a glider crash against a wall positioned perpendicularly to the object moving at a speed of 77 km/h. Computations allowed for determination of time intervals of the signals that are required to assess the behavior of the cockpit and pilot's body - accelerations and displacements in selected points of the glider's structure and loads applied to the pilot's body: head and chest accelerations, forces at femur, lumbar spine and safety belts. Computational results were compared with the results of a previous experimental test that had been designed to verify the numerical model. The glider's cockpit was completely destroyed in the crash and the loads transferred to the pilot's body were very substantial - way over the permitted levels. Since modeling results are fairly consistent with the experimental test, the numerical model can be used for simulation of plane crashes in the future.
The aim of this work is to develop a numerical model capable of predicting the grain density in the Mg-based matrix phase of an AZ91/SiC composite, as a function of the total mass fraction of the embedded SiC particles. Based on earlier work in a range of alloy systems, we assume an exponential relationship between the grain density and the maximum supercooling during solidification. Analysis of data from cast samples with different thicknesses, and mass fractions of added SiCp, permits conclusions to be drawn on the role of SiCp in increasing grain density. By fitting the data, an empirical nucleation law is derived that can be used in a micro model. Numerical simulation based on the model can predict the grain density of magnesium alloys containing SiC particles, using the mass fraction of the particles as inputs. These predictions are compared with measured data.
Experiments of filling the model moulds cavity of various inner shapes inserted in rectangular cavity of the casting die (dimensions: 280 mm (height) x 190 mm (width) x 10 mm (depth) by applying model liquids of various density and viscosity are presented in the paper. Influence of die venting as well as inlet system area and inlet velocity on the volumetric rate of filling of the model liquid – achieved by means of filming the process in the system of a cold-chamber casting die was tested. Experiments compared with the results of simulation performed by means of the calculation module Novacast (Novaflow&Solid) for the selected various casting conditions – are also presented in the paper.
The paper presents a simulation model of the hybrid magnetic bearing dedicated to simulations of transient state. The proposed field-circuit model is composed of two components. The first part constitutes a set of ordinary differential equations that describes electrical circuits and mechanics. The second part of the simulation model consists of parameters such as magnetic forces, dynamic inductances and velocity-induced voltages obtained from the 3D finite element analysis. The MATLAB/Simulnik softwarewas used to implement the simulation model with the required control system. The proposed field-circuit model was validated by comparison of time responses with the prototype of the hybrid magnetic bearing.
The paper presents a conception of power electronics voltage controlled current source (VCCS) which is able much more precise mapping of its an output current in a reference signal, compared to a typical converter solution. It can be achieved by means of such interconnection of two separate converters that one of them corrects a total output current towards a reference signal. An output power of auxiliary converter is much smaller than an output power of main one. Thanks to continuous work of this converter also pulse modulation components in this current are minimized. These benefits are paid for by a relatively small increase in the complexity and the cost of the system. This conception of a converter has been called the double-converter topology (DCT). In the author opinion presented solution of the system can find application in many power electronics equipment and, therefore, will be developed. In the paper DCT basics, simulation experiments, and possible practical arrangement of the DCT are presented.
The results of bearing capacity, deformability and fracture toughness of reinforced concrete beams with the external reinforcement in the form of steel cut and stretchy sheet, obtained due to the conducting of the experiment and mathematical simulation which were made of concrete of C40/50 class are given in the article. Mathematical simulation of beam structures is done on the basis of the deformation model which allows to conduct calculations of the unified methodological positions of different elements with diverse configuration of cross section and reinforcement as well as take into consideration elastic and plastic properties of concrete and reinforcement, assessing the actual stress-strain state of sections of reinforced concrete elements at different loading levels, including ultimate one. The deformation model is based on the actual diagrams use of concrete and reinforcement materials deformation and conditions of efforts balance in the normal section and hypothesis of flat sections. The theoretical value of bearing capacity and deformability, obtained as a result of the mathematical simulation was compared to the experimental data. The satisfactory coincidence of the mathematical calculation of bearing capacity, deformability, fracture toughness and experimental data gives an opportunity to use the algorithm not only for beam structures with bar reinforcement but also for beam structures with the external reinforcement in the form of steel cut and stretchy sheet.
This paper presents a three-dimensional model of an airbag located outside of a small city car at the front bumper, which is intended to protect the vehicle against the effects of road traffic collisions. Results of numerical simulations of airbag operation in case of collision with two types of obstacles are presented: a flat, vertical wall and a circular pillar with a diameter of 200 mm. The paper presents the physical model, which is the subject of simulation, along with its mathematical description and the numerical calculation scheme used.
In this work we investigate the present capabilities of computational fluid dynamics for wall boiling. The computational model used combines the Euler/Euler two-phase flow description with heat flux partitioning. This kind of modeling was previously applied to boiling water under high pressure conditions relevant to nuclear power systems. Similar conditions in terms of the relevant non-dimensional numbers have been realized in the DEBORA tests using dichlorodifluoromethane (R12) as the working fluid. This facilitated measurements of radial profiles for gas volume fraction, gas velocity, bubble size and liquid temperature as well as axial profiles of wall temperature. After reviewing the theoretical and experimental basis of correlations used in the ANSYS CFX model used for the calculations, we give a careful assessment of the necessary recalibrations to describe the DEBORA tests. The basic CFX model is validated by a detailed comparison to the experimental data for two selected test cases. Simulations with a single set of calibrated parameters are found to give reasonable quantitative agreement with the data for several tests within a certain range of conditions and reproduce the observed tendencies correctly. Several model refinements are then presented each of which is designed to improve one of the remaining deviations between simulation and measurements. Specifically we consider a homogeneous MUSIG model for the bubble size, modified bubble forces, a wall function for turbulent boiling flow and a partial slip boundary condition for the liquid phase. Finally, needs for further model developments are identified and promising directions discussed.
The paper considers the modeling and estimation of the stochastic frontier model where the error components are assumed to be correlated and the inefficiency error is assumed to be autocorrelated. The multivariate Farlie-Gumble-Morgenstern (FGM) and normal copula are used to capture both the contemporaneous and the temporal dependence between, and among, the noise and the inefficiency components. The intractable multiple integrals that appear in the likelihood function of the model are evaluated using the Halton sequence based Monte Carlo (MC) simulation technique. The consistency and the asymptotic efficiency of the resulting simulated maximum likelihood (SML) estimators of the present model parameters are established. Finally, the application of model using the SML method to the real life US airline data shows significant noise-inefficiency dependence and temporal dependence of inefficiency.
The article is devoted to the development of technogenic risk management models and formalization of the process of support in making decision in the sphere of industrial safety. The structural, informative and mathematical models, used to process information in the technological risks management, as well as a formal model of the process of support of making decision in achieving an acceptable level of technical risk are presented and analyzed.
The paper deals with application of the Gumbel model to evaluation of the environmental loads. According to recommendations of Eurocodes, the conventional method of determining return period and characteristic values of loads utilizes the theory of extremes and implicitly assumes that the cumulative distribution function of the annual or other basic period extremes is the Gumbel distribution. However, the extreme value theory shows that the distribution of extremes asymptotically approaches the Gumbel distribution when the number of independent observations in each observation period from which the maximum is abstracted increases to infinity. Results of calculations based on simulation show that in practice the rate of convergence is very slow and significantly depends on the type of parent results distribution, values of coefficient of variation, and number of observation periods. In this connection, a straightforward purely empirical method based on fitting a curve to the observed extremes is suggested.