Questions tagged [constraints]
For questions about computationally solving a problem subject to some constraints on the solution. For optimizations, apply the more specific [constrained-optimization] tag instead.
55
questions
2
votes
1
answer
156
views
How to add damped constraint force to constrained dynamics simulation?
I have implemented a constraint dynamics physics simulation as proposed by Andrew Witkin et al 1990, but I cannot get the initial constraint "snapping" correctly.
I implemented
$$ JWJ^{T} \...
1
vote
1
answer
191
views
Solving linear system of equations with constraints on unknowns
I would like to solve a system of linear equations $y=Uh$ for an unknown vector $h$, where I have a few constraints on some of the elements of $h$. The matrix $U$ is composed of a vector $u$ (length $...
0
votes
0
answers
34
views
Order of Error - Confusion: Clarifying Constraints on Constants and Determining Order of Error
I'm struggling to determine the order of error when considering the error value denoted by $\text{err}$ in relation to the variable $h$. Specifically, I aim to ascertain the value of $x$ in the ...
0
votes
0
answers
158
views
Damping not working in Verlet simulation with absolute constraint
I'm trying to build a simple physic simulation with Verlet intergration, which I'm currently implementing as such:
...
2
votes
1
answer
69
views
Solving constrained odes's using inbuilt solvers in Matlab/Octave
I would like to solve a set of coupled second order differential equations using inbuilt Matlab/Octave subroutines. These equations arise when trying to model sliding of mass ($m_2$) over a wedge of ...
5
votes
1
answer
130
views
What are the various methods in adding an additional constraint to the quadratic spline interpolation problem
I am taking a class on numerical analysis. While the professor was deriving the theory behind quadratic splines, the professor mentioned that a quadratic spline function has the form:
$$
p_{i}(x)=a_{i}...
1
vote
0
answers
299
views
SLSQP solver scipy with linear subset constraints
I have been trying to solve a least squares problem of the following form:
$$
\begin{equation}
\min_{\vec{x}} \frac{1}{2} \lVert f(\vec{x}) - f_{\text{target}} \rVert_{2}^2 + \alpha\Big( \frac{1-\rho}{...
2
votes
0
answers
90
views
How "kinematic" rigid bodies are implemented in physics engines
In most physics engines there's this concept of "static" bodies, which act as rigid bodies with infinite mass. Then there are "kinematic" bodies that act as static bodies, but ...
0
votes
0
answers
49
views
The condtion for Augmented Langrangian Multiplier
I am currently learning the usage of Augmented Lagrangian Multiplier to achieve my equality constraint. I have learnt from the https://en.wikipedia.org/wiki/Augmented_Lagrangian_method that I have two ...
3
votes
0
answers
229
views
Hanging nodes in deal.ii tutorials: how is the continuity constraint imposed?
While looking at step6 of deal.ii tutorials, I decided to try to understand how the constraints coming from hanging nodes are imposed. So I started by watching video lecture 16 by prof. Bangerth
As ...
2
votes
0
answers
37
views
Differences in using Clausius-Duhem inequality vs Principle of Virtual Work/Power in derriving constitutive equations?
I am a novice getting my toes wet in continuum mechanics and nonlinear elasticity. I have seen papers that use both approaches in developing constitutive connections to compliment balance equations of ...
1
vote
1
answer
546
views
How to best code a problem with scipy, cvxpy or Convex.jl with given generated data
I have a curve fitting problem of the form:
$$
\textbf{y} = f(\textbf{x}, a,b,c,d) + \varepsilon
$$
$$
f(x, a,b,c,d) = \frac{b}{e^{x\cdot a}+c}+d
$$
with the constraint
\begin{equation}
\begin{aligned}...
2
votes
1
answer
2k
views
Is there a way to bound the values of a variable when using scipy.integrate.solve_ivp in python?
I want to solve an IVP in python with two variables, x and u, but I need the values of u to be between 0 and 1. Right now it is giving me a solution with negative values for u. Here is the code I have....
1
vote
1
answer
81
views
Constraint programming problem with conditional constraints and some unknown indicator variables
I have an interesting little problem that I believe can be formulated in terms of optimization or constraint programming. I have a few dozen variables $a$, $b$, $c$ ... and a set of constraints that ...
3
votes
1
answer
1k
views
Setting up optimization problem in GEKKO
I have the following dynamical system,
$\frac{d \phi}{dt} = -M^TDM\phi \tag{1}\label{1}$
$\frac{d \hat\phi}{dt} = -M^T\tilde{D}M\hat \phi \tag{2} \label{2}$
$\eqref{1}$ represents the exact ...