skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: On the Use of Half-Implicit Numerical Integration in Multibody Dynamics
Abstract This work highlights the use of half-implicit numerical integration in the context of the index three differential algebraic equations (DAEs) of multibody dynamics. Although half-implicit numerical integration is well established for ordinary differential equations problems, to the best of our knowledge, no formal discussion covers its use in the context of index three DAEs of multibody dynamics. We wish to address this since when compared to fully implicit methods, half-implicit integration has two attractive features: (i) the solution method does not require the computation of the Jacobian associated with the constraint, friction, contact, or user-defined applied forces; and (ii) the solution is simpler to implement. Moreover, for nonstiff problems, half-implicit numerical integration yields a faster solution. Herein, we outline the numerical method and demonstrate it in conjunction with three mechanisms. We report on convergence order behavior and solution speed. The Python software developed to generate the results reported is available as open in a public repository for reproducibility studies.  more » « less
Award ID(s):
2209791
PAR ID:
10427418
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Computational and Nonlinear Dynamics
Volume:
18
Issue:
1
ISSN:
1555-1415
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, IEEE International Conference on Communications, Control, and Computing Technologies for Smart Grids (SmartGridComm))
    null (Ed.)
    Transient analysis is vital to the planning and operation of electric power systems. Traditional transient analysis utilizes numerical methods to solve the differential-algebraic equations (DAEs) to compute the trajectories of quantities in the grid. For this, various numerical integration methods have been developed and used for decades. On the other hand, solving the DAEs for a relatively large system such as power grids is computationally intensive and is particularly challenging to perform online. In this paper, a novel machine learning (ML) based approach is proposed and developed to predict post-contingency trajectories of a generator in the time domain. The training data are generated by using an off-line simulation platform considering random disturbance occurrences and clearing times. As a proof-of-concept study, the proposed ML-based approach is applied to a single generator. A Long Short Term Memory (LSTM) network representation of the selected generator is successfully trained to capture the dependencies of its dynamics across a sufficiently long time span. In the online assessment stage, the LSTM network predicts the entire post-contingency transient trajectories given initial conditions of the power system triggered by system changes due to fault scenarios. Numerical experiments in the New York/New England 16-machine 86-bus power system show that the trained LSTM network accurately predicts the generator’s transient trajectories. Compared to existing numerical integration methods, the post-disturbance trajectories of generator’s dynamic states are computed much faster using the trained predictor, offering great promises for significantly accelerating both offline and online transient studies. 
    more » « less
  2. ; (, Journal of Mechanisms and Robotics)
    Abstract This paper presents the design, dynamic modeling, and integration of a single degree of freedom (DOF) robotic leg mechanism intended for tailed quadruped locomotion. The design employs a lightweight six-bar linkage that couples the hip and knee flexion/extension joints mechanically, requiring only a single degree of actuation. By utilizing a parametric optimization, a unique topological arrangement is achieved that results in a foot trajectory that is well suited for dynamic gaits including trot-running, bounding, and galloping. Furthermore, a singular perturbation is introduced to the hybrid dynamic framework to address the lack of robust methods that provide a solution for the differential algebraic equations (DAEs) that characterize closed kinematic chain (CKC) structures as well as the hybrid nature of legged locomotion. By approximating the system dynamics as ordinary differential equations (ODEs) and asymptotically driving the constraint error to zero, CKCs can adopt existing real-time model-based/model-predictive/hybrid-control frameworks. The dynamic model is verified through simulations and the foot trajectory was experimentally validated. Preliminary open-loop planar running demonstrated speeds up to 3.2 m/s. These advantages, accompanied by low-integration costs, warrant this leg as a robust, effective platform for future tailed quadruped research. 
    more » « less
  3. ; ; (, IMA Journal of Numerical Analysis)
    null (Ed.)
    Abstract We implement the unified transform method of Fokas as a numerical method to solve linear evolution partial differential equations on the half-line. The method computes the solution at any $$x$$ and $$t$$ without spatial discretization or time stepping. With the help of contour deformations and oscillatory integration techniques, the method’s complexity does not increase for large $x,t$ and the method is more accurate as $x,t$ increase (absolute errors are smaller, relative errors are bounded). Our goal is to make no assumptions on the functional form of the initial or boundary functions beyond some decay and smoothness, while maintaining high accuracy in a large region of the $(x,t)$ plane. 
    more » « less
  4. ; (, Numerical Methods for Partial Differential Equations)
    Abstract This work focuses on modeling the interaction between an incompressible, viscous fluid and a poroviscoelastic material. The fluid flow is described using the time‐dependent Stokes equations, and the poroelastic material using the Biot model. The viscoelasticity is incorporated in the equations using a linear Kelvin–Voigt model. We introduce two novel, noniterative, partitioned numerical schemes for the coupled problem. The first method uses the second‐order backward differentiation formula (BDF2) for implicit integration, while treating the interface terms explicitly using a second‐order extrapolation formula. The second method is the Crank–Nicolson and Leap‐Frog (CNLF) method, where the Crank–Nicolson method is used to implicitly advance the solution in time, while the coupling terms are explicitly approximated by the Leap‐Frog integration. We show that the BDF2 method is unconditionally stable and uniformly stable in time, while the CNLF method is stable under a CFL condition. Both schemes are validated using numerical simulations. Second‐order convergence in time is observed for both methods. Simulations over a longer period of time show that the errors in the solution remain bounded. Cases when the structure is poroviscoelastic and poroelastic are included in numerical examples. 
    more » « less
  5. ; (, International Journal for Numerical Methods in Engineering)
    Implicit-explicit (IMEX) time integration schemes are well suited for non-linear structural dynamics because of their low computational cost and high accuracy. However, the stability of IMEX schemes cannot be guaranteed for general non-linear problems. In this article, we present a scalar auxiliary variable (SAV) stabilization of high-order IMEX time integration schemes that leads to unconditional stability. The proposed IMEX-BDFk-SAV schemes treat linear terms implicitly using kth-order backward difference formulas (BDFk) and non-linear terms explicitly. This eliminates the need for iterations in non-linear problems and leads to low computational costs. Truncation error analysis of the proposed IMEX-BDFk-SAV schemes confirms that up to kth-order accuracy can be achieved and this is verified through a series of convergence tests. Unlike existing SAV schemes for first-order ordinary differential equations (ODEs), we introduce a novel SAV for the proposed schemes that allows direct solution of the second-order ODEs without transforming them into a system of first-order ODEs. Finally, we demonstrate the performance of the proposed schemes by solving several non-linear problems in structural dynamics and show that the proposed schemes can achieve high accuracy at a low computational cost while maintaining unconditional stability. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email