Layerwise mixed models for analysis of multilayered piezoelectric composite plates using least-squares formulation

This work provides an assessment of layerwise mixed models using least-squares formulation for the coupled electromechanical static analysis of multilayered plates. In agreement with three-dimensional (3D) exact solutions, due to compatibility and equilibrium conditions at the layers interfaces, certain mechanical and electrical variables must fulfill interlaminar C-0 continuity, namely: displacements, in-plane strains, transverse stresses, electric potential, in-plane electric field components and transverse electric displacement (if no potential is imposed between layers). Hence, two layerwise mixed least-squares models are here investigated, with two different sets of chosen independent variables: Model A, developed earlier, fulfills a priori the interiaminar C-0 continuity of all those aforementioned variables, taken as independent variables; Model B, here newly developed, rather reduces the number of independent variables, but also fulfills a priori the interlaminar C-0 continuity of displacements, transverse stresses, electric potential and transverse electric displacement, taken as independent variables. The predictive capabilities of both models are assessed by comparison with 3D exact solutions, considering multilayered piezoelectric composite plates of different aspect ratios, under an applied transverse load or surface potential. It is shown that both models are able to predict an accurate quasi-3D description of the static electromechanical analysis of multilayered plates for all aspect ratios.

Multiobjective design of viscoelastic laminated composite sandwich panels

The optimal design of laminated sandwich panels with viscoelastic core is addressed in this paper, with the objective of simultaneously minimizing weight and material cost and maximizing modal damping. The design variables are the number of layers in the laminated sandwich panel, the layer constituent materials and orientation angles and the viscoelastic layer thickness. The problem is solved using the Direct MultiSearch (DMS) solver for multiobjective optimization problems which does not use any derivatives of the objective functions. A finite element model for sandwich plates with transversely compressible viscoelastic core and anisotropic laminated face layers is used. Trade-off Pareto optimal fronts are obtained and the results are analyzed and discussed.

Critical Velocity obtained using Simplified Models of the Railway Track: Viability and Applicability

Increased demands on the capacity of the railway network gave rise to new issues related to the dynamic response of railway tracks subjected to moving vehicles. Thus, it becomes important to evaluate the applicability of traditionally used simplified models which have a closed form solution. Regarding simplified models, transversal vibrations of a beam on a visco-elastic foundation subjected to a moving load are considered. Governing equations are obtained by Hamilton’s principle. Shear distortion, rotary inertia and effect of axial force are accounted for. The load is introduced as a time varying force moving at a constant velocity. Transversal vibrations induced by the load are solved by the normal-mode analysis. Reflected waves at the extremities of the full beam are avoided by introduction of semi-infinite elements. Firstly, the critical velocity obtained from this model is compared with results of an undamped Euler- Bernoulli formulation with zero axial force. Secondly, a finite element model in ABAQUS is examined. The new contribution lies in the introduction of semi- infinite elements and in the first step to a systematic comparison, which have not been published so fa

Avaliação do comportamento de uma ponte metálica centenária sujeita à ação de tráfego ferroviário ligeiro

A presente dissertação tem como objetivo principal a análise numérica do comportamento dinâmico de uma ponte ferroviária, sob ação de tráfego ligeiro ferroviário. Neste contexto são apresentados alguns fundamentos teóricos a ter em conta nestes domínios, visando uma melhor compreensão dos fenómenos existentes no comportamento dinâmico de pontes ferroviárias quando sujeitas ao tráfego. O caso de estudo teve como foco a ponte Luiz I, uma ponte metálica situada sobre o rio Douro, que liga as cidades do Porto e Vila Nova de Gaia, sob ação de tráfego ligeiro ferroviário no seu tabuleiro superior para a condição anterior aos trabalhos de reabilitação e reforço realizados entre 2004 e 2005. Para o efeito foi desenvolvido um modelo numérico de elementos finitos da ponte realizado com recurso ao programa ANSYS, assim como um modelo numérico do veículo do Metro de Lisboa. Com base nestes modelos foram obtidos os parâmetros modais, nomeadamente as frequências naturais e os modos de vibração de toda a estrutura e do veículo. O estudo do comportamento dinâmico da ponte foi realizado por intermédio de uma metodologia de cargas móveis e de interação veículo-estrutura, através da ferramenta computacional Train-Bridge Interaction (TBI). As análises dinâmicas foram efetuadas para a passagem dos veículos de passageiros das redes de Metros do Porto e Lisboa. Nestas análises é estudada a resposta da estrutura em função da variabilidade ao nível da secção transversal, dependência do tramo, influência do veículo, da sua velocidade de circulação e impacto das frequências de vibração estimadas pelo modelo numérico.

Experimental and numerical approaches for structural assessment in new footbridge designs (SFRSCC–GFPR hybrid structure)

Within the civil engineering field, the use of the Finite Element Method has acquired a significant importance, since numerical simulations have been employed in a broad field, which encloses the design, analysis and prediction of the structural behaviour of constructions and infrastructures. Nevertheless, these mathematical simulations can only be useful if all the mechanical properties of the materials, boundary conditions and damages are properly modelled. Therefore, it is required not only experimental data (static and/or dynamic tests) to provide references parameters, but also robust calibration methods able to model damage or other special structural conditions. The present paper addresses the model calibration of a footbridge bridge tested with static loads and ambient vibrations. Damage assessment was also carried out based on a hybrid numerical procedure, which combines discrete damage functions with sets of piecewise linear damage functions. Results from the model calibration shows that the model reproduces with good accuracy the experimental behaviour of the bridge.

Structural assessment of a masonry vault in Portugal

This paper reports on a structural safety assessment and performance evaluation of the upper choir of the Santa Maria de Belém Church in the Jerónimos monastery, Lisbon, one of the most important cultural heritage buildings in Portugal. The possibility of adding a new 20 t organ to the upper choir and its effects on the church structure's response are presented. A refined and a simplified finite-element model is developed to investigate the structure's performance under self-weight and seismic actions. A sensitivity analysis is performed to investigate the effect of masonry mechanical properties and rib cross-sections on the structural response, given the difficulty in accurately obtaining this information. The results show that the safety level of the structure is acceptable, even in the case of adding a heavy new organ.

Analysis of the seismic performance of a two storey log house

The dearth of knowledge on the load resistance mechanisms of log houses and the need for developing numerical models that are capable of simulating the actual behaviour of these structures has pushed efforts to research the relatively unexplored aspects of log house construction. The aim of the research that is presented in this paper is to build a working model of a log house that will contribute toward understanding the behaviour of these structures under seismic loading. The paper presents the results of a series of shaking table tests conducted on a log house and goes on to develop a numerical model of the tested house. The finite element model has been created in SAP2000 and validated against the experimental results. The modelling assumptions and the difficulties involved in the process have been described and, finally, a discussion on the effects of the variation of different physical and material parameters on the results yielded by the model has been drawn up.

Blast loading of masonry infills: testing and simulation

This work intends to present a newly developed test setup for dynamic out-of-plane loading using underWater Blast Wave Generators (WBWG) as loading source. Underwater blasting operations have been, during the last decades, subject of research and development of maritime blasting operations (including torpedo studies), aquarium tests for the measurement of blasting energy of industrial explosives and confined underwater blast wave generators. WBWG allow a wide range for the produced blast impulse and surface area distribution. It also avoids the generation of high velocity fragments and reduces atmospheric sound wave. A first objective of this work is to study the behavior of masonry infill walls subjected to blast loading. Three different masonry walls are to be studied, namely unreinforced masonry infill walls and two different reinforcement solutions. These solutions have been studied previously for seismic action mitigation. Subsequently, the walls will be simulated using an explicit finite element code for validation and parametric studies. Finally, a tool to help designers to make informed decisions on the use of infills under blast loading will be presented.

Otimização das características de um chassis tubular de automóvel de competição com tração a duas rodas motrizes

Dissertação de mestrado integrado em Engenharia Mecânica

Estudo numérico sobre o reforço de pilares de betão armado com cantoneiras e presilhas metálicas: análise paramétrica e modelo de cálculo

Dissertação de mestrado integrado em Engenharia Civil

Anatomic shoulder arthroplasty: consequences of humeral head malposition

Purpose Third generation anatomic total shoulder prostheses offer a wide range of adaptability (size, thickness, retroversion and offset of the humeral head, cervico-diaphyseal angle) in order to reproduce anatomy and biomechanics of the shoulder as normal as possible. The large variability of the implants may also induce malposition. Our goal was to analyse the consequences of a humeral head malposition, which is one of the most frequent placement errors. Material and Methods A 3D finite element model of the glenohumeral joint, including the rotator cuff muscles and the deltoid, was used with the Aequalis anatomic prosthesis. Active abduction was simulated. Three humeral head placements were compared : anatomic positioning (A), 5 mm inferior positioning (B), 5 mm superior positioning (C). The effect of humeral head malposition was evaluated through the following quantities : the range of motion free of impingements, the glenohumeral contact pattern, and the stress within the polyethylene and the cement. Results Inferior positioning (B) of the humeral head produced a superior impingement before 90° of abduction, an inferior eccentric contact point on the glenoid, and 165% increase of cement stress. Superior positioning (C) of the humeral head produced a postero-superior eccentric contact point on the glenoid, 300% increase of glenohumeral contact pressure, 450% increase of polyethylene stress, and 207% increase of cement stress. Conclusion Malposition of the humeral head of anatomic prostheses induces biomechanical consequences that may preclude the glenoid survival. Particular attention must be paid to reproduce the humeral anatomy as normal as possible.

Simulated joint and muscle forces in reversed and anatomic shoulder prostheses.

Reversed shoulder prostheses are increasingly being used for the treatment of glenohumeral arthropathy associated with a deficient rotator cuff. These non-anatomical implants attempt to balance the joint forces by means of a semi-constrained articular surface and a medialised centre of rotation. A finite element model was used to compare a reversed prosthesis with an anatomical implant. Active abduction was simulated from 0 degrees to 150 degrees of elevation. With the anatomical prosthesis, the joint force almost reached the equivalence of body weight. The joint force was half this for the reversed prosthesis. The direction of force was much more vertically aligned for the reverse prosthesis, in the first 90 degrees of abduction. With the reversed prosthesis, abduction was possible without rotator cuff muscles and required 20% less deltoid force to achieve it. This force analysis confirms the potential mechanical advantage of reversed prostheses when rotator cuff muscles are deficient.

Benefit of Cup Medialization in Total Hip Arthroplasty is Associated With Femoral Anatomy.

BACKGROUND: Medialization of the cup with a respective increase in femoral offset has been proposed in THA to increase abductor moment arms. Insofar as there are potential disadvantages to cup medialization, it is important to ascertain whether the purported biomechanical benefits of cup medialization are large enough to warrant the downsides; to date, studies regarding this question have disagreed. QUESTIONS/PURPOSES: The purpose of this study was to quantify the effect of cup medialization with a compensatory increase in femoral offset compared with anatomic reconstruction for patients undergoing THA. We tested the hypothesis that there is a (linear) correlation between preoperative anatomic parameters and muscle moment arm increase caused by cup medialization. METHODS: Fifteen patients undergoing THA were selected, covering a typical range of preoperative femoral offsets. For each patient, a finite element model was built based on a preoperative CT scan. The model included the pelvis, femur, gluteus minimus, medius, and maximus. Two reconstructions were compared: (1) anatomic position of the acetabular center of rotation, and (2) cup medialization compensated by an increase in the femoral offset. Passive abduction-adduction and flexion-extension were simulated in the range of normal gait. Muscle moment arms were evaluated and correlated to preoperative femoral offset, acetabular offset, height of the greater trochanter (relative to femoral center of rotation), and femoral antetorsion angle. RESULTS: The increase of muscle moment arms caused by cup medialization varied among patients. Muscle moment arms increase by 10% to 85% of the amount of cup medialization for abduction-adduction and from -35% (decrease) to 50% for flexion-extension. The change in moment arm was inversely correlated (R(2) = 0.588, p = 0.001) to femoral antetorsion (anteversion), such that patients with less femoral antetorsion gained more in terms of hip muscle moments. No linear correlation was observed between changes in moment arm and other preoperative parameters in this series. CONCLUSIONS: The benefit of cup medialization is variable and depends on the individual anatomy. CLINICAL RELEVANCE: Cup medialization with compensatory increase of the femoral offset may be particularly effective in patients with less femoral antetorsion. However, cup medialization must be balanced against its tradeoffs, including the additional loss of medial acetabular bone stock, and eventual proprioceptive implications of the nonanatomic center of rotation and perhaps joint reaction forces. Clinical studies should better determine the relevance of small changes of moment arms on function and joint reaction forces.

Effect of tibial tray inclination on the femoro-tibial contact pattern of a mobile total knee arthroplasty

Introduction: The posterior inclination of the tibial component is an important factor that can affect the success of total knee arthroplasty. It can reduce the posterior impingement and thus increase the range of flexion, but it may also induce instability in flexion, anterior impingement between the polyethylene of postero-stabilizing knee prosthesis, and anterior conflict with the cortical bone and the stem. Although the problem is identified, there is still a debate on the ideal inclination angle and the surgical technique to avoid an excessive posterior inclination. The aim of this study was to predict the effect of a posterior inclination of the tibial component on the contact pattern on the tibial insert, using a numerical musculoskeletal model of the knee joint. Methods: A 3D finite element model of the knee joint was developed to simulate an active and loaded squat movement after total knee arthroplasty. Flexion was actively controlled by the quadriceps muscle and muscle activations were estimated from EMG data and were synchronized by a feedback algorithm. Two inclinations of the tibial tray were considered: a posterior inclination of 0° or 10°. During the entire range of flexion, the following quantities were calculated: the tibiofemoral and patello-femoral contact force, and the contact pattern on polyethylene insert. The antero-posterior displacement of the contact pattern was also measured. Abaqus 6.7 was used for all analyses. Results: The tibio-femoral and patello-femoral contact forces increased during flexion and reached respectively 4 and 7 BW (bodyweight) at 90° of flexion. They were slightly affected by the inclination of the tibial tray. Without posterior inclination, the contact pattern on the tibial insert remained centered. The contact pressure was lower than 5 MPa below 60° of flexion, but exceeded 20 MPa at 90° of flexion. The posterior inclination displaced the contact point posteriorly by 2 to 4 mm. Conclusion: The inclination of the tibial tray displaced the contactpattern towards the posterior border of the tibial insert. However, even for 10° of inclination, the contact center remained far from the posterior border (12 mm). There was no instability predicted for this movement.

Early stage disc degeneration does not have an appreciable affect on stiffness and load transfer following vertebroplasty and kyphoplasty.

Vertebroplasty and kyphoplasty have been reported to alter the mechanical behavior of the treated and adjacent-level segments, and have been suggested to increase the risk for adjacent-level fractures. The intervertebral disc (IVD) plays an important role in the mechanical behavior of vertebral motion segments. Comparisons between normal and degenerative IVD motion segments following cement augmentation have yet to be reported. A microstructural finite element model of a degenerative IVD motion segment was constructed from micro-CT images. Microdamage within the vertebral body trabecular structure was used to simulate a slightly (I = 83.5% of intact stiffness), moderately (II = 57.8% of intact stiffness), and severely (III = 16.0% of intact stiffness) damaged motion segment. Six variable geometry single-segment cement repair strategies (models A-F) were studied at each damage level (I-III). IVD and bone stresses, and motion segment stiffness, were compared with the intact and baseline damage models (untreated), as well as, previous findings using normal IVD models with the same repair strategies. Overall, small differences were observed in motion segment stiffness and average stresses between the degenerative and normal disc repair models. We did however observe a reduction in endplate bulge and a redistribution in the microstructural tissue level stresses across both endplates and in the treated segment following early stage IVD degeneration. The cement augmentation strategy placing bone cement along the periphery of the vertebra (model E) proved to be the most advantageous in treating the degenerative IVD models by showing larger reductions in the average bone stresses (vertebral and endplate) as compared to the normal IVD models. Furthermore, only this repair strategy, and the complete cement fill strategy (model F), were able to restore the slightly damaged (I) motion segment stiffness above pre-damaged (intact) levels. Early stage IVD degeneration does not have an appreciable effect in motion segment stiffness and average stresses in the treated and adjacent-level segments following vertebroplasty and kyphoplasty. Placing bone cement in the periphery of the damaged vertebra in a degenerative IVD motion segment, minimizes load transfer, and may reduce the likelihood of adjacent-level fractures.