Formulation of Locally Reacting Surfaces in FDTD/K-DWM Modelling of Acoustic Spaces

In this paper, we present new methods for constructing and analysing formulations of locally reacting surfaces that can be used in finite difference time domain (FDTD) simulations of acoustic spaces. Novel FDTD formulations of frequency-independent and simple frequency-dependent impedance boundaries are proposed for 2D and 3D acoustic systems, including a full treatment of corners and boundary edges. The proposed boundary formulations are designed for virtual acoustics applications using the standard leapfrog scheme based on a rectilinear grid, and apply to FDTD as well as Kirchhoff variable digital waveguide mesh (K-DWM) methods. In addition, new analytic evaluation methods that accurately predict the reflectance of numerical boundary formulations are proposed. numerical experiments and numerical boundary analysis (NBA) are analysed in time and frequency domains in terms of the pressure wave reflectance for different angles of incidence and various impedances. The results show that the proposed boundary formulations structurally adhere well to the theoretical reflectance. In particular, both reflectance magnitude and phase are closely approximated even at high angles of incidence and low impedances. Furthermore, excellent agreement was found between the numerical boundary analysis and the experimental results, validating both as tools for researching FDTD boundary formulations.

Automated mixed dimensional modelling from 2D and 3D CAD models

The motivation for this paper is to present procedures for automatically creating idealised finite element models from the 3D CAD solid geometry of a component. The procedures produce an accurate and efficient analysis model with little effort on the part of the user. The technique is applicable to thin walled components with local complex features and automatically creates analysis models where 3D elements representing the complex regions in the component are embedded in an efficient shell mesh representing the mid-faces of the thin sheet regions. As the resulting models contain elements of more than one dimension, they are referred to as mixed dimensional models. Although these models are computationally more expensive than some of the idealisation techniques currently employed in industry, they do allow the structural behaviour of the model to be analysed more accurately, which is essential if appropriate design decisions are to be made. Also, using these procedures, analysis models can be created automatically whereas the current idealisation techniques are mostly manual, have long preparation times, and are based on engineering judgement. In the paper the idealisation approach is first applied to 2D models that are used to approximate axisymmetric components for analysis. For these models 2D elements representing the complex regions are embedded in a 1D mesh representing the midline of the cross section of the thin sheet regions. Also discussed is the coupling, which is necessary to link the elements of different dimensionality together. Analysis results from a 3D mixed dimensional model created using the techniques in this paper are compared to those from a stiffened shell model and a 3D solid model to demonstrate the improved accuracy of the new approach. At the end of the paper a quantitative analysis of the reduction in computational cost due to shell meshing thin sheet regions demonstrates that the reduction in degrees of freedom is proportional to the square of the aspect ratio of the region, and for long slender solids, the reduction can be proportional to the aspect ratio of the region if appropriate meshing algorithms are used.

Dealing with non-linearity in shape modelling of articulated objects

We address the problem of non-linearity in 2D Shape modelling of a particular articulated object: the human body. This issue is partially resolved by applying a different Point Distribution Model (PDM) depending on the viewpoint. The remaining non-linearity is solved by using Gaussian Mixture Models (GMM). A dynamic-based clustering is proposed and carried out in the Pose Eigenspace. A fundamental question when clustering is to determine the optimal number of clusters. From our point of view, the main aspect to be evaluated is the mean gaussianity. This partitioning is then used to fit a GMM to each one of the view-based PDM, derived from a database of Silhouettes and Skeletons. Dynamic correspondences are then obtained between gaussian models of the 4 mixtures. Finally, we compare this approach with other two methods we previously developed to cope with non-linearity: Nearest Neighbor (NN) Classifier and Independent Component Analysis (ICA).

Efficient modelling and optimisation of hybrid multilayered plates subject to ballistic impact

Among the key challenges present in the modelling and optimisation of composite structures against impact is the computational expense involved in setting up accurate simulations of the impact event and then performing the iterations required to optimise the designs. It is of more interest to find good designs given the limitations of the resources and time available rather than the best possible design. In this paper, low cost but sufficiently accurate finite element (FE) models were generated in LS Dyna for several experimentally characterised materials by semi-automating the modelling process and using existing material models. These models were then used by an optimisation algorithm to generate new hybrid offspring, leading to minimum weight and/or cost designs from a selection of isotropic metals, polymers and orthotropic fibre-reinforced laminates that countered a specified impact threat. Experimental validation of the optimal designs thus identified was then successfully carried out using a single stage gas gun. With sufficient computational hardware, the techniques developed in this pilot study can further utilise fine meshes, equations of state and sophisticated material models, so that optimal hybrid systems can be identified from a wide range of materials, designs and threats.

Element-Free Galerkin modelling of composite damage

Delamination and matrix cracking are routine damage mechanisms, observed by post-mortem analysis of laminated structures containing geometrical features such as notches or bolts. Current finite element tools cannot explicitly model an intralaminar matrix microcrack, except if the location of the damage is specified a priori. In this work, a meshless technique, the Element-Free Galerkin (EFG) method, is utilized for the first time to simulate delamination (interlaminar) and intralaminar matrix microcracking in composite laminates.

The Behavior of compressively loaded stiffener runout specimens - Part II: Finite element analysis

The termination of stiffeners in composite aircraft structures give rise to regions of high interlaminar shear and peel stresses as the load in the stiffener is diffused into the skin. This is of particular concern in co-cured composite stiffened structures where there is a relatively low resistance to through-thickness stress components at the skin-stiffener interface. In Part I, experimental results of tested specimens highlighted the influence of local design parameters on their structural response. Indeed some of the observed behavior was unexpected. There is a need to be able to analyse a range of changes in geometry rapidly to allow the analysis to form an integral part of the structural design process.

This work presents the development of a finite element methodology for modelling the failure process of these critical regions. An efficient thick shell element formulation is presented and this element is used in conjuction with the Virtual Crack Closure Technique (VCCT) to predict the crack growth characteristics of the modelled specimens. Three specimens were modelled and the qualitative aspects of crack growth were captured successfully. The shortcomings in the quantitative correlation between the predicted and observed failure loads are discussed. There was evidence to suggest that high through-thickness compressive stresses enhanced the fracture toughness in these critical regions.

Modelling the chromosphere and transition region of ε Eri (K2 V)

Measurements of ultraviolet line fluxes from Space Telescope Imaging Spectrograph and Far-Ultraviolet Spectroscopic Explorer spectra of the K2-dwarf e Eri are reported. These are used to develop new emission measure distributions and semi-empirical atmospheric models for the chromosphere and lower transition region of the star. These models are the most detailed constructed to date for a main-sequence star other than the Sun. New ionization balance calculations, which account for the effect of finite density on dielectronic recombination rates, are presented for carbon, nitrogen, oxygen and silicon. The results of these calculations are significantly different from the standard Arnaud & Rothenflug ion balance, particularly for alkali-like ions. The new atmospheric models are used to place constraints on possible first ionization potential (FIP)-related abundance variations in the lower atmosphere and to discuss limitations of single-component models for the interpretation of certain optically thick line fluxes. © 2005 RAS.

Regularized Online Mixture of Gaussians for Background Subtraction

Mixture of Gaussians (MoG) modelling [13] is a popular approach to background subtraction in video sequences. Although the algorithm shows good empirical performance, it lacks theoretical justification. In this paper, we give a justification for it from an online stochastic expectation maximization (EM) viewpoint and extend it to a general framework of regularized online classification EM for MoG with guaranteed convergence. By choosing a special regularization function, l1 norm, we derived a new set of updating equations for l1 regularized online MoG. It is shown empirically that l1 regularized online MoG converge faster than the original online MoG .

Modelling of poly(ethylene terephthalate) in injection stretch-blow moulding

Simulations of the injection stretch-blow moulding process have been developed for the manufacture of poly(ethylene terephthalate) bottles using the commercial finite element package ABAQUS/standard. Initially a simulation of the manufacture of a 330 mL bottle was developed with three different material models (hyperelastic, creep, and a non-linear viscoelastic model (Buckley model)) to ascertain their suitability for modelling poly(ethylene terephthalate). The Buckley model was found to give results for the sidewall thickness that matched best with those measured from bottles off the production line. Following the investigation of the material models, the Buckley model was chosen to conduct a three-dimensional simulation of the manufacture of a 2 L bottle. It was found that the model was also capable of predicting the wall thickness distribution accurately for this bottle. In the development of the three-dimensional simulation a novel approach, which uses an axisymmetric model until the material reaches the petaloid base, was developed. This resulted in substantial savings in computing time. © 2000 IoM Communication Ltd.

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.

NONLINEAR NUMERICAL MODELLING OF LIGHTENING STRIKE EFFECT ON COMPOSITE PANELS WITH TEMPERATURE DEPENDANT MATERIAL PROPERTIES

Lightning strike is one of the challenges that the aerospace industry is facing in an effort to increase the percentage of composite materials used in aircraft structures. Lightning strike damage is due to high orthotropic electric resistivity of the composite panels, which leads to high thermal loads that cause decomposition of the epoxy and delimitations of the laminates. Yet, experimental testing of lightning strike on aircraft panels is expensive due to the large number of design parameters that can control the inflicted damage. A coupled thermal-electrical finite element analysis is used to investigate the design variables space that can affect lightning strike damage on epoxy/graphite composite panels. The contribution of this study is modeling the composite panels’ material properties as temperature dependent, which was excluded by other researchers. A number of practical solutions to minimize the damage effect are proposed. Two set of experimental results are used to verify the numerical ones. One experimental set for plain composite panel, and second one for composite panels with joints

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.

Simulation of three-dimension cutting force and tool deflection in the end milling operation based on finite element method

This paper presents a 3D simulation system which is employed in order to predict cutting forces and tool deflection during end-milling operation. In order to verify the accuracy of 3D simulation, results (cutting forces and tool deflection) were compared with those based on the theoretical relationships, in terms of agreement with experiments. The results obtained indicate that the simulation is capable of predicting the cutting forces and tool deflection.

Integrating geophysical data and upscaling techniques to support regional groundwater flow modelling: A practical example of the Lagan Valley, Northern Ireland

When studying heterogeneous aquifer systems, especially at regional scale, a degree of generalization is anticipated. This can be due to sparse sampling regimes, complex depositional environments or lack of accessibility to measure the subsurface. This can lead to an inaccurate conceptualization which can be detrimental when applied to groundwater flow models. It is important that numerical models are based on observed and accurate geological information and do not rely on the distribution of artificial aquifer properties. This can still be problematic as data will be modelled at a different scale to which it was collected. It is proposed here that integrating geophysics and upscaling techniques can assist in a more realistic and deterministic groundwater flow model. In this study, the sedimentary aquifer of the Lagan Valley in Northern Ireland is chosen due to intruding sub-vertical dolerite dykes. These dykes are of a lower permeability than the sandstone aquifer. The use of airborne magnetics allows the delineation of heterogeneities, confirmed by field analysis. Permeability measured at the field scale is then upscaled to different levels using a correlation with the geophysical data, creating equivalent parameters that can be directly imported into numerical groundwater flow models. These parameters include directional equivalent permeabilities and anisotropy. Several stages of upscaling are modelled in finite element. Initial modelling is providing promising results, especially at the intermediate scale, suggesting an accurate distribution of aquifer properties. This deterministic based methodology is being expanded to include stochastic methods of obtaining heterogeneity location based on airborne geophysical data. This is through the Direct Sample method of Multiple-Point Statistics (MPS). This method uses the magnetics as a training image to computationally determine a probabilistic occurrence of heterogeneity. There is also a need to apply the method to alternate geological contexts where the heterogeneity is of a higher permeability than the host rock.