Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinematics and dynamics”

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 577 records · Page 32

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.↗

A new approach to global control of redundant manipulators

A new and simple approach to configuration control of redundant manipulators is presented. In this approach, the redundancy is utilized to control the manipulator configuration directly in task space, where the task will be performed. A number of kinematic functions are defined to reflect the desirable configuration that will be achieved for a given end-effector position. The user-defined kinematic functions and the end-effector Cartesian coordinates are combined to form a set of task-related configuration variables as generalized coordinates for the manipulator. An adaptive scheme is then utilized to globally control the configuration variables so as to achieve tracking of some desired reference trajectories. This accomplishes the basic task of desired end-effector motion, while utilizing the redundancy to achieve any additional task through the desired time variation of the kinematic functions. The control law is simple and computationally very fast, and does not require the complex manipulator dynamic model.

Seraji, Homayoun↗

How to push a block along a wall

Some robot tasks require manipulation of objects that may be touching other fixed objects. The effects of friction and kinematic constraint must be anticipated, and may even be exploited to accomplish the task. An example task, a dynamic analysis, and appropriate effector motions are presented. The goal is to move a rectangular block along a wall, so that one side of the block maintains contact with the wall. Two solutions that push the block along the wall are discussed.

Mason, Matthew T.↗

The dynamic effects of internal robots on Space Station Freedom

Many of the planned experiments of the Space Station Freedom (SSF) will require acceleration levels to be no greater than microgravity (10 exp -6 g) levels for long periods of time. Studies have demonstrated that without adequate control, routine operations may cause disturbances which are large enough to affect on-board experiments. One way to both minimize disturbances and make the SSF more autonomous is to utilize robots instead of astronauts for some operations. The present study addresses the feasibility of using robots for microgravity manipulation. Two methods for minimizing the dynamic disturbances resulting from the robot motions are evaluated. The first method is to use a robot with kinematic redundancy (redundant links). The second method involves the use of a vibration isolation device between the robot and the SSF laboratory module. The results from these methods are presented along with simulations of robots without disturbance control.

Miller, Jeffrey H.↗

The dynamic effects of internal robots on Space Station Freedom

Many of the planned experiments of the Space Station Freedom (SSF) will require acceleration levels to be no greater than microgravity (10 exp -6 g) levels for long periods of time. Studies have demonstrated that without adequate control, routine operations may cause disturbances which are large enough to affect on-board experiments. One way to both minimize disturbances and make the SSF more autonomous is to utilize robots instead of astronauts for some operations. The present study addresses the feasibility of using robots for microgravity manipulation. Two methods for minimizing the dynamic disturbances resulting from the robot motions are evaluated. The first method is to use a robot with kinematic redundancy (redundant links). The second method involves the use of a vibration isolation device between the robot and the SSF laboratory module. The results from these methods are presented along with simulations of robots without disturbance control.

Miller, Jeffrey H.↗

Characterizing the GD-1 Stream with DESI DR2 Data: Thin Stream and Hot Cocoon

GD-1 is among the longest, coldest stellar streams in the Milky Way, making it an ideal target for probing dark matter substructure through dynamical heating. We present a catalog of 608 spectroscopically confirmed GD-1 members from the first three years of Dark Energy Spectroscopic Instrument (DESI) observations. This constitutes the largest homogeneous spectroscopic sample of GD-1, doubling the number of members previously available only through heterogeneous compilations combining multiple surveys with different systematics. Using these data, we derive updated stream tracks in sky position, proper motion, and radial velocity that extend over $100^\circ$ of the stream. We apply a Gaussian mixture model to decompose the stream into a dynamically cold thin component ($σ_V = 2.49\pm 0.28$ km s$^{-1}$, width $= 0.23\pm0.01^\circ$) and a kinematically hot cocoon ($σ_V = 6.13\pm0.75$ km s$^{-1}$, width $= 2.18\pm0.17^\circ$). The cocoon contains $\sim30\%$ of members and its velocity dispersion is consistent with $\sim11$ Gyr of heating by cold dark matter subhalos. We also detect a large proper motion dispersion ($41.36\pm4.98$ km s$^{-1}$) along the stream direction in the cocoon component. This feature indicates a significant line-of-sight distance spread in the cocoon, and its origin will be further explored in a forthcoming paper. These measurements demonstrate the power of DESI spectroscopy for characterizing the multi-component phase-space structure of stellar streams and constraining small-scale dark matter substructure.

Jarvis, Emma [Toronto U.] (ORCID:0009000656127336)↗

Short Range Correlations and Few Body Precision Measurements

This talk will discuss the exciting experimental program studying short range correlations and few body dynamics at Jefferson Lab in ?=3 nuclei. Quasi-elastic electron scattering on these light nuclei have readily calculable ground states but can already exhibit complex nuclear behavior including nucleon-nucleon (NN) short range correlations (SRCs). 3H and 3He, being mirror nuclei, exploit the maximum available isospin asymmetry and can be compared directly with ab initio calculations. This talk will summarize the recent Hall A 3H and 3He results. The first measurement of SRC nucleons in 3He using the CLAS detector from Hall B and identifying lead neutrons in ?? pairs at high missing momentum will be shown. These fascinating results will be extended in an approved experiment using the CLAS12 detector to measure high-statistics and absolute cross sections on 2H, 3He, and 3H with a wide range of kinematical coverage. This talk will discuss the recent and proposed measurements and their implications on nuclear models for few body nuclear dynamics and short range structure.

Szumila-Vance, Holly↗

Dynamic analysis of astronaut motions in microgravity: Applications for Extravehicular Activity (EVA)

Simulations of astronaut motions during extravehicular activity (EVA) tasks were performed using computational multibody dynamics methods. The application of computational dynamic simulation to EVA was prompted by the realization that physical microgravity simulators have inherent limitations: viscosity in neutral buoyancy tanks; friction in air bearing floors; short duration for parabolic aircraft; and inertia and friction in suspension mechanisms. These limitations can mask critical dynamic effects that later cause problems during actual EVA's performed in space. Methods of formulating dynamic equations of motion for multibody systems are discussed with emphasis on Kane's method, which forms the basis of the simulations presented herein. Formulation of the equations of motion for a two degree of freedom arm is presented as an explicit example. The four basic steps in creating the computational simulations were: system description, in which the geometry, mass properties, and interconnection of system bodies are input to the computer; equation formulation based on the system description; inverse kinematics, in which the angles, velocities, and accelerations of joints are calculated for prescribed motion of the endpoint (hand) of the arm; and inverse dynamics, in which joint torques are calculated for a prescribed motion. A graphical animation and data plotting program, EVADS (EVA Dynamics Simulation), was developed and used to analyze the results of the simulations that were performed on a Silicon Graphics Indigo2 computer. EVA tasks involving manipulation of the Spartan 204 free flying astronomy payload, as performed during Space Shuttle mission STS-63 (February 1995), served as the subject for two dynamic simulations. An EVA crewmember was modeled as a seven segment system with an eighth segment representing the massive payload attached to the hand. For both simulations, the initial configuration of the lower body (trunk, upper leg, and lower leg) was a neutral microgravity posture. In the first simulation, the payload was manipulated around a circular trajectory of 0.15 m radius in 10 seconds. It was found that the wrist joint theoretically exceeded its ulnal deviation limit by as much as 49. 8 deg and was required to exert torques as high as 26 N-m to accomplish the task, well in excess of the wrist physiological limit of 12 N-m. The largest torque in the first simulation, 52 N-m, occurred in the ankle joint. To avoid these problems, the second simulation placed the arm in a more comfortable initial position and the radius and speed of the circular trajectory were reduced by half. As a result, the joint angles and torques were reduced to values well within their physiological limits. In particular, the maximum wrist torque for the second simulation was only 3 N-m and the maximum ankle torque was only 6 N-m.

Newman, Dava J.↗

Investigation into the instantaneous centre of rotation for enhanced design of floating offshore wind turbines

The dynamic behaviour of floating offshore wind turbines (FOWTs) involves complex interactions of multivariate loads from wind, waves, and currents, which result in complex motion characteristics. Although methods for analysing global motion responses are well-established, the time- and location-dependent kinematics remain underexplored. This paper investigates the instantaneous centre of rotation (ICR), a point of zero velocity at a time instance of general plane motion. Understanding and strategically positioning the ICR can reduce the dynamic motion in critical structural locations, enhancing the performance and structural robustness of FOWTs. The paper presents a method for computing the ICR using time-domain simulation results and proposes a statistical analysis approach suitable for design studies. Building on prior research, it examines the sensitivity of the ICR to external loading and design features, providing insights into how these factors influence motion response and how the motion response influences the statistics of the ICR, structural loads, and other performance metrics of interest. The study explores two FOWT configurations, a spar and a semisubmersible, identifying design variables that most effectively control the ICR statistics and identifying the ICR statistics most correlated with the responses of interest. Finally, through two case studies, we demonstrate how to apply these new insights in a practical design scenario. By adjusting the design variables most correlated with the ICR (fairlead vertical position and centre of mass for the spar and mooring line length and offset column diameter for the semisubmersible), we successfully modified the designs of the floating support structures to reduce the loads in the mooring lines, tower base, and blade roots, improving the ultimate strength and fatigue characteristics compared to the original designs.

17 WIND ENERGY↗

Probing nuclear effects with neutrino-induced charged-current neutral pion production

We study neutrino-induced charged-current (CC) $\pi^0$ production on carbon nuclei using events with fully imaged final-state proton-$\pi^0$ systems. Novel use of final-state correlations based on transverse kinematic imbalance enable the first measurements of the struck nucleon's Fermi motion, of the intranuclear momentum transfer (IMT) dynamics, and of the final-state hadronic momentum configuration in neutrino pion production. Event distributions are presented for i) the momenta of neutrino-struck neutrons below the Fermi surface, ii) the direction of missing transverse momentum characterizing the strength of IMT, and iii) proton-pion momentum imbalance with respect to the lepton scattering plane. The observed Fermi motion and IMT strength are compared to the previous MINERvA measurement of neutrino CC quasielastic-like production. The measured shapes and absolute rates of these distributions, as well as the cross-section asymmetries show tensions with predictions from current neutrino generator models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Motion models in attitude estimation

Attitude estimator use observations from different times to reduce the effects of noise. If the vehicle is rotating, the attitude at one time needs to be propagated to that at another time. If the vehicle measures its angular velocity, attitude propagating entails integrating a rotational kinematics equation only. If a measured angular velocity is not available, torques can be computed and an additional rotational dynamics equation integrated to give the angular velocity. Initial conditions for either of these integrations come from the estimation process. Sometimes additional quantities, such as gyro and torque parameters, are also solved for. Although the partial derivatives of attitude with respect to initial attitude and gyro parameters are well known, the corresponding partial derivatives with respect to initial angular velocity and torque parameters are less familiar. They can be derived and computed numerically in a way that is analogous to that used for the initial attitude and gyro parameters. Previous papers have demonstrated the feasibility of using dynamics models for attitude estimation but have not provided details of how each angular velocity and torque parameters can be estimated. This tutorial paper provides some of that detail, notably how to compute the state transition matrix when closed form expressions are not available. It also attempts to put dynamics estimation in perspective by showing the progression from constant to gyro-propagated to dynamics-propagated attitude motion models. Readers not already familiar with attitude estimation will find this paper an introduction to the subject, and attitude specialists may appreciate the collection of heretofore scattered results brought together in a single place.

Chu, D.↗

DEVELOPMENT OF AN INFLIGHT COUNTERMEASURE TO MITIGATE POSTFLIGHT GAIT DYSFUNCTION

Following spaceflight crewmembers experience gait and postural instabilities due to inflight adaptive alterations in sensorimotor function. These changes can pose a risk to crew safety if nominal or emergency vehicle egress is required immediately following long-duration spaceflight. At present, no operational countermeasure is available to mitigate postflight locomotor disturbances. Therefore, the goal of this study is to develop an inflight training regimen that facilitates the recovery of locomotor function after long-duration spaceflight. The countermeasure we are developing is based on the concept of variable practice. During this type of training the subject gains experience producing the appropriate adaptive motor behavior under a variety of sensory conditions and response constraints. This countermeasure is built around current ISS treadmill exercise activities. Crewmembers will conduct their nominal inflight treadmill exercise while being exposed to variations in visual flow patterns. These variations will challenge the postural and locomotor systems repeatedly, thereby promoting adaptive reorganization in locomotor behavior. As a result of this training a subject learns to solve a class of motor problems, rather than a specific motor solution to one problem, Le., the subject learns response generalizability or the ability to "learn to learn" under a variety of environmental constraints. We anticipate that this training will accelerate recovery of postural and locomotor function during readaptation to gravitational environments following spaceflight facilitating neural adaptation to unit (Earth) and partial (Mars) gravity after long-duration spaceflight. The study calls for one group of subjects to perform the inflight treadmill training regimen while a control group of subjects performs only the nominal exercise procedures. Locomotor function in both groups is assessed before and after spaceflight using two tests of gait function: The Integrated Treadmill Locomotion Test (ITLT) and the Functional Mobility Test (FMT). The ITLT characterizes alterations in the integrated function of multiple sensorimotor subsystems responsible for the control of locomotion. This test calls for subjects to walk on a motorized treadmill while we assess changes in dynamic postural stability, head-trunk coordination, short-latency head stabilization responses, dynamic visual acuity, lower limb coordination strategies and gait cycle timing. To make these assessments we measure the following parameters while subjects walk on the treadmill: 1) full body 3-dimensional kinematics using a motion capture system (Motion Analysis Corp., Santa Rosa, CA); 2) the shock-wave transmitted from heel-strike to the head using triaxial accelerometers placed on the tibia and head (Entran, Fairfield, NJ); 3) vertical forces using an instrumented treadmill (Kistler Instrument Corp., Amherst, NY); 4) Dynamic visual acuity using Landolt Cs presented on a laptop computer located 4m from the eyes and 5) Gait cycle timing using foot-switches (Motion Lab Systems, Inc., Baton Rouge, LA) attached to the plantar surface of each shoe at the heel and toe. The FMT evaluates a subject's ability to perform challenging locomotor maneuvers similar to those encountered during an egress from a space vehicle. Subjects step over and duck under obstacles along with negotiating a series of pylons set up on a base of 10 cm thick medium density foam. The dependent measures for the FMT are time to complete the course and the number of obstacles touched. To date, we have collected pre and postflight locomotion data from Expeditions 5-9 who will serve as part of the control group for this study. Preliminary results comparing the recovery rates in gait control sub-systems obtained from the ITLT and FMT performance showed two recovery patterns: 1) a concordant recovery trend between gait control parameters and FMT performance indicating a restitution pattern of recovery and 2) gait controecovery that lagged recovery in FMT performance suggesting that improvement in locomotor function was attained through a pattern of substitution. These data suggest that recovery of postflight locomotor function may occur through adaptive mechanisms that lead to either restitution or substitution of function. Understanding the modes of postflight readaptation has implications for countermeasure development and testing and in astronaut postflight rehabilitation.

Bloomberg, J. J.↗

Development of an Inflight Countermeasure to Mitigate Postflight Gait Dysfunction

Following spaceflight crewmembers experience gait and postural instabilities due to inflight adaptive alterations in sensorimotor function. These changes can pose a risk to crew safety if nominal or emergency vehicle egress is required immediately following long-duration spaceflight. At present, no operational countermeasure is available to mitigate postflight locomotor disturbances. Therefore, the goal of this study is to develop an inflight training regimen that facilitates the recovery of locomotor function after long-duration spaceflight. The countermeasure we are developing is based on the concept of variable practice. During this type of training the subject gains experience producing the appropriate adaptive motor behavior under a variety of sensory conditions and response constraints. This countermeasure is built around current ISS treadmill exercise activities. Crewmembers will conduct their nominal inflight treadmill exercise while being exposed to variations in visual flow patterns. These variations will challenge the postural and locomotor systems repeatedly, thereby promoting adaptive reorganization in locomotor behavior. As a result of this training a subject learns to solve a class of motor problems, rather than a specific motor solution to one problem, Le., the subject learns response generalizability or the ability to "learn to learn" under a variety of environmental constraints. We anticipate that this training will accelerate recovery of postural and locomotor function during readaptation to gravitational environments following spaceflight facilitating neural adaptation to unit (Earth) and partial (Mars) gravity after long-duration spaceflight. The study calls for one group of subjects to perform the inflight treadmill training regimen while a control group of subjects performs only the nominal exercise procedures. Locomotor function in both groups is assessed before and after spaceflight using two tests of gait function: The Integrated Treadmill Locomotion Test (ITLT) and the Functional Mobility Test (FMT). The ITLT characterizes alterations in the integrated function of multiple sensorimotor subsystems responsible for the control of locomotion. This test calls for subjects to walk on a motorized treadmill while we assess changes in dynamic postural stability, head-trunk coordination, short-latency head stabilization responses, dynamic visual acuity, lower limb coordination strategies and gait cycle timing. To make these assessments we measure the following parameters while subjects walk on the treadmill: 1) full body 3-dimensional kinematics using a motion capture system (Motion Analysis Corp., Santa Rosa, CA); 2) the shock-wave transmitted from heel-strike to the head using triaxial accelerometers placed on the tibia and head (Entran, Fairfield, NJ); 3) vertical forces using an instumented treadmill (Kistler Instrument Corp., Amherst, NY); 4) Dynamic visual acuity using Landolt Cs presented on a laptop computer located 4m from the eyes and 5) Gait cycle timing using foot-switches (Motion Lab Systems, Inc., Baton Rouge, LA) attached to the plantar surface of each shoe at the heel and toe. The FMT evaluates s. subject's ability to perform challenging locomotor maneuvers similar to those encountered during an egress from a space vehicle. Subjects step over and duck under obstacles along with negotiating a series of pylons set up on a base of 10 cm thick medium density foam. The dependent measures for the FMT are time to complete the course and the number of obstacles touched. To date, we have collected pre and postflight locomotion data from Expeditions 5-9 who will serve as part of the control group for this study. Preliminary results comparing the recovery rates in gait control sub-systems obtained from the ITLT and FMT performance showed two recovery patterns: 1) a concordant recovery trend between gait control parameters and FMT performance indicating a restitution pattern of recovery and 2) gait controecovery that lagged recovery in FMT performance suggesting that improvement in locomotor function was attained through a pattern of substitution. These data suggest that recovery of postflight locomotor function may occur through adaptive mechanisms that lead to either restitution or substitution of function. Understanding the modes of postflight readaptation has implications for countermeasure development and testing and in astronaut postflight rehabilitation.

Bloomberg, J. J.↗

Modeling and Simulation Efforts to Support Improved Comfort in ARGOS

BACKGROUND: The Active Response Gravity Offload System (ARGOS) provides an analog environment for extravehicular activity (EVA) testing and training. Discomfort has been observed during longer suited test sessions. While the subject’s core is offloaded during surface EVA evaluations, his/her arms experience full Earth gravity and can become overly fatigued, especially during suited tests which involve reaching and prolonged arm extensions. A device (ARGOS Negation of Gravitational Effects on the Limbs: ANGEL) to offload the weight of the arms and suit sleeves is being developed by JSC’s Flight Systems Branch, and here we present preliminary modeling of that device using the open-source biomechanical tool OpenSim [1,2] with an in-house developed plugin. We analyze a series of motions performed by a single shirt-sleeved subject with goals of characterizing the device, validating the model, and predicting whether reduced gravity conditions (i.e., lunar gravity (Lg) or Martian gravity (Mg)) can be accurately simulated with the device, as well as providing comfort to the ARGOS user. METHODS AND RESULTS: To model the offload device, we augment the OpenSim human model topology with the offload mechanism components and joints, using CAD models to represent the mechanism graphically. The joint angles of the device are either obtained from (1) inverse kinematics (IK) using motion capture markers on the various components of the device or (2) calculated in the OpenSim plugin by modeling how the components configure themselves under the offloading spring tension given a particular IK-derived arm position. Given the joint angles of the device, the resulting force on the arm is computed by the plugin and applied as an external load in inverse dynamics (ID) in order to enable study of overall shoulder joint torques as well as offload achieved. We verify the calculated joint angles by using the inverse kinematic data and the forces from manual measurements of the spring both independently and integrated within the device. We found that calculated joint angles generally represent the angles measured and computed with IK, supporting a possible analysis workflow inputting human motion data and observing system behavior under varied design parameters. In two different device configurations in which the maximum applied force was 131 N, our current model accurately captured force with a difference of 2-3 N from measured loads. Though our initial test was performed with a shirt-sleeve subject, arm weights were added to emulate the weight of the suit sleeve and the subject was positioned in a test stand with a Mark-III Hard Upper Torso (HUT) and Portable Life Support System (PLSS) mockup. Arm range of motion tasks were performed outside of the HUT, inside the HUT, and inside the HUT while using the device. A variety of other upper body tasks were completed as well. In summary, we have developed a model to investigate and verify an upper limb offload device currently in development. We believe this model will be a valuable tool not only for device characterization but also to predict proper configurations to simulate Lg or Mg conditions, investigate range of motion concerns, predict limitations such as internal collisions and contacts, and inform future design improvements.

L B Nilsson↗

An adaptive space-time phase field formulation for dynamic fracture of brittle shells based on LR NURBS

Here, we present an adaptive space-time phase field formulation for dynamic fracture of brittle shells. Their deformation is characterized by the Kirchhoff–Love thin shell theory using a curvilinear surface description. All kinematical objects are defined on the shell’s mid-plane. The evolution equation for the phase field is determined by the minimization of an energy functional based on Griffith’s theory of brittle fracture. Membrane and bending contributions to the fracture process are modeled separately and a thickness integration is established for the latter. The coupled system consists of two nonlinear fourth-order PDEs and all quantities are defined on an evolving two-dimensional manifold. Since the weak form requires C 1 -continuity, isogeometric shape functions are used. The mesh is adaptively refined based on the phase field using Locally Refinable (LR) NURBS. Time is discretized based on a generalized-α method using adaptive time-stepping, and the discretized coupled system is solved with a monolithic Newton–Raphson scheme. Finally, the interaction between surface deformation and crack evolution is demonstrated by several numerical examples showing dynamic crack propagation and branching.

42 ENGINEERING↗

On Practical Aspects of Variational Consistency in Contact Dynamics

Usage of contact mechanics methodologies is a pervasive modeling requirement in dynamic simulations. While for some trivial problems, solutions taken from analytical geometry are available, use of a finite element framework is common to achieve formulation generality. This work explores two dynamic contact formulations: one based on the traditional node-to-segment (NTS) approach, and a variationally consistent segment-to-segment (STS) mortar formulation. The NTS formulation employed here enforces the constraints kinematically (i.e., the interpenetration is enforced to the solver tolerance), whereas the mortar approach uses Lagrange multipliers to enforce the contact constraints. Both approaches are implemented in the open-source finite element framework Multiphysics Object-Oriented Simulation Environment (MOOSE). The results highlight two relevant contact-interface-related dynamic phenomena in finite element simulations. First, stabilization of contact constraints is discussed, taking into account the evolution of the total energy in a benchmark problem. Second, the influence of finite element discretization on both of the aforementioned contact formulations is analyzed by exercising a large-deformation example with continuous relative sliding. Variationally consistent contact approaches such as the mortar formulation lead to improved energy preservation and avoid spurious excitation of the system's frequencies. This is especially relevant in settings where inertia and vibrations are of importance.

42 ENGINEERING↗

Nonlinear Interface Reduction for Time-Domain Analysis of Hurty/Craig-Bampton Superelements with Frictional Contact

Virtual prototyping in engineering design relies today on modern numerical models of contacting structures with accurate resolution of interface mechanics, which strongly affect the system-level stiffness and energy dissipation due to frictional losses. High-fidelity modeling within the localized interfaces is required to resolve local quantities of interest that may drive design decisions. The high-resolution finite element meshes necessary to resolve inter-component stresses tend to be computationally expensive, particularly when the analyst is interested in response time histories. The Hurty/Craig-Bampton (HCB) transformation is a widely used method in structural dynamics for reducing the interior portion of a finite element model while having the ability to retain all nonlinear contact degrees of freedom (DOF) in physical coordinates. These models may still require many DOF to adequately resolve the kinematics of the interface, leading to inadequate reduction and computational savings. This study proposes a novel interface reduction method to overcome these challenges by means of system-level characteristic constraint (SCC) modes and properly orthogonal interface modal derivatives (POIMDs) for transient dynamic analyses. Both SCC modes and POIMDs are computed using the reduced HCB mass and stiffness matrices, which can be directly computed from many commercial finite element analysis software. Comparison of time history responses to an impulse-type load in a mechanical beam assembly indicate that the interface-reduced model correlates well with the HCB truth model. Localized features like slip and contact area are well-represented in the time domain when the beam assembly is loaded with a broadband excitation. The proposed method also yields reduced-order models with greater critical timestep lengths for explicit integration schemes.

42 ENGINEERING↗

Asynchronicity in opposed-piston RCMs: Does it matter?

Rapid Compression Machines (RCMs) are widely utilized to study combustion phenomena at engine-relevant conditions, and significant efforts are typically made to create a quiescent environment, particularly for investigations of autoignition chemistry. Opposed-piston configurations can be advantageous due to shorter compression times and reduced surface area to volume ratios. Each side must be actuated simultaneously, but this can be challenging in practice. These devices, like most RCMs, utilize hydraulics for actuation, speed control and arrestation of the piston at the end of the stroke; there is no mechanical control or linkage of the two piston trajectories. To quantify the magnitudes and effects of piston asynchronous behavior, this work employs both detailed experimental measurements and, for the first time, high-fidelity, Direct Numerical Simulation (DNS). The boundary conditions are carefully considered applying insight from high-resolution linear variable differential transformer (LVDT) measurements of the piston trajectory and a zero-dimensional kinematics model of the piston-shaft assembly. Sufficient resolution in the piston crevice region is used. The complicated fluid dynamical behavior that can evolve during piston compression and the ensuing delay processes due to offset timings from t offset = 0-10 ms is elucidated. It is found that near t offset = 6 ms and beyond, the boundary layer on the face of the first-seating piston can be sufficiently perturbed, due initially to reemergence of gas from the crevice of the firstseating piston, so that the adiabatic core can become degraded at long ignition delay times. Substantial mixing of colder gas into the interior of the reaction chamber can alter the measurements, similar to effects previously observed for improper piston crevice configuration. In conclusion, experimental techniques to mitigate asynchronous behavior are discussed and demonstrated.

33 ADVANCED PROPULSION SYSTEMS↗