Numerical study on the surface depression of the concrete runway pavement under impact load

Airport runway pavement always subjected to huge impact loading due to the hard landing of aircraft on the pavement surface. Therefore runway pavements should have sufficient impact resistance capability to avoid damage causing by hard impact like surface deflection in downward or penetration since the repair works is cumbersome within the operating condition of airport and also increases the service life cost of the pavement structure. Several research works have been carried out on airport runway pavement to measure the present condition of pavement and also to predict future performance of it. However, most of the works are confined by pavement response under moving aircraft loading. Nevertheless, no comprehensive research work is yet conducted to identify the controlling factors which might have significant effect in changing the common pavements damage like surface penetration depth under impact of aircraft. Therefore, a 3D FE study is conducted to determine some effective factors in controlling the top surface penetration depth of runway pavement. Among the exterior factors, mass of the impactor, velocity of the impactor, impact angle and boundary conditions are selected and as interior factors, thickness of the runway pavement, compressive strength and density of materials used in the runway pavement are selected.

Computation of turbulent flow in an IFRF quarl burner

A computational model for isothermal axisymmetric turbulent flow in a quarl burner is set up using the CFD package FLUENT, and numerical solutions obtained from the model are compared with available experimental data. A standard k-e model and and two versions of the RNG k-e model are used to model the turbulence. One of the aims of the computational study is to investigate whether the RNG based k-e turbulence models are capable of yielding improved flow predictions compared with the standard k-e turbulence model. A difficulty is that the flow considered here features a confined vortex breakdown which can be highly sensitive to flow behaviour both upstream and downstream of the breakdown zone. Nevertheless, the relatively simple confining geometry allows us to undertake a systematic study so that both grid-independent and domain-independent results can be reported. The systematic study includes a detailed investigation of the effects of upstream and downstream conditions on the predictions, in addition to grid refinement and other tests to ensure that numerical error is not significant. Another important aim is to determine to what extent the turbulence model predictions can provide us with new insights into the physics of confined vortex breakdown flows. To this end, the computations are discussed in detail with reference to known vortex breakdown phenomena and existing theories. A major conclusion is that one of the RNG k-e models investigated here is able to correctly capture the complex forward flow region inside the recirculating breakdown zone. This apparently pathological result is in stark contrast to the findings of previous studies, most of which have concluded that either algebraic or differential Reynolds stress modelling is needed to correctly predict the observed flow features. Arguments are given as to why an isotropic eddy-viscosity turbulence model may well be able to capture the complex flow structure within the recirculating zone for this flow setup. With regard to the flow physics, a major finding is that the results obtained here are more consistent with the view that confined vortex breakdown is a type of axisymmetric boundary layer separation, rather than a manifestation of a subcritical flow state.

A numerical method for separation of stresses in photo-orthotropic elasticity

A numerical method is suggested for separation of stresses in photo-orthotropic elasticity using the numerical solution of compatibility equation for orthotropic case. The compatibility equation is written in terms of a stress parameter S analogous to the sum of principal stresses in two-dimensional isotropic case. The solution of this equation provides a relation between the normal stresses. The photoelastic data give the shear stress and another relation between the two normal stresses. The accuracy of the numerical method and its application to practical problems are illustrated with examples.

An absolute digital positioning numerical control system

A module containing all the functional components required for a digital absolute positioning process of one axis of a machine tool has been designed and constructed. Circuit realization makes use of integrated circuit elements.

Numerical study of heat transfer in laminar film boiling by the finite-difference method

The finite-difference form of the basic conservation equations in laminar film boiling have been solved by the false-transient method. By a judicious choice of the coordinate system the vapour-liquid interface is fitted to the grid system. Central differencing is used for diffusion terms, upwind differencing for convection terms, and explicit differencing for transient terms. Since an explicit method is used the time step used in the false-transient method is constrained by numerical instability. In the present problem the limits on the time step are imposed by conditions in the vapour region. On the other hand the rate of convergence of finite-difference equations is dependent on the conditions in the liquid region. The rate of convergence was accelerated by using the over-relaxation technique in the liquid region. The results obtained compare well with previous work and experimental data available in the literature.

Vibrations of Multifiber Composite Shells-Some Numerical Results

A study of vibrations of multifiber composite shells is presented. Special attention is paid to the effect of composition of different fibers on the frequency spectrum of a freely vibrating cylindrical shell. The numerical results indicate clustering of frequency spectrum of a freely vibrating cylindrical composite shell as compared with the isotropic shell, and the spectrum varies considerably with the composition of the constituent materials.

Experimental and numerical studies of hollow flange channel beams subject to web crippling under ETF and ITF load cases

This paper presents the details of experimental and numerical studies on the web crippling behaviour of hollow flange channel beams, known as LiteSteel beams (LSB). The LSB has a unique shape of a channel beam with two rectangular hollow flanges, made using a unique manufacturing process. Experimental and numerical studies have been carried out to evaluate the behaviour and design of LSBs subject to pure bending actions, predominant shear actions and combined actions. To date, however, no investigation has been conducted into the web crippling behaviour and strength of LSB sections under ETF and ITF load conditions. Hence experimental studies consisting of 28 tests were first conducted in this research to assess the web crippling behaviour and strengths of LSBs under two flange load cases (ETF and ITF). Experimental web crippling capacity results were then compared with the predictions from AS/NZS 4600 and AISI S100 design rules, which showed that AS/NZS 4600 and AISI S100 design equations are very unconservative for LSBs under ETF and ITF load cases. Hence improved equations were proposed to determine the web crippling capacities of LSBs. Finite element models of the tested LSBs were then developed, and used to determine the elastic buckling loads of LSBs under ETF and ITF load cases. New equations were proposed to determine the corresponding elastic buckling coefficients of LSBs. Finally suitable design rules were also developed under the Direct Strength Method format using the test results and buckling analysis results from finite element analyses.

Numerical modelling of rivet fastened rectangular hollow flange channel beams subject to shear

The intermittently rivet fastened Rectangular Hollow Flange Channel Beam (RHFCB) is a new cold-formed hollow section proposed as an alternative to welded hollow flange beams. Many experimental and numerical studies have been carried out in the past to investigate the shear behaviour of lipped channel beams. However, no research has been undertaken on the shear behaviour of rivet fastened RHFCBs. Therefore experimental and numerical studies were undertaken to investigate the shear behaviour and strength of rivet fastened RHFCBs. In this research finite element models of rivet fastened RHFCBs were developed to investigate their nonlinear shear behaviour including their buckling characteristics and ultimate shear strength. This paper presents the details of the finite element models of rivet fastened RHFCBs and the results. Both finite element analysis and experimental results showed that the current design rules are very conservative for the shear design of rivet fastened RHFCBs. Appropriate improvements have been proposed for the design rules of shear strength of rivet fastened RHFCBs within the Direct Strength Method format.

Thermal performance of magnesium oxide wall board using numerical modelling

Fire resistance of light-gauge steel frame (LSF) walls can be enhanced by lining them with single or multiple layers of wall boards. This research is focused on the thermal per-formance of Magnesium Oxide (MgO) wall boards in comparison to the conventional gypsum plasterboards exposed to standard fire on one side. Thermal properties of MgO board and gypsum plasterboard were measured first and then used in the finite element heat transfer models of the two types of panels. The measured thermal property results show that MgO board will perform better than the gypsum plasterboards due to its higher specific heat values at elevated temperatures. However, MgO board loses 50% of its ini-tial mass at about 500 °C compared to 16% for gypsum plasterboard. The developed finite element models were validated using the fire test results of gypsum plasterboards and then used to study the thermal performance of MgO board panels. Finite element analysis re-sults show that when MgO board panels are exposed to standard fire on one side the rate of temperature rise on the ambient side is significantly reduced compared to gypsum plas-terboard. This has the potential to improve the overall thermal performance of MgO board lined LSF walls and their fire resistance levels (FRL). However, full scale fire tests are needed to confirm this. This paper presents the details of this investigation and the results.

Material characterisation and numerical modelling of gypsum plasterboards in fire

The fire resistance characteristic of LSF wall systems mainly depends on the protective linings in use, commonly gypsum plasterboards. However, unclassified boards with varying composition and more notably with ambiguous thermal properties are increasingly becoming available in the market. Therefore a study was undertaken with an aim to set minimum standards for fire protective boards used in LSF wall applications. This paper presents the details of this study based on material characterisation and finite element thermal modelling of the most commonly used fire protective board, gypsum plasterboards, to address these critical issues related to fire safety design. In the material characterisation phase of this study, thermal properties of three different gypsum plasterboards manufactured in Australia were measured, analysed and compared. Subsequently, it proposes a thermal property based “k-factor” capable of giving an overall measure of the fire performance of boards, so that it can be used in appropriately classifying fire protective boards. As it is not known how this factor relates to the overall fire performance of LSF wall systems, numerical models were also developed and used to simulate the performance of LSF walls exposed to the standard fire. Finally, a correlation between time-temperature profiles from numerical analyses and calculated k-factors was established.

Effect of soil variability on the bearing capacity of clay and in slope stability problems

Site-specific geotechnical data are always random and variable in space. In the present study, a procedure for quantifying the variability in geotechnical characterization and design parameters is discussed using the site-specific cone tip resistance data (qc) obtained from static cone penetration test (SCPT). The parameters for the spatial variability modeling of geotechnical parameters i.e. (i) existing trend function in the in situ qc data; (ii) second moment statistics i.e. analysis of mean, variance, and auto-correlation structure of the soil strength and stiffness parameters; and (iii) inputs from the spatial correlation analysis, are utilized in the numerical modeling procedures using the finite difference numerical code FLAC 5.0. The influence of consideration of spatially variable soil parameters on the reliability-based geotechnical deign is studied for the two cases i.e. (a) bearing capacity analysis of a shallow foundation resting on a clayey soil, and (b) analysis of stability and deformation pattern of a cohesive-frictional soil slope. The study highlights the procedure for conducting a site-specific study using field test data such as SCPT in geotechnical analysis and demonstrates that a few additional computations involving soil variability provide a better insight into the role of variability in designs.

Numerical integration over arbitrary polygonal domains based on Schwarz-Christoffel conformal mapping

This paper presents a new numerical integration technique oil arbitrary polygonal domains. The polygonal domain is mapped conformally to the unit disk using Schwarz-Christoffel mapping and a midpoint quadrature rule defined oil this unit disk is used. This method eliminates the need for a two-level isoparametric mapping Usually required. Moreover, the positivity of the Jacobian is guaranteed. Numerical results presented for a few benchmark problems in the context of polygonal finite elements show that the proposed method yields accurate results.

Stationary crack tip fields in elastic-plastic solids: an overview of recent numerical simulations

In this paper, an overview of some recent numerical simulations of stationary crack tip fields in elastic-plastic solids is presented. First, asymptotic analyses carried out within the framework of 2D plane strain or plane stress conditions in both pressure insensitive and pressure sensitive plastic solids are reviewed. This is followed by discussion of salient results obtained from recent computational studies. These pertain to 3D characteristics of elastic-plastic near-front fields under mixed mode loading, mechanics of fracture and simulation of near-tip shear banding process of amorphous alloys and influence of crack tip constraint on the structure of near-tip fields in ductile single crystals. These results serve to illustrate several important features associated with stress and strain distributions near the crack tip and provide the foundation for understanding the operative failure mechanisms. The paper concludes by highlighting some of the future prospects for this field of study.

A semi-analytical particle filter for identification of nonlinear oscillators

Particle filters find important applications in the problems of state and parameter estimations of dynamical systems of engineering interest. Since a typical filtering algorithm involves Monte Carlo simulations of the process equations, sample variance of the estimator is inversely proportional to the number of particles. The sample variance may be reduced if one uses a Rao-Blackwell marginalization of states and performs analytical computations as much as possible. In this work, we propose a semi-analytical particle filter, requiring no Rao-Blackwell marginalization, for state and parameter estimations of nonlinear dynamical systems with additively Gaussian process/observation noises. Through local linearizations of the nonlinear drift fields in the process/observation equations via explicit Ito-Taylor expansions, the given nonlinear system is transformed into an ensemble of locally linearized systems. Using the most recent observation, conditionally Gaussian posterior density functions of the linearized systems are analytically obtained through the Kalman filter. This information is further exploited within the particle filter algorithm for obtaining samples from the optimal posterior density of the states. The potential of the method in state/parameter estimations is demonstrated through numerical illustrations for a few nonlinear oscillators. The proposed filter is found to yield estimates with reduced sample variance and improved accuracy vis-a-vis results from a form of sequential importance sampling filter.

Numerical Study of Heat Transfer From Pin-Fin Heat Sink Using Steady and Pulsated Impinging Jets

This paper investigates numerically the heat transfer characteristics of confined slot jet impingement on a pin-fin heat sink. A variety of pin-fin heat sinks is investigated, and the resulting enhancement of heat transfer studied. The distribution of heat transfer coefficient on the top surface of the base plate and that along the fin height are examined. Both steady and pulsated jets are studied. It is observed that for a steady jet impingement on a pin-fin heat sink, the effective heat transfer coefficient increases with fin height, leading to a corresponding decrease in base plate temperature for the same heat flux. In the case of pulsated jets, the influence of pulse frequency and the Reynolds number is examined, and their effect on the effective heat transfer coefficient is studied.