Numerical simulation of a Savonius turbine above an infinite-width forward facing step

The Savonius turbine, although simple in construction, typically has a maximum power coefficient (cP) of about 0.2. This is significantly lower than the cP of the axial flow propeller-type turbine which typically can be as high as 0.5. However, a simple means to improve the cP of a Savonius turbine is to install it above a forward facing step, for example, a cliff or a building. In this work, prior experimental results of the tow testing of a Savonius turbine installed above a finite-width bluff body were used to validate computational fluid dynamics simulation of the same experimental conditions. The validated simulation settings were then used to obtain the maximum cP of a similar turbine of finite width but installed above an infinite-width forward facing step over a range of installation positions above and behind the step.

Simulation of the fastener manufacturing process

Cold forging of steel has been modelled using a physical simulation and the finite element method (FEM) to predict the conditions under which failure occurs. FEM was performed using FORGE2 software, which allows analysis of two-dimensional, axi-symmetric forging processes. Various cold forming processes were modelled. The material model has been validated through wire drawing and bulging experiments. The results will be applied to cold forging of industrial products such as fasteners.

Numerical modeling of the dynamic tensile behavior of irregular fibers

Most fibers are irregular, and they are often subjected to rapid straining during mechanical processing and end-use applications. In this paper, the effect of fiber dimensional irregularities on the dynamic tensile behavior of irregular fibers is examined, using the finite element method (FEM). Fiber dimensional irregularities are simulated with sine waves of different magnitude (10%, 30% and 50% level of diameter variation). The tensile behavior of irregular fibers is examined at different strain rates (333%/sec, 3,333%/sec and 30,000%/sec). The breaking load and breaking extension of irregular fibers at different strain rates are then calculated from the finite element model. The results indicate that strain rate has a significant effect on the dynamic tensile behavior of an irregular fiber, and that the position of the thinnest segment along the fiber affects the simulation results markedly. Under dynamic conditions, an irregular fiber does not necessarily break at the thinnest segment, which is different from the quasi-static results.

Numerical study of the mode I delamination toughness of Z-pinned laminates

A finite element (FE) model is developed to investigate mode I delamination toughness of z-pin reinforced composite laminates. The z-pin pullout process is simulated by the deformation of a set of non-linear springs. A critical crack opening displacement (COD) criterion is used to simulate crack growth in a double-cantilever-beam (DCB) made of z-pinned laminates. The toughness of the structure is quantified by the energy release rate, which is calculated using the contour integral method. The FE model is verified for both unpinned and z-pinned laminates. Predicted loading forces from FE analysis are compared to available test data. Good agreement is achieved. Our numerical results indicate that z-pins can greatly increase the mode I delamination toughness of the composite laminates. The influence of design parameters on the toughness enhancement of z-pinned laminates is also investigated, which provides important information to optimise and improve the z-pinning technique.

Crushing simulation of foam-filled aluminium tubes

A numerical study is presented in this paper to investigate the energy absorption of foam-filled aluminium tubes during crushing. The post-buckling mode of the foam-tube structures has been successfully simulated. The predicted compressive load-displacement is in a good agreement with experimental results. The energy absorption ability of the composite structure due to plastic deformation in a crushing process is evaluated by comparison with the tube structure without foam. The results indicate that the energy absorption of a foam-filled tube structure is superior to the tube without foam. The influences of the friction and the geometric parameters of the structure on the energy absorption have also been investigated. Results from this study will assist automotive industry to design crashworthy components based on foam-filled tubes.

Dynamic simulation of collisions of heavy high-speed trucks with concrete barriers

Real vehicle collision experiments on full-scale road safety barriers are important to determine the outcome of a vehicle versus barrier impact accident. However, such experiments require large investment of time and money. Numerical simulation has therefore been imperative as an alternative method for testing concrete barriers. In this research, spring subgrade models were first developed to simulate the ground boundary of concrete barriers. Both heavy trucks and concrete barriers were modeled using finite element methods (FEM) to simulate dynamic collision performances. Comparison of the results generated from computer simulations and on-site full-scale experiments demonstrated that the developed models could be applied to simulate the collision of heavy trucks with concrete barriers to provide the data to design new road safety barriers and analyze existing ones.

Accuracy and precision assesment of stochastic simulation tools for springback variation

The sheet forming industry is plagued by inherent variations in its many input variables, making quality control and improvements a major hurdle. This is particularly poignant for Advanced High Strength Steels (AHSS), which exhibit a large degree of property variability. Current FE-based simulation packages are successful at predicting the manufacturability of a particular sheet metal components, however, due to their numerical deterministic nature are inherently unable to predict the performance of a real-life production process. Though they are now beginning to incorporate the stochastic nature of production in their codes. This work investigates the accuracy and precision of a current stochastic simulation package, AutoForm Sigma v4.1, by developing an experimental data set where all main sources of variation are captured through precise measurements and standard tensile tests. Using a Dual Phase 600Mpa grade steel a series of semi-cylindrical channels are formed at two Blank Holder Pressure levels where the response metric is the variation in springback determined by the flange angle. The process is replicated in AutoForm Sigma and an assessment of accuracy and precision of the predictions are performed. Results indicate a very good correspondence to the experimental trials, with mean springback response predicted to within 1 ° of the flange angle and the interquartile spread of results to within 0.22°.

Computational aspects of the numerical solution of SDEs

In the last 30 to 40 years, many researchers have combined to build the knowledge base of theory and solution techniques that can be applied to the case of differential equations which include the effects of noise. This class of ``noisy'' differential equations is now known as stochastic differential equations (SDEs). Markov diffusion processes are included within the field of SDEs through the drift and diffusion components of the Itô form of an SDE. When these drift and diffusion components are moderately smooth functions, then the processes' transition probability densities satisfy the Fokker-Planck-Kolmogorov (FPK) equation -- an ordinary partial differential equation (PDE). Thus there is a mathematical inter-relationship that allows solutions of SDEs to be determined from the solution of a noise free differential equation which has been extensively studied since the 1920s. The main numerical solution technique employed to solve the FPK equation is the classical Finite Element Method (FEM). The FEM is of particular importance to engineers when used to solve FPK systems that describe noisy oscillators. The FEM is a powerful tool but is limited in that it is cumbersome when applied to multidimensional systems and can lead to large and complex matrix systems with their inherent solution and storage problems. I show in this thesis that the stochastic Taylor series (TS) based time discretisation approach to the solution of SDEs is an efficient and accurate technique that provides transition and steady state solutions to the associated FPK equation. The TS approach to the solution of SDEs has certain advantages over the classical techniques. These advantages include their ability to effectively tackle stiff systems, their simplicity of derivation and their ease of implementation and re-use. Unlike the FEM approach, which is difficult to apply in even only two dimensions, the simplicity of the TS approach is independant of the dimension of the system under investigation. Their main disadvantage, that of requiring a large number of simulations and the associated CPU requirements, is countered by their underlying structure which makes them perfectly suited for use on the now prevalent parallel or distributed processing systems. In summary, l will compare the TS solution of SDEs to the solution of the associated FPK equations using the classical FEM technique. One, two and three dimensional FPK systems that describe noisy oscillators have been chosen for the analysis. As higher dimensional FPK systems are rarely mentioned in the literature, the TS approach will be extended to essentially infinite dimensional systems through the solution of stochastic PDEs. In making these comparisons, the advantages of modern computing tools such as computer algebra systems and simulation software, when used as an adjunct to the solution of SDEs or their associated FPK equations, are demonstrated.

Numerical analysis of heat transfer and the optimisation of regenerators

To optimize a regenerator’s structure and its operation parameters and, consequently, to increase the efficiency of heat recovery and to save energy, a computational approach is used to study the unsteady three-dimensional flow and heat transfer. The simulation is performed in two steps. In the first step, the gas flow and heat transfer in a typical sphere-bed unit is simulated to deduce a dimensionless equation of heat transfer between gas and sphere. In the second step, a model is developed to simplify the prototype and to simulate the gas flow and heat transfer in the whole regenerator. The heat exchange process in regenerators and the effects of the regenerator’s structure and operation parameters, such as gas mass flux, reversal time, regenerator height, sphere diameter, and thermophysical properties of the spheres, are studied with the model to determine efficiency of heat recovery.

Transient fluid–structure coupling for simulation of a trileaflet heart valve using weak coupling

In this article, a three-dimensional transient numerical approach coupled with fluid–structure interaction for the modeling of an aortic trileaflet heart valve at the initial opening stage is presented. An arbitrary Lagrangian–Eulerian kinematical description together with an appropriate fluid grid was used for the coupling strategy with the structural domain. The fluid dynamics and the structure aspects of the problem were analyzed for various Reynolds numbers and times. The fluid flow predictions indicated that at the initial leaflet opening stage a circulation zone was formed immediately downstream of the leaflet tip and propagated outward as time increased. Moreover, the maximum wall shear stress in the vertical direction of the leaflet was found to be located near the bottom of the leaflet, and its value decreased sharply toward the tip. In the horizontal cross section of the leaflet, the maximum wall shear stresses were found to be located near the sides of the leaflet.

Potential applications of discrete-event simulation and fuzzy rule-based systems to structural reliability and availability

This chapter discusses and illustrates some potential applications of discrete-event simulation (DES) techniques in structural reliability and availability analysis, emphasizing the convenience of using probabilistic approaches in modern building and civil engineering practices. After reviewing existing literature on the topic, some advantages of probabilistic techniques over analytical ones are highlighted. Then, we introduce a general framework for performing structural reliability and availability analysis through DES. Our methodology proposes the use of statistical distributions and techniques – such as survival analysis – to model component-level reliability. Then, using failure- and repair-time distributions and information about the structural logical topology (which allows determination of the structural state from their components’ state), structural reliability, and availability information can be inferred. Two numerical examples illustrate some potential applications of the proposed methodology to achieving more reliable and structural designs. Finally, an alternative approach to model uncertainty at component level is also introduced as ongoing work. This new approach is based on the use of fuzzy rule-based systems and it allows the introduction of experts’ opinions and evaluations in our methodology.

Numerical simulations on warm forming of stainless steel with trip-effect

In steels with TRIP-effect, a phase transformation from the retained-austenite to martensite occurs during forming, and it significantly affects hardening behaviours. Such an effect is sensitive to the amount of strain as well as the temperature variation. For materials with a strong TRIP-effect, new forming techniques are needed to develop that can lead to lighter and stronger components in automotive industry. This paper presents a coupled thermo-mechanical finite element modelling and simulation of a warm deep drawing of austenitic stainless steel (including a TRIP-effect) using LS-DYNA and temperature effect on forming process of such materials is investigated.

A full body musculoskeletal model based on flexible multibody simulation approach utilised in bone strain analysis during human locomotion

Load-induced strains applied to bone can stimulate its development and adaptation. In order to quantify the incident strains within the skeleton, in vivo implementation of strain gauges on the surfaces of bone is typically used. However, in vivo strain measurements require invasive methodology that is challenging and limited to certain regions of superficial bones only such as the anterior surface of the tibia. Based on our previous study [Al Nazer et al. (2008) J Biomech. 41:1036–1043], an alternative numerical approach to analyse in vivo strains based on the flexible multibody simulation approach was proposed. The purpose of this study was to extend the idea of using the flexible multibody approach in the analysis of bone strains during physical activity through integrating the magnetic resonance imaging (MRI) technique within the framework. In order to investigate the reliability and validity of the proposed approach, a three-dimensional full body musculoskeletal model with a flexible tibia was used as a demonstration example. The model was used in a forward dynamics simulation in order to predict the tibial strains during walking on a level exercise. The flexible tibial model was developed using the actual geometry of human tibia, which was obtained from three-dimensional reconstruction of MRI. Motion capture data obtained from walking at constant velocity were used to drive the model during the inverse dynamics simulation in order to teach the muscles to reproduce the motion in the forward dynamics simulation. Based on the agreement between the literature-based in vivo strain measurements and the simulated strain results, it can be concluded that the flexible multibody approach enables reasonable predictions of bone strain in response to dynamic loading. The information obtained from the present approach can be useful in clinical applications including devising exercises to prevent bone fragility or to accelerate fracture healing.

Integrated fluid-thermal-structural numerical analysis for the quenching of metallic components

The quenching of a metal component with a channel section in a water tank is numerically simulated. Computational fluid dynamics (CFD) is used to model the multiphase flow and the heat transfer in film boiling, nucleate boiling and convective cooling processes to calculate the difference in heat transfer rate around the component and then combining with the thermal simulation and structure analysis of the component to study the effect of heat transfer rate on the distortion of the U-channel component. A model is also established to calculate the residual stress produced by quenching. The coupling fluid-thermal-structural simulation provides an insight into the deformation of the component and can be used to perform parameter analysis to reduce the distortion of the component. © 2011 Shanghai Jiaotong University and Springer-Verlag Berlin Heidelberg.