The 3-Second Throw That Couldn’t Happen — a legal science story

How I used physics, PyBullet, and a hollow steel door to test a contested accusation against the laws of motion.

When the laws of physics become your defence

Most people who get falsely accused of something call a lawyer. I started a Python project.

The accusation was specific. Three seconds in a corridor. I was supposed to have pulled the alleged victim back, swapped places with her, thrown her back-first into the elevator door, and swapped places again. All of this in front of a third party — a court-appointed social curator — who, by her own testimony, had her back turned at the precise moment of the alleged impact.

Three seconds. Two doors, two metres apart. A 70 kg adult moving as a rigid object. A 90 kg actor doing the moving.

I am not a lawyer. I am the kind of person who reaches for scipy.stats and compute_impact() when faced with an emotional problem. So I built henryk-simulations.

This article walks through the deconstruction. The headline finding: at the most charitable possible interpretation of the accusation, the impact would have delivered 18 kilonewtons of peak force, 26 g of deceleration, and a peak sound pressure level of 124 dB at the phone microphone that was recording the entire visit. The medical examination afterwards documented one bruise on the right shoulder. The recording contains neither a clipping spike nor any panel-ringing acoustic signature. The third-party witness reports no loud noise of any kind.

You can decide for yourself what that means.

The corridor is two metres wide from the apartment door (north wall) to the elevator door (south wall). The elevator door is a standard hollow steel design: two 2 mm plates with a 3 cm air gap, plus a small glass window. The apartment door is about 1 m wide and opens into the corridor.

I stood with my back pressed flat against the elevator door (maximum retreat — I had been asked to step aside). The alleged victim stood at the apartment door. The third-party witness stood to the west, in the narrower entrance segment of the corridor, with a clear line of sight.

The accusation required three things to happen inside this 2 m by 3 s envelope:

  1. I would pull the victim ~1.5 m south, toward the elevator door
  2. I would somehow swap positions with her and throw her back-first into that door (the swap is non-trivial — we would be passing through each other or each performing a half-rotation)
  3. I would swap positions again, ending up myself facing the door

The third-party witness corroborates one part of this last claim. When she turned around, the victim was indeed in front of me — but she was leaning against me front-first, slid down the door front-first, then crawled forward on all fours. The scream started at the moment the witness turned, not at the alleged moment of impact.

I’ll come back to these details. First, the physics.

A minimum-phase reconstruction

The framing I landed on early was this. Don’t try to simulate what happened. Define an impossible-movement test.

The idea is this: decompose the accusation into the smallest set of phases that still admits the verbatim claim, and give each phase the maximum time the 3 s budget allows. That produces the lowest possible peak velocity, acceleration, and force any motion compatible with the accusation could have.

If that minimum is already extreme — if even the most charitable reconstruction demands superhuman acceleration or biomechanically-impossible forces — then any richer reconstruction, which would compress phases into shorter intervals, is strictly worse.

The decomposition came out to three phases of one second each:

This is formally an ELBO-style lower bound on the required demand. The Evidence Lower BOund (ELBO) comes from variational inference: when the true posterior is intractable, you pick a tractable surrogate distribution and lean on Jensen’s inequality to drop a lower bound underneath the log-evidence. Jensen’s says that for a concave function the expectation of the function is at most the function of the expectation:

Apply that to the log-evidence with surrogate q(z) standing in for the intractable posterior p(z|x)

Equivalently the gap is exactly a Kullback-Leibler divergence:

KL is the distance between our story and the true one — and our lower bound undershoots reality by exactly that much. Tighten the story, shrink the distance, and the bound climbs toward the truth. Jensen is the engine, the concavity of logis the chassis.

Here we are doing the same move one storey up. The intractable object is the true physical demand D(M_true); the tractable family is the set of minimum-phase reconstructions, optimised within it to give a surrogate D_min(q*). By construction:

A violation at the lower bound is a violation at every richer decomposition.

That is — I am cheating on the prosecution’s behalf. They get the friendliest possible reading of the accusation.

The numbers

The kinematics are not subtle. With a 1.5 m translation in 1.0 s starting from rest, peak acceleration during the pull phase is 6.0 m/s² — about z = 3.75 against a recreational-sprint reference (extreme band). The 180° rotations in each of the other two phases require 6.28 rad/s of peak angular velocity — implausible against published standing-pivot rates of 3.5 ± 1.0 rad/s.

But that’s the motion side. The interesting numbers are on the impact side.

The continuous-velocity model treats swap-throw as a continuation of the pull. The actor does not release the victim at the moment of door contact — they keep applying force — so the victim hits the door at her peak velocity, not at the velocity she had at the end of the pull. This works out to:

That is 11.6 kph. Modest, against a typical body. But the elevator door is not a typical body. The steel panel is effectively rigid — it barely deforms — so almost all of the deceleration distance comes from the body itself: the ~2 cm of soft-tissue and pre-yield rib-cage compliance before bone fails (Lobdell 1973; Kroell 1971). Plug in d = 2 cmtoa = v²/2dfor a 70 kg torso:

Eighteen kilonewtons is the kind of load you get from a small SUV nudging a wall at ten kilometres an hour. Twenty-six g is in the race-car-crash range. The kinetic energy puts you in the “serious thoracic injury, organ contusion” band per Sturdivan 2004. The peak force is in AIS 5+ — flail chest, serious organ injury, life-threatening — per Viano 1989.

You will notice an apparent contradiction in that figure. The impact velocity (11.6 kph) reads “moderate, bruising” while the force, acceleration, and kinetic energy all read “serious” or “critical”. This is not a contradiction. The published velocity bands assume a soft target — chest deformation of 5–10 cm against a compliant surface. A steel door is not a compliant surface. The same 11.6 kph absorbed over 2 cm rigid instead of 10 cm soft multiplies the peak force by roughly five times. The four metrics are four views of the same impact. The velocity is the cause; the force, the g, and the deposited energy are the outcomes.

The forensic angle: predicted acoustics

Here is where the project goes from “biomechanics simulation” to “forensic test you can actually run”.

If 70 kg of human torso slams into a 2 m × 1 m steel panel at 11.6 kph, the panel should make a noise. A predictable noise. Kirchhoff thin-plate theory gives the flexural modal frequencies of a simply supported rectangular plate:

with flexural rigidity D = Eh³/(12(1-v²)and areal σ = ρh density . For a 2 m × 1 m × 2 mm steel plate, the first few modes land at 6, 9.6, 15.6, 20.4, 24, 24 Hz. The 3 cm sealed air cavity between the two steel plates of the hollow door gives a half-wave resonance c /2dat 5717 Hz — the bright “clang” you hear when something hits a hollow metal door.

That is the spectrum. The amplitude is the other half.

Acoustic power:

For a struck steel panel, typically falls between 0.001 (lossy mounting) and 0.05 (well-coupled).

Intensity at the listener:

Across that radiation-efficiency bracket, here is what the predicted impact should sound like:

At the phone microphone, two metres from the door: 114 to 131 dB SPL peak. Consumer phone microphones clip at around 120 dB. The audio recording of the visit, which was running throughout, should therefore contain at least four things:

  1. A clipping spike at the moment of impact
  2. Sub-bass panel ringing in the 6–24 Hz band, decaying over hundreds of milliseconds
  3. A bright cavity-mode resonance around 5.7 kHz
  4. An audible reaction from the third-party witness — she was 1.5 m away, predicted SPL 126 dB, jackhammer-at-one-metre territory

The recording contains none of these. The witness reports no bang. By the time the scream begins, on her own testimony, she had already turned back around to look.

Predicted vs. documented

Here is the side-by-side I find hardest to read with a straight face:

Of ten predicted observables, zero match. Seven are absent or directly contradicted. One is temporally decoupled. Two are pending direct audio spectral analysis. The single right-shoulder bruise is geometrically and energetically inconsistent with a back-first whole-torso impact at 18 kN.

One more thing about that bruise. The defendant’s position is that it was self-inflicted by the alleged victim after the visit, to seed a corroborating mark in the medical record. The defendant has chosen not to dispute the bruise’s origin directly. The argument above runs on physics, not on contested medical interpretation. The mechanical impossibility of producing a single localised shoulder bruise from the alleged whole-torso back-first impact at 18 kN is enough, on its own, to falsify the verbatim claim. Newton was kind enough to do the talking.

The takeaway

I am not the first person to use physics in a courtroom. Crash reconstruction is a recognised forensic discipline, and so is gunshot-residue analysis, and so is blood-pattern analysis. What is unusual here is the level the test runs at: not “did this specific impact pattern occur” but “could this story have happened at all, given Newton’s second law and the geometry of a Polish elevator door”.

The maths is unsubtle. There are only a few free knobs — phase durations, stopping distance, body masses. The references are published, peer-reviewed, and consistent across sources. The Kirchhoff plate equation has been solved since the 1850s. None of this is novel science.

What is novel, maybe, is that all of it now fits in a few hundred lines of Python in a repository anyone can clone. The lower-bound argument lives in compute_scenario. The injury threshold mapping is a dict literal. The acoustic prediction is one call to predict_signature() against a DEFAULT_ELEVATOR_DOOR. Every number in this article is reproducible by running one notebook.

There is a serious version of this story and a less serious version. The serious version is that two and a half years of contested family-court litigation can be substantially clarified, if not resolved, by twenty minutes of biomechanics. The less serious version is the README of the project, which states the same finding in the form a friend would tell you at a bar.

Both versions agree on the headline: the laws of physics would have to take an unscheduled coffee break for the events to unfold as described.

Limitations

The model is honest about what it does not establish:

  • Whether the alleged motion actually occurred (only kinematic feasibility under the most charitable model)
  • Specific injury patterns beyond rough AIS-band mapping
  • Whether the audio recording’s apparent absence of an impact signature is itself confirmation, vs other explanations (compression artefacts, microphone position, ambient noise floor)
  • Soft-tissue injury biomechanics in any detail
  • Anything outside the physics of the alleged motion — emotional dynamics, motivations, the wider parental-alienation context that prompted the false accusation in the first place

This is just one tool in a larger picture. But it is a tool that does not lie, does not have an interest, and does not get tired.

If you want to look

The repository is public at github.com/stellarshenson/henryk-simulations. The rendered simulation is on YouTube. The full forensic deconstruction in long form lives in docs/incident_analysis.md. Names of everyone involved are anonymised.

Clone it. Change the masses. Change the corridor width. Change the phase durations. See if you can find a configuration that makes the impact pattern reasonable.

I could not.

Argue with me. Please.

This is a public repository, an executable notebook, and a single nested PARAMS dict. The whole pipeline takes about twelve seconds to re-run. If you think the physics is wrong, the methodology is wrong, the assumptions are uncharitable, or the conclusion is overstated, break it open and tell me why. Issues and pull requests on the repo. Alternative reconstructions, better stopping distances, different radiation efficiencies, alternative readings of the testimony - all welcome, all for the benefit of science.

If you find a configuration that makes the impact plausible and consistent with a single localised shoulder bruise, I owe you a coffee. If you find one that also explains the witness reporting no acoustic event at 1.5 m from a steel door, I owe you two coffees and we publish a joint follow-up.

Full bibliography. Each entry: author / year, full title, venue, key value used in this and other documents in the project. Have fun.

  1. Daams, B.J. (1994).
    Human force exertion in user-product interaction.
    Delft University Press, Delft.
    Key value: two-arm peak push force standing 800 ± 200 N; single-arm dominant 400 ± 100 N.
  2. Mital, A. & Kumar, S. (1995).
    Human muscle strength definitions, measurement, and usage. Part I — Guidelines for the practitioner.
    International Journal of Industrial Ergonomics, 16(4), 237–256.
    Key value: confirms Daams push-force budget; sustained 30 s force ≈ 60% of peak.
  3. Mero, A., Komi, P.V. & Gregor, R.J. (1992).
    Biomechanics of sprint running. A review.
    Sports Medicine, 13(6), 376–392.
    Key value: peak horizontal CoM acceleration 3.0 ± 0.8 m/s² (recreational male) and 5.0 ± 0.5 m/s² (elite sprinter) within first 0.2 s of a standing start.
  4. Cross, R. (2004).
    Physics of overarm throwing.
    American Journal of Physics, 72(3), 305–312.
    Key value: 5 kg overhand-throw release velocity ~8 ± 2.5 m/s; release kinetic energy ~160 ± 80 J.
  5. Hodgson, A.J., Lewis, J. & Drury, C.G. (2008).
    A turning-while-walking task with cognitive load.
    Applied Ergonomics, 39(3), 386–396.
    Key value: peak yaw angular velocity during voluntary 180° standing pivot 3.5 ± 1.0 rad/s.
  6. Plagenhoef, S., Evans, F.G. & Abdelnour, T. (1983).
    Anatomical data for analyzing human motion.
    Research Quarterly for Exercise and Sport, 54(2), 169–178.
    Key value: whole-body yaw moment of inertia about CoM for a 75 kg standing male 1.5 ± 0.4 kg·m².
  7. Marteniuk, R.G., MacKenzie, C.L. & Leavitt, J.L. (1990).
    Functional relationships between grasp and transport components in a prehension task.
    Human Movement Science, 9(2), 149–176.
    Key value: hand peak velocity during goal-directed reach 2.5 ± 0.8 m/s (arm extension 0.4–0.7 m).
  8. Viano, D.C. (1989).
    Biomechanical responses and injuries in blunt lateral impact.
    SAE Transactions 98, Paper 892432, 1690–1719.
    Key value: Abbreviated Injury Scale (AIS) thoracic mapping; peak impact force ≥ 12 kN corresponds to AIS 5+ (critical / life-threatening: flail chest, serious organ injury).
  9. Cavanaugh, J.M. (1989).
    The biomechanics of thoracic trauma.
    In Nahum, A.M. & Melvin, J.W. (eds.), Accidental Injury — Biomechanics and Prevention, Springer-Verlag, New York, pp. 362–390.
    Key value: force, deflection, and viscous-criterion thresholds for rib fracture and pulmonary contusion.
  10. Stapp, J.P. (1971).
    Voluntary human tolerance levels.
    In Gurdjian, E.S., Lange, W.A., Patrick, L.M. & Thomas, L.M. (eds.), Impact Injury and Crash Protection, Charles C. Thomas, Springfield IL, pp. 308–349.
    Key value: whole-body deceleration tolerance bands; documented survival at ~45 g; ≥ 100 g typically lethal.
  11. Eiband, A.M. (1959).
    Human tolerance to rapidly applied accelerations: a summary of the literature.
    NASA Memorandum 5–19–59E, National Aeronautics and Space Administration, Washington DC.
    Key value: duration-vs-magnitude survival envelopes for whole-body acceleration (the “Eiband curve”).
  12. Sturdivan, L.M., Viano, D.C. & Champion, H.R. (2004).
    Analysis of injury criteria to assess chest and abdominal injury risks in blunt and ballistic impacts.
    Journal of Trauma, 56(3), 651–663.
    Key value: kinetic-energy → injury correlation for blunt thoracic impact; ≥ 500 J → serious thoracic injury, organ contusion.
  13. Viano, D.C. & Lau, I.V. (1985).
    Thoracic impact: a viscous tolerance criterion.
    SAE Paper 851687, Tenth Experimental Safety Vehicle Conference, Oxford, England.
    Key value: velocity-injury bands assuming compliant chest deformation (soft-target geometry, 5–10 cm chest compression).
  14. Kemper, A.R. et al. Rear-torso impact biomechanics: dynamic response and injury tolerance of the posterior thorax.
    Journal of Biomechanics (Elsevier), article PII S0021929015003772.
    https://www.sciencedirect.com/science/article/abs/pii/S0021929015003772
    (Publisher page returned HTTP 403 to automated fetch; full citation per database record.)
    Key value: experimental rear-torso loads of 6.9–10.5 kN associated with costo-vertebral and costo-transverse joint injuries, thoracic spine ligament damage, and rib fractures.
  15. Brown, R. & Lefferdo, J. Thorax Injury Biomechanics.
    Chapter in Grant’s Atlas of Anatomy (13th ed.), Lippincott Williams & Wilkins, 2013.
    https://musculoskeletalkey.com/thorax-injury-biomechanics/
    Key value: lateral thorax impacts of 1.6–1.9 kN at ~4.3 m/s (15.5 kph) producing 4–13 rib fractures; chest deformation 51–66 mm associated with severe osseo-ligamentous injury.
  16. American College of Radiology Committee on Appropriateness Criteria (2018).
    ACR Appropriateness Criteria® — Rib Fractures.
    https://acsearch.acr.org/docs/69450/Narrative/
    Key value: clinical-guideline reference confirming multiple rib fractures as the typical consequence of high-energy blunt thoracic trauma.
  17. Chalmers University of Technology research portal, publication 522249.
    Rib fracture probability as a function of impact velocity in blunt thoracic trauma.
    https://research.chalmers.se/publication/522249/file/522249_Fulltext.pdf
    (Underlying PDF binary; title inferred from content described in ./impact_analysis.md source list.)
    Key value: significant rise in multiple-rib-fracture probability at impact velocities ~20 kph.

Konrad Jelen is a data scientist and CTO specialising in AI solutions for manufacturing, finance and market research, and the alienated father of his son Henry


The 3-Second Throw That Couldn’t Happen — a legal science story was originally published in Towards AI on Medium, where people are continuing the conversation by highlighting and responding to this story.

Leave a Comment

Your email address will not be published. Required fields are marked *

Scroll to Top