Contacts that only push

Free to turn and unable to

An ellipse in a pocket the size of its own bounding box has four contacts whose rows span two dimensions, so the cone of permitted twists is a whole line and the first-order answer is that it spins both ways. It cannot turn by any amount whatever: the penetration grows as the square of the angle, with a fitted exponent of 1.9944, and a circle in the same pocket turns for ever.

Assumes The test is a program, not a rank and It moves to first order and not at all.

Everything in this field so far is a statement about a linearisation. A contact contributes a linear inequality in the twist, the cone is an intersection of half-spaces, and the margin is a distance in the space those rows live in. It is worth finding out what that costs before trusting any of it.

The case that shows it is not a degenerate drawing and not a mechanism at a singularity. It is an ordinary part in an ordinary pocket.

Free to spin, and it cannot turn at allAn ellipse of semi-axes 1.4 and 0.9 inside four flat walls that touch it at the ends of its own axes. Every normal points at the centre, so every row's moment is nought and the four rows span two dimensions rather than three: the cone of permitted twists is the whole spin axis, a **line** through the origin, and the first-order answer is that the part is free to turn either way. Drag the angle and watch what happens. The ellipse's reach in the direction of the top and bottom walls is √(a²sin²θ + b²cos²θ), which is smallest at θ = 0 and grows from there, so any rotation whatever drives it into both of them — by 0.016 at this angle. **A nullity is a candidate and not a motion**, and this is the shape of case the fields before this one could not produce: not a mechanism at a singularity, but an ordinary part in an ordinary pocket. positioned by solving, not by drawing.rank 2 of 3 · the cone is a linepenetration 0.016 at 0.16 rad
Fig. 1 An ellipse in a pocket the size of its own bounding box. Drag the angle. The first-order test says it may turn either way.

Rank two, and what that means here

Take an ellipse of semi-axes aa and bb and put four flat walls against it, touching at the ends of its own axes. Four contacts, each with its normal pointing at the ellipse’s centre.

A normal through the centre has (p×n)z=0(p \times n)_z = 0: the moment entry of its row is exactly nought, not nearly so. So all four rows lie in the plane ω=0\omega = 0 of the twist space, the rank is two rather than three, and the cone {t:At0}\{t : At \ge 0\} contains the whole spin axis — a line through the origin rather than any number of rays.

Read as a region of centres, that is a uniformly shaded plane: every point, in both senses, permitted.

The rank decision, and how close it is. Every singular value of every arrangement's rows, as a fraction of that arrangement's largest, on a logarithmic scale. Five of the seven have three ordinary values and a rank of three; two of them have two ordinary values and a third at the arithmetic's own floor, sixteen orders of magnitude down. The decision is not close on any row, which is what makes a rank quoted here a measurement rather than an opinion — and it is worth saying because the two rows with a deficient rank are the two the field's most surprising arguments are about, so a reader is entitled to ask whether the deficiency is real. It is: every normal of a circular or elliptical part passes through the part's centre, so the moment entry is exactly nought rather than nearly so.
Fig. 2 The rank decision on every arrangement in the ledger. Two of the seven have a third singular value at the arithmetic’s own floor, sixteen orders down, so the deficiency is a fact rather than a threshold.

The first-order answer is therefore that the part spins freely, both ways, without limit. The answer is wrong, and it is wrong for a reason no amount of care with the rank decision could have caught: the decision is not close.

The support function does the work

Rotate the ellipse by θ\theta about its centre and ask how far it now reaches toward the top wall. That is its support function in the direction (0,1)(0, 1), evaluated on the rotated shape:

h(θ)  =  a2sin2θ+b2cos2θ  .h(\theta) \;=\; \sqrt{a^2 \sin^2\theta + b^2 \cos^2\theta}\;.

At θ=0\theta = 0 that is bb, which is where the wall is. And bb is the minimum of hh, since a>ba > b. So any rotation whatever — either sense — increases the reach and drives the ellipse into both the top and the bottom wall at once. Expanding,

h(θ)b  =  a2b22bθ2+O(θ4),h(\theta) - b \;=\; \frac{a^2 - b^2}{2b}\,\theta^2 + O(\theta^4),

quadratic, with a strictly positive coefficient, and with no linear term at all. The missing linear term is the first-order analysis being right: to first order nothing happens, which is exactly what a rank of two said.

Quadratic, which is why no angle is small enough. The deepest penetration against the angle asked for, on log axes, for the ellipse in its bounding box and for a circle in the same box. The ellipse's line has a fitted slope of 1.9944 against a predicted 2 — the reach grows as b + (a² − b²)θ²/2b, so the depth is quadratic — and the circle's points are not on the chart at all, because its penetration is exactly nought at every angle and a logarithm has nothing to say about that. The two together are the whole of the rung: the first-order analysis returns the same answer for both, a rank of two and a spin left over, and the second-order behaviour is the entire difference between a part that turns freely for ever and one that cannot turn by a thousandth of a radian.
Fig. 3 The depth against the angle asked for, on log axes, with a fitted slope of 1.9944 against a predicted 2. A circle in the same pocket is not on the chart at all.

Measured rather than expanded: at a hundredth of a radian the deepest penetration is 6.4×1056.4 \times 10^{-5}; at a sixteenth, 1.6×1021.6 \times 10^{-2}; and a straight-line fit through five decades of angle gives a slope of 1.9944 against a predicted 2.

The two side walls do the opposite. The reach toward (1,0)(1,0) after a rotation is a2cos2θ+b2sin2θ\sqrt{a^2\cos^2\theta + b^2\sin^2\theta}, which decreases from aa, so the ellipse pulls away from them. Only two of the four contacts do the blocking, and which two is decided by which axis is shorter.

The circle, which is the whole argument

The same computation on a circle in the same pocket returns exactly nought at every angle, because a circle’s support function is constant.

Where the first-order answer is the answer, and where it is not. The ledger with a column added. On five of the seven rows the linear answer is the whole answer: what the cone says the part can do, the part does, and following the escape takes it as far as it likes. On the last two the rank is the same number for opposite reasons — both leave a line, both are the spin of a part whose normals all point at its centre — and the disc turns for ever while the ellipse cannot turn at all. Nothing in the rank distinguishes them, and the thing that does is a second derivative. This is the same asymmetry the network field's obstruction has, arrived at from a different direction: a cone is an upper bound on what a part can do, exact when there is nothing to check and a candidate list when there is.
Fig. 4 The ledger with a column added. Two rows have the same rank and opposite answers, and nothing in the first-order analysis distinguishes them.

So here are two arrangements with:

  • the same number of contacts, four;
  • the same normals, all pointing at the centre;
  • the same rank, two;
  • the same escape cone, the whole spin axis;
  • the same margin, nought;

and opposite behaviour. One turns for ever and the other cannot turn by a thousandth of a radian. The entire difference is a second derivative, and every instrument in the four rungs before this one returns identical output for both.

That is a stronger statement than the corresponding one in the network field, and worth comparing. There, a flex the rank permits may be blocked, and the obstruction is a number computed from the dependencies among the constraints — a framework with no dependencies has no obstruction, which covers every mechanism in eighteen of this site’s fields. Here the obstruction has nothing to do with dependencies: the rows of the ellipse’s arrangement are independent in pairs and their deficiency is a rank deficiency, and the thing that blocks the motion is the curvature of the part, which is not in the constraint matrix at all.

Where the curvature went

That is worth saying twice, because it is the sharpest thing in the rung.

A contact’s row is built from a point and a normal. Nothing about the surface’s shape enters it — not its curvature, not its extent, not whether it is a flat or an arc. Two parts with the same contact points and the same normals produce identical matrices and identical cones, however differently they are shaped between those points.

What one contact forbids, drawn as a placeA single contact on one edge of a square, and the whole plane coloured by what it permits. A rotation about a point is a twist, and a twist is permitted when it does not drive the part into the obstacle; because a rotation about (x, y) is affine in the point, the condition is a **half-plane** and the boundary is a straight line — the line through the contact along its own surface. On one side of it only anticlockwise rotations are permitted, on the other only clockwise, and the two together are the whole plane bar the line itself. So one contact rules out exactly half of what the part could do and leaves the other half untouched, which is why the count of contacts a hold needs is one more than the dimension rather than equal to it: the first 1 of them cannot leave nothing over. The picture is exact — the regions are clipped polygons, not a sampled grid.anticlockwise permittedclockwise permittedone contact · one half-planethe boundary is the surface's own line
Fig. 5 The half-plane one contact leaves. It is decided by the tangent line and by nothing else about the surface, at first order.

So the first-order analysis is not merely a linearisation of the geometry; it is an analysis that has thrown the geometry away except at the contact points. Most of the time that is harmless, because the linear answer is already trivial and there is nothing left for curvature to do. It stops being harmless exactly when the linear answer is not trivial and the contacts are curved, which is what a rank deficiency on a smooth part looks like.

Reading it the other way: the field’s whole machinery works because most parts are polyhedral. A polygon has flat faces, a flat face has zero curvature, and a flat contact’s second-order term vanishes — so for a polygonal part the first-order answer is the answer.

Measured rather than argued

The expansion above is exact and the site’s habit is to have a second route, so the numbers here come from a walk rather than from a series.

The part is rotated about the twist’s own centre by a finite angle — a rigid motion, not a step of one — and the depth is how far past each obstacle plane the part now reaches, computed from its support function in the direction that plane is being pushed. A genuine motion gives a depth at the arithmetic’s floor at every angle and no exponent worth fitting; a blocked one gives a depth growing as θ2\theta^2, and the exponent is the evidence rather than the assumption.

That is the same instrument the network field walks a flex with, and it is here for the same reason: a cone is a candidate list and a walk is what settles it.

There is one difference and it is in the field’s favour. A framework’s walk needs a Newton projection at every step, because the constraints are equations and a step off the manifold has to be corrected. Here the motion is a rigid rotation and there is nothing to correct — the part is rigid, the walk is exact, and the only question is whether the moved part is inside the obstacle. So the second-order test in this field costs a support function evaluation and is exact at any angle, which is why the figures can put a slider on it.

Second-order form closure, and what it is not

An arrangement that is not form-closed at first order and cannot move at all is called second-order form closure, and the name is worth handling carefully because it invites two wrong readings.

It is not a weaker kind of hold. The ellipse in its box cannot rotate by any amount, at all, ever; there is no motion, not a small one. What is weaker is the test, not the restraint.

And it is not friction. Nothing above needs a coefficient, a normal force or a material. The blocking is done by the shapes: the pocket is the ellipse’s bounding box and the ellipse does not fit into it turned.

How far inside the hull the origin actually is. The same seven arrangements with their margins drawn rather than tabulated, because the shape of this chart is the argument: the quantity is not a probability and not a percentage, it is a distance — how far the origin sits from the nearest face of the hull of the contact rows, with every row a unit vector so the number is comparable across arrangements. The two that hold come in at 0.211 and 0.091; the five that do not come in at exactly nought, and they are drawn at nought rather than left off. A margin that falls smoothly to nothing is what makes this a measurement: an arrangement approaching one that lets go says so before it does.
Fig. 6 The margin, which reports nought for the ellipse in its box. That number is right about what it measures and is not the answer to whether the part moves.

What is genuinely weaker is the robustness. A first-order hold survives a perturbation: move the contacts a little and the margin changes a little and stays positive. A second-order hold does not, in one direction — give the pocket a clearance and the ellipse can turn, by an angle that goes as the square root of the clearance, since the penetration is quadratic. A clearance of a micron on a part 10 mm across buys about 0.03 radians of rotation, which is nearly two degrees. That is the same square-root relationship the network field measured on a blocked flex, fitted there at 0.512 against a half, and it arrives here for the same reason: a quadratic obstruction crossed at a given tolerance is crossed at the square root of it.

So a second-order hold is a hold that a tolerance destroys much faster than a first-order one, and a designer who has one should know which kind they have.

The pocket nobody would draw, and the one everybody does

The ellipse in its bounding box reads as a contrived example, and it is worth saying why the situation is not.

Every part with a smooth boundary held in a pocket cut to fit is this arrangement. A bearing race in a housing, a lens in a cell, a cam follower roller between its cheeks, a shaft in a bush — each is a curved part against a curved or flat obstacle, touching where the two are tangent, with every normal running along the line of centres. The moment entries are nought or nearly so, the rank is short, and whether the part turns is decided by the relative curvature and not by anything the contact rows contain.

For a circular part in a circular pocket that spin is real and is the point of the arrangement: a bearing is exactly a second-order failure to hold, deliberately, and it is the same observation the counting field makes about a pin from the other side. For anything not circular it is not, and the difference is invisible to the first-order test.

The three points, and the circles that stand in for them. The same three points as the log-log measurement, drawn where they sit on the coupler with the circle each one's path is momentarily on. They are ordinary-looking points and their circles are ordinary-looking circles; nothing in the picture distinguishes 3.00, 3.94, 4.97 orders of contact. That is the argument for measuring rather than drawing: the difference between these three is entirely a difference in how long the agreement lasts. positioned by solving, not by drawing.
Fig. 7 Two curves in contact to a stated order, from the curvature field. The order of contact is the quantity this rung turns on, and it is not in any contact row.

Which suggests a way to read the whole rung. The first-order test asks whether the contacts surround the part in direction; the second-order test asks whether they surround it in shape. A polygon has no shape between its contacts and the first question is the whole question. A smooth part has shape everywhere, and the arrangements where that matters are exactly the ones where the normals line up.

Where a reader has met this before

Two places on this site, and they are opposite cases, which is why neither of them prepared anybody for this one.

A four-bar at a toggle has a momentarily degenerate Jacobian and a perfectly good motion straight through it. The rank is short and the flex is real. A parallel platform at a direct singularity gains a freedom with every actuator locked, and that freedom is real too — the platform genuinely moves.

In both, a rank deficiency was a warning that the mechanism was about to do something. Here it is a rank deficiency that means nothing at all, on a part that is going nowhere, and the sign of it is indistinguishable.

Two rank deficiencies, and how to tell them apart

Since the first-order test cannot separate the two cases, it is worth saying what can, cheaply, before anybody reaches for a walk.

The null direction of a rank-two contact arrangement in the plane is always a pure rotation about some point — that is what a deficient moment column means. So the test is: rotate the part about that point and ask whether its support function in each blocked contact’s normal direction is at a minimum there. If it is, any rotation increases the reach and the part is blocked. If it is at a maximum the part pulls away and the contact was never doing anything. If it is stationary to all orders, the part is a circle about that point and the spin is real.

That is a statement about one scalar function of one variable and it costs nothing. It also explains the ledger’s two smooth rows in one sentence each: the ellipse’s support function toward its walls is at a minimum, so it is blocked; the disc’s is constant, so it turns.

The reason the field does not build that test into the routine is that it is planar and specific, and the walk is neither. A walk works in space, works for any twist, and returns a number with an exponent attached rather than a verdict — which is what the site’s habit asks for.

The coefficient, in closed form

The penetration is measured as quadratic with a fitted exponent of 1.9944, and the coefficient has a closed form worth writing down, because it says exactly what the second-order hold depends on and when it disappears.

The ellipse’s support function in the direction of the top wall, after a rotation by θ\theta, is a2sin2θ+b2cos2θ\sqrt{a^2\sin^2\theta + b^2\cos^2\theta}. At θ=0\theta = 0 that is bb, which is where the wall is. Expanding,

h(θ)b+a2b22bθ2,h(\theta) \approx b + \frac{a^2 - b^2}{2b}\,\theta^2 ,

so the penetration is 12κθ2\tfrac{1}{2}\kappa\theta^2 with κ=(a2b2)/b\kappa = (a^2-b^2)/b — no linear term, a strictly positive coefficient, and everything in it a function of the two semi-axes alone.

Check it against the measurement. With a=1.4a = 1.4 and b=0.9b = 0.9 the coefficient (a2b2)/2b(a^2-b^2)/2b is 1.15/1.8=0.63891.15/1.8 = 0.6389, and at a hundredth of a radian that gives a penetration of 6.39×1056.39\times10^{-5} — against the 6.4×1056.4\times10^{-5} the sweep reports. The closed form and the rotated-and-measured penetration are the same number, which is the second route this claim was owed.

The expression settles the circle in one line as well. At a=ba = b the numerator is zero, the coefficient vanishes identically, and the penetration is nought at every angle rather than small — which is the circle’s free spin, arriving as a special case rather than as a separate computation.

And it says how strong the hold is, which the exponent alone cannot. The coefficient is (a2b2)/2b(a^2-b^2)/2b, so it grows with how far the part is from circular and collapses as it approaches a circle. A part with semi-axes differing by one per cent has a coefficient two orders of magnitude below this ellipse’s, so the angle at which it penetrates by a given amount is ten times larger. A second-order hold on a nearly round part is nearly no hold at all, and the arithmetic says precisely how nearly.

That gives the robustness objection its number. A clearance cc is defeated when the penetration reaches it, at θ=2c/κ\theta = \sqrt{2c/\kappa} — so the permitted rotation goes as the square root of the clearance and the inverse square root of the coefficient. Ten times the clearance buys about three times the rotation; a part ten times rounder buys about three times the rotation as well. Both are gentle dependences, which is exactly why a second-order hold degrades so much faster than a first-order one, where the permitted motion is linear in the clearance and vanishes with it.

What is not settled

Second order is not the end. An arrangement that clears the second-order test can still be blocked at third, and the general question of whether a first-order escape extends to a finite motion is not decided by any finite number of derivatives. The walk is the instrument that settles it and a walk that goes a long way is evidence rather than proof — the same caution a search that has stopped finding new answers earns on this site.

And nothing here is a general test. What is computed is the penetration of one particular part along one particular twist, from that part’s own support function. That works because every obstacle in this field is a plane and every part is convex. A concave part against a concave obstacle has a relative curvature that can go either way, the blocking condition involves both surfaces, and this field does not build one — which is a statement about scope rather than about difficulty.

A point, its pole, and the centre it is turning aboutThe tracing point is on the coupler at (0.45, 0.5) of its length. The cross is the pole, the faint curve is the path the point traces over a whole turn, and the circle is the one that path is momentarily on — centre marked, radius 0.267. **The point, the pole and the centre are collinear**, which is not an accident of this position: a point's centre of curvature always lies on its own ray from the pole, and Euler and Savary's relation says where on it. Here that relation puts the centre 5.6e-16 of a unit from where differentiating the loop equation three times puts it. positioned by solving, not by drawing.the pointits centretwo routes agree to 5.6e-16positioned by solving, not by drawing
Fig. 8 Two curved surfaces in contact, from the meshing field, where the relative geometry of the pair is the whole subject. That is where a general second-order condition would have to come from.

And the sign of the second-order term is the whole answer. For the ellipse in its bounding box it is positive at both blocking contacts and there is nothing to trade against, which is why the block is unconditional. An arrangement whose second-order terms have mixed signs permits some finite motion and refuses others, and the boundary between them is a curve rather than a cone. That case exists, it is not in the ledger, and the honest reason is that nothing this field is about produces one.

About the same objects

Not linked from either essay — found by the objects both name.

What links here

Essays that link to this one from their own argument.

The objects this essay names

Each one links to every other essay that touches it.

ConditioningEscape coneForm closureInfinitesimal flexNull spacePath curvatureRankSecond-orderSupport functionUnilateral constraint