Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Matrix stiffness”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4Linked to original sources

Investigation of scattering of elastic waves by cylinders in 1-3 piezocomposites.

The scattering behavior of P-waves in piezoelectric composites with 1-3 connectivity is studied. The method of wave function expansion is adopted for the theoretical derivations. Analytical expressions are obtained for the distributions of mechanical displacement in z-direction along the circumferences of piezoelectric cylinders. These solutions are used to study the influence of each element of the stiffness matrix and the piezoelectric matrix on the various resonant modes of vibration. Numerical results obtained indicate that perturbations of the elements c44 and e15 significantly affect resonant frequencies and amplitudes, perturbations of c11 and c12 have pronounced effects on resonant modes of high frequencies also. However, the resonant modes are not so sensitive to the perturbations of c13, e31 and e33. The dynamic characteristics of 1-3 connectivity piezoelectric composites exposed here are meaningful for the design and manufacture of sensor/actuator elements by this kind of composites as well as the on-line health monitoring of the mechanical properties variations of the composites itself.

Ceramics↗

A study on the dispersions of piezoelectric plate with laser ultrasound measurement and theoretical modeling.

This research is focused in the measurement and modeling for the dispersion relations of Lamb waves propagating in LiNbO(3) and LiTaO(3) plates. A theoretical model based on a stiffness matrix method with recursion computation algorism is used to provide numerical calculations for the dispersion relations. Also, a dry, noncontact laser ultrasound technique is used to measured dispersion multi-mode dispersion relations. For all the experiments, the measured dispersion curves show good agreement with the theoretical calculation, indicating the reliabilities in the measurement and modeling.

Journal Article↗

Elastic properties of cancellous bone derived from finite element models of parameterized microstructure cells.

Evaluation of constitutive properties of cancellous bone and their relationships to microstructural parameters is a crucial issue in analysis of stresses and strains in bone tissues and simulation of their remodelling. Known limitations of experimental methods as well as of the micro-FE techniques make the analysis and homogenization of 'equivalent' trabecular microstructures an advantageous tool for this task. In this study, parameterized orthotropic constitutive models of cancellous bone are derived from finite element analysis of repeatable microstructure cells. Two cell types are analysed: cube- and prism-based. The models are fully three-dimensional, have realistic curvilinear shapes and are parameterized with three shape parameters. Variation of the parameters allows to imitate most of the typical microstructure patterns observed in real bones, along with variety of intermediate geometries. Finite element models of cells are generated by a special-purpose structured mesh generator for any arbitrary set of shape parameter values. Six static numerical tests are performed for an exhaustive number of parameter value sets (microstructure instances). Multi-point boundary conditions imposed on the models ensure mutual fitting of deformed neighbouring cells. Values of computed stresses allow to determine all coefficients of elastic orthotropic stiffness matrix. Results have a form of tabularized functions of elastic constants versus the shape parameters. Comparison of the results with micro-FE data obtained for a large set of cancellous bone specimens proves a good agreement, though evidently better in the case of the prism-based cell model.

Animals↗

Anisotropic properties of human tibial cortical bone as measured by nanoindentation.

The purpose of this study was to investigate the effects of elastic anisotropy on nanoindentation measurements in human tibial cortical bone. Nanoindentation was conducted in 12 different directions in three principal planes for both osteonic and interstitial lamellae. The experimental indentation modulus was found to vary with indentation direction and showed obvious anisotropy (one-way analysis of variance test, P < 0.0001). Because experimental indentation modulus in a specific direction is determined by all of the elastic constants of cortical bone, a complex theoretical model is required to analyze the experimental results. A recently developed analysis of indentation for the properties of anisotropic materials was used to quantitatively predict indentation modulus by using the stiffness matrix of human tibial cortical bone, which was obtained from previous ultrasound studies. After allowing for the effects of specimen preparation (dehydrated specimens in nanoindentation tests vs. moist specimens in ultrasound tests) and the structural properties of bone (different microcomponents with different mechanical properties), there were no statistically significant differences between the corrected experimental indentation modulus (Mexp) values and corresponding predicted indentation modulus (Mpre) values (two-tailed unpaired t-test, P > 0.5). The variation of Mpre values was found to exhibit the same trends as the corrected Mexp data. These results show that the effects of anisotropy on nanoindentation measurements can be quantitatively evaluated.

Anisotropy↗

A new variational shape-from-orientation approach to correcting intensity inhomogeneities in magnetic resonance images.

A new intensity inhomogeneity correction algorithm based on a variational shape-from-orientation formulation is presented. Unlike most previous methods, the proposed algorithm is fully automatic, widely applicable and very efficient. Since no prior classification knowledge about the image is assumed in the proposed algorithm, it can be applied to correct intensity inhomogeneities for a wide variety of medical images. In this paper, a finite-element method is used to model the smooth bias-field function. Orientation constraints for the bias-field function are computed at the nodal locations of the regular discretization grid away from the boundary between different class regions. The selection of reliable orientation constraints is facilitated by the goodness of fit of a first-order polynomial model to the neighborhood of each nodal location. The automatically selected orientation constraints are integrated in a regularization framework, which leads to minimization of a convex and quadratic energy function. This energy minimization is accomplished by solving a linear system with a large, sparse, symmetric and positive semi-definite stiffness matrix. We employ an adaptive preconditioned conjugate-gradient algorithm to solve the linear system very efficiently. Experimental results on a variety of magnetic resonance images are given to demonstrate the effectiveness and efficiency of the proposed algorithm.

Algorithms↗

Orthotropic properties of cancellous bone modelled as parameterized cellular material.

Analysis of stresses and strains in bone tissues and simulation of their adaptive remodelling require exhaustive information about distribution of constitutive properties of cancellous bone and their relationships to microstructural parameters. Homogenization of "equivalent" trabecular microstructures appears to be an advantageous tool for this task. In this study, parameterized orthotropic constitutive models of cancellous bone are derived from finite element analysis of repeatable microstructure cells. The models, based on a space-filling dodecahedron, are fully three-dimensional and are parameterized with four shape parameters. Variation of the parameters allows to imitate most of typical microstructure patterns observed in real bones, along with a variety of intermediate geometries. Finite element models of cells are generated by a special-purpose structured mesh generator for any arbitrary set of shape parameter values. Static numerical tests are performed for an exhaustive number of parameter value sets (microstructure instances). Coefficients of elastic orthotropic stiffness matrix are determined as tabularized functions of elastic constants versus the shape parameters. Additionally, they are correlated to apparent density and principal fabric tensor values. Comparison of the results with micro-FE data obtained for a large set of cancellous bone specimens proves a good agreement.

Animals↗

Ultrasound critical-angle reflectometry (UCR): a new modality for functional elastometric imaging.

This paper discusses the measurement of velocity in a solid based on the analysis of the amplitude and phase of ultrasound waves reflected by a solid, a technique called ultrasound critical-angle reflectometry (UCR). To this end, the complete formulation of ultrasound wave reflection and refraction from a liquid-solid interface is described. Differences between this formulation and previously published ones are briefly discussed. Based on this analysis it is in particular possible to measure by this technique not only pressure but also, for the first time in such studies, shear wave velocities, an experimentally confirmed result. The measurement of the complete stiffness matrix of a transversely isotropic solid, specifically cortical bone, by applying UCR elastometry to any point on the solid's surface is demonstrated. Finally this method is extended to functional elastometric imaging. The techniques presented in this paper offer new opportunities for applications of UCR imaging to the assessment of bone metabolism, formation and disease and also the analysis of composite materials in general.

Anisotropy↗

A second-order finite element algorithm for solving the three-dimensional EEG forward problem.

A finite element algorithm has been developed to solve the electroencephalogram (EEG) forward problem. A new computationally efficient approach to calculate the stiffness matrix of second-order tetrahedral elements has been developed for second-order tetrahedral finite element models. The present algorithm has been evaluated by means of computer simulations, by comparing with analytic solutions in a multi-spheres concentric head model. The developed finite element method (FEM) algorithm has also been applied to address questions of interest in the EEG forward problem. The present simulation study indicates that the second-order FEM provides substantially enhanced numerical accuracy and computational efficiency, as compared with the first-order FEM for comparable numbers of tetrahedral elements. The anisotropic conductivity distribution of the head tissue can be taken into account in the present FEM algorithm. The effects of dipole eccentricity, size of finite elements and local mesh refinement on solution accuracy are also addressed in the present simulation study.

Algorithms↗

Indeterminacy and the onset of motion in a simple granular packing.

We examine the relatively simple problem of a disk placed in a symmetric V-shaped channel, and subjected to gravity and a torque. We obtain analytic predictions of the contact forces and disk motion using two different models. In the first model, the disk is assumed to be perfectly rigid, leading to force indeterminacy. In the second model, the disk is assumed to interact with the walls of the channel via linear springs, leading to a unique solution for the contact forces. The results of these two models are compared. It is shown that there are two possible ways motion can occur--through the appearance of a null eigenvector of the stiffness matrix, or through an instability. When motion occurs through an instability, the first model cannot predict when the disk will rotate; it is necessary to know the undetermined forces in order to predict the motion of the disk. It is also shown how indeterminacy in the first model is linked to memory in the second. The analytical results are also compared with numerical simulations using two different methods, each related to one of the models.

Journal Article↗

Finite element modeling of transient ultrasonic waves in linear viscoelastic media.

Linear viscoelasticity offers a minimal framework within which to construct a causal model for wave propagation in absorptive media. Viscoelastic media are often described as media with 'fading memory,' that is, the present state of stress is dependent on the present strain and the complete time history of strain convolved with appropriate time-dependent shear and bulk stress relaxation moduli. An axisymmetric, displacement-based finite element method for modeling pulsed ultrasonic waves in linear, homogeneous, and isotropic (LHI) viscoelastic media is developed that does not require storage of the complete time history of displacement at every node. This is accomplished by modeling stress relaxation moduli as discrete or continuous spectra of decaying exponentials and relaxation times. Details of the construction and computation of the time-dependent stiffness matrix are presented. As an application of the finite element method, a finite number of exponentials (amplitudes and relaxation times) are employed to represent a typical model for a continuous relaxation spectrum. It is demonstrated that a small number of discrete exponentials are required to model ultrasonic wave propagation of a typical band-limited pulse in a model material accurately. Previous work has shown this model to be consistent with other analytic models for wave propagation in viscoelastic media [1].

Journal Article↗

Interactive deformation of soft tissues with haptic feedback for medical learning.

An effective deformable model based on a successive force propagation process is proposed. It avoids the laborious stiffness matrix formulation and is scalable simply by controlling the penetration depth. Mechanical tests are performed to evaluate its feasibility for modeling real tissues. An interactive system is developed using a commercial haptic device.

Computer Simulation↗

Time-resolved line focus acoustic microscopy of layered anisotropic media: application to composites.

This paper presents theoretical and experimental studies of the time-domain response of line focus acoustic microscopy from a layered anisotropic medium. A method for elastic constant reconstruction from acoustic microscopy signatures also is presented. The microscopy response is complicated by multiple reflections in the layers and by the anisotropic nature of the material. The model is based on a new, stable recursive stiffness matrix algorithm developed for a multilayered anisotropic medium, which is applied to the interpretation of the time-resolved acoustic microscopy signature. Specific examples are given for unidirectional and multidirectional graphite epoxy composites. It is shown that the fluid load has a significant effect on the leaky surface waves in these composites, increasing surface wave speed above that for the slow transverse wave. This results in its absence from the microscopy signature of the surface wave. The theoretical results are compared with experiments carried out using a line focus PVDF transducer developed at National Institute of Standards and Technology (NIST). Time-resolved acoustic microscopy has been applied to the determination of elastic constants of a unidirectional composite or of one lamina in a cross-ply composite. The lateral waves and multiple reflections of bulk waves appearing in the microscopy signatures are used for the elastic properties reconstruction. The reconstruction results are compared to data obtained by the self-reference double-through-transmission ultrasonic bulk wave method.

Journal Article↗

A bidomain model based BEM-FEM coupling formulation for anisotropic cardiac tissue.

A hybrid boundary element method (BEM)/finite element method (FEM) approach is proposed in order to properly consider the anisotropic properties of the cardiac muscle in the magneto- and electrocardiographic forward problem. Within the anisotropic myocardium a bidomain model based FEM formulation is applied. In the surrounding isotropic volume conductor the BEM is adopted. Coupling is enabled by requesting continuity of the electric potential and the normal of the current density across the boundary of the heart. Here, the BEM part is coupled as an equivalent finite element to the finite element stiffness matrix, thus preserving in part its sparse property. First, continuous convergence of the coupling scheme is shown for a spherical model comparing the computed results to an analytic reference solution. Then, the method is extended to the depolarization phase in a fibrous model of a dog ventricle. A precomputed activation sequence obtained using a fine mesh of the heart was downsampled and used to calculate body surface potentials and extracorporal magnetic fields considering the anisotropic bidomain conductivities. Results are compared to those obtained by neglecting in part or totally (oblique or uniform dipole layer model) anisotropic properties. The relatively large errors computed indicate that the cardiac muscle is one of the major torso inhomogeneities.

Animals↗

The influence of the fixed negative charges on mechanical and electrical behaviors of articular cartilage under unconfined compression.

Unconfined compression test has been frequently used to study the mechanical behaviors of articular cartilage, both theoretically and experimentally. It has also been used in explant and gel-cell-complex studies in tissue engineering. In biphasic and poroelastic theories, the effect of charges fixed on the proteoglycan macromolecules in articular cartilage is embodied in the apparent compressive Young's modulus and the apparent Poisson's ratio of the tissue, and the fluid pressure is considered to be the portion above the osmotic pressure. In order to understand how proteoglycan fixed charges might affect the mechanical behaviors of articular cartilage, and in order to predict the osmotic pressure and electric fields inside the tissue in this experimental configuration, it is necessary to use a model that explicitly takes into account the charged nature of the tissue and the flow of ions within its porous interstices. In this paper, we used a finite element model based on the triphasic theory to study how fixed charges in the porous-permeable soft tissue can modulate its mechanical and electrochemical responses under a step displacement in unconfined compression. The results from finite element calculations showed that: 1) A charged tissue always supports a larger load than an uncharged tissue of the same intrinsic elastic moduli. 2) The apparent Young's modulus (the ratio of the equilibrium axial stress to the axial strain) is always greater than the intrinsic Young's modulus of an uncharged tissue. 3) The apparent Poisson's ratio (the negative ratio of the lateral strain to the axial strain) is always larger than the intrinsic Poisson's ratio of an uncharged tissue. 4) Load support derives from three sources: intrinsic matrix stiffness, hydraulic pressure and osmotic pressure. Under the unconfined compression, the Donnan osmotic pressure can constitute between 13%-22% of the total load support at equilibrium. 5) During the stress-relaxation process following the initial instant of loading, the diffusion potential (due to the gradient of the fixed charge density and the associated gradient of ion concentrations) and the streaming potential (due to fluid convection) compete against each other. Within the physiological range of material parameters, the polarity of the electric potential depends on both the mechanical properties and the fixed charge density (FCD) of the tissue. For softer tissues, the diffusion effects dominate the electromechanical response, while for stiffer tissues, the streaming potential dominates this response. 6) Fixed charges do not affect the instantaneous strain field relative to the initial equilibrium state. However, there is a sudden increase in the fluid pressure above the initial equilibrium osmotic pressure. These new findings are relevant and necessary for the understanding of cartilage mechanics, cartilage biosynthesis, electromechanical signal transduction by chondrocytes, and tissue engineering.

Animals↗

A finite element based method to determine the properties of planar soft tissue.

A finite element based method to determine the incremental elastic material properties of planar membranes was developed and evaluated. The method is applicable to tissues that exhibit inhomogeneity, geometric and material nonlinearity, and anisotropy. Markers are placed on the tissue to form a four-node quadrilateral element. The specimen is loaded to an initial reference state, then three incremental loading sets are applied and the nodal displacements recorded. One of these loadings must include shear. These data are used to solve an over-determined system of equations for the tangent stiffness matrix. The method was first verified using analytical data. Next, data obtained from a latex rubber sheet were used to evaluate experimental procedures. Finally, experiments conducted on preconditioned rat skin revealed nonlinear orthotropic behavior. The vector norm comparing the applied and calculated nodal force vectors was used to evaluate the accuracy of the solutions.

Algorithms↗

Incorporation of spinal flexibility measurements into finite element analysis.

This technical note demonstrates two methods of incorporating the experimental stiffness of spinal motion segments into a finite element analysis of the spine. The first method is to incorporate the experimental data directly as a stiffness matrix. The second method approximates the experimental data as a beam element.

Elasticity↗

Kinematic analysis and simulation of vertebral motion under static load-part II: simulation study.

Based on the earlier work by the author on the development of a theoretical motion simulation model of an open-loop, kinematic chain, the present paper examines the motion characteristics of a human spine motion-segment. Based on the experimental data collected from the human spine cadaver specimens of lumbar spine, a static simulation model of a lumbar segment is proposed. The equivalent motion characteristic of the intervertebral joint is described by a spherical pair which is located for an intervertebral joint by fitting in a least-square sense a cone to a set of motion data describing the axode characteristics of the joint. The intervertebral disk and the ligaments are represented by an equivalent stiffness matrix whose elements are assumed to be constant over the entire range of mobility. The proposed model showed a good agreement with the observed motion data of an intervertebral joint.

Humans↗

Finite element prediction of wave motion in structural waveguides.

A method is presented by which the wavenumbers for a one-dimensional waveguide can be predicted from a finite element (FE) model. The method involves postprocessing a conventional, but low order, FE model, the mass and stiffness matrices of which are typically found using a conventional FE package. This is in contrast to the most popular previous waveguide/FE approach, sometimes termed the spectral finite element approach, which requires new spectral element matrices to be developed. In the approach described here, a section of the waveguide is modeled using conventional FE software and the dynamic stiffness matrix formed. A periodicity condition is applied, the wavenumbers following from the eigensolution of the resulting transfer matrix. The method is described, estimation of wavenumbers, energy, and group velocity discussed, and numerical examples presented. These concern wave propagation in a beam and a simply supported plate strip, for which analytical solutions exist, and the more complex case of a viscoelastic laminate, which involves postprocessing an ANSYS FE model. The method is seen to yield accurate results for the wavenumbers and group velocities of both propagating and evanescent waves.

Acoustics↗