Comparison of three-dimensional lower extremity running kinematics of young adult and elderly runners

The objective of this study was to compare the three-dimensional lower extremity running kinematics of young adult runners and elderly runners. Seventeen elderly adults (age 67-73 years) and 17 young adults (age 26-36 years) ran at 3.1ms-1 on a treadmill while the movements of the lower extremity during the stance phase were recorded at 120Hz using three-dimensional video. The three-dimensional kinematics of the lower limb segments and of the ankle and knee joints were determined, and selected variables were calculated to describe the movement. Our results suggest that elderly runners have a different movement pattern of the lower extremity from that of young adults during the stance phase of running. Compared with the young adults, the elderly runners had a substantial decrease in stride length (1.97 vs. 2.23m; P=0.01), an increase in stride frequency (1.58 vs. 1.37Hz; P=0.002), less knee flexion/extension range of motion (26 vs. 33; P=0.002), less tibial internal/external rotation range of motion (9 vs. 12; P0.001), larger external rotation angle of the foot segment (toe-out angle) at the heel strike (-5.8 vs. -1.0; P=0.009), and greater asynchronies between the ankle and knee movements during running. These results may help to explain why elderly individuals could be more susceptible to running-related injuries.

Aerodynamic study of three-dimensional larynx models using finite element methods

The airflow velocities and pressures are calculated from a three-dimensional model of the human larynx by using the finite element method. The laryngeal airflow is assumed to be incompressible, isothermal, steady, and created by fixed pressure drops. The influence of different laryngeal profiles (convergent, parallel, and divergent), glottal area, and dimensions of false vocal folds in the airflow are investigated. The results indicate that vertical and horizontal phase differences in the laryngeal tissue movements are influenced by the nonlinear pressure distribution across the glottal channel, and the glottal entrance shape influences the air pressure distribution inside the glottis. Additionally, the false vocal folds increase the glottal duct pressure drop by creating a new constricted channel in the larynx, and alter the airflow vortexes formed after the true vocal folds. (C) 2007 Elsevier Ltd. All rights reserved.

Analyzing static three-dimensional elastic domains with a new infinite boundary element formulation

A new two-dimensionally mapped infinite boundary element (IBE) is presented. The formulation is based on a triangular boundary element (BE) with linear shape functions instead of the quadrilateral IBEs usually found in the literature. The infinite solids analyzed are assumed to be three-dimensional, linear-elastic and isotropic, and Kelvin fundamental solutions are employed. One advantage of the proposed formulation over quadratic or higher order elements is that no additional degrees of freedom are added to the original BE mesh by the presence of the IBEs. Thus, the IBEs allow the mesh to be reduced without compromising the accuracy of the result. Two examples are presented, in which the numerical results show good agreement with authors using quadrilateral IBEs and analytical solutions. (C) 2010 Elsevier Ltd. All rights reserved.

An alternative multi-region BEM technique for three-dimensional elastic problems

The main objective of this work is to present an alternative boundary element method (BEM) formulation for the static analysis of three-dimensional non-homogeneous isotropic solids. These problems can be solved using the classical boundary element formulation, analyzing each subregion separately and then joining them together by introducing equilibrium and displacements compatibility. Establishing relations between the displacement fundamental solutions of the different domains, the alternative technique proposed in this paper allows analyzing all the domains as one unique solid, not requiring equilibrium or compatibility equations. This formulation also leads to a smaller system of equations when compared to the usual subregion technique, and the results obtained are even more accurate. (C) 2008 Elsevier Ltd. All rights reserved.

GENERALIZED FINITE ELEMENT METHOD FOR NONLINEAR THREE-DIMENSIONAL ANALYSIS OF SOLIDS

The Generalized Finite Element Method (GFEM) is employed in this paper for the numerical analysis of three-dimensional solids tinder nonlinear behavior. A brief summary of the GFEM as well as a description of the formulation of the hexahedral element based oil the proposed enrichment strategy are initially presented. Next, in order to introduce the nonlinear analysis of solids, two constitutive models are briefly reviewed: Lemaitre`s model, in which damage and plasticity are coupled, and Mazars`s damage model suitable for concrete tinder increased loading. Both models are employed in the framework of a nonlocal approach to ensure solution objectivity. In the numerical analyses carried out, a selective enrichment of approximation at regions of concern in the domain (mainly those with high strain and damage gradients) is exploited. Such a possibility makes the three-dimensional analysis less expensive and practicable since re-meshing resources, characteristic of h-adaptivity, can be minimized. Moreover, a combination of three-dimensional analysis and the selective enrichment presents a valuable good tool for a better description of both damage and plastic strain scatterings.

A hybrid-mixed finite element formulation for the geometrically exact analysis of three-dimensional framed structures

This paper addresses the development of a hybrid-mixed finite element formulation for the quasi-static geometrically exact analysis of three-dimensional framed structures with linear elastic behavior. The formulation is based on a modified principle of stationary total complementary energy, involving, as independent variables, the generalized vectors of stress-resultants and displacements and, in addition, a set of Lagrange multipliers defined on the element boundaries. The finite element discretization scheme adopted within the framework of the proposed formulation leads to numerical solutions that strongly satisfy the equilibrium differential equations in the elements, as well as the equilibrium boundary conditions. This formulation consists, therefore, in a true equilibrium formulation for large displacements and rotations in space. Furthermore, this formulation is objective, as it ensures invariance of the strain measures under superposed rigid body rotations, and is not affected by the so-called shear-locking phenomenon. Also, the proposed formulation produces numerical solutions which are independent of the path of deformation. To validate and assess the accuracy of the proposed formulation, some benchmark problems are analyzed and their solutions compared with those obtained using the standard two-node displacement/ rotation-based formulation.

Hybrid and multi-field variational principles for geometrically exact three-dimensional beams

This paper addresses the development of several alternative novel hybrid/multi-field variational formulations of the geometrically exact three-dimensional elastostatic beam boundary-value problem. In the framework of the complementary energy-based formulations, a Legendre transformation is used to introduce the complementary energy density in the variational statements as a function of stresses only. The corresponding variational principles are shown to feature stationarity within the framework of the boundary-value problem. Both weak and linearized weak forms of the principles are presented. The main features of the principles are highlighted, giving special emphasis to their relationships from both theoretical and computational standpoints. (C) 2010 Elsevier Ltd. All rights reserved.

A transient three-dimensional heat transfer model of the human body

The objective of this work is to develop an improved model of the human thermal system. The features included are important to solve real problems: 3D heat conduction, the use of elliptical cylinders to adequately approximate body geometry, the careful representation of tissues and important organs, and the flexibility of the computational implementation. Focus is on the passive system, which is composed by 15 cylindrical elements and it includes heat transfer between large arteries and veins. The results of thermal neutrality and transient simulations are in excellent agreement with experimental data, indicating that the model represents adequately the behavior of the human thermal system. (C) 2009 Elsevier Ltd. All rights reserved.

Swarm Intelligence applied in synthesis of hunting strategies in a three-dimensional environment

Systems of distributed artificial intelligence can be powerful tools in a wide variety of practical applications. Its most surprising characteristic, the emergent behavior, is also the most answerable for the difficulty in. projecting these systems. This work proposes a tool capable to beget individual strategies for the elements of a multi-agent system and thereof providing to the group means on obtaining wanted results, working in a coordinated and cooperative manner as well. As an application example, a problem was taken as a basis where a predators` group must catch a prey in a three-dimensional continuous ambient. A synthesis of system strategies was implemented of which internal mechanism involves the integration between simulators by Particle Swarm Optimization algorithm (PSO), a Swarm Intelligence technique. The system had been tested in several simulation settings and it was capable to synthesize automatically successful hunting strategies, substantiating that the developed tool can provide, as long as it works with well-elaborated patterns, satisfactory solutions for problems of complex nature, of difficult resolution starting from analytical approaches. (c) 2007 Elsevier Ltd. All rights reserved.

Three-dimensional analysis of fracture, corrosion and wear surfaces

A brief look at the history of fractography has shown a recent trend in the quantification of topographic parameters through the use of three-dimensional reconstruction techniques, which associate SEM stereoscopy and stereophotogrammetry software, allowing the calculation of the elevation measurement at numerous points of the topography due to the parallax that takes place during the tilting of the sample along the microscope eucentric plane. Several investigators have used reconstruction techniques to correlate some fractographic parameters, such as fractal dimension and fractured to projected area ratio, to the mechanical properties of materials, such as fracture toughness and tensile strength. So far, the search for a clear relationship between the fracture topography and mechanical properties has provided ambiguous results. The present work applied a surface metrology software to reconstruct three-dimensionally fracture surfaces (transgranular cleavage, intergranular and dimple fracture), corrosion pits and tribo-surfaces in order to explore the potential of this stereophotogrammetry technique. The existence of a variation in the calculated topographic parameters with the conditions of SEM image acquisition reinforces the importance of both good image acquisition and accurate calibration methods in order to validate this 3D reconstruction technique in metrological terms. Preliminary results did not indicate the existence of a clear relationship between either the true to project area ratio and CVN absorbed energy or the fractal dimension and CVN absorbed energy. It is likely that each fracture mechanism presents a proper relationship between the fractographic parameters and mechanical properties. (C) 2009 Elsevier Ltd. All rights reserved.

Three-Dimensional Stress Distribution in the Human Periodontal Ligament in Masticatory, Parafunctional, and Trauma Loads: Finite Element Analysis

Background: The presence of the periodontal ligament (PDL) makes it possible to absorb and distribute loads produced during masticatory function and other tooth contacts into the alveolar process via the alveolar bone proper. However, several factors affect the integrity of periodontal structures causing the destruction of the connective matrix and cells, the loss of fibrous attachment, and the resorption of alveolar bone. Methods: The purpose of this study was to evaluate the stress distribution by finite element analysis in a PDL in three-dimensional models of the upper central incisor under three different load conditions: 100 N occlusal loading at 45 degrees (model 1: masticatory load); 500 N at the incisal edge at 45 degrees (model 2: parafunctional habit); and 800 N at the buccal surface at 90 degrees (model 3: trauma case). The models were built from computed tomography scans. Results: The stress distribution was quite different among the models. The most significant values (harmful) of tensile and compressive stresses were observed in models 2 and 3, with similarly distinct patterns of stress distributions along the PDL. Tensile stresses were observed along the internal and external aspects of the PDL, mostly at the cervical and middle thirds. Conclusions: The stress generation in these models may affect the integrity of periodontal structures. A better understanding of the biomechanical behavior of the PDL under physiologic and traumatic loading conditions might enhance the understanding of the biologic reaction of the PDL in health and disease. J Periodontol 2009;80:1859-1867.

Fast Structure-Based Assignment of 15N HSQC Spectra of Selectively 15N-Labeled Paramagnetic Proteins

A novel strategy for fast NMR resonance assignment of N-15 HSQC spectra of proteins is presented. It requires the structure coordinates of the protein, a paramagnetic center, and one or more residue-selectively N-15-labeled samples. Comparison of sensitive undecoupled N-15 HSQC spectra recorded of paramagnetic and diamagnetic samples yields data for every cross-peak on pseudocontact shift, paramagnetic relaxation enhancement, cross-correlation between Curie-spin and dipole-dipole relaxation, and residual dipolar coupling. Comparison of these four different paramagnetic quantities with predictions from the three-dimensional structure simultaneously yields the resonance assignment and the anisotropy of the susceptibility tensor of the paramagnetic center. The method is demonstrated with the 30 kDa complex between the N-terminal domain of the epsilon subunit and the theta subunit of Escherichia Coll DNA polymerase III. The program PLATYPUS was developed to perform the assignment, provide a measure of reliability of the assignment, and determine the susceptibility tensor anisotropy.

Solution structure and proposed binding mechanism of a novel potassium channel toxin kappa-conotoxin PVIIA

Background: kappa-PVIIA is a 27-residue polypeptide isolated from the venom of Conus purpurascens and is the first member of a new class of conotoxins that block potassium channels. By comparison to other ion channels of eukaryotic cell membranes, voltage-sensitive potassium channels are relatively simple and methodology has been developed for mapping their interactions with small-peptide toxins, PVIIA, therefore, is a valuable new probe of potassium channel structure. This study of the solution structure and mode of channel binding of PVIIA forms the basis for mapping the interacting residues at the conotoxin-ion channel interface. Results: The three-dimensional structure of PVIIA resembles the triple-stranded beta sheet/cystine-knot motif formed by a number of toxic and inhibitory peptides. Subtle structural differences, predominantly in loops 2 and 4, are observed between PVIIA and other conotoxins with similar structural frameworks, however. Electrophysiological binding data suggest that PVIIA blocks channel currents by binding in a voltage-sensitive manner to the external vestibule and occluding the pore, Comparison of the electrostatic surface of PVIIA with that of the well-characterised potassium channel blocker charybdotoxin suggests a likely binding orientation for PVIIA, Conclusions: Although the structure of PVIIA is considerably different to that of the alpha K scorpion toxins, it has a similar mechanism of channel blockade. On the basis of a comparison of the structures of PVIIA and charybdotoxin, we suggest that Lys19 of PVIIA is the residue which is responsible for physically occluding the pore of the potassium channel.

Changes in three-dimensional architecture of microfilaments in cultured vascular smooth muscle cells during phenotypic modulation

To investigate changes in the three-dimensional microfilament architecture of vascular smooth muscle cells (SMC) during the process of phenotypic modulation, rabbit aortic SMCs cultured under different conditions and at different time points were either labelled with fluorescein-conjugated probes to cytoskeletal and contractile proteins for observation by confocal laser scanning microscopy, or extracted with Triton X-100 for scanning electron microscopy. Densely seeded SMCs in primary culture, which maintain a contractile phenotype, display prominent linear myofilament bundles (stress fibres) that are present throughout the cytoplasm with alpha-actin filaments predominant in the central part and beta-actin filaments in the periphery of the cell. Intermediate filaments form a meshed network interconnecting the stress fibres and linking directly to the nucleus. Moderately and sparsely seeded SMCs, which modulate toward the synthetic phenotype during the first 5 days of culture, undergo a gradual redistribution of intermediate filaments from the perinuclear region toward the peripheral cytoplasm and a partial disassembly of stress fibres in the central part of the upper cortex of the cytoplasm, with an obvious decrease in alpha-actin and myosin staining. These changes are reversed in moderately seeded SMCs by day 8 of culture when they have reached confluence. The results reveal two changes in microfilament architecture in SMCs as they undergo a change in phenotype: the redistribution of intermediate filaments probably due to an increase in synthetic organelles in the perinuclear area, and the partial disassembly of stress fibres which may reflect a degradation of contractile components.

Evolution of three-dimensional folds for a non-Newtonian plate in a viscous medium

We derive analytical solutions for the three-dimensional time-dependent buckling of a non-Newtonian viscous plate in a less viscous medium. For the plate we assume a power-law rheology. The principal, axes of the stretching D-ij in the homogeneously deformed ground state are parallel and orthogonal to the bounding surfaces of the plate in the flat state. In the model formulation the action of the less viscous medium is replaced by equivalent reaction forces. The reaction forces are assumed to be parallel to the normal vector of the deformed plate surfaces. As a consequence, the buckling process is driven by the differences between the in-plane stresses and out of plane stress, and not by the in-plane stresses alone as assumed in previous models. The governing differential equation is essentially an orthotropic plate equation for rate dependent material, under biaxial pre-stress, supported by a viscous medium. The differential problem is solved by means of Fourier transformation and largest growth coefficients and corresponding wavenumbers are evaluated. We discuss in detail fold evolutions for isotropic in-plane stretching (D-11 = D-22), uniaxial plane straining (D-22 = 0) and in-plane flattening (D-11 = -2D(22)). Three-dimensional plots illustrate the stages of fold evolution for random initial perturbations or initial embryonic folds with axes non-parallel to the maximum compression axis. For all situations, one dominant set of folds develops normal to D-11, although the dominant wavelength differs from the Biot dominant wavelength except when the plate has a purely Newtonian viscosity. However, in the direction parallel to D-22, there exist infinitely many modes in the vicinity of the dominant wavelength which grow only marginally slower than the one corresponding to the dominant wavelength. This means that, except for very special initial conditions, the appearance of a three-dimensional fold will always be governed by at least two wavelengths. The wavelength in the direction parallel to D-11 is the dominant wavelength, and the wavelength(s) in the direction parallel to D-22 is determined essentially by the statistics of the initial state. A comparable sensitivity to the initial geometry does not exist in the classic two-dimensional folding models. In conformity with tradition we have applied Kirchhoff's hypothesis to constrain the cross-sectional rotations of the plate. We investigate the validity of this hypothesis within the framework of Reissner's plate theory. We also include a discussion of the effects of adding elasticity into the constitutive relations and show that there exist critical ratios of the relaxation times of the plate and the embedding medium for which two dominant wavelengths develop, one at ca. 2.5 of the classical Biot dominant wavelength and the other at ca. 0.45 of this wavelength. We propose that herein lies the origin of parasitic folds well known in natural examples.