Search NASA⌕ Search

SEARCH · Search NASA

Results for “joint coordinates method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 19 records

Systematic generation of multibody equations of motion suitable for recursive and parallel manipulation

The formulation of a method known as the joint coordinate method for automatic generation of the equations of motion for multibody systems is summarized. For systems containing open or closed kinematic loops, the equations of motion can be reduced systematically to a minimum number of second order differential equations. The application of recursive and nonrecursive algorithms to this formulation, computational considerations and the feasibility of implementing this formulation on multiprocessor computers are discussed.

Nikravesh, Parviz E.↗

Multibody for Everybody (M4E): A Symbolic Dynamics Modeling Tool with Applications in Simulation, Control, and Optimization

Developing the analytical model of a multibody system is often the initial step in control and optimization. The analytical model (equations of motion) describes a system’s time evolution under specified forcing conditions. Although developing these equations is easy for simple systems, this process becomes more complex for systems composed of multiple bodies. Deriving equations of motion for complex multibody systems requires specialized expertise in multibody dynamics, is time-consuming, and is susceptible to error. To address this issue, this paper presents an open-source, easy-to-use, systematic framework to derive symbolic equations of motion in both Python and MATLAB using the joint coordinate formulation. This formulation results in a set of ordinary differential equations that use the minimum set of coordinates needed to model a system. The symbolic representation provides better insight into the influence of design parameters on system performance, facilitates sensitivity analysis and parameter studies, and supports direct implementation of control and optimization routines. The tool enables numerical simulation for specified parameter sets, is modular for straightforward integration with other tools and libraries, and allows incorporation of hydrodynamics, mooring, and other external forces. The result is a reproducible, extensible pipeline for modeling, simulation, and design of complex multibody systems. The proposed tool is versatile and can be applied to domains such as robotics, control, and design. In addition, we integrated external libraries that provide capabilities for modeling offshore systems such as underwater robots and marine energy converters.

16 TIDAL AND WAVE POWER↗

Augmented Kalman filter with a reduced mechanical model to estimate tower loads on a land-based wind turbine: a step towards digital-twin simulations

This article presents an application of the Kalman filtering technique to estimate loads on a wind turbine. The approach combines a mechanical model and a set of measurements to estimate signals that are not available in the measurements, such as wind speed, thrust, tower position, and tower loads. The model is severalfold faster than real time and is intended to be run online, for instance, to evaluate real-time fatigue life consumption of a field turbine using a digital twin, perform condition monitoring, or assess loads for dedicated control strategies. The mechanical model is built using a Rayleigh–Ritz approach and a set of joint coordinates. We present a general method and illustrate it using a 2-degrees-of-freedom (DOF) model of a wind turbine and using rotor speed, generator torque, pitch, and tower-top acceleration as measurement signals. The different components of the model are tested individually. The overall method is evaluated by computing the errors in estimated tower-bottom-equivalent moment from a set of simulations. From this preliminary study, it appears that the tower-bottom-equivalent moment is obtained with about 10 % accuracy. The limitation of the model and the required steps forward are discussed.

17 WIND ENERGY↗

Computational Simplification of Robot-Arm Dynamics

Report presents two general methods for reducing mathematical complexities of state equations of robot-arm dynamics. Robot arms contain both rotary and linear joints. Both methods start with homogenous coordinates and Lagrangian formulation of mechanics briefly summarized in report. First method uses matrix-analysis techniques; second, vector-analysis techniques. Vector-analysis method includes new differential-vector representation of centripetal and Coriolis forces.

Bejczy, A. K.↗

Using Analog Field Tests To Link and Prepare Science and In-Situ Resource Utilization for Future Space Missions

A major goal of NASA s human exploration program is to learn how to use the resources of space, known as In-Situ Resource Utilization (ISRU), to lower the cost and risk of human space exploration. Successful implementation of ISRU requires detailed knowledge of surface and subsurface materials, minerals, and volatiles that may be present. This same information is required to better understand the physical and geologic composition, structure, origin, and evolution of the Moon, Mars, and other extraterrestrial bodies of interest. It is also important to recognize that while ISRU and science objectives may be similar, the desired method or hardware to achieve the information desired may be drastically different. One method to promote understanding, coordination, and joint development of instruments and operations between Science and ISRU is the use of analog field demonstrations.

Sanders, Gerald B.↗

Kinematic control of robot with degenerate wrist

Kinematic resolved rate equations allow an operator with visual feedback to dynamically control a robot hand. When the robot wrist is degenerate, the computed joint angle rates exceed operational limits, and unwanted hand movements can result. The generalized matrix inverse solution can also produce unwanted responses. A method is introduced to control the robot hand in the region of the degenerate robot wrist. The method uses a coordinated movement of the first and third joints of the robot wrist to locate the second wrist joint axis for movement of the robot hand in the commanded direction. The method does not entail infinite joint angle rates.

Barker, L. K.↗

A feasibility study of hand kinematics for EVA analysis using magnetic resonance imaging

A new method for analyzing the kinematics of joint motion using magnetic resonance imaging (MRI) is described. The reconstruction of the metacarpalphalangeal joint of the left index finger into a 3D graphic display is shown. From the reconstructed volumetric images, measurements of the angles of movement of the applicable bones are obtained and processed by analyzing the screw motion of the joint. Landmark positions are chosen at distinctive locations of the joint at fixed image threshold intensity levels to ensure repeatability. The primarily 2D planar motion of this joint is then studied using a method of constructing coordinate systems using three or more points. A transformation matrix based on a world coordinate system describes the location and orientation of the local target coordinate system. The findings show the applicability of MRI to joint kinematics for gaining further knowledge of the hand-glove design for EVA.

Dickenson, Reuben D.↗

A feasibility study of hand kinematics for EVA analysis using magnetic resonance imaging

A new method of analyzing the kinematics of joint motion is developed. Magnetic Resonance Imaging (MRI) offers several distinct advantages. Past methods of studying anatomic joint motion have usually centered on four approaches. These methods are x-ray projection, goniometric linkage analysis, sonic digitization, and landmark measurement of photogrammetry. Of these four, only x-ray is applicable for in vivo studies. The remaining three methods utilize other types of projections of inter-joint measurements, which can cause various types of error. MRI offers accuracy in measurement due to its tomographic nature (as opposed to projection) without the problems associated with x-ray dosage. Once the data acquisition of MR images was complete, the images were processed using a 3D volume rendering workstation. The metacarpalphalangeal (MCP) joint of the left index finger was selected and reconstructed into a three-dimensional graphic display. From the reconstructed volumetric images, measurements of the angles of movement of the applicable bones were obtained and processed by analyzing the screw motion of the MCP joint. Landmark positions were chosen at distinctive locations of the joint at fixed image threshold intensity levels to ensure repeatability. The primarily two dimensional planar motion of this joint was then studied using a method of constructing coordinate systems using three (or more) points. A transformation matrix based on a world coordinate system described the location and orientation of a local target coordinate system. Future research involving volume rendering of MRI data focusing on the internal kinematics of the hand's individual ligaments, cartilage, tendons, etc. will follow. Its findings will show the applicability of MRI to joint kinematics for gaining further knowledge of the hand-glove (power assisted) design for extravehicular activity (EVA).

Dickenson, Rueben D.↗

A complete analytical solution for the inverse instantaneous kinematics of a spherical-revolute-spherical (7R) redundant manipulator

Using a method based upon resolving joint velocities using reciprocal screw quantities, compact analytical expressions are generated for the inverse solution of the joint rates of a seven revolute (spherical-revolute-spherical) manipulator. The method uses a sequential decomposition of screw coordinates to identify reciprocal screw quantities used in the resolution of a particular joint rate solution, and also to identify a Jacobian null-space basis used for the direct solution of optimal joint rates. The results of the screw decomposition are used to study special configurations of the manipulator, generating expressions for the inverse velocity solution for all non-singular configurations of the manipulator, and identifying singular configurations and their characteristics. Two functions are therefore served: a new general method for the solution of the inverse velocity problem is presented; and complete analytical expressions are derived for the resolution of the joint rates of a seven degree of freedom manipulator useful for telerobotic and industrial robotic application.

Podhorodeski, R. P.↗

The coupling effects of kinematics and flexibility on the Lagrangian dynamic formulation of open chain deformable links

A nonlinear Lagrangian formulation for the spatial kinematic and dynamic analysis of open chain deformable links consisting of cylindrical joints that connect pairs of flexible links is developed. The special cases of revolute or prismatic joint can also be obtained from the kinematic equations. The kinematic equations are described using a 4x4 matrix method. The configuration of each deformable link in the open loop kinematic chain is identified using a coupled set of relative joint variables, constant geometric parameters, and elastic coordinates. The elastic coordinates define the link deformation with respect to a selected joint coordinate system that is consistent with the kinematic constraints on the boundary of the deformable link. These coordinates can be introduced using approximation techniques such as Rayleigh-Ritz method, finite element technique or any other desired approach. The large relative motion between two neighboring links are defined by a set of joint coordinates which describes the large relative translational and rotational motion between two neighboring joint coordinate systems. The origin of these coordinate systems are rigidly attached to the neighboring links at the joint definition points along the axis of motion.

Changizi, Koorosh↗

A finite element approach for large motion dynamic analysis of multibody structures in space

A three-dimensional finite element formulation for modeling the transient dynamics of constrained multibody space sructures with truss-like configurations is presented. Convected coordinate systems are used to define rigid-body motion of individual elements in the system. These systems are located at one end of each element and are oriented such that one axis passes through the other end of the element. Deformation of each element, relative to its convected coordinate system, is defined by cubic flexural shape functions as used in finite element methods of structural analysis. The formulation is oriented toward joint dominated structures and places the generalized coordinates at the joint. A transformation matrix is derived to integrate joint degree-of-freedom into the equations of motion of the element. Based on the derivation, a general-purpose code LATDYN (Large Angle Transient DYNamics) was developed. Two examples are presented to illustrate the application of the code. For the spin-up of a flexible beam, results are compared with existing solutions available in the literature. For the deployment of one bay of a deployable space truss (the Minimast), results are verified by the geometric knowledge of the system and converged solution of a successively refined model.

Chang, Che-Wei↗

Spherical Coordinate Systems for Streamlining Suited Mobility Analysis

Introduction: When describing human motion, biomechanists generally report joint angles in terms of Euler angle rotation sequences. However, there are known limitations in using this method to describe complex motions such as the shoulder joint during a baseball pitch. Euler angle notation uses a series of three rotations about an axis where each rotation is dependent upon the preceding rotation. As such, the Euler angles need to be regarded as a set to get accurate angle information. Unfortunately, it is often difficult to visualize and understand these complex motion representations. It has been shown that using a spherical coordinate system allows Anthropometry and Biomechanics Facility (ABF) personnel to increase their ability to transmit important human mobility data to engineers, in a format that is readily understandable and directly translatable to their design efforts. Objectives: The goal of this project was to use innovative analysis and visualization techniques to aid in the examination and comprehension of complex motions. Methods: This project consisted of a series of small sub‐projects, meant to validate and verify a new method before it was implemented in the ABF's data analysis practices. A mechanical test rig was built and tracked in 3D using an optical motion capture system. Its position and orientation were reported in both Euler and spherical reference systems. In the second phase of the project, the ABF estimated the error inherent in a spherical coordinate system, and evaluated how this error would vary within the reference frame. This stage also involved expanding a kinematic model of the shoulder to include the rest of the joints of the body. The third stage of the project involved creating visualization methods to assist in interpreting motion in a spherical frame. These visualization methods will be incorporated in a tool to evaluate a database of suited mobility data, which is currently in development. Results: Initial results demonstrated that a spherical coordinate system is helpful in describing and visualizing the motion of a space suit. The system is particularly useful in describing the motion of the shoulder, where multiple degrees of freedom can lead to very complex motion paths.

Benson, Elizabeth↗

Connection forces in deformable multibody dynamics

In the dynamic formulation of holonomic and nonholonomic systems based on D'Alembert-Lagrange equation, the forces of constraints are maintained in the dynamic equations by introducing auxiliary variables, called Lagrange multipliers. This approach introduces a set of generalized reaction forces associated with the system generalized coordinates. Different sets of variables can be used as generalized coordinates and accordingly, the generalized reactions associated with these generalized coordinates may not be the actual reaction forces at the joints. In rigid body dynamics, the generalized reaction forces and the actual reaction forces at the joints represent equipollent systems of forces since they produce the same total forces and moments at and about any point on the rigid body. This is not, however, the case in deformable body analyses wherein the generalized reaction forces depend on the system generalized reference and elastic coordinates. In this paper, a method for determining the actual reaction forces at the joints from the generalized reaction forces in deformable multibody systems is presented.

Shabana, A. A.↗

Force control of a multi-arm robot system

A force-compensated control method for multiple manipulators is presented that allows coordinated manipulation of a jointly grasped object. In this scheme, each arm independently carries out the motions required to realize the desired motion of a prescribed point on the manipulated object. The approach has been implemented and demonstrated on a laboratory system consisting of two industrial, computer-controller manipulators.

Alberts, Thomas E.↗

On the determination of certain astronomical, selenodesic, and gravitational parameters of the moon

A method was examined for joint construction of a selenocentric fundamental system which can be realized by a coordinate catalog of reference contour points uniformly positioned over the entire lunar surface, and determination of the parameters characterizing the gravitational field, rotation, and orbital motion of the moon. Characteristic of the problem formulation is the introduction of a new complex of inconometric measurements which can be made using pictures obtained from an artificial lunar satellite. The proposed method can be used to solve similar problems on any other planet for which surface images can be obtained from a spacecraft. Characteristic of the proposed technique for solving the problem is the joint statistical analysis of all forms of measurements: orbital iconometric, earth-based trajectory, and also a priori information on the parameters in question which is known from earth-based astronomical studies.

Aleksashin, Y. P.↗

Comparing Theoretically Scaled Biomechanical Models

BACKGROUND An investigation into incorporating space suit aspects required a 50th-percentile male model which was not part of the existing dataset of models. Previously, theoretical models for 5th-percentile female and 95th-percentile male were created utilizing a scaled test subject as close as possible to the target height and weight. Scaling factors were created by taking the height ratio and applying it uniformly to the model body segments, then fine tuning to ensure the theoretical model’s height is as expected. This study was initiated to examine the existing method of isometrically scaling in OpenSim [1,2] and to create alternative methods which do not rely on an existing subject being close in height and weight to the theoretical model of interest. Generating a theoretical biomechanical model provides additional abilities without the reliance on available OpenSim models or real subjects. Gained abilities include creating different percentile models and attaining representative anthropometry for any target height/weight, such as targeting specific crew populations or gaps within the current dataset. METHODS AND RESULTS This study includes 3 methods for generating the theoretical model; utilizing the modified unscaled OpenSim Full Body Rajagopal Model (FBRM) [3,4], a dataset of scaled OpenSim models, and existing scaling factors from Dumas et al. [5]. All methods use the Anthropometric Survey of US Army Personnel (ANSUR II) [6] as an input of necessary anthropometric measurements. The following measurements are retrieved from the ANSUR II collection: mass, stature, cervicale height, acromial height, axilla height, waist height, trochanterion height, lateral femoral epicondyle height, lateral malleolus height, acromion-radiale length, radiale-stylion length, palm length, ball of foot length, bicristal breadth, bimalleolar breadth. Some of the measurements are direct segment lengths and others are utilized to derive segment lengths. Height, weight, and age ranges are the required inputs to parse the collection for a mean value of the measurements. A simple iterative process may be necessary for the output of mean mass and stature consistent with the theoretical model of interest. In some cases, one may specify exact values instead of a range but that is dependent upon whether those exact values pertain to a single subject from the ANSUR II collection. Equations from Dumas et al. were used to calculate scaling factors for the two methods that utilize OpenSim scaled and unscaled models. In the third method, Dumas’ scaling factors were used directly instead of generating our own. After scaling factors are achieved, they are applied along with the mass and segment lengths retrieved from ANSUR II in order to calculate the segment mass, Center of Mass (CoM), mass Moment of Inertia (MoI), and the joint location in the parent frame. The theoretical models from the different methods were compared to the previous existing method of isometrically scaling in OpenSim, as well as the standards found in NASA STD-3000 [7] and the NASA Human Integration Design Handbook [8]. The whole-body center of mass is the main comparison performed between the models and the standards. The body segment mass properties were also compared when possible. In some cases, the standards do not provide complete data, or the center of mass location is not provided with respect to the appropriate joint coordinate system. For a 50th-percentile male subject, results from isometric scaling in OpenSim and the method utilizing the unscaled FBRM were compared to the standards in NASA STD-3000. According to the standard, the whole-body vertical CoM location was taken from the head vertex to the CoM. The standards provide 80.2 cm for the CoM vertical component whereas isometric scaling in OpenSim and utilizing the unscaled FBRM provide 82.2 and 82.3 cm, respectively. Through definition, the whole-body CoM in the medial/lateral direction agreed across the models. Investigation is on-going to clearly identify the reference point pertaining to dorsal/ventral whole-body CoM position in standards and relate that to the theoretical models. REFERENCES [1] Delp, S.L., et al., “OpenSim: Open-source Software to Create and Analyze Dynamic Simulations of Movement”, IEEE Transactions on Biomedical Engineering, (2007). [2] Seth, A., Hicks J.L., Uchida, T.K., Habib, A., Dembia, C.L., Dunne, J.J., Ong, C.F., DeMers, M.S., Rajagopal, A., Millard, M., Hamner, S.R., Arnold, E.M., Yong, J.R., Lakshmikanth, S.K., Sherman,n M.A., Delp, S.L. OpenSim: Simulating musculoskeletal dynamics and neuromuscular control to study human and animal movement. Plos Computational Biology, 14(7). (2018) [3] R. K. Huffman, W. K. Thompson, C. A. Gallo, L. J. Quiocho, “Improvement of Scaling and Inverse Kinematic Results with Additional Upper Body Joints Added to Opensim Rajagopal Model”, NASA Human Research Program Investigator’s Workshop, (2019). [4] Huffman R.K., Thompson W., Gallo C., “Modified OpenSim Rajagopal Full Body Model”, New Technology Report, MSC-26872-1. (August 2020). [5] Dumas R., Cheze L. and Verriest J., 2007. “Adjustments to McConville et al. and Young et al. body segment inertial parameters”. Journal of Biomechanics 40, pp. 543-553 [6] Gordon C.C., 2012 “Anthropometric Survey of U.S. Army Personnel: Methods and Summary Statistics”, US Army Natick Soldier RD&E Center. [7] NASA STD-3000/REV-B, 1995. “The Man-System Integration Standards” [8] NASA/SP-2010-3407/REV-1, 2014. “The Human Integration Design Handbook (HIDH)”

F N Matari↗

A Finite-Volume "Shaving" Method for Interfacing NASA/DAO''s Physical Space Statistical Analysis System to the Finite-Volume GCM with a Lagrangian Control-Volume Vertical Coordinate

Toward the development of a finite-volume Data Assimilation System (fvDAS), a consistent finite-volume methodology is developed for interfacing the NASA/DAO's Physical Space Statistical Analysis System (PSAS) to the joint NASA/NCAR finite volume CCM3 (fvCCM3). To take advantage of the Lagrangian control-volume vertical coordinate of the fvCCM3, a novel "shaving" method is applied to the lowest few model layers to reflect the surface pressure changes as implied by the final analysis. Analysis increments (from PSAS) to the upper air variables are then consistently put onto the Lagrangian layers as adjustments to the volume-mean quantities during the analysis cycle. This approach is demonstrated to be superior to the conventional method of using independently computed "tendency terms" for surface pressure and upper air prognostic variables.

Lin, Shian-Jiann↗

A finite element method for time varying geometry in multibody structures

A three-dimensional finite element formulation using convected coordinates is presented for the multibody dynamics of truss-like configurations. Unlike existing formulations, the present one does not superimpose nonlinear rigid body kinematics with linear structural mode shapes, an approach that has recently been shown to be grossly inaccurate under certain conditions. Instead, the finite element method is extended to treat large motions/deformations. The formulation is oriented toward joint dominated structures and places the generalized coordinates at the joints. For the planar spin-up of a flexible beam, results are compared with those derived from a commercially available computer program. The two programs predict nearly identical results.

Housner, J. M.↗