{"componentChunkName":"component---src-templates-blog-post-js","path":"/blog/making-a-verlet-physics-engine-in-javascript","result":{"data":{"markdownRemark":{"excerpt":"Have you ever wondered if you can make your own physics engine in JavaScript? If so, you have come to the right place. We are going to build…","html":"

Have you ever wondered if you can make your own physics engine in JavaScript? If so, you have come to the right place. We are going to build a Physics engine from scratch in JavaScript.

\n

Before we start, I should mention that this tutorial assumes you have a good understanding of Vectors.\nDon’t worry if you do not yet have this understanding — Vectors are simple: get the Vector.js.

\n

What Is Verlet Physics?

\n

According to Wikipedia:

\n
\n

Verlet integration is a numerical method used to integrate Newton’s equations of motion. It is frequently used to calculate trajectories of particles in molecular dynamics simulations and computer graphics. The algorithm was first used in 1791 by Delambre and has been rediscovered many times since then, most recently by Loup Verlet in the 1960s for use in molecular dynamics.

\n
\n

In simple terms, Verlet physics is just a system of connected dots.

\n

In a Verlet system, we have 2 main components:

\n
    \n
  1. Points (dots)
    2. \n
  2. Constraints (sticks)
    3. \n
\n

\n \n \n \n \n \n \n
A Box. Connected With Dots And Sticks
\n

\n

Dots

\n

The dots have basic physics applied to them.\nWe have to keep track of the dots’ current and old positions to add the physics behavior to them — you'll see this when we implement it.

\n
// Dot.js\nclass Dot {\n  constructor(x, y) {\n    this.pos = new Vector(x, y)    this.oldpos = new Vector(x, y)    this.friction = 0.97\n    this.groundFriction = 0.7\n    this.gravity = new Vector(0, 1)\n    this.radius = 5\n    this.color = \"#e62a4f\"\n    this.mass = 1\n  }\n}
\n

We have the basic setup, now let's render the dots and make them move.

\n
// Dot.js\nupdate() {\n  let vel = Vector.sub(this.pos, this.oldpos);  vel.mult(this.friction);\n  // if the point touches the ground set groundFriction\n  if (this.pos.y >= CANVAS_HEIGHT - this.radius && vel.magSq() > 0.000001) {\n    var m = vel.mag();\n    vel.x /= m;\n    vel.y /= m;\n    vel.mult(m * this.groundFriction);\n  }\n  this.oldpos.setXY(this.pos.x, this.pos.y);  this.pos.add(vel);  this.pos.add(this.gravity);\n}
\n

The update function will update the position and handle the physics of the dot.

\n

To do Verlet integration, we have to calculate velocity based on dots old position.

\n

In the first line, we are subtracting the current position from the old position to get the desired velocity. After calculating the velocity, we apply friction to the dots, so they come to rest instead of sliding forever.

\n

Then, we update the old position by saying this.oldpos.setXY(this.pos.x, this.pos.y) and add the velocity and gravity to the position.

\n

We also want to make them stay inside the canvas, so we have to add some checks. We will also add another function: constrain():

\n
// Dot.js\nconstrain() {\n  if (this.pos.x > CANVAS_WIDTH - this.radius) {\n    this.pos.x = CANVAS_WIDTH - this.radius;\n  }\n  if (this.pos.x < this.radius) {\n    this.pos.x = this.radius;\n  }\n  if (this.pos.y > CANVAS_HEIGHT - this.radius) {\n    this.pos.y = CANVAS_HEIGHT - this.radius;\n  }\n  if (this.pos.y < this.radius) {\n    this.pos.y = this.radius;\n  }\n};
\n

Let’s add the render method too:

\n
// Dot.js\nrender() {\n  ctx.beginPath();\n  ctx.fillStyle = this.color;\n  ctx.arc(this.pos.x, this.pos.y, this.radius, 0, Math.PI * 2);\n  ctx.fill();\n  ctx.closePath();\n}
\n

Setting up:

\n
// index.js | setup\nlet canvas = document.getElementById(\"c\")\nlet ctx = canvas.getContext(\"2d\")\nlet CANVAS_WIDTH = window.innerWidth\nlet CANVAS_HEIGHT = window.innerHeight\ncanvas.width = CANVAS_WIDTH\ncanvas.height = CANVAS_HEIGHT\n\nlet dots = []\n\nfor (let i = 0; i < 50; i++) {\n  dots.push(\n    new Dot(Math.random() * CANVAS_WIDTH, Math.random() * CANVAS_HEIGHT)  )\n}\n\nfunction animate() {\n  ctx.clearRect(0, 0, CANVAS_WIDTH, CANVAS_HEIGHT)\n\n  for (let d of dots) {    d.update()    d.constrain()    d.render(ctx)  }\n  requestAnimationFrame(animate)\n}\nanimate()
\n

We added a lot more dots in random positions and then update() constrain() render() them.

\n

Let's see what that looks like:

\n

\n

Excellent — it's just a start, but finally, we have something.

\n

Now we are going to add the sticks!

\n

Sticks

\n

Sticks are the core of Verlet physics. Sticks make sure that the dots don’t get too far or too close to each other — they constrain dots to a certain distance.

\n

\n \n \n \n \n \n \n
https://slsdo.github.io/blob-family/#constraint
\n

\n

Stick.js class is relatively simple. It will take two Dots as an argument and a length. But if the length is not defined, we will calculate the length based on the dot’s position.

\n
// Stick.js\nclass Stick {\n  constructor(p1, p2, length) {\n    this.startPoint = p1    this.endPoint = p2    this.stiffness = 2\n    this.color = \"#f5476a\"\n\n    // if the length is not given then calculate the distance based on the position\n    if (!length) {\n      this.length = this.startPoint.pos.dist(this.endPoint.pos)    } else {\n      this.length = length\n    }\n  }\n}
\n

Now let's add the actual core of the algorithm, This resolves and updates the dot’s position based on the stick’s distance, ultimately constraining it to a certain distance from all other dots.

\n
// Stick.js\nupdate() {\n  // calculate the distance between two dots\n  let dx = this.endPoint.pos.x - this.startPoint.pos.x;\n  let dy = this.endPoint.pos.y - this.startPoint.pos.y;\n  // pythagoras theorem\n  let dist = Math.sqrt(dx * dx + dy * dy);\n  // calculate the resting distance betwen the dots\n  let diff = (this.length - dist) / dist * this.stiffness;\n\n  // getting the offset of the points\n  let offsetx = dx * diff * 0.5;\n  let offsety = dy * diff * 0.5;\n\n  // calculate mass\n  let m1 = this.startPoint.mass + this.endPoint.mass;\n  let m2 = this.startPoint.mass / m1;\n  m1 = this.endPoint.mass / m1;\n\n  // and finally apply the offset with calculated mass\n  if (!this.startPoint.pinned) {\n    this.startPoint.pos.x -= offsetx * m1;\n    this.startPoint.pos.y -= offsety * m1;\n  }\n  if (!this.endPoint.pinned) {\n    this.endPoint.pos.x += offsetx * m2;\n    this.endPoint.pos.y += offsety * m2;\n  }\n}
\n

Okay, I think we are good to go! Let’s add the render function:

\n
// Stick.js\nrender(ctx) {\n  ctx.beginPath();\n  ctx.strokeStyle = this.color;\n  ctx.moveTo(this.startPoint.pos.x, this.startPoint.pos.y);\n  ctx.lineTo(this.endPoint.pos.x, this.endPoint.pos.y);\n  ctx.stroke();\n  ctx.closePath();\n}
\n

Setting up:

\n
// index.js | setup\nlet canvas = document.getElementById(\"c\");\nlet ctx = canvas.getContext(\"2d\");\nlet CANVAS_WIDTH = window.innerWidth;\nlet CANVAS_HEIGHT = window.innerHeight;\ncanvas.width = CANVAS_WIDTH;\ncanvas.height = CANVAS_HEIGHT;\n\n// arrays\nlet dots = [];\nlet sticks = [];\n// forming a BOX\ndots.push(new Dot(100, 100, (Math.random() - 0.5) * 100.0); // x, y, vx, vy\ndots.push(new Dot(200, 100));\ndots.push(new Dot(200, 200));\ndots.push(new Dot(100, 200));\n\n// sticks\nsticks.push(new Stick(dots[0], dots[1]))sticks.push(new Stick(dots[1], dots[2]))sticks.push(new Stick(dots[2], dots[3]))sticks.push(new Stick(dots[3], dots[0]))sticks.push(new Stick(dots[3], dots[1]))\nfunction animate() {\n  ctx.clearRect(0, 0, CANVAS_WIDTH, CANVAS_HEIGHT);\n\n  for (let d of dots) {\n    d.update();\n    d.constrain();\n    d.render(ctx);\n  }\n  for (let s of sticks) {    s.update();    s.render(ctx);  }\n  requestAnimationFrame(animate);\n}\nanimate();
\n

Let’s see the results:

\n

\n

As you can see, we have a box! Well, a box that acts like it’s made of Jello, That’s because, in a single frame, one update per stick is not enough to make the dots rest at their length. We have to iterate the simulation as many times as we can — the more the iterations, the more the rigid box will be.

\n

Adding these lines to the existing code will make the box rigid:

\n
// index.js | setup\nconst ITERATION = 100 // max physics iterations per frame\nfunction animate() {\n  ctx.clearRect(0, 0, CANVAS_WIDTH, CANVAS_HEIGHT)\n\n  // iterate over the simulation\n  for (let i = 0; i < ITERATION; i++) {    for (let d of dots) {\n      d.constrain()\n    }\n    for (let s of sticks) {\n      s.update()\n    }\n  }\n\n  for (let d of dots) {\n    d.update()\n    d.render(ctx)\n  }\n  for (let s of sticks) {\n    s.update()\n    s.render(ctx)\n  }\n\n  requestAnimationFrame(animate)\n}\nanimate()
\n

\n

Looks amazing, isn’t it?

\n

Entity.js: managing dots and sticks in one place

\n

Okay, now we have Dot.js and Stick.js. Both are working nicely, but the problem is we don't have much control over how we use them.

\n

We will create an Entity Class that will easily handle the Updates and Renders of the Verlet Object. I'm going to paste the whole code here — it’s nothing special:

\n
// Entity.js\nclass Entity {\n  constructor(iterations) {\n    this.dots = []\n    this.sticks = []\n    this.iterations = iterations || 16\n  }\n\n  addDot(x, y, vx, vy) {\n    this.dots.push(new Dot(x, y, vx, vy))\n  }\n\n  addStick(p1, p2, length) {\n    this.sticks.push(new Stick(this.dots[p1], this.dots[p2], length))\n  }\n\n  updatePoints() {\n    for (let i = 0; i < this.dots.length; i++) {\n      this.dots[i].update()\n    }\n  }\n\n  updateSticks() {\n    for (let i = 0; i < this.sticks.length; i++) {\n      this.sticks[i].update()\n    }\n  }\n\n  updateContrains() {\n    for (let i = 0; i < this.dots.length; i++) {\n      this.dots[i].constrain()\n    }\n  }\n\n  renderPoints(ctx) {\n    for (let i = 0; i < this.dots.length; i++) {\n      this.dots[i].render(ctx)\n    }\n  }\n  renderSticks(ctx) {\n    for (let i = 0; i < this.sticks.length; i++) {\n      this.sticks[i].render(ctx)\n    }\n  }\n\n  update(ctx) {\n    this.updatePoints()\n    for (let j = 0; j < this.iterations; j++) {\n      this.updateSticks()\n      this.updateContrains()\n    }\n    this.renderPoints(ctx)\n    this.renderSticks(ctx)\n  }\n}
\n

Using the Entity Class:

\n
// index.js | setup\n...\nlet box = new Entity(100)\n// x, y, vx, vy added random velocity in first dot to make the box rotate a little bit\nbox.addDot(100, 100, (Math.random() - 0.5) * 100.0)\nbox.addDot(200, 100)\nbox.addDot(200, 200)\nbox.addDot(100, 200)\n\nbox.addStick(0, 1)\nbox.addStick(1, 2)\nbox.addStick(2, 3)\nbox.addStick(3, 0)\nbox.addStick(3, 1)\n...
\n

Yes, this is looking pretty clean and manageable!

\n

Let's take a look at the final result:

\n

\n

Now we can create lots of fun Verlet Shapes with this Entity.js class!

\n

Verly.js: A physics engine I wrote

\n

Verly.js is a Robust Verlet physics engine I wrote. It has many cool features:

\n
    \n
  1. Attraction — Repulsion behavior
    2. \n
  2. Basic Shapes
    3. \n
  3. Box, Hexagon, Cloth, Rope, Ragdoll
    4. \n
  4. Cloth Tearing
    5. \n
  5. Typography and Text
    6. \n
\n

With Verly.js you can create a tearable cloth in just 15 lines of code:

\n
// Verly.js | demo\nlet canvas = document.getElementById(\"c\")\nlet ctx = canvas.getContext(\"2d\")\nlet width = (canvas.width = 800)\nlet height = (canvas.height = 600)\n\nlet verly = new Verly(16)let cloth = verly.createCloth(150, 150, 250, 250, 15, 2)\nfunction animate() {\n  ctx.clearRect(0, 0, width, height)\n  verly.update()  cloth.tear(100) // tear threshold  requestAnimationFrame(animate)\n}\nanimate()
\n

Demo: https://anuraghazra.github.io/Verly.js/

\n

Source-code: https://github.com/anuraghazra/Verly.js

\n

Examples: https://anuraghazra.github.io/Verly.js/examples

\n

Take a look at some of the examples in my Verly.js Physics Engine, which I created some time ago. I added almost everything in that engine.

\n

\n \n \n \n \n \n \n
Spherical Shapes anuraghazra.github.io
\n
\n
\n \n \n \n \n \n \n
RagDoll anuraghazra.github.io
\n
\n
\n \n \n \n \n \n \n
Cloth anuraghazra.github.io
\n

\n

Codepen examples:

\n

\n

Verly.js’s API is easy to use and flexible because of its Entity Component Structure.

\n

Thanks for reading — I hope you learned something!

\n
\n

Other places to learn about Verlet physics:

\n
\n

Keith Peter’s CodingMath Verlet Physics Videos.

\n
\n
\n

Amazing Article At DataGenetic

\n
","timeToRead":9,"id":"5e71fd2a-4c9c-58ca-9a5d-dd1700b6cf67","frontmatter":{"date":"August 04, 2019","title":"Making a Verlet Physics Engine in Javascript"}}},"pageContext":{"slug":"/blog/making-a-verlet-physics-engine-in-javascript"}},"staticQueryHashes":["1033876704","3435786681","550521386"]}