A method for planar biaxial mechanical testing that includes in-plane shear.

A limitation in virtually all planar biaxial studies of soft tissues has been the inability to include the effects of in-plane shear. This is due to the inability of current mechanical testing devices to induce a state of in-plane shear, due to the added cost and complexity. In the current study, a straightforward method is presented for planar biaxial testing that induces a combined state of in-plane shear and normal strains. The method relies on rotation of the test specimen's material axes with respect to the device axes and on rotating carriages to allow the specimen to undergo in-plane shear freely. To demonstrate the method, five glutaraldehyde treated bovine pericardium specimens were prepared with their preferred fiber directions (defining the material axes) oriented at 45 deg to the device axes to induce a maximum shear state. The test protocol included a wide range of biaxial strain states, and the resulting biaxial data re-expressed in material axes coordinate system. The resulting biaxial data was then fit to the following strain energy function W: [equation: see text] where E'ij is the Green's strain tensor in the material axes coordinate system and c and Ai are constants. While W was able to fit the data very well, the constants A5 and A6 were found not to contribute significantly to the fit and were considered unnecessary to model the shear strain response. In conclusion, while not able to control the amount of shear strain independently or induce a state of pure shear, the method presented readily produces a state of simultaneous in-plane shear and normal strains. Further, the method is very general and can be applied to any anisotropic planar tissue that has identifiable material axes.     

Fertile and sterile frond segments of the lyginopterid seed fern Feraxotheca

Fertile and sterile frond segments of the lyginopterid seed fern Feraxotheca culcitaus are described. In both the fertile and sterile specimen three orders of frond axes are borne alternately and in the same plane. Antepenultimate and penultimate frond axes are characterized by a C-shaped vascular bundle and numerous canals containing a golden-colored material in the fossilized condition. Synangia with only four sporangia demonstrate radial symmetry; bilateral symmetry is present in larger synangia. Pinnules possess dichotomous venation, with a single vein entering each lobe of the deeply incised lamina. A single canal occurs below each vein, and vein ends have an amplified vascular bundle. Morphologically Feraxotheca compares most closely with the compression genus Crossotheca. The morphology of the vascular bundle in the largest pinna is similar to a Lower Carboniferous specimen of Rhodea.     

Determination of 3D spinal kinematics without defining a local vertebral coordinate system

In this paper a method is presented to calculate Euler's angles of rotation of a body segment during locomotion without a priori defining the location of the center of rotation, and without defining a local vertebral coordinate system. The method was applied to in vivo spinal kinematics. In this method, the orientation of each segment is identified by a set of three markers. The orientation of the axes of rotation is calculated based on the average position of the markers during one stride cycle. Some restrictions and assumptions should be made. The approach is viable only when the average orientation of the anatomical axes of rotation of each spinal segment during a stride cycle coincides with the three axes of the laboratory coordinate system. Furthermore, the rotations should be symmetrical with respect to both sides of the plane of symmetry of the spinal segment, and the subject should move parallel to one axis of the laboratory coordinate system. Since in experimental conditions these assumptions will only be met approximately, errors will be introduced in the calculated angles of rotation. The magnitude of the introduced errors was investigated in a computer simulation experiment. Since the maximal errors did not exceed 0.7° in a range of misalignments up to 10° between the two coordinate systems, the approach proved to be a valid method for the estimation of spinal kinematics.     

A scapular coordinate frame for clinical and kinematic analyses

The aim of this study was to define a body-fixed coordinate frame for the scapula that minimises axes variability and is closely related to the clinical frame of reference. Medical images of 21 scapulae were used to quantify 14 different axes from identifiable landmarks. The plane of the blade of the scapula was defined. The orientations of the quantified axes were calculated. The angular relationships between axes were quantified and applied to grade the sensitivity of each axis to inter-scapular variations in the others. The volume of data required to define an axis was noted for its dependency on pathology and the three criteria were weighted according to relative importance. The two axes with the highest weighting were applied to define a body-fixed Cartesian coordinate frame for the scapula. A least square medio-lateral line through the centre of the spine root was the most optimal axis. The plane formed by the spine root line and a least square line through the centre of the lateral border ridge was the most optimal scapular plane. This body-fixed Cartesian coordinate frame is closely aligned to the cardinal planes in the anatomical position and thus is a clinically applicable, specimen invariant coordinate frame that can be used in patient-specific kinematics modelling.     

Effects of damage on the orthotropic material symmetry of bovine tibial trabecular bone

The macroscopic mechanical properties of trabecular bone can be predicted by its architecture using theoretical relationships between the elastic and architectural properties. Microdamage caused by overloading or fatigue decreases the apparent elastic moduli of trabecular bone requiring these relationships to be modified to predict the damaged elastic properties. In the case of isotropic damage, the apparent level elastic properties could be determined by multiplying all of the elastic constants by a single scalar factor. If the damage is anisotropic, the elastic constants may change by differing factors and the material coordinate system could become misaligned with the fabric coordinate system. High-resolution finite element models were used to simulate damage overloading on seven trabecular bone specimens subjected to pure shear strain in two planes. Comparison of the apparent elastic moduli of the specimens before and after damage showed that the reduction of the elastic moduli was anisotropic. This suggests that the microdamage within the specimens was inhomogeneous. However, after damage the specimens exhibited nearly orthotropic material symmetry as they did before damage. Changes in the orientation of the orthotropic material coordinate system were also small and occurred primarily in the transverse plane. Thus, while damage in trabecular bone is anisotropic, the material coordinate system remains aligned with the fabric tensor.     

Bilateral Symmetry Detection on the Basis of Scale Invariant Feature Transform

The automatic detection of bilateral symmetry is a challenging task in computer vision and pattern recognition. This paper presents an approach for the detection of bilateral symmetry in digital single object images. Our method relies on the extraction of Scale Invariant Feature Transform (SIFT) based feature points, which serves as the basis for the ascertainment of the centroid of the object; the latter being the origin under the Cartesian coordinate system to be converted to the polar coordinate system in order to facilitate the selection symmetric coordinate pairs. This is followed by comparing the gradient magnitude and orientation of the corresponding points to evaluate the amount of symmetry exhibited by each pair of points. The experimental results show that our approach draw the symmetry line accurately, provided that the observed centroid point is true.     

A kinematic model of the human ankle

The spatial gross motion of the foot with respect to the shank is modelled as rotations about two fixed ankle axes: the upper ankle rotation axis (plantarflexion/dorsiflexion) and the subtalar rotation axis (inversion/eversion). The positions of the axes are determined by externally visible bony landmarks of the lower leg and are measured for a living subject. The model input data are the plantarflexion/dorsiflexion and inversion/eversion rotation angles; the model output is a 4 × 4 transformation matrix which quantitatively describes the relative position of a foot coordinate system with respect to a shank coordinate system.     

A continuous wave technique for the measurement of the elastic properties of cortical bone

A continuous wave technique is described for measuring the nine independent orthotropic elastic coefficients from a single cubic specimen. The side dimensions of this cubic specimen are on the order of 5 mm. Because of the small size of the specimen, the spatial resolution of material inhomogeneity using this technique is quite good. Although it is possible to apply this technique to any elastic material such as woods or metals, the elastic properties of human and canine cortical femora are presented here. The orthotropic elastic coefficients and the variation of these coefficients are presented as a function of anatomical position.     

The sophisticated visual system of a tiny Cambrian crustacean: analysis of a stalked fossil compound eye

Fossilized compound eyes from the Cambrian, isolated and three-dimensionally preserved, provide remarkable insights into the lifestyle and habitat of their owners. The tiny stalked compound eyes described here probably possessed too few facets to form a proper image, but they represent a sophisticated system for detecting moving objects. The eyes are preserved as almost solid, mace-shaped blocks of phosphate, in which the original positions of the rhabdoms in one specimen are retained as deep cavities. Analysis of the optical axes reveals four visual areas, each with different properties in acuity of vision. They are surveyed by lenses directed forwards, laterally, backwards and inwards, respectively. The most intriguing of these is the putatively inwardly orientated zone, where the optical axes, like those orientated to the front, interfere with axes of the other eye of the contralateral side. The result is a three-dimensional visual net that covers not only the front, but extends also far laterally to either side. Thus, a moving object could be perceived by a two-dimensional coordinate (which is formed by two axes of those facets, one of the left and one of the right eye, which are orientated towards the moving object) in a wide three-dimensional space. This compound eye system enables small arthropods equipped with an eye of low acuity to estimate velocity, size or distance of possible food items efficiently. The eyes are interpreted as having been derived from individuals of the early crustacean Henningsmoenicaris scutula pointing to the existence of highly efficiently developed eyes in the early evolutionary lineage leading towards the modern Crustacea.     

Determining the optimal flexion-extension axis of the elbow in vivo - a study of interobserver and intraobserver reliability

In the current study the interobserver and intraobserver reliability of a recently developed method to obtain the position and orientation vectors of the flexion-extension axis of the elbow in vivo is determined. The method uses the Flock of Birds six degrees-of-freedom electromagnetic tracking device. Ten subjects performed three trials comprising five flexion and extension cycles. The movements of the forearm with respect to the upper arm were recorded. Observer A measured two trials and observer B measured one trial. Optimal instantaneous helical axes were calculated in a humeral coordinate system for each trial. Intraclass correlation coefficients and 99% confidence intervals were computed to compare the three measurements. Zero was in the range of all the narrow confidence intervals, which is strong indication for resemblance. Interobserver intraclass correlation coefficients values for orientation vectors were good to excellent and intraobserver values were fair to good. The intraclass correlation coefficients values for position vectors were lower, probably due to the lack of variance between subjects. It is concluded that the method is reliable and can be used in certain clinical settings.     

Alignment of 3D structures of macromolecular assemblies

MOTIVATION: A number of macromolecular assemblies are being reconstructed in 3D from electron micrographs. The analysis yields a 3D matrix representing the protein density map. In reconstruction processes and in comparing the results of different experiments, it is often necessary to obtain all models oriented the same way in three dimensions. The problem is not trivial since there exist no 3D counterpart of correlation analysis used for 2D images. It is usually solved by time consuming trial and error algorithms. RESULTS: 3D density distributions can be brought to a 'canonical' orientation. The tensor of inertia of the distribution is determined and its eigenvectors are oriented along the coordinate axes. The method is fast and essentially free of reference. It is suitable for structures whose inertial axes do not completely degenerate as they do in icosahedral viruses or if symmetry is cubic. Applications are presented for asymmetric objects and for molecules possessing symmetry axes higher than twofold. IMPLEMENTATION: The implementation simply requires the accumulation of the inertial tensor and its diagonalisation. Volume data rotation has been already illustrated in this journal by the authors.     

Symmetry types, systems and their multiplicity in the structure of adenovirus capsid. II. Rotational facet groups of five-, three- and two-fold symmetry axes

The icosahedral adenovirus capsid has three rotational symmetry axes of different types. The six five-fold, ten three-fold and the fifteen two-fold axes have two superficial points each, altogether 62. The axes determine the number and location of the identical rotational facet groups and that during the different rotational phases which other regular facets and with what multiplicity shall be covered by them. The number of rotational facets of the five-, three- and two-fold rotational symmetry axes is 4, 6.66 and 10, respectively. In all the three cases, there are two kinds of possible arrangements of the facets. During the rotation--when the facets of the facet group placed on one by one to the neighbouring identical facet groups--at the five-fold axes, the facets of the rotational facet group get into cover position 12 times with all the 20 regular capsid facets, 20 times at the three-fold axes, and 30 times at the two-fold axes in a way that a different facet combination (facet hit) falls to every facet, and the original symmetry is not disturbed. After all, this means 240, 400 and 600 facet combinations, i.e. multiplicity in case of five-, three- and two-fold symmetry axes respectively, and these numbers correspond with that of the theoretically possible variations. The same results can be calculated by multiplying the number of real rotations of the capsid bringing the body into itself, i.e. the number 60 with the number of facets contributing to the five-, three- and two-fold rotational phases. The other way of the determination of multiplicity takes into account that all the facet groups of the capsid rotate simultaneously during all the rotational phases, and this multiplies the number of multiplicity with the number of the rotational types five-, three- and two-fold which result in one and the same multiplicity number in the case of five-, three- and two-fold symmetry, alike 1200. Perpendicular to the five-fold symmetry axes with the line of intersection drawn horizontally in the middle along the 6 geodetic ribbon like motifs a regular decagonal intersection forms and the capsid can be cut into two equal parts, in which the polypeptides show a 72 degree rotation from each other, but with a proper rotation the polypeptides get into a congruent position, which means 300 or 600 specific facet combinations. The capsid similar to the icosahedron has also 15 virtual mirror planes which divide the capsid into two, identically arranged halves, forming six right angle triangles on each facet, altogether 120 smaller rectangular so-called Mobius-triangles on the surface. In the three-fold symmetry axis of the facets, these triangles in two separate groups of three can be rotated symmetrically with 120 degrees according to the orientation of the polypeptide subunits in a way that the hexon and other polypeptides too nearly cover each other. Consequently, the adenovirus capsid is a symmetrically arranged body in which several various symmetry types and symmetry systems can be found and their structural symmetry elements exist simultaneously and covering each other. The icosahedral symmetry types and systems are valid and functional simultaneously and in parallel with great multiplicity, but the existence of more than 1500 structural elements in several depth levels, their order of location and distribution make the symmetry of the capsid richer and more complex.     

Errors induced by off-axis measurement of the elastic properties of bone

Misalignment between the axes of measurement and the material symmetry axes of bone causes error in anisotropic elastic property measurements. Measurements of Poisson's ratio were strongly affected by misalignment errors. The mean errors in the measured Young's moduli were 9.5 and 1.3 percent for cancellous and cortical bone, respectively, at a misalignment angle of 10 degrees. Mean errors of 1.1 and 5.0 percent in the measured shear moduli for cancellous and cortical bone, respectively, were found at a misalignment angle of 10 degrees. Although, cancellous bone tissue was assumed to have orthotropic elastic symmetry, the possibility of the greater symmetry of transverse isotropy was investigated. When the nine orthotropic elastic constants were forced to approximate the five transverse isotropic elastic constants, errors of over 60 percent were introduced. Therefore, it was concluded that cancellous bone is truly orthotropic and not transversely isotropic. A similar but less strong result for cortical bone tissue was obtained.     

Comparison of the hindfoot axes of a multi-segment foot model to the underlying bony anatomy

Musculoskeletal models used in gait analysis require coordinate systems to be identified for the body segments of interest. It is not obvious how hindfoot (or rearfoot) axes defined by skin-mounted markers relate to the anatomy of the underlying bones. The aim of this study was to compare the marker-based axes of the hindfoot in a multi-segment foot model to the orientations of the talus and calcaneus as characterized by their principal axes of inertia. Twenty adult females with no known foot deformities had radio-opaque markers placed on their feet and ankles at the foot model marker locations. CT images of the feet were acquired as the participants lay supine with their feet in a semi-weight bearing posture. The spatial coordinates of the markers were obtained from the images and used to define the foot model axes. Segmented masks of the tali and calcanei were used to create 3D bone models, from which the principal axes of the bones were obtained. The orientations of the principal axes were either within the range of typical values reported in the imaging literature or differed in ways that could be explained by variations in how the angles were defined. The model hindfoot axis orientations relative to the principal axes of the bones had little bias but were highly variable. Consideration of coronal plane hindfoot alignment as measured clinically and radiographically suggested that the model hindfoot coordinate system represents the posterior calcaneal tuberosity, rather than the calcaneus as a whole.     

Structural analysis of the S-Layer of Lampropedia hyalina

The two-dimensional (2D) structure of the regularly structured surface layer (S-layer) of the gram-negative eubacterium Lampropedia hyalina has been determined at the molecular level to a nominal resolution of 2.1 nm by transmission electron microscopy and digital image processing. The inner, or “perforate,” layer consists of dimeric block-shaped units located at two-fold symmetry axes. These morphological dimers associate around three-fold symmetry axes to form a continuous layer with p6 symmetry and a lattice constant of 14.6 ± 0.4 nm. Scanning transmission electron microscopy (STEM) yields a mass-per-area (MPA) value for the perforate layer of 3.5 kDa/nm². The outer, or “punctate,” layer is composed of long, roughly cylindrical units centered on six-fold symmetry axes, which are connected by six fine linking arms joining at the three-fold symmetry axes to create a hexagonal layer with a lattice constant of 25.6 ± 0.5 nm. The MPA of the “composite”-i.e., perforate plus punctate—layer is 10.2 kDa/nm².     

Tilt aftereffects generated by symmetrical dot patterns with two or four axes of symmetry

This paper follows from studies by Joung, van der Zwan and Latimer (2000) in which symmetrical dot patterns with one axis of symmetry were used to produce tilt aftereffects (TAEs). The present paper investigates TAE functions produced by symmetrical dot patterns with multiple axes of symmetry. In Experiments 1 and 2, TAE functions produced by dot patterns with two axes of symmetry were compared with TAE functions produced by line stimuli arranged in the same orientation and location as the axes of symmetry in the dot patterns. Similar functions were found. In Experiments 3 and 4, functions produced by dot patterns with four axes of symmetry were compared with functions produced by line stimuli arranged in the same orientation and location as the four axes of symmetry. Again, similar functions were found. These experiments demonstrate that line stimuli and dot stimuli produce similar TAE functions. The implications of these results are discussed.     

A procedure for determining angular positional data relative to the principal axes of the human body

This paper describes a simple computational procedure for determining angular displacement-time histories of human motion from three-dimensional cine data. The method is based on algebraic transformations of coordinates and coordinate axes. Through these coordinate transformations data was acquired for a multi-axial tumbling skill to illustrate angular displacement-time data relative to the moving coordinate system described by the human body through space.     

A joint coordinate system proposal for the study of the trapeziometacarpal joint kinematics

The International Society of Biomechanics (ISB) has recommended a standardisation for the motion reporting of almost all human joints. This study proposes an adaptation for the trapeziometacarpal joint. The definition of the segment coordinate system of both trapezium and first metacarpal is based on functional anatomy. The definition of the joint coordinate system (JCS) is guided by the two degrees of freedom of the joint, i.e. flexion-extension about a trapezium axis and abduction-adduction about a first metacarpal axis. The rotations obtained using three methods are compared on the same data: the fixed axes sequence proposed by Cooney et al., the mobile axes sequence proposed by the ISB and our alternative mobile axes sequence. The rotation amplitudes show a difference of 9 degrees in flexion-extension, 2 degrees in abduction-adduction and 13 degrees in internal-external rotation. This study emphasizes the importance of adapting the JCS to the functional anatomy of each particular joint.     

Evaluation of a 3D object registration method for analysis of humeral kinematics

In 3D image-based studies of joint kinematics, 3D registration methods should be automatic, insensitive to segmentation inconsistencies and use coordinate systems that have clinically relevant orientations and locations because this is important for analyzing rotation angles and translation directions. We developed and evaluated a registration method, which is based on the cylindrical geometry of the humerus shaft and an analysis of the inertia moments of the humerus head, in order to consistently and automatically orient the humerus coordinate system according to its anatomy. Registration techniques must be thoroughly evaluated. In this study we used a well-detectable marker as reference, from which coordinate system determination errors of a 3D object could be measured. This allowed us to quantify by means of unique error analysis the translational and rotational errors in terms of how much and about/along which humeral axis errors occurred. The evaluation experiments were performed using virtual rotations of 3D humeral binary image, a humerus model and a 3D image of a volunteer's shoulder. They indicated that the humeral coordinate system determination errors primarily originated from segmentation inconsistencies, which influenced mostly the humeral transverse axes orientation. The error analysis revealed that the developed registration method reduced the effect of manual segmentation inconsistencies on the orientation of the humeral transverse axes up to 37%, in comparison to the commonly used inertia registration.     

A new non-orthogonal decomposition method to determine effective torques for three-dimensional joint rotation

This paper describes a new non-orthogonal decomposition method to determine effective torques for three-dimensional (3D) joint rotation. A rotation about a joint coordinate axis (e.g. shoulder internal/external rotation) cannot be explained only by the torque about the joint coordinate axis because the joint coordinate axes usually deviate from the principal axes of inertia of the entire kinematic chain distal to the joint. Instead of decomposing torques into three orthogonal joint coordinate axes, our new method decomposes torques into three "non-orthogonal effective axes" that are determined in such a way that a torque about each effective axis produces a joint rotation only about one of the joint coordinate axes. To demonstrate the validity of this new method, a simple internal/external rotation of the upper arm with the elbow flexed at 90 degrees was analyzed by both orthogonal and non-orthogonal decomposition methods. The results showed that only the non-orthogonal decomposition method could explain the cause-effect mechanism whereby three angular accelerations at the shoulder joint are produced by the gravity torque, resultant joint torque, and interaction torque. The proposed method would be helpful for biomechanics and motor control researchers to investigate the manner in which the central nervous system coordinates the gravity torque, resultant joint torque, and interaction torque to control 3D joint rotations.