pylinkage 0.2.2

Build and optimize planar linkages using PSO

pylinkage

A linkage builder written in Python. This package is made to create planar linkages and optimize them kinematically thanks to Particle Swarm Optimization. It is still an early work, so it should receive great changes in the future.

Installation

Using pip

This package is in the PyPi as pylinkage, and can be installed using: pip install pylinkage

It is the recommended way downloading it, since the release versions are synced.

Setting up Virtual Environment

We provide an environment.yml file for conda. Use conda env update --file environment.yml --name pylinkage-env to install the requirements in a separate environment.

Usage

As of today, the code is segmented in three parts:

• geometry.py that module handles geometric primitives, such as circle intersections, distance claculation. It works in Euclidian space only. Aside from dist and sqr_dist functions, you might not use it directly.
• linkage.py this module describes joints and linkages
  • Due to the geometric approach, joints (instances of Joint object) are defined without links.
  • The Linkage class that will make your code shorter.
• optimizer.py proposes three optimizations based on three techniques:
  • The "exhaustive" optimization (exhaustive_optimization function) is a dumb optimization method, consisting or trying sequencially all positions. It is here for demonstration purposes only, and you should not use it if you are looking for an efficient technique.
  • The built-in Particle Swarm Optimizer (PSO). I started with it, so it offers a large set of useful options for linkage optimization. However, it is here for legacy purposes, and is much short than the PySwarms module.
  • PSO using PySwarms. We provide a wrapper function to PySwarm from ljvmiranda921, that will progressively be extended.
• visualizer.py can make graphic illustrations of your linkage using matplotlib.
  • It is also used to visualise your n-dimensional swarm, which is not supported by PySwarms.

Requirements

Python 3, numpy for calculation, matplotlib for drawing, and standard libraries. You will also need PySwarms for the optimization since the built-in PSO is deprecated and will be removed soon.

Example

Let's start with a crank-rocker four-bar linkage, as classic of mechanics.

Joints definition

The first thing you need if to define points belonging to the frame:

import pylinkage.linkage as pl

frame_first = pl.Static(0, 0)
frame_second = pl.Static(3, 0)

• O, O: "x, y" position of the first point.
• 3, 0: same for the second point.

Then we have to define at least one crank because we do a kinematic simulation.

crank = pl.Crank(0, 1, joint0=frame_first, angle=0.31, distance=1)

Here you need some explanations:

• 0, 1: x and y initial coordinates of the tail of the crank link.
• joint0: the parent Joint to link with. The pin will be created on the position of the parent, which is the head of the crank link.
• angle: the crank will rotate with this angle, in radians, at each iteration.
• distance: distance to keep constant between crank link tail and head.

Now we add a pin joint to close the kinematic loop.

pin = pl.Pivot(3, 2, joint0=crank, joint1=frame_second, distance0=3, distance1=1)

Here comes some help:

• joint0, joint1: first and second Joints you want to link to, the order is not important.
• distance0, distance1: distance to keep constant between this joint and his two parents.

And here comes the pain: Why do we specify initial coordinates 3, 2? Moreover, they seem incompatible with distance to parents/parents' positions!

• This explanation is simple: mathematically a pin joint the intersection of two circles. The intersection is often two points. To choose the strating point, we calculate both intersection (when possible), then we keep the intersection closer to the previous position as the solution.

Wait! A linkage with a single motor and only one pin joint? That doesn't make sense! :Behind the curtain, many joints are created on the fly. When you define a Crank joint, it creates a motor and a pin joint on the crank's link head. For a Pivot joint, it creates 3 pin joints: one on the position of each of its parents, and one its position, forming a reshapable triangle. This is why pylinkage is so short to write.

Linkage definition and simulation

Once your linkage is finished, you can either can the reload method of each Joint in a loop, or put everything in a Linkage that will handle this can of thinks for you.

Linkage definition is simple:

my_linkage = pl.Linkage(joints=(frame_first, frame_second, crank, pin))

That's all!

Now we want to simulate it and to get the locus of pin. Just use the step method of Linkage to make a complete rotation.

locus = my_linkage.step()

You can also specify the number of steps with the iteration argument, or subdivisions of each iteration withdt.

Let's recape.

import pylinkage.linkage as pl

# Static points in space, belonging to the frame
frame_first = pl.Static(0, 0)
frame_second = pl.Static(3, 0)
# Main motor
crank = pl.Crank(0, 1, joint0=frame_first, angle=.31, distance=1)
# Close the loop
pin = pl.Pivot(3, 2, joint0=crank, joint1=frame_second, 
               distance0=3, distance1=1)

my_linkage = pl.Linkage(joints=(frame_first, frame_second, crank, pin))

locus = my_linkage.step()

Visualization

Firsting first, you made a cool linkage, but only you know what it is. Let's add friendly names to joints, so the communication is simplified.

frame_first.name = "A"
crank.name = "B"
pin.name = "C"
frame_second.name = "D"
# Linkage can also have names
my_linkage.name = "Four-bar linkage"

Then you can view your linkage!

import pylinkage.visualizer as visu

visu.show_linkage(my_linkage)

Last recap, rearranging names:

import pylinkage.linkage as pl
import pylinkage.visualizer as visu

# Static points in space, belonging to the frame
frame_first = pl.Static(0, 0, name="A")
frame_second = pl.Static(3, 0, name="D")
# Main motor
crank = pl.Crank(0, 1, joint0=frame_first, angle=.31, distance=1, name="B")
# Close the loop
pin = pl.Pivot(3, 2, joint0=crank, joint1=frame_second,
               distance0=3, distance1=1, name="C")

# Linkage definition
my_linkage = pl.Linkage(
    joints=(frame_first, frame_second, crank, pin),
    order=(frame_first, frame_second, crank, pin),
    name="My four-bar linkage")

# Visualization
visu.show_linkage(my_linkage)

Optimization

Now, we want automatic optimization of our linkage, using a certain criterion. Let's find the four-bar linkage with the minimal distance between two points as output.

Our objective function, often called the fitness function, is the following:

def fitness_func(linkage, params):
    """Return some stuff."""
    linkage.set_constraints(*params)
    return 0
    """