ABSTRACTAs the industry becomes more efficient all the time and the budgets narrower the demand forsafety equipment has also increased to minimize downtime and injuries. With the efficienciesgrowing in Swedish industry the supplier of this project has also seen an increased demand forsafety equipment. One product that has increased is their “Industryprotection” a productcreated to protect forklifts from colliding with walls in industry spaces.The purpose of this project is to generate a concept for a new design of the“Industryprotection”. It also involves creating a static FEM analysis of a collision between aforklift and the “Industryprotection”. The goal for the analysis is to create a steady modelwhich easily could be applied to different concepts. A study of the different calculationsmethods will also be executed.The project is split into two phases one conceptual phase and one calculation phase, the laststep of the project is to analyze the recommended concept versus the current one using theanalysis created.The work is limited to only go as far as a recommended concept is created, the dimensioning,final drawing, manufacturing preparation and validating is given to the supplier of the projectThe theory for this project is gathered from studying literature from market evaluation,computer aided calculation and design process. After theory for each phase of the project hasbeen gathered they have been modified to suit the phase. Also a large study of everyparameter in the FEM analysis has been executed to secure a stable model for the analysis.The output of the model became a simple and steady model which easily could be applied todifferent concepts. The study of different calculation methods of the model generated a fivetime increase of the amount of momentum the “Industryprotection” could endure. / SAMMANFATTNINGDå industrin hela tiden utvecklas mot att bli mer effektiv och budgeterna blir allt snävareställs också högre krav på att minimera driftstopp och skador. I takt med att arbetet aveffektiviseringar inom svensk industri har ökat har uppdragsgivaren EAB, sett en störreefterfrågan av säkerhetsutrustning .En kategori som har ökat är deras industriskydd som ärskapat för skydda väggar i industrilokaler.Syftet med projektet är att generera ett rekommenderat koncept över hur ett nytt industriskyddkan konstrueras. Vidare skall även en statisk analysmodell skapas över ett kollision scenariomellan truck och industriskyddet. Målet med analysen är att skapa en stabil modell som enkeltkan appliceras på olika koncept. Vidare skall även en studie kring olika beräkningsmetoderutföras.Arbetet är uppdelat i två olika faser en konceptfas samt en beräkningsfas. Som sista steg iprojektet används den skapade beräkningsmodellen för att utvärdera det rekommenderadekonceptetProjektet är begränsat till det stadie där ett rekommenderat koncept är framtaget. Slutligdimensionering, detaljkonstruktion, tillverkningsberedning och verifiering överlämnas tilluppdragsgivaren.Teorin kring projektet har samlats in genom att studera relevant litteratur inommarknadsundersökningar, datorstödberäkningshjälp samt konstruktionsprocessen. Efterstudien har teori anpassats och applicerats för detta specifika projekt. Under beräkningsfasenhar en omfattande studie kring samtliga parametrar utförts för att skapa en stabilberäkningsmodell.Utfallet av modellen blev en enkel och stabil beräkningsmodell som kan appliceras på olikakoncept med goda resultat. Utvärdering av beräkningsmetoder genererade en ökning av denabsorberade energin med en faktor på fem.

In this thesis, the possibility to numerically describing the behaviour that signifies shear type cracking in concrete is studied. Different means for describing cracking are evaluated where both methods proposed in design codes based on experiments and advanced finite element analyses with a non-linear material description are evaluated. It is shown that there is a large difference in the estimation of the crack width based on the calculation methods in design codes. The large difference occurs due to several of these methods do not account for shear friction in the crack face. The finite element method is an important tool for analysing the non-linear behaviour caused by cracking. It is especially of importance when combined with experimental investigations for evaluating load bearing capacity or establishing the structural health. It is shown that non-linear continuum material models can successfully be used to accurately describe the shear type cracking in concrete. A method based on plasticity and damage theory was shown to provide accurate estimations of the behaviour. The methods based on fracture mechanics with or without inclusion of damage theory, overestimated the stiffness after crack initiation considerably. The rotated crack approach of these methods gave less accurate descriptions of the crack pattern and underestimated the crack widths. After verification of the material model, realistic finite element models based on plasticity and damage theory are developed to analyse the cause for cracking in two large concrete structures. The Storfinnforsen hydropower buttress dam is evaluated where the seasonal temperature variation in combination with the water pressure have resulted in cracking. With the numerical model the cause for cracking can be explained and the crack pattern found in-situ is accurately simulated. The model is verified against measurements of variation in crest displacement and crack width with close agreement. The construction process of a balanced cantilever bridge, Gröndal Bridge, is numerically simulated and a rational explanation of the cause for cracking is presented. It is shown that large stresses and micro-cracks develop in the webs during construction, especially after tensioning the continuing tendons in the bottom flange. Further loads from temperature variation cause cracking in the webs that is in close agreement with the cracking found in-situ. The effect of strengthening performed on this bridge is also evaluated where the vertical Dywidag tendons so far seem to have been successful in stopping further crack propagation. / QC 20100730

Este trabalho apresenta a aplicação de técnicas de processamento paralelo na resolução de equações algébricas oriundas do Método de Elementos Finitos aplicado ao Eletromagnetismo, nos regimes estático e harmônico. As técnicas de programação paralelas utilizadas foram OpenMP, CUDA e GPUDirect, sendo esta última para as plataformas do tipo Multi-GPU. Os métodos iterativos abordados incluem aqueles do subespaço Krylov: Gradientes Conjugados, Gradientes Biconjugados, Conjugado Residual, Gradientes Biconjugados Estabilizados, Gradientes Conjugados para equações normais (CGNE e CGNR) e Gradientes Conjugados ao Quadrado. Todas as implementações fizeram uso das bibliotecas CUSP, CUSPARSE e CUBLAS. Para problemas estáticos, os seguintes pré-condicionadores foram adotados, todos eles com implementações paralelizadas e executadas na GPU: Decomposições Incompletas LU e de Cholesky, Multigrid Algébrico, Diagonal e Inversa Aproximada. Para os problemas harmônicos, apenas os dois primeiros pré-condicionadores foram utilizados, porém na sua versão sequencial, com execução na CPU, resultando em uma implementação híbrida CPU-GPU. As ferramentas computacionais desenvolvidas foram testadas na simulação de problemas de aterramento elétrico. No caso do regime harmônico, em que o fenômeno é regido pela Equação de Onda completa com perdas e não homogênea, a formulação adotada foi aquela em dois potenciais, A-V aresta-nodal. Em todas as situações, os aplicativos desenvolvidos para GPU apresentaram speedups apreciáveis, demonstrando a potencialidade dessa tecnologia para a simulação de problemas de larga escala na Engenharia Elétrica, com excelente relação custo-benefício. / This work presents the use of parallel processing techniques in Graphics Processing Units (GPU) for the solution of algebraic equations arising from the Finite Element modeling of electromagnetic phenomena, both in steadystate and time-harmonic regime. The techniques used were parallel programming OpenMP, CUDA and GPUDirect, the latter for those platforms of type Multi-GPU. The iterative methods discussed include those of the Krylov subspace: Conjugate Gradients, Bi-conjugate Gradients, Conjugate Residual, Bi-conjugate Gradients Stabilized, Conjugate Gradients for Normal Equations (CGNE and CGNR) and Conjugate Gradients Squared. All implementations have made use of CUSP, CUSPARSE and CUBLAS libraries. For the static problems, the following pre-conditioners were adopted, all with parallelized implementations and executed on the GPU: Incomplete decompositions, both LU and Cholesky, Algebraic Multigrid, Diagonal and Approximate Inverse. For the time-harmonic varying problems, only the first two pre-conditioners were used, but in their sequential version and running in the CPU, which yielded a hybrid CPU-GPU implementation. The developed computational tools were tested in the simulation of electrical grounding systems. In the case of the harmonic regime, in which the phenomenon is governed by the driven, lossy wave equation, the formulation adopted was that in two potential, the ungauged edge A-V formulation. In all cases, the developed GPU-based tools showed considerable speedups, showing that this is a promising technology for the simulation of large-scale Electrical Engineering problems, with excellent cost-benefit.

Este trabalho apresenta a aplicação de técnicas de processamento paralelo na resolução de equações algébricas oriundas do Método de Elementos Finitos aplicado ao Eletromagnetismo, nos regimes estático e harmônico. As técnicas de programação paralelas utilizadas foram OpenMP, CUDA e GPUDirect, sendo esta última para as plataformas do tipo Multi-GPU. Os métodos iterativos abordados incluem aqueles do subespaço Krylov: Gradientes Conjugados, Gradientes Biconjugados, Conjugado Residual, Gradientes Biconjugados Estabilizados, Gradientes Conjugados para equações normais (CGNE e CGNR) e Gradientes Conjugados ao Quadrado. Todas as implementações fizeram uso das bibliotecas CUSP, CUSPARSE e CUBLAS. Para problemas estáticos, os seguintes pré-condicionadores foram adotados, todos eles com implementações paralelizadas e executadas na GPU: Decomposições Incompletas LU e de Cholesky, Multigrid Algébrico, Diagonal e Inversa Aproximada. Para os problemas harmônicos, apenas os dois primeiros pré-condicionadores foram utilizados, porém na sua versão sequencial, com execução na CPU, resultando em uma implementação híbrida CPU-GPU. As ferramentas computacionais desenvolvidas foram testadas na simulação de problemas de aterramento elétrico. No caso do regime harmônico, em que o fenômeno é regido pela Equação de Onda completa com perdas e não homogênea, a formulação adotada foi aquela em dois potenciais, A-V aresta-nodal. Em todas as situações, os aplicativos desenvolvidos para GPU apresentaram speedups apreciáveis, demonstrando a potencialidade dessa tecnologia para a simulação de problemas de larga escala na Engenharia Elétrica, com excelente relação custo-benefício. / This work presents the use of parallel processing techniques in Graphics Processing Units (GPU) for the solution of algebraic equations arising from the Finite Element modeling of electromagnetic phenomena, both in steadystate and time-harmonic regime. The techniques used were parallel programming OpenMP, CUDA and GPUDirect, the latter for those platforms of type Multi-GPU. The iterative methods discussed include those of the Krylov subspace: Conjugate Gradients, Bi-conjugate Gradients, Conjugate Residual, Bi-conjugate Gradients Stabilized, Conjugate Gradients for Normal Equations (CGNE and CGNR) and Conjugate Gradients Squared. All implementations have made use of CUSP, CUSPARSE and CUBLAS libraries. For the static problems, the following pre-conditioners were adopted, all with parallelized implementations and executed on the GPU: Incomplete decompositions, both LU and Cholesky, Algebraic Multigrid, Diagonal and Approximate Inverse. For the time-harmonic varying problems, only the first two pre-conditioners were used, but in their sequential version and running in the CPU, which yielded a hybrid CPU-GPU implementation. The developed computational tools were tested in the simulation of electrical grounding systems. In the case of the harmonic regime, in which the phenomenon is governed by the driven, lossy wave equation, the formulation adopted was that in two potential, the ungauged edge A-V formulation. In all cases, the developed GPU-based tools showed considerable speedups, showing that this is a promising technology for the simulation of large-scale Electrical Engineering problems, with excellent cost-benefit.

In this dissertation, the accuracy and reliability of non-linear finite element computations in application to surgical simulation is evaluated. The evaluation is performed through comparison between the experiment and finite element analysis of indentation of soft tissue phantom and human brain phantom. The evaluation is done in terms of the forces acting on the cylindrical Aluminium indenter and deformation of the phantoms due to these forces. The deformation of the phantoms is measured by tracking 3D motions of X-ray opaque markers implanted in the direct neighbourhood under the indenter using a custom-made biplane X-ray image intensifiers (XRII) system. The phantoms are made of Sylgard® 527 gel to simulate the hyperelastic constitutive behaviour of the brain tissue. The phantoms are prepared layer by layer to facilitate the implantation of the X-ray opaque markers. The modelling of soft tissue phantom indentation and human brain phantom indentation is performed using the ABAQUSTM/Standard finite element solver. Realistic geometry model of the human brain phantom obtained from Magnetic Resonance images is used. Specific constitutive properties of the phantom layers determined through uniaxial compression tests are used in the model. The models accurately predict the indentation force-displacement relations and marker displacements in both soft tissue phantom indentation and human brain phantom indentation. Good agreement between the experimental and modelling results verifies the reliability and accuracy of the finite element analysis techniques used in this study and confirms the predictive power of these techniques in application to surgical simulation.

Cette thèse porte sur l’étude du comportement mécanique et vibratoire d’un composite biosourcé incorporant un matériau viscoélastique. Les matériaux étudiés sont des stratifiés en composite lin/greenpoxy et des stratifiés viscoélastiques composés d'un noyau viscoélastique en caoutchouc naturel confiné entre deux composites. La première partie du travail a été consacrée à l’étude de l’influence de l’intégration de la couche viscoélastique sur le comportement mécanique des composites. L’analyse des résultats expérimentaux et l’observation des signaux d’émission acoustique obtenus dans ces composites soumis à différentes sollicitations mécaniques en statique et en fatigue ont permis d’identifier les signatures acoustiques des mécanismes d’endommagement prépondérants dans les deux matériaux. Dans un deuxième temps, les propriétés dynamiques de ces composites ont été déterminées à partir des essais de vibration. Les résultats obtenus ont montré que la couche viscoélastique a joué un rôle majeur dans l'amortissement et la dissipation d'énergie des composites. Suite à cette analyse, nous avons mis en place une procédure, utilisant la méthode des éléments finis, pour calculer l’amortissement de ces matériaux. Dans le but de mettre en évidence l’influence des caractéristiques de la couche viscoélastique, une étude paramétrique a été menée sur le composite viscoélastique, permettant d’optimiser l’amortissement de ce matériau en faisant varier divers paramètres. Enfin, le comportement visqueux des composites a été caractérisé par la méthode de résonance non linéaire en faisant varier l’amplitude d’excitation. / This thesis focuses on the study of the mechanical and vibration behaviour of a flax fibre reinforced composites with and without an interleaved natural viscoelastic layer. The composite materials have been characterized experimentally using different mechanical and vibrational tests. First, both types of composites were studied using uni-axial tensile and three-points bending tests. Acoustic emission (AE) has been often used for the identification and characterization of micro failure mechanisms in composites. The results showed that these composites have very high specific characteristics. It can be used for applications currently using composites reinforced with synthetic fibres such glass, carbon…. Next, experimental and finite element vibration analyses were carried out on the composites with and without an interleaved natural viscoelastic layer. A good agreement between the two methods was obtained. It has been shown that the viscoelastic layer plays a major role in damping because it has a high level of energy dissipation. Therefore, it improves with a significant way the modal properties of the composite. Finally, nonlinear resonance tests were performed on the composites. It has been shown that the viscoelastic layer generates a nonlinear behaviour in the material. The linear and nonlinear, elastic and dissipative parameters have been calculated to deduce finally that nonlinear parameters are more sensitive to heterogeneities than those derived from linear vibration tests.