Nonlinear numerical modelling of lightning strike effect on composite panels with temperature dependent material properties

This paper presents a physics based modelling procedure to predict the thermal damage of composite material when struck by lightning. The procedure uses the Finite Element Method with non-linear material models to represent the extreme thermal material behaviour of the composite material (carbon/epoxy) and an embedded copper mesh protection system. Simulation predictions are compared against published experimental data, illustrating the potential accuracy and computational cost of virtual lightning strike tests and the requirement for temperature dependent material modelling. The modelling procedure is then used to examine and explain a number of practical solutions to minimize thermal material damage. © 2013 Elsevier Ltd.

Modelling and characterization of ultrasonic consolidation process of Aluminium alloys

Ultrasonic consolidation process is a rapid manufacturing process used to join thin layers of metal at low temperatures and low energy consumption. In this work, finite element method has been used to simulate the ultrasonic consolidation of Aluminium alloys 6061 (AA-6061) and 3003 (AA-3003). A thermomechanical material model has been developed in the framework of continuum cyclic plasticity theory which takes into account both volume (acoustic softening) and surface (thermal softening due to friction) effects. A friction model based on experimental studies has been developed, which takes into account the dependence of coefficient of friction upon contact pressure, amount of slip, temperature and number of cycles. Using the developed material and friction model ultrasonic consolidation (UC) process has been simulated for various combinations of process parameters involved. Experimental observations are explained on the basis of the results obtained in the present study. The current research provides the opportunity to explain the differences of the behaviour of AA-6061 and AA-3003 during the ultrasonic consolidation process. Finally, trends of the experimentally measured fracture energies of the bonded specimen are compared to the predicted friction work at the weld interface resulted from the simulation at similar process condition. Similarity of the trends indicates the validity of the developed model in its predictive capability of the process. © 2008 Materials Research Society.

Implementation of a non-orthogonal model for the finite element simulation of textile composite draping

In the pursuit of producing high quality, low-cost composite aircraft structures, out-of-autoclave manufacturing processes for textile reinforcements are being simulated with increasing accuracy. This paper focuses on the continuum-based, finite element modelling of textile composites as they deform during the draping process. A non-orthogonal constitutive model tracks yarn orientations within a material subroutine developed for Abaqus/Explicit, resulting in the realistic determination of fabric shearing and material draw-in. Supplementary material characterisation was experimentally performed in order to define the tensile and non-linear shear behaviour accurately. The validity of the finite element model has been studied through comparison with similar research in the field and the experimental lay-up of carbon fibre textile reinforcement over a tool with double curvature geometry, showing good agreement.

Numerical Simulation of Tanpura String Vibrations

A numerical model of a tanpura string is presented, based on a recently developed, stability-preserving way of incorporating the non-smooth forces involved in the impactive distributed contact between the string and the bridge. By defining and modelling the string-bridge contact over the full length of the bridge, the simulated vibrations can be monitored through the force signals at both the bridge and the nut. As such it offers a reference model for both measurements and sound synthesis. Simulations starting from different types of initial conditions demonstrate that the model reproduces the main characteristic feature of the tanpura, namely the sustained appearance of a precursor in the force waveforms, carrying a band of overtones which decrease in frequency as the string vibrations decay. Results obtained with the numerical model are used to examine, through comparison, the effect of the bridge and of the thread on the vibrations.

Defining Simulation Intent

Defining Simulation Intent involves capturing high level modelling and idealisation decisions in order to create an efficient and fit-for-purpose analysis. These decisions are recorded as attributes of the decomposed design space.

An approach to defining Simulation Intent is described utilising three known technologies: Cellular Modelling, the subdivision of space into volumes of simulation significance (structures, gas paths, internal and external airflows etc.); Equivalencing, maintaining a consistent and coherent description
of the equivalent representations of the spatial cells in different analysis models; and Virtual Topology, which offers tools for partitioning and de-partitioning the model without disturbing the manufacturing oriented design geometry. The end result is a convenient framework to which high level analysis attributes can be applied, and from which detailed analysis models can be generated
with a high degree of controllability, repeatability and automation. There are multiple novel aspects to the approach, including its reusability, robustness to changes in model topology and the inherent links created between analysis models at different levels of fidelity and physics.

By utilising Simulation Intent, CAD modelling for simulation can be fully exploited and simulation work-flows can be more readily automated, reducing many repetitive manual tasks (e.g. the definition of appropriate coupling between elements of different types and the application of boundary conditions). The approach has been implemented and tested with practical examples, and
significant benefits are demonstrated.

Multiphysics Simulation of Thermal Plasma Produced by Lightning Strike on Aircraft Structures

A novel numerical technique is proposed to model thermal plasma of microseconds/milliseconds time-scale effect. Modelling thermal plasma due to lightning strike will allow the estimation of electric current density, plasma pressure, and heat flux at the surface of the aircraft structure. These input data can then be used for better estimation of the mechanical/thermal induced damage on the aircraft structures for better protection systems design. Thermal plasma generated during laser cutting, electric (laser) welding and other plasma processing techniques have been the focus of many researchers. Thermal plasma is a gaseous state that consists from a mixture of electrons, ions, and natural particles. Thermal plasma can be assumed to be in local thermodynamic equilibrium, which means the electrons and the heavy species have equal temperature. Different numerical techniques have been developed using a coupled Navier Stokes – Heat transfer – Electromagnetic equations based on the assumption that the thermal plasma is a single laminar gas flow. These previous efforts focused on generating thermal plasma of time-scale in the range of seconds. Lighting strike on aircraft structures generates thermal plasma of time-scale of milliseconds/microseconds, which makes the previous physics used not applicable. The difficulty comes from the Navier-Stokes equations as the fluid is simulated under shock load, this introducing significant changes in the density and temperature of the fluid.

Energy conserving schemes for the simulation of musical instrument contact dynamics

Collisions are an innate part of the function of many musical instruments. Due to the nonlinear nature of contact forces, special care has to be taken in the construction of numerical schemes for simulation and sound synthesis. Finite difference schemes and other time-stepping algorithms used for musical instrument modelling purposes are normally arrived at by discretising a Newtonian description of the system. However because impact forces are non-analytic functions of the phase space variables, algorithm stability can rarely be established this way. This paper presents a systematic approach to deriving energy conserving schemes for frictionless impact modelling. The proposed numerical formulations follow from discretising Hamilton׳s equations of motion, generally leading to an implicit system of nonlinear equations that can be solved with Newton׳s method. The approach is first outlined for point mass collisions and then extended to distributed settings, such as vibrating strings and beams colliding with rigid obstacles. Stability and other relevant properties of the proposed approach are discussed and further demonstrated with simulation examples. The methodology is exemplified through a case study on tanpura string vibration, with the results confirming the main findings of previous studies on the role of the bridge in sound generation with this type of string instrument.

The performance of concrete exposed to marine environments: predictive modelling and use of laboratory/on site test methods

This paper reports an approach by which laboratory based testing and numerical modelling can be combined to predict the long term performance of a range of concretes exposed to marine environments. Firstly, a critical review of the test methods for assessing the chloride penetration resistance of concrete is given. The repeatability of the different test results is also included. In addition to the test methods, a numerical simulation model is used to explore the test data further to obtain long-term chloride ingress trends. The combined use of testing and modelling is validated with the help of long-term chloride ingress data from a North Sea exposure site. In summary, the paper outlines a methodology for determining the long term performance of concrete in marine environments.

Advanced Modelling and Optimization of Infared Oven in Injection Stretch Blow-moulding for Energy Saving

In the production process of polyethylene terephthalate (PET) bottles, the initial temperature of preforms plays a central role on the final thickness, intensity and other structural properties of the bottles. Also, the difference between inside and outside temperature profiles could make a significant impact on the final product quality. The preforms are preheated by infrared heating oven system which is often an open loop system and relies heavily on trial and error approach to adjust the lamp power settings. In this paper, a radial basis function (RBF) neural network model, optimized by a two-stage selection (TSS) algorithm combined with partial swarm optimization (PSO), is developed to model the nonlinear relations between the lamp power settings and the output temperature profile of PET bottles. Then an improved PSO method for lamp setting adjustment using the above model is presented. Simulation results based on experimental data confirm the effectiveness of the modelling and optimization method.

FE Simulation of Viscoplastic Consistency Model

Finite element (FE) simulations provides an inexpensive alternative for material testingof new metal alloys. Carrying out experimental testing is expensive. Nanoindentation is particularly costly due to the equipment needed to work on such a scale. FE simulations provide an inexpensive means of material testing if accurately carried out. This paper will demonstrate the applicability and accuracy of using FE modelling for basic material tests and will propose that the viscoplastic model may be used for nanoindentation testing. The simulations will test the Young’s modulus of materials during analysis when an Abaqus VUMAT is used. The viscoplastic model is incorporated into a subroutine and is tested at the macroscopic scale against previous published results.

Measuring and modelling air mass flow rate in the injection stretch blow moulding process

The injection stretch blow moulding process involves the inflation and stretching of a hot preform into a mould to form bottles. A critical process variable and an essential input for process simulations is the rate of pressure increase within the preform during forming, which is regulated by an air flow restrictor valve. The paper describes a set of experiments for measuring the air flow rate within an industrial ISBM machine and the subsequent modelling of it with the FEA package ABAQUS. Two rigid containers were inserted into a Sidel SBO1 blow moulding machine and subjected to different supply pressures and air flow restrictor settings. The pressure and air temperature were recorded for each experiment enabling the mass flow rate of air to be determined along with an important machine characteristic known as the ‘dead volume’. The experimental setup was simulated within the commercial FEA package ABAQUS/Explicit using a combination of structural, fluid and fluid link elements that idealize the air flowing through an orifice behaving as an ideal gas under isothermal conditions. Results between experiment and simulation are compared and show a good correlation.

Modelling of Organic Rankine Cycle for Waste Heat Recovery Process in Supercritical Condition

Organic Rankine Cycle (ORC) is the most commonly used method for recovering energy from small sources of heat. The investigation of the ORC in supercritical condition is a new research area as it has a potential to generate high power and thermal efficiency in a waste heat recovery system. This paper presents a steady state ORC model in supercritical condition and its simulations with a real engine’s exhaust data. The key component of ORC, evaporator, is modelled using finite volume method, modelling of all other components of the waste heat recovery system such as pump, expander and condenser are also presented. The aim of this paper is to investigate the effects of mass flow rate and evaporator outlet temperature on the efficiency of the waste heat recovery process. Additionally, the necessity of maintaining an optimum evaporator outlet temperature is also investigated. Simulation results show that modification of mass flow rate is the key to changing the operating temperature at the evaporator outlet.

Finite element modelling of composite structures under crushing load

This paper details the theory and implementation of a composite damage model, addressing damage within a ply (intralaminar) and delamination (interlaminar), for the simulation of crushing of laminated composite structures. It includes a more accurate determination of the characteristic length to achieve mesh objectivity in capturing intralaminar damage consisting of matrix cracking and fibre failure, a load-history dependent material response, an isotropic hardening nonlinear matrix response, as well as a more physically-based interactive matrix-dominated damage mechanism. The developed damage model requires a set of material parameters obtained from a combination of standard and non-standard material characterisation tests. The fidelity of the model mitigates the need to manipulate, or "calibrate", the input data to achieve good agreement with experimental results. The intralaminar damage model was implemented as a VUMAT subroutine, and used in conjunction with an existing interlaminar damage model, in Abaqus/Explicit. This approach was validated through the simulation of the crushing of a cross-ply composite tube with a tulip-shaped trigger, loaded in uniaxial compression. Despite the complexity of the chosen geometry, excellent correlation was achieved with experimental results.

Multi-agents modelling of EV purchase willingness based on questionnaires

Traditional experimental economics methods often consume enormous resources of qualified human participants, and the inconsistence of a participant’s decisions among repeated trials prevents investigation from sensitivity analyses. The problem can be solved if computer agents are capable of generating similar behaviors as the given participants in experiments. An experimental economics based analysis method is presented to extract deep information from questionnaire data and emulate any number of participants. Taking the customers’ willingness to purchase electric vehicles (EVs) as an example, multi-layer correlation information is extracted from a limited number of questionnaires. Multi-agents mimicking the inquired potential customers are modelled through matching the probabilistic distributions of their willingness embedded in the questionnaires. The authenticity of both the model and the algorithm is validated by comparing the agent-based Monte Carlo simulation results with the questionnaire-based deduction results. With the aid of agent models, the effects of minority agents with specific preferences on the results are also discussed.