Numerical analyses of crack spacing due to shrinkage in reinforced concrete slabs

Shrinkage cracking is commonly observed in concrete flat structures such as highway pavements, slabs, and bridge decks. Crack spacing due to shrinkage has received considerable attention for many years [1-3]. However, some aspects concerning the mechanism of crack spacing still remain un-clear. Though it is well known that the interval of the cracks generally falls with a range, no satisfactory explanation has been put forward as to why the minimum spacing exists.

Numerical Analyses of Discrete Gust Response for an Aircraft

Based on Navier-Stokes equations and structural and flight dynamic equations of motion, dynamic responses in vertical discrete gust flow perturbation are investigated for a supersonic transport model. A tightly coupled method was developed by subiterations between aerodynamic equations and dynamic equations of motion. First, under the assumption of rigid-body and single freedom of motion in the vertical plunging, the results of a direct-coupling method are compared with the results of quasi-steady model method. Then, gust responses for the one-minus-cosine gust profile arc analyzed with two freedoms of motion in plunging and pitching for the airplane configurations with and without the consideration of structural deformation.

Numerical analyses for the material bifurcation in plane sheet under tension

The various patterns (shear banding, surface wrinkling and necking) of material bifurcation in plane sheet under tension are investigated in this paper by means of a numerical method. It is found that numerical analysis can provide better ground for searching for the lowest critical loads. The inhomogeneity caused by void damage and the nonuniformity in the stress distribution across sheet thickness are proved to have detrimental effects on the material bifurcation. Nevertheless, material stability can be promoted by any means of depressing void damage or alleviating stress, even locally across the thickness. Besides, the peculiar behaviour of material bifurcation under slight biaxiality state is demonstrated. Copyright (C) 1996 Elsevier Science Ltd

Numerical Analyses of Transonic Fluid/Structure interaction in Aircraft Engineering

Through the coupling between aerodynamic and structural governing equations, a fully implicit multiblock aeroelastic solver was developed for transonic fluid/stricture interaction. The Navier-Stokes fluid equations are solved based on LU-SGS (lower-upper symmetric Gauss-Seidel) Time-marching subiteration scheme and HLLEW (Harten-Lax-van Leer-Einfeldt-Wada) spacing discretization scheme and the same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Transfinite interpolation (TFI) is used for the grid deformation of blocks neighboring the flexible surfaces. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between fluid and structure. The developed code was fort validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. In the subsonic and transonic range, the calculated flutter speeds and frequencies agree well with experimental data, however, in the supersonic range, the present calculation overpredicts the experimental flutter points similar to other computations. Then the flutter character of a complete aircraft configuration is analyzed through the calculation of the change of structural stiffness. Finally, the phenomenon of aileron buzz is simulated for the weakened model of a supersonic transport wing/body model at Mach numbers of 0.98 and l.05. The calculated unsteady flow shows, on the upper surface, the shock wave becomes stronger as the aileron deflects downward, and the flow behaves just contrary on the lower surface of the wing. Corresponding to general theoretical analysis, the flow instability referred to as aileron buzz is induced by a stronger shock alternately moving on the upper and lower surfaces of wing. For the rigid structural model, the flow is stable at all calculated Mach numbers as observed in experiment

Numerical analyses of tunnel-induced settlement damage to a masonry wall

Settlements due to underground construction represent a risk for the architectural heritage, especially in The Netherlands, because of the combination of soft soil, fragile pile foundation and brittle, un-reinforced masonry façade. Modelling of soil-structure interaction is fundamental to assess the risk of building damage due to tunnelling. This paper presents results of finite element analyses carried out with different models for a simple masonry wall. Focus is paid on the comparison between coupled, uncoupled and semi-coupled analyses, in which the soil-structure interaction is represented in different ways. In particular, the implementation of a soil-structure interface model in the numerical analyses is analysed, in order to asses its validity. The aim of the research project is the development of a damage classification system for different building typologies.

Numerical analyses of Runge–Kutta implicit–explicit schemes for horizontally explicit, vertically implicit solutions of atmospheric models

With the prospect of exascale computing, computational methods requiring only local data become especially attractive. Consequently, the typical domain decomposition of atmospheric models means horizontally-explicit vertically-implicit (HEVI) time-stepping schemes warrant further attention. In this analysis, Runge-Kutta implicit-explicit schemes from the literature are analysed for their stability and accuracy using a von Neumann stability analysis of two linear systems. Attention is paid to the numerical phase to indicate the behaviour of phase and group velocities. Where the analysis is tractable, analytically derived expressions are considered. For more complicated cases, amplification factors have been numerically generated and the associated amplitudes and phase diagnosed. Analysis of a system describing acoustic waves has necessitated attributing the three resultant eigenvalues to the three physical modes of the system. To do so, a series of algorithms has been devised to track the eigenvalues across the frequency space. The result enables analysis of whether the schemes exactly preserve the non-divergent mode; and whether there is evidence of spurious reversal in the direction of group velocities or asymmetry in the damping for the pair of acoustic modes. Frequency ranges that span next-generation high-resolution weather models to coarse-resolution climate models are considered; and a comparison is made of errors accumulated from multiple stability-constrained shorter time-steps from the HEVI scheme with a single integration from a fully implicit scheme over the same time interval. Two schemes, “Trap2(2,3,2)” and “UJ3(1,3,2)”, both already used in atmospheric models, are identified as offering consistently good stability and representation of phase across all the analyses. Furthermore, according to a simple measure of computational cost, “Trap2(2,3,2)” is the least expensive.

Experimental and numerical analyses about the efficiency of flow through devices for the sediment controll in urban runoff

As land is developed, the impervious surfaces that are created increase the amount of runoff during rainfall events, disrupting the natural hydrologic cycle, with an increment in volume of runoff and in pollutant loadings. Pollutants deposited or derived from an activity on the land surface will likely end up in stormwater runoff in some concentration, such as nutrients, sediment, heavy metals, hydrocarbons, gasoline additives, pathogens, deicers, herbicides and pesticides. Several of these pollutants are particulate-bound, so it appears clear that sediment removal can provide significant water-quality improvements and it appears to be important the knowledge of the ability of stromwater treatment devices to retain particulate matter. For this reason three different units which remove sediments have been tested through laboratory. In particular a roadside gully pot has been tested under steady hydraulic conditions, varying the characteristics of the influent solids (diameter, particle size distribution and specific gravity). The efficiency in terms of particles retained has been evaluated as a function of influent flow rate and particles characteristics; results have been compared to efficiency evaluated applying an overflow rate model. Furthermore the role of particles settling velocity in efficiency determination has been investigated. After the experimental runs on the gully pot, a standard full-scale model of an hydrodynamic separator (HS) has been tested under unsteady influent flow rate condition, and constant solid concentration at the input. The results presented in this study illustrate that particle separation efficiency of the unit is predominately influenced by operating flow rate, which strongly affects the particles and hydraulic residence time of the system. The efficiency data have been compared to results obtained from a modified overflow rate model; moreover the residence time distribution has been experimentally determined through tracer analyses for several steady flow rates. Finally three testing experiments have been performed for two different configurations of a full-scale model of a clarifier (linear and crenulated) under unsteady influent flow rate condition, and constant solid concentration at the input. The results illustrate that particle separation efficiency of the unit is predominately influenced by the configuration of the unit itself. Turbidity measures have been used to compare turbidity with the suspended sediments concentration, in order to find a correlation between these two values, which can allow to have a measure of the sediments concentration simply installing a turbidity probe.

Study of wheel-rail impact at insulated rail joint through experimental and numerical methods

As part of an ongoing research on the development of a longer life insulated rail joint (IRJ), this paper reports a field experiment and a simplified 2D numerical modelling for the purpose of investigating the behaviour of rail web in the vicinity of endpost in an insulated rail joint (IRJ) due to wheel passages. A simplified 2D plane stress finite element model is used to simulate the wheel-rail rolling contact impact at IRJ. This model is validated using data from a strain gauged IRJ that was installed in a heavy haul network; data in terms of the vertical and shear strains at specific positions of the IRJ during train passing were captured and compared with the results of the FE model. The comparison indicates a satisfactory agreement between the FE model and the field testing. Furthermore, it demonstrates that the experimental and numerical analyses reported in this paper provide a valuable datum for developing further insight into the behaviour of IRJ under wheel impacts.

Numerical modelling and design of LiteSteel Beams subject to lateral buckling

The LiteSteel Beam (LSB) is a new hollow flange channel section developed using a patented dual electric resistance welding and cold-forming process. It has a unique geometry consisting of torsionally rigid rectangular hollow flanges and a slender web, and is commonly used as flexural members. However, the LSB flexural members are subjected to a relatively new lateral distortional buckling mode, which reduces their moment capacities. Unlike lateral torsional buckling, the lateral distortional buckling of LSBs is characterised by simultaneous lateral deflection, twist and cross sectional change due to web distortion. Therefore a detailed investigation into the lateral buckling behaviour of LSB flexural members was undertaken using experiments and finite element analyses. This paper presents the details of suitable finite element models developed to simulate the behaviour and capacity of LSB flexural members subject to lateral buckling. The models included all significant effects that influence the ultimate moment capacities of such members, including material inelasticity, lateral distortional buckling deformations, web distortion, residual stresses, and geometric imperfections. Comparison of elastic buckling and ultimate moment capacity results with predictions from other numerical analyses and available buckling moment equations, and experimental results showed that the developed finite element models accurately predict the behaviour and moment capacities of LSBs. The validated model was then used in a detailed parametric study that produced accurate moment capacity data for all the LSB sections and improved design rules for LSB flexural members subject to lateral distortional buckling.

Numerical and experimental studies of cold-formed steel floor systems under standard fire conditions

Light gauge cold-formed steel frame (LSF) structures are increasingly used in industrial, commercial and residential buildings because of their non-combustibility, dimensional stability, and ease of installation. A floor-ceiling system is an example of its applications. LSF floor-ceiling systems must be designed to serve as fire compartment boundaries and provide adequate fire resistance. Fire rated floor-ceiling assemblies formed with new materials and construction methodologies have been increasingly used in buildings. However, limited research has been undertaken in the past and hence a thorough understanding of their fire resistance behaviour is not available. Recently a new composite panel in which an external insulation layer is used between two plasterboards has been developed at QUT to provide a higher fire rating to LSF floors under standard fire conditions. But its increased fire rating could not be determined using the currently available design methods. Research on LSF floor systems under fire conditions is relatively recent and the behaviour of floor joists and other components in the systems is not fully understood. The present design methods thus require the use of expensive fire protection materials to protect them from excessive heat increase during a fire. This leads to uneconomical and conservative designs. Fire rating of these floor systems is provided simply by adding more plasterboard sheets to the steel joists and such an approach is totally inefficient. Hence a detailed fire research study was undertaken into the structural and thermal performance of LSF floor systems including those protected by the new composite panel system using full scale fire tests and extensive numerical studies. Experimental study included both the conventional and the new steel floor-ceiling systems under structural and fire loads using a gas furnace designed to deliver heat in accordance with the standard time- temperature curve in AS 1530.4 (SA, 2005). Fire tests included the behavioural and deflection characteristics of LSF floor joists until failure as well as related time-temperature measurements across the section and along the length of all the specimens. Full scale fire tests have shown that the structural and thermal performance of externally insulated LSF floor system was superior than traditional LSF floors with or without cavity insulation. Therefore this research recommends the use of the new composite panel system for cold-formed LSF floor-ceiling systems. The numerical analyses of LSF floor joists were undertaken using the finite element program ABAQUS based on the measured time-temperature profiles obtained from fire tests under both steady state and transient state conditions. Mechanical properties at elevated temperatures were considered based on the equations proposed by Dolamune Kankanamge and Mahendran (2011). Finite element models were calibrated using the full scale test results and used to further provide a detailed understanding of the structural fire behaviour of the LSF floor-ceiling systems. The models also confirmed the superior performance of the new composite panel system. The validated model was then used in a detailed parametric study. Fire tests and the numerical studies showed that plasterboards provided sufficient lateral restraint to LSF floor joists until their failure. Hence only the section moment capacity of LSF floor joists subjected to local buckling effects was considered in this research. To predict the section moment capacity at elevated temperatures, the effective section modulus of joists at ambient temperature is generally considered adequate. However, this research has shown that it leads to considerable over- estimation of the local buckling capacity of joist subject to non-uniform temperature distributions under fire conditions. Therefore new simplified fire design rules were proposed for LSF floor joist to determine the section moment capacity at elevated temperature based on AS/NZS 4600 (SA, 2005), NAS (AISI, 2007) and Eurocode 3 Part 1.3 (ECS, 2006). The accuracy of the proposed fire design rules was verified with finite element analysis results. A spread sheet based design tool was also developed based on these design rules to predict the failure load ratio versus time, moment capacity versus time and temperature for various LSF floor configurations. Idealised time-temperature profiles of LSF floor joists were developed based on fire test measurements. They were used in the detailed parametric study to fully understand the structural and fire behaviour of LSF floor panels. Simple design rules were also proposed to predict both critical average joist temperatures and failure times (fire rating) of LSF floor systems with various floor configurations and structural parameters under any given load ratio. Findings from this research have led to a comprehensive understanding of the structural and fire behaviour of LSF floor systems including those protected by the new composite panel, and simple design methods. These design rules were proposed within the guidelines of the Australian/New Zealand, American and European cold- formed steel structures standard codes of practice. These may also lead to further improvements to fire resistance through suitable modifications to the current composite panel system.

Experimental and numerical investigation of the behaviour of CFRP strengthened CHS beams subjected to bending

This paper presents the results of an experimental and numerical program to investigate the circular hollow section (CHS) beams, strengthened using Carbon Fibre Reinforced Polymer (CFRP) sheets. The circular hollow shaped steel beams bonded with different CFRP layer orientations were tested under four-point bending. The mid-span deflection, service load and failure load were recorded. The LHL (where L, first inner longitudinal layer, H, second hoop layer and L, third outer longitudinal layer) and LLH (where L, first inner longitudinal layer, L, second longitudinal layer and H, third outer hoop layer) layer oriented strengthened beams perform slightly better than HHL (where H, first inner hoop layer, H, second hoop layer and L, third outer longitudinal layer) layer oriented strengthened beams. The LHL and LLH layer oriented treated beams showed very similar structural behaviour. Numerical analyses were then conducted on the CFRP strengthened steel CHS beams. The validity of the models has been assessed by comparing the failure loads and mid-span deflections. The effects of various parameters such as bond length, section types, tensile modulus of CFRP, adhesive layer thickness and adhesive types have been studied.

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.

Numerical Analysis of LEC Growth of GaAs with an Axial Magnetic Field

A set of numerical analyses for momentum and heat transfer For a 3 in. (0.075 m) diameter Liquid Encapsulant Czochralski (LEC) growth of single-crystal GaAs with or without all axial magnetic field was carried Out using the finite-element method. The analyses assume a pseudosteady axisymmetric state with laminar floats. Convective and conductive heat transfers. radiative heat transfer between diffuse surfaces and the Navier-Stokes equations for both melt and encapsulant and electric current stream function equations Cor melt and crystal Lire considered together and solved simultaneously. The effect of the thickness of encapsulant. the imposed magnetic field strength as well as the rotation rate of crystal and crucible on the flow and heat transfer were investigated. (C) 2002 Published by Elsevier Science Ltd.

Numerical analysis of LEC growth of GaAs with an axial magnetic field

A set of numerical analyses for momentum and heat transfer For a 3 in. (0.075 m) diameter Liquid Encapsulant Czochralski (LEC) growth of single-crystal GaAs with or without all axial magnetic field was carried Out using the finite-element method. The analyses assume a pseudosteady axisymmetric state with laminar floats. Convective and conductive heat transfers. radiative heat transfer between diffuse surfaces and the Navier-Stokes equations for both melt and encapsulant and electric current stream function equations Cor melt and crystal Lire considered together and solved simultaneously. The effect of the thickness of encapsulant. the imposed magnetic field strength as well as the rotation rate of crystal and crucible on the flow and heat transfer were investigated. (C) 2002 Published by Elsevier Science Ltd.

Numerical Simulation of Waves Generated by Seafloor Movements

Waves generated by vertical seafloor movements are simulated by use of a fully nonlinear two-dimensional numerical wave tank. In the source region, the seafloor lifts to a designated height by a generation function. The numerical tests show that file linear theory is only valid for estimating the wave behaviors induced by the seafloor movements with a small amplitude, and the fully nonlinear numerical model should be adopted in the simulation of the wave generation by the large amplitude seafloor movements. Without the background surface waves, many numerical tests on the stable maximum elevations eta(max)(0) are carried out by both the linear theory and the fully nonlinear model. The results of two models are compared and analyzed. For the fully nonlinear model, the influences of the amplitudes and the horizontal lengths on eta(max)(0) are stronger than that of the characteristic duration times. Furthermore, results reveal that there are significant differences between the linear theory and the fully nonlinear model. When the influences of the background surface waves are considered, the corresponding numerical analyses reveal that with the fully nonlinear model the eta(max)(0) near-linearly varies with the wave amplitudes of the surface waves, and the eta(max)(0) has significant dependences on the wave lengths and the wave phases of the surface waves. In addition, the differences between the linear theory and the fully nonlinear model are still obvious, aid these differences are significantly affected by The wave parameters of the background surface waves, such as the wave amplitude, the wave length and the wave phase.