Underwater explosions are often modeled with Arbitrary Lagrangian-Eulerian (ALE) Finite Element Methods. The objective of this thesis is to validate the simulation method, with respect to the propagating shock wave. A two-dimensional axisymmetric model of a spherical TNT charge submerged in water is simulated using LS-DYNA. The explosive is modeled with the Burn Fraction technique and the Jones-Wilkins-Lee equation of state. Water is modeled as a non-viscous fluid, with the Grüneisen equation of state. The convergence for different mesh resolutions, the effect of different advection methods, and varied constants in the artificial viscosity are examined. Generally, the simulations agree well with empirical results, but the maximum pressure diminishes more rapidly with distance compared to experiments. The excessive dampening is most notable in the early stages of the propagation. Also, unexpected oscillations are observed near the discontinuity. The choice of advection scheme and constants in the artificial viscosity do not resolve the issues suggesting that other numerical techniques for treating the discontinuity should be considered. / Undervattensexplosioner simuleras ofta med ALE-baserade finita elementmetoder. Detta examensarbete avser att validera simuleringsmetoden med hänsyn till stötvågens utbredning i vattnet. En tvådimensionell axisymmetrisk modell av en sfärisk TNT-laddning nedsänkt i vatten simuleras i LS-DYNA. Laddningen modelleras med hjälp av brinnfraktioner och Jones-Wilkins-Lee tillståndsekvation. Vattnet modelleras som en inviskös fluid tillsammans med Grüneisens tillståndsekvation. Nätkonvergens, val av advektionsmetod och ändring av konstanter i den artificiella viskositeten studeras. Övergripande resultat stämmer väl överens med empirisk data, men stötvågens topptryck avtar fortare än väntat. Denna dämpning är tydligast i utredningens tidiga skeden. Dessutom observeras oväntade oscillationer kring stötvågens diskontinuerliga tryckprofil. Val av advektionsmetod och konstanter tillhörande artificiella viskositeten verkar ha liten betydelse för resultaten. En alternativ numerisk metod för behandling av stötvågens diskontinuitet bör implementeras.

This thesis comprehends a series of nonlinear numerical studies with the finite element software's LS-Dyna and Impetus AFEA. The main focus lies on a comparative crash analysis of an aluminium beam profile which the company Sapa technology has used during their crash analysis. The aluminium profile has the characteristic of having different thickness over span ratios within the profile. This characteristic provided the opportunity to conduct a performance investigation of shell and solid elements with finite element analysis. Numerical comparisons were made between shell and solid elements where measurable parameters such as internal energy, simulation times, buckling patterns and material failures were compared to physical tests conducted prior to this thesis by Sapa technology. The performance investigation of shell and solid elements was initiated by creating models of the aluminium profile for general visualization and to facilitate the meshing of surfaces. The meshing procedure was considered to be an important factor of the analysis. The mesh quality and element orientations were carefully monitored in order to achieve acceptable results when the models were compared to physical tests. Preliminary simulations were further conducted in order to obtain a clear understanding of software parameters when performing crash simulations in LS-Dyna and Impetus AFEA. The investigated parameters were element formulations and material models. A general parameter understanding facilitated in the selection of parameters for actual simulations, where material failure and damage models were used. In conclusion, LS-Dyna was observed to provide a bigger internal energy absorption during the crushing of the beam with longer simulation times for solid elements when compared to shell elements. Impetus AFEA did on the other hand provide results close to physical test data with acceptable simulation times when compared to physical tests. The result difference obtained from the FE-software's in relation to physical crash experiments were considered to be varied but did indicate that shell elements were efficient enough for the specific profile during simulations with LS-Dyna. Impetus AFEA proved that the same time to be numerically efficient for energy approximations with solid elements refined with the third polynomial.

The basis for performing a thermo-mechanical staggered coupled heat source analysis of a welding simulation is implemented into LS-DYNA. In this report, three methods for initiating the heat source’s mechanical motion during girth and tubular joint welding are developed and evaluated. The first method is a reformulation of the equations used at Det Norske Veritas, the second is an incorporation of the equations into excel and the third is a standalone third party software. The most efficient of the developed methods turned out to be the software which creates k-files which are implemented into the main k-file using LS-PrePost. All methods have been visually and numerically evaluated using Excel, LS-DYNA and LS-PrePost.

This study investigates the mechanical and geometrical attributes of egg–box energy absorbing structures as crash safety barriers in the automotive industry. The research herein was originated from the earlier work of Prof. Shirvani, patented and further investigated by Cellbond Composites Ltd. who has invested in further research, for developing an analytical tool for geometric optimisation as an enhanced resolution to various shapes and materials. Energy absorption in egg-box occurs through plastic deformation of cell walls, examined through non–linear finite element simulations using ANSYS® and ANSYS/LS–DYNA® FE packages. Experimental dynamic crash tests have been designed to verify the validity of the FE simulations. Geometrical models are defined as 3D graphical representations, outlined in detail. Further, the impact behaviour of commercially pure aluminium egg-box energy absorbers is studied to identify the optimum design parameters describing the geometry of the structure. A simulation-based multi-objective optimisation strategy is employed to find a set of Pareto-optimal solutions where each solution represents a trade-off point with respect to the two conflicting objectives: the maximum impact force and the energy absorption capacity of the structure. The aim is to simultaneously minimise the former and maximise the latter, in the attempt to find purpose–specific optimal egg–box geometries. In light of the associated outcomes, it is shown that egg–box geometries with < ω ), thin walls (t < 1mm), short inter–peak distances and small peak diameters. M – < ω ), thin walls (t < 1mm), lengthy inter–peak distances and smaller peak diameters. It is concluded that, egg–box structures combined in the form of sandwich panels can be designed per application to act as optimised energy absorbers. Results of the proposed optimised sandwich structure are verified using analytical techniques.

Nowadays FE-simulations for sheet metal forming process are used to reduce the tryout phase in automotive industries. But the complex forming simulation processes are very challenging. One of the challenges is to create an FE-model which can be used to analyze the effects of cushion motion on the forming process. This thesis is focused on creating an FE model for two dies single cushion sheet metal forming press which can be used to analyze the effects of cushion motion on the contact pressure between the stamping tools and blanks. Using LS-PrePost a model with rigid stamping tools and cushion was created where the two blanks were of different thicknesses. After the model creation LS-DYNA was used for the simulations. The results showed that the cushion is moving in all DOFs and due to this movements, uneven contact pressure distribution is seen on the blanks and stamping tools.

A numerical simulation is conducted using LS-DYNA to simulate hard projectile impact on a friction stir welded (FSW) plate. As the hard projectile has a wide range of velocity, mass and shape, when referring to AMC 25.963(e) of CS-25, ―Fuel Tank Access Cover‖, the hard projectile can be defined as 9.5 mm cubic-shaped steel engine debris with an initial impact velocity of 213.4 m/s (700 ft/s). This preliminary study was to evaluate whether the fuel tank adjacent skin panel joined by FSW would pass the regulation. First, the geometry and Johnson-Cook material model of the FSW joint were developed based on previous experimental research and validated by comparison with the tensile test on the FSW specimen. Then the impact on an Aluminium Alloy 2024 (AA 2024) plate without FSW was modelled. The minimum thickness of a homogeneous AA 2024 plate which could withstand the impact from engine debris is 3 mm. Finally the impact on 3 mm thick AA 2024 FSW plate was simulated. The welding induced residual stress was implemented in the plate model. The impact centre was changed from the nugget zone to the thermo-mechanically affected zone, heat-affected zone and base material zone of the FSW joint. Penetration only occurred in the model with impact centre on the nugget zone. Additional simulation indicated that increasing the thickness of the FSW plate to 3.6 mm could prevent the penetration.

The bump at the end of the railway bridge is a result of differential movement between the bridge deck and the approach embankment. The movement can have the form of a bump or a dip. Either defect in the track geometry can cause significant problems in track performance. The current state of practice was evaluated by conducting a literature review and an industry survey. According to the survey, approximately half of all railway bridges are affected by the bump/dip. The total annual cost for repairing these bridge transition problems is estimated at $26 million. This does not take into account the considerable cost resulting from speed reductions that railroads must place on trains at these locations. In addition to the increased maintenance costs, the bump/dip leads to higher impact loads, uncomfortable rides and possible safety hazards. The track response due to the bump at the end of the bridge was evaluated by creating a 4-D finite element model of the train, track structure and track substructure. The motion of the train model across a bridge/approach transition, with and without a bump/dip, was then simulated using LS-DYNA. It was found that a track modulus differential alone (no bump/dip) at a bridge/approach location leads to impact forces as well as increased ballast and subgrade pressures on the approach. This instigates the formation of a bump or dip in the track. The track response is increased when a bump/dip is present in the track profile. A parametric study looking at the influence of train direction, train speed, bump/dip size, approach embankment soil modulus, approach tie material, bridge tie material, bridge deck type, ballast thickness and approach tie length on the magnitude of impact forces, track deflection, ballast and subgrade pressures was also performed with the model. Finally, a design solution to minimize the bump at the end of the bridge is proposed. The solution involves installing varying length steel bars into a soft subgrade approach embankment. The solution addresses both the settlement and track modulus differential between the bridge and the embankment. A full-scale field test of this prototype solution is underway.

In this thesis large deformation behavior of shell structures under high intensity transient loading conditions is investigated by means of finite element method. For this purpose an explicit finite element program is developed with interactive user interface. The developed program deals with geometric and material nonlinearities which stem from large deformation elastic - plastic behavior. Results of the developed code are compared with the experimental data taken from the literature and simulation results of the commercial finite element program Ls-Dyna. Moreover, sensitivity study is carried out for mesh size, element type and material model parameters. After the comparison and verification of the obtained results, it is concluded that converged and reasonable results are achieved.

The demand for artificial heart valve replacements is increasing as a result of birth defects, ageing and disease. Collaboration between engineers, biologists and mathematicians is necessary to handle problems related to biocompatibility and fluid dynamics. As a result of the increased demand for artificial heart valves, many new designs have been developed recently. A method to test those designs is to use mathematical modeling. This method has a relatively low-cost and can be used as a preliminary tool before expensive prototypes are created. This research analyzes the use of the numerical modeling software LS-DYNA for large-displacement fluid-structure interaction. It is a preliminary study aimed at the analysis of heart valve dynamics. In particular, a channel with flap model is created in LS-DYNA. The model's physics, convergence and ability to handle large deformations is investigated.

A popular way to reduce weight in industrial products without compromising the strength or stiffness is to replace components made of metal by plastics that have been reinforced by glass fibers. When fibers are introduced in a plastic, the resulting composite usually becomes anisotropic, which makes it much more complex to work with in simulation software. This thesis looks at modeling of such a composite using the multi-scale material modeling tool Digimat. An injection molding simulation of a brush cutter casing made of a short fiber reinforced plastic has been performed in order to obtain information about the glass fiber orientations, and thus the anisotropy, in each material point. That information has then been transferred over from the injection mesh to the structural mesh via a mapping routine. An elasto-viscoplastic material model with failure has been employed and calibrated against experimental data to find the corresponding material parameters. Lastly, a finite element analysis simulating a drop test has been performed. The results from the analysis have been compared with a physical drop test in order to evaluate the accuracy of the methodology used. The outcome has been discussed, conclusions have been drawn and suggestions for further studies have been presented.