31:T58f,forceSim = { const gravity = 0.00005; // Regulates mechanics speed const planet = { id: 'planet', initialState: { y: -height / 4, vx: 0.5 * height * Math.sqrt(gravity) } }; const moon = { id: 'moon', initialState: { y: -height * 3/8, vx: 0.25 * height * Math.sqrt(gravity * 100) } }; const bodies = [{ id: 'sun', initialState: { x: 0, y: 0 } }, planet, moon], links = [ { source: 'planet', target: 'sun-anchor', attraction: gravity }, { source: 'moon', target: 'planet-anchor', attraction: gravity * 100 } // Higher attraction = more revolutions ]; // Generate body anchors const anchors = bodies.map(body => { return { id: body.id + '-anchor', body: body } }); const forceSim = d3.forceSimulation() .alphaDecay(0) .velocityDecay(0) .nodes([...bodies, ...anchors]) .force("gravitate-around", d3.forceLink(links) .id(d => d.id) .distance(0) .strength(link => link.attraction || 0.1) ) .on('tick', () => { mutable planetPos = [planet.x, planet.y]; mutable moonPos = [moon.x, moon.y]; fixAnchors(); }); // Set initial states bodies.forEach(body => { for (let k in body.initialState) { body[k] = body.initialState[k]; } }); function fixAnchors() { anchors.forEach(anchor => { anchor.fx = anchor.body.x; anchor.fy = anchor.body.y; }); } return forceSim; }32:[] 2f:{"title":"ForceLink based Hierarchical Orbits","cells":[{"id":0,"value":"# ForceLink based Hierarchical Orbits\n\nAttempt at simulation of multiple particles orbitting around each other in a hierarchical system (e.g. sun > planet > moon).\n\nUses [D3 force-directed](https://github.com/d3/d3-force) physics engine's [forceLink](https://github.com/d3/d3-force#links) to connect particles to their moving gravitational centers via fixed position anchor nodes.\n\nThe system fails to reach equilibrium due to the \"spring\" nature of forceLink, in which contrary to gravitational forces the intensity decreases with proximity, causing bodies to inevitably fall onto their attraction center. The [forceMagnetic](https://github.com/vasturiano/d3-force-magnetic) force is a better fit for modelling this type of motion.\n\nCompare with the [forceMagnetic version](forcemagnetic-based-hierarchical-orbits).","pinned":false,"hidden":false,"mode":"md"},{"id":10,"value":"\n${trailsCanvas}\n","pinned":false,"hidden":false,"mode":"html"},{"id":20,"value":"mutable planetPos = [0, 0]","pinned":false,"hidden":false,"mode":"ojs"},{"id":23,"value":"mutable moonPos = [0, 0]","pinned":false,"hidden":false,"mode":"ojs"},{"id":29,"value":"trailsCanvas = html``","pinned":false,"hidden":false,"mode":"ojs"},{"id":37,"value":"addTrailDot(planetPos)","pinned":false,"hidden":false,"mode":"ojs"},{"id":110,"value":"addTrailDot(moonPos)","pinned":false,"hidden":false,"mode":"ojs"},{"id":113,"value":"addTrailDot = pos => {\n d3.select(trailsCanvas).append('circle')\n .attr('id', 'trail')\n .attr('r', 1.5)\n .attr('cx', pos[0])\n .attr('cy', pos[1])\n .transition().duration(5000)\n .style('opacity', 0)\n .remove();\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":18,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":13,"value":"**Settings**","pinned":false,"hidden":false,"mode":"md"},{"id":15,"value":"height = 400","pinned":false,"hidden":false,"mode":"ojs"},{"id":5,"value":"**Dependencies**","pinned":false,"hidden":false,"mode":"md"},{"id":7,"value":"d3 = require('d3@7')","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W32","stdlib":"1","titling":"auto"}