An official website of the United States government
A finite element based computational scheme is developed and employed to assess a duality based variational approach to the solution of the linear heat and transport PDE in one space dimension and time, and the nonlinear system of ODEs of Euler for the rotation of a rigid body about a fixed point. The formulation turns initial-(boundary) value problems into degenerate elliptic boundary value problems in (space)-time domains representing the Euler-Lagrange equations of suitably designed dual functionals in each of the above problems. We demonstrate reasonable success in approximating solutions of this range of parabolic, hyperbolic, and ODE primal problems, which includes energy dissipation as well as conservation, by a unified dual strategy lending itself to a variational formulation. The scheme naturally associates a family of dual solutions to a unique primal solution; such ‘gauge invariance’ is demonstrated in our computed solutions of the heat and transport equations, including the case of a transient dual solution corresponding to a steady primal solution of the heat equation. Primal evolution problems with causality are shown to be correctly approximated by noncausal dual problems.
more »
« less
This content will become publicly available on October 1, 2025
A Hidden Convexity of Nonlinear Elasticity
Abstract A technique for developing convex dual variational principles for the governing PDE of nonlinear elastostatics and elastodynamics is presented. This allows the definition of notions of a variational dual solution and a dual solution corresponding to the PDEs of nonlinear elasticity, even when the latter arise as formal Euler–Lagrange equations corresponding to non-quasiconvex elastic energy functionals whose energy minimizers do not exist. This is demonstrated rigorously in the case of elastostatics for the Saint-Venant Kirchhoff material (in all dimensions), where the existence of variational dual solutions is also proven. The existence of a variational dual solution for the incompressible neo-Hookean material in 2-d is also shown. Stressed and unstressed elastostatic and elastodynamic solutions in 1 space dimension corresponding to a non-convex, double-well energy are computed using the dual methodology. In particular, we show the stability of a dual elastodynamic equilibrium solution for which there are regions of non-vanishing length with negative elastic stiffness, i.e. non-hyperbolic regions, for which the corresponding primal problem is ill-posed and demonstrates an explosive ‘Hadamard instability;’ this appears to have implications for the modeling of physically observed softening behavior in macroscopic mechanical response.
more »
« less
- Award ID(s):
- 2021019
- PAR ID:
- 10548670
- Publisher / Repository:
- Springer
- Date Published:
- Journal Name:
- Journal of Elasticity
- Volume:
- 156
- Issue:
- 3
- ISSN:
- 0374-3535
- Page Range / eLocation ID:
- 975 to 1014
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Natural dynamics, nonlinear optimization, and, more recently, convex optimization are available methods for stiffness design of energy-efficient series elastic actuators. Natural dynamics and general nonlinear optimization only work for a limited set of load kinetics and kinematics, cannot guarantee convergence to a global optimum, or depend on initial conditions to the numerical solver. Convex programs alleviate these limitations and allow a global solution in polynomial time, which is useful when the space of optimization variables grows (e.g., when designing optimal nonlinear springs or co-designing spring, controller, and reference trajectories). Our previous work introduced the stiffness design of series elastic actuators via convex optimization when the transmission dynamics are negligible, which is an assumption that applies mostly in theory or when the actuator uses a direct or quasi-direct drive. In this work, we extend our analysis to include friction at the transmission. Coulomb friction at the transmission results in a non-convex expression for the energy dissipated as heat, but we illustrate a convex approximation for stiffness design. We experimentally validated our framework using a series elastic actuator with specifications similar to the knee joint of the Open Source Leg, an open-source robotic knee-ankle prosthesis.more » « less
-
Compared to rigid actuators, Series Elastic Actuators (SEAs) offer a potential reduction of motor energy consumption and peak power, though these benefits are highly dependent on the design of the torque-elongation profile of the elastic element. In the case of linear springs, natural dynamics is a traditional method for this design, but it has two major limitations: arbitrary load trajectories are difficult or impossible to analyze and it does not consider actuator constraints. Parametric optimization is also a popular design method that addresses these limitations, but solutions are only optimal within the space of the parameters. To overcome these limitations, we propose a non-parametric convex optimization program for the design of the nonlinear elastic element that minimizes energy consumption and peak power for an arbitrary periodic reference trajectory. To obtain convexity, we introduce a convex approximation to the expression of peak power; energy consumption is shown to be convex without approximation. The combination of peak power and energy consumption in the cost function leads to a multiobjective convex optimization framework that comprises the main contribution of this paper. As a case study, we recover the elongation-torque profile of a cubic spring, given its natural oscillation as the reference load. We then design nonlinear SEAs for an ankle prosthesis that minimize energy consumption and peak power for different trajectories and extend the range of achievable tasks when subject to actuator constraints.more » « less
-
Abstract In this work, we investigate the existence and uniqueness properties of a composite structure (multilayered)–fluid interaction PDE system which arises in multi-physics problems, and particularly in biofluidic applications related to the mammalian blood transportation process. The PDE system under consideration consists of the interactive coupling of 3D Stokes flow and 3D elastic dynamics which gives rise to an additional 2D elastic equation on the boundary interface between these 3D PDE systems. By means of a nonstandard mixed variational formulation, we show that the PDE system generates a$$C_0$$ -semigroup on the associated finite energy space of data. In this work, the presence of the pressure term in the 3D Stokes equation adds a great challenge to our analysis. To overcome this difficulty, we follow a methodology which is based on the necessarily non-Leray-based elimination of the associated pressure term, via appropriate nonlocal operators. Moreover, while we express the fluid solution variable via decoupling of the Stokes equation, we construct the elastic solution variables by solving a mixed variational formulation via a Babuska–Brezzi approach.more » « less
