Saturday, 5 October 2013

Mendel90 design improvements

I temporarily suspended kit production for a few weeks in order to have some time to make some improvements to the design. The following improvements apply to kits currently being shipped.

Hot End

It came to light during the hot summer that the top of the J-Head can get too hot for extruding PLA. The problem seems to manifest itself when extruding slowly because it takes longer for the filament to pass through the hex grub screw (which is the Achilles heel of the J-Head because it is wider than the rest of the filament path) and so has enough time to soften and buckle. Ironically it is also more of a problem when extruding at lower temperatures because the filament gets harder to push and so buckles at a lower temperature.

On investigation I found what I had suspected for a long time: that the heat at the top is mainly due to convection rather than conduction. This is easily demonstrated by running the extruder on its side. A lot of people add a second fan but a simple solution is to insulate the heater block. I now wrap it in two turns of silicone "self amalgamating" or "self fusing" tape, which is rated for 260°C. There are several brands but I use E-Z Fuse Tape which is 25mm wide. I stretch 110mm of it around the bottom of the PEEK, just below the bottom vent slots. I then trim it flush with the bottom of the heater block with a pair of scissors.


This makes enough difference to allow it to extrude very slowly without jamming. To cool it further I added a small hole to the fan duct.

A small airflow makes a big difference because it disrupts the convection but it is only active after the first layer when the fan is on.

It is possible to drill the hole to upgrade an old fan duct but straight drilling will spit the layers apart as it is a very thin wall. The way I did it was to melt a hole using a soldering iron with a conical bit set to 200°C. At that temperature you can wipe the ABS off the bit before it has chance to burn and make a mess of it. I then clean up the hole by drilling the melted ring it with a 4.5mm drill and then trim it flush with a penknife.

Another mod I made was to apply a bit more pressure to the idler by adding washers to each end of the springs. The same effect can be achieved by tightening the screws so they go about 1.5mm past the back of the bearing block but I like them flush with the block for an easier reference point.


These three little changes made a big difference and I can now extrude PLA at 185°C, which gives better print quality. I have updated the android gcode I include on the SD card and I have put the new Skeinforge settings on Github.

Settings

Other improvements I made to the settings are: -
The PLA settings were erroneously under the ABS profile, sorry for the confusion. The un-configured profiles have been removed and two new ones: PLA0.3 and PLA0.2 have been added. The android sample is sliced with the 0.3mm layer profile because it is easier to get the Z calibration right with the highest layer height possible with PLA. The 0.2mm profile gives a much nicer print.
The ooze free unattended start was not very reliable with PLA. Lowering the temperature helped but I also improved the extruder priming by extruding a very fat line along the front edge of the bed instead of a blob.
I have come to realise that you get better prints if you do the outlines at the same speed as the infill and looking back I think that is what I have seen commercial machines do. The old profiles had the outlines at 25mm/s and infill at 50mm/s. I now do both at 40mm/s.
The reason it works better is because when the extruder rewinds at the end of one run and then fast forwards the same distance to start another it will restore the same nozzle pressure as it had on the previous run. That means the initial flow rate will be the same, so if the speed has changed it will be too fast or too slow. This leads to blobs on the outline after doing faster infill and gaps at the start of infill after doing slow outlines. 
Another advantage is with a constant flow rate the plug of plastic that forms at the hot end transition zone will be a constant length. If the speed increases the length reduces, which means less of the barrel has plastic expanded to the barrel's diameter and more has it at the original filament diameter. This implies there is less plastic in the barrel, so an excess must have been extruded while the plug was shortening. Conversely when the extrusion rate drops the plug extends so too little plastic will be extruded during that time.
Yet another advantage of using constant flow rate is that the die swell will be constant, so the plastic is always stretched by the same amount. When outlines are done slower the die swell is less, so they are stretched less, and so don't span gaps as well as the infill.
I have set the outline flow rates slightly lower than the infill. Skeinforge's volumetric calculations fill a rectangle equal to the filament width and layer height. This is OK for the infill but the outline does not have a rectangular profile, it has rounded sides, so it comes out a little wider than intended. The error is less when the filament path is wide compared to its height. The correct value is 1 + (π/4 - 1) / (w / h). Using this formula I get dimensionally correct prints instead of outlines that are too big and holes that are too small.
PLA prints better when cooled by the fan but you don't want it cooling the print surface when the first layer is being laid down. Ironically print cooling seems to be more important when extruding at 185°C than at 220°C for reasons I can't yet explain. To be able to automatically turn the fan on at a specified layer I had to incorporate an update to the Cool plugin by DeuxVis and Joost-b plus my own little bug fix. I have put the patched version of Skeinforge50 that I include on the SD card onto Github.

Firmware

I made three small tweaks to the version of Marlin I distribute: -
I reduced the XY acceleration from 4000 to 2000 as the Y axis can lose steps due to resonance in high frequency infill otherwise.
I found that the Y axis would overshoot when homing if the X axis hits the end stop first. This was because the axis moves 40% faster on it own than when moving diagonally. I reduced the homing speeds but multiplied them by 1.4 when moving diagonally. The net effect is diagonal homing is the same speed as before but it now doesn't speed up when one axis stops.
I updated the PID constants for J-Head MK5B as I hadn't changed them from MK4B. It didn't seem to make much difference.

X ends

A few people started reporting the bar clamps on the X motor bracket and the X idler bracket had cracked. I am not sure why this became a problem after hundreds of kits. Some reported they could not get enough tension to pull the belt tight, so perhaps the belt got stiffer, or the plastic parts weaker. I built a kit myself to investigate using a problem belt that one of my customers returned. It was inconclusive because I managed to get the machine printing well without breaking the clamps but I think some belts are stiffer than they were in the first kits, despite coming from the same manufacturer and I also think the plastic parts might not be as strong as they used to be due to ABS supply changing. I tackled the problem from both sides by immediately switching to more flexible belt and then redesigning the X end clamps.

I switched to Brecoflex® belt which has more thinner wires in it and also boasts that it has a better arrangement of strands within each wire. This makes it more flexible without being less strong. The wires are also parallel, whereas cheaper belts have them in a spiral so they meet the edge every so often.

The clamps now only have plastic in compression so are much stronger. They do use a lot more fasteners but because the screws are accessible from above they no longer need to have hex heads, so that is one less type of screw required.


The nut trap for the lead screw is now at the top (don't be fooled by the hexagonal hole at the bottom, that is just clearance to allow the nut to be fitted to the leadscrew first). This greatly increases the distance of the nut from the Z coupling when the Z axis is at the bottom. When the nut was at the bottom it got to within about 10mm so any eccentricity was magnified by about 20 at the top of the leadscrew.

Because holes can't be printed in mid air there are one layer thick support membranes at the bottom and the top of the nut trap. This creates a completely enclosed void which revealed a bug in OpenScad. It generates an STL with the internal faces facing the wrong way, i.e. outwards instead of inwards. Fortunately Skeinforge does not care and slices it correctly but Slic3r doesn't.

Another little tweak to the design was to angle the hole for the idler axle. It is a clearance hole so the belt tension pulls the front of the bolt to the right, which tends to make the belt run to the front of the idler pulley, requiring a shim under the back washer to correct it. The new design angles the bolt hole slightly to correct for this so doesn't require a shim. Here is a cross section showing that when the axle bolt is orthogonal the front rests against the right side of the hole and the nut at the back rests on the left side of its trap.


The X motor bracket has been simplified and made much easier to print. The previous design fastened the motor with three screws at the front and used a box construction to be stiff enough not to bend under the belt tension. The new design has one screw in the back of the motor so that only two are needed at the front leading to a much smaller bracket.


The length of stepper motors has a poor tolerance of ±1mm because it is made from a stack of laminations pressed together. For this reason I couldn't use a normal screw through the plastic. Instead the screw head is packed out with as many washers needed to fill the gap between the motor and the bracket. This will be one if the motor is top tolerance and five if it is bottom tolerance. The hole in the bracket is a snug fit for the head of the screw so it acts as a dowel.


Like the idler, the clearance screw holes for the motor would allow the front to be pulled to the left causing the belt to run to the front of the pulley. The motor holes have been offset by half their clearance to avoid this.

Z couplings

The original Z couplings worked well on all of my machines until I built another one that produced objects with some Z banding. Investigation with a dial gauge showed that the Z position was not increasing linearly leading to uneven layer heights. I found that the lead screw was eccentric leaving the coupling due to the bore of the PVC tubing not being central. Also PVC is more plastic than elastic, so when it is squeezed by the coupling it deforms permanently and does not give much flexibility. The couplings where designed with a gap so it was possible to over tighten them leading to too much constraint. The top needed to be fully closed because the studding has a smaller outside diameter than a normal M6 bolt would have.

To address these issues I switched to neoprene rubber tubing that seems to have the bore more central and is also a lot more flexible. I changed the plastic parts so they can be closed completely to give the correct pressure rather than relying on a gap being maintained. I also made the screw positions diagonally opposite which makes them a bit easier to tighten up.


The only downside to the neoprene is that it doesn't give enough grip to be able to turn the motor manually when it is enabled. To allow bed levelling I changed the Pronterface buttons for back left and back right to disable the motors. The new config file is on Github.

X carriage fan bracket

The length of the J-Head MK5B is incorrect in the Mendel90 model. This makes the fan duct too low and it hits the clips on the bed. A simple solution is to make the slots in the fan bracket deeper so it can be fitted a bit higher. I also added strengthening ribs as the old bracket allowed the duct to droop.

Tube end caps

I reversed the direction of the screws that hold the aluminium tubes onto the base. The nut is now trapped inside the tube by the end cap. This means the tube no longer needs the large holes for the screw head, which greatly reduces the amount of machining I need to do to make them.


The two extra washers are needed at the front to avoid using an additional screw length. The two rear fixing blocks are now different from the others because they don't have a nut trap in the middle position.

Y carriage heat shield

Instead of using aluminium foil to reflect heat I found that simply filling the gap under the bed with multiple layers of corrugated cardboard allows the bed to reach higher temperatures.

Spool Holders

These are now assembled with the large washer on the outside of the spool rather than the inside, requiring four more small washers. This is to accommodate spools that originate from Orbitech that are slightly wider than the ones I used to get.

Kits have been shipping with all these mods for a while now but the time out taken to make the changes led to a backlog and a large waiting list which we have been working flat out to clear and hope to get back to next day dispatch this week.

Friday, 7 December 2012

Mendel90 updates

It has been a long road but I am now in a position to sell complete kits for Mendel90. My original plan was to use laser cut acrylic but the companies I got quotes from could not guarantee the holes sizes would be accurate enough to be tapped. I also found out that 3mm Dibond is stiffer than 6mm acrylic as well as being a lot lighter and cheaper. The downside is that the polyethylene core is too soft and the aluminium wall too thin to tap a thread into reliably. That meant replacing all the screws with nuts and bolts, which made the machine fiddly to put together. To fix that I reversed all the bolts and put nut traps in the plastic parts. The design still supports the MDF and acrylic variants with screws, but I only sell the plastic parts for these.

Picture courtesy of  Alzibiff, a Mendel90 owner who is also a photographer: thefullpicture.co.uk 
Because the underside of the base now has screw heads protruding, the machine has to have some form of feet to raise it. I added two aluminium square tubes, which also provide stiffness as they are directly under the stays. An added  bonus is that the wires for the Y axis can now run underneath the base to improve the appearance. Another advantage is that the printed cable clips can now be replaced by zip ties that go through holes in the panels.

I changed the default electronics from Sanguinololu to Melzi because it includes the SD card interface and fan drive on board and uses screw connectors rather than friction fit. These are the only connectors I have found to be reliable in the long term and it will allow me to offer a no solder version of the kit. After a couple of false starts I now have a reliable source of these and they are 100% functionally tested using the same motors as I provide in the kit.


I also changed the PSU to an ATX500 because I can buy them in the UK with CE approval and they come with integral mains inlet and switch, obviating the need for the kit builder to do their own mains wiring.  It is also shorter, allowing room for the Melzi (which is very long and thin). The downside is it needs a pair of load resistors on the 5V and 3.3V to get the 12V rail to be close to 12V under load. These are provided in the kit.


To stop the deeper power supply sticking out beyond the base I moved the right hand stay inwards. It is now directly behind the  Z bar that it braces, which is better from a structural point of view and gives maximum space for the electronics bay. A couple of plastic brackets hold the ATX PSU in place.


The build height is increased to 200mm as there is very little downside to doing so and it makes room for the Melzi electronics.

I added a spool holder which suspends the spool between the stays on 608 bearings using a couple of triangular brackets.



In order to turn the filament through 90 degrees I use a PTFE Bowden tube. That also has the benefit of removing any drag on the head when the extruder pulls filament from the reel. The only force on the head is that required to bend the loop of filament and the tube, which has only a thin wall.

The end of the tube is terminated by a printed connector with a flange that just sits on top of the extruder. That allows the extruder to reverse without having to push the filament back up the tube.

There is dust wiper on entry to the tube which consists of a block of foam squeezed into a smaller box with the filament running through a slot in the middle. It prevents dust that settles on the spool being dragged into the extruder.


I made a few tweaks to the extruder. I now use J-Head nozzles so I integrated a groove mount into the Wade's block. The nozzle is a tight fit and I press it in with printed jig and a vice with rubber jaws.


After the first five kits I found I could no longer get the J-Head MK4B and had to switch to the MK5B. That isn't long enough for a machine like Mendel90 because the bottom of the carriage ends up too close to the bed. I had to extend the bottom of the Wades' block so that it now protrudes below the carriage. In order to do that I had to make it a little slimmer to fit through the hole and that had the side effect of making the bearing housing symmetrical about the filament path. I think the reason it wasn't before is because there isn't room for the bolt heads on the motor plate side. I solved that by turning the bolts around and using captive nuts in the bearing block. That also makes the springs much easier to fit as the nut engages before the spring has to be compressed. The hob position is now further up the bolt at 25mm meaning more of the shoulder enters the bearing.

Other minor changes were that I lowered the motor a few mm so that it traps the head of the mounting screw  making it easier to fit to the carriage. I also replaced the spring that I used to retain the nut on the hobbed bolt with a lock nut and a star washer. They only need to be finger tight.

In order to be able to offer a solder-less version of the kit I made a tiny break out PCB for the extruder motor and heater connections.


I also made the socket on the end of the X cable an IDC version. That meant the pin to wire connections had to be 1:1, so I had to increase the number of pins from 9 to 15. The heater wires were previously doubled up to handle 2 Amp heaters but now I have enough pins for three wires giving about 4 Amps. There are also wires for a fan and a Z probe.

I simplified the X ribbon cable layout. Previously I had a pair of grounded wires acting as guard between the X limit switch and the noisy motor wires. Instead I simply moved the limit switch wires to the other side of the cable, where the quiet signals are. It actually makes the X end wiring neater and gives me the extra two wires for the extruder heater so that the cable remains 20 way, which is a standard size.


On the subject of ribbon cables: I increased the one for the bed from 24 way to 26 way making that a standard size also and two extra wires for the bed heater is not a bad thing.

In order to be able to print small items and items with steep overhangs in PLA I took a leaf out of Richard Gain's book and added a ducted fan to the carriage. I couldn't use his design directly because he has a longer nozzle mounted in a different orientation. This is my compressed version.


 The exit of the duct produces a ring of air directed inwards. The idea is to direct the air onto the part close to where the new plastic is being laid down without cooling the nozzle.


Because the Melzi only has a single fan output I removed the bed cooling fan but left the hole and fan guard for it. It is something I use to speed up production by cooling the bed rapidly at the end of a build, but the 80mm high airflow fans are expensive. They used to be cheap when they were used in PCs but they have all switched to quiet ones nowadays with less flow. If you want to add one you can hack a MOSFET onto the expansion port of the Melzi and control it with M42 in Marlin.

The OpenScad model now includes everything in the kit, which is everything needed to build the machine apart from some sticky tape used to secure the PTFE tubing. Most things are visible in the rendering but a few generate BOM entries only. These are things like wires and ribbon cables, which are hard to draw. Resistors, thermistors, sleeving and heat shrink are all drawn in places close to where they are used.

Having everything in the model is the only way to keep the BOM 100% accurate. I also re-structured the sub assemblies so that they reflect the order things are assembled rather than their placement in the machine. For example the Z lead nuts were in the X end assemblies because that is where they end up, but when assembling they are fitted to the leadscrews and then inserted into the X ends, so they need to be part of the Z axis assembly.

I created a detailed build manual for the kit in OpenOffice format. It links to a lot of pictures generated by the model so that they automatically update. The manual is also checked into GitHub, so there is a version matching each revision of the machine. I had started with instructions in the Reprap wiki but that soon become impractical as it can only represent a single version and all the images need to be manually updated. It can't link to images on GitHub for instance. Although the manual is for the kit version of the machine, it will also be useful to people making the other versions. The main differences are that different fasteners are used. The PDF version is here.

The first five kits were commissioned for a build a weekend hosted by the GIST lab in Sheffield during the university's Festival of the Mind. Due to various supply problems I only just managed to get all the parts together in time and then only with a lot of help from a couple of friends. I ran out of time completely to do the instructions so I stayed up all night and hacked some together. Despite that the build weekend went well.

Two of the teams completed their machines on day one and were extruding plastic. By the end of day two all the machines were completed, two teams were printing items downloaded from Thingiverse, another was extruding. Unfortunately two were held up by faulty Melzis. I hadn't tested them beforehand because I had five spares, but it turned out one of the originals had a small fault and none of the spares worked at all. Needless to say I won't be using those suppliers again!

The Derby Makers team won the "Golden Spanner Award" awarded for "the best demonstration of the Craft and Magic of Technology". Here they are printing.

Picture courtesy of Derby Maker Glyn Smith.
More pictures of the weekend here.

The build weekend was very helpful for ironing out snags. This lead to tweaks in the design and the instructions being greatly expanded.

Kit details:

Specification

Build volume 200mm x 200mm x 200mm.
Filament 3mm.
Nozzle size 0.4mm.
Footprint 465mm x 419mm.
Height 400mm, with spool 609mm.

Contents

The kit contains everything needed to get printing including ~ 50m of Faberdashery PLA. You will need a PC and some tools to put it together and run it, details at the start of the manual. Some things included in the kit are: -

  • Melzi electronics with Atmega1284P, ROSC shorted for correct microstepping, tested and programmed with Marlin firmware configured for Mendel90
  • Micro SD card for pause free printing, comes pre-loaded with all the software needed to print on Windows (Linux equivalents can be downloaded) and a USB to SD adaptor for direct connection to a PC.
  • 0.4mm J-Head MK5B hot end for 3mm filament, pre-assembled with thermistor, resistor, wires and sleeving.
  • ATX500 PSU with quiet fan and dummy load resistors.
  • USB cable.
  • Prusa MK2 heated bed with 2mm glass sheet.
  • 3mm aluminium composite panels CNC cut and drilled.
  • 5 NEMA17 43Ncm stepper motors.
  • 8mm hardened high carbon rods with h6 tolerance as required by the linear bearings.
  • T2.5 metal pulleys and polyurethane belts with steel reinforcement.
  • CNC hobbed bolt.
  • Extruder break out PCB.
  • 1% thermistors to avoid the need for temperature calibration.
The complete BOM is here.

The big advantage of the design is that the axes need no alignment to ensure they are orthogonal. The only calibration required is bed levelling, Z height and extruder flow rate. Details in the manual.

The price of the kit is £499, plus VAT in the EU, plus shipping. For availability check the forum post here.

Sunday, 25 November 2012

More accurate thermistor tables

A couple of weeks ago I wanted to create some thermistor tables for Marlin. At that time it had a copy of createTemperatureLookup.py, which I think was written by Zach Smith based on my article "Measuring temperature the easy Way". It uses the simple two constant thermistor equation based on the resistance at 25°C and beta.

The two constant formula was adequate for my own software because I have separate constants for each thermistor that I use and I calibrate them against a thermocouple at the working temperature and room temperature. Marlin however has a thermistor table for each type of thermistor, so if you use the same thermistor for the bed and the hot end they share a table. The problem with the simple equation is that beta is not very constant and depends on temperature. There are several values given on the datasheet for different temperature ranges
, none of them very applicable to our application. If you have beta correct for the hot end at say 250°C it is about 7°C out at the bed temperature, say 130°C.

I decided to make a new script which uses the three constant Steinhart–Hart equation. The graph shows the difference between the two equations over a large temperature range: -
The two constant equation is only accurate around the two temperatures the constants are calculated at (in this case 25°C and 256°C). When these are at opposite ends of the thermistors range the error in the middle is quite large.

The script I made is MakeTempTable.py. Its parameters are resistances at three temperatures. The tables is makes look like this: -

    {     344,       300     }, // r=   101 adc=  21.47
    {     369,       295     }, // r=   108 adc=  23.08
    {     397,       290     }, // r=   117 adc=  24.83
    {     428,       285     }, // r=   126 adc=  26.75
    {     461,       280     }, // r=   136 adc=  28.84
    {     498,       275     }, // r=   147 adc=  31.12
    {     538,       270     }, // r=   160 adc=  33.63
    {     582,       265     }, // r=   173 adc=  36.37
    {     630,       260     }, // r=   188 adc=  39.38
    {     683,       255     }, // r=   205 adc=  42.69
    {     741,       250     }, // r=   223 adc=  46.32
    {     805,       245     }, // r=   243 adc=  50.31
    {     875,       240     }, // r=   266 adc=  54.71
    {     953,       235     }, // r=   290 adc=  59.55
    {    1038,       230     }, // r=   318 adc=  64.88
    {    1132,       225     }, // r=   349 adc=  70.77
    {    1236,       220     }, // r=   384 adc=  77.26
    {    1351,       215     }, // r=   423 adc=  84.42
    {    1477,       210     }, // r=   466 adc=  92.32
    {    1617,       205     }, // r=   515 adc= 101.05
    {    1771,       200     }, // r=   570 adc= 110.68
    {    1941,       195     }, // r=   632 adc= 121.30
    {    2128,       190     }, // r=   702 adc= 133.01
    {    2335,       185     }, // r=   782 adc= 145.91
    {    2562,       180     }, // r=   872 adc= 160.11
    {    2811,       175     }, // r=   975 adc= 175.70
    {    3085,       170     }, // r=  1092 adc= 192.81
    {    3384,       165     }, // r=  1225 adc= 211.53
    {    3711,       160     }, // r=  1378 adc= 231.95
    {    4066,       155     }, // r=  1554 adc= 254.15
    {    4451,       150     }, // r=  1756 adc= 278.21
    {    4866,       145     }, // r=  1989 adc= 304.15
    {    5312,       140     }, // r=  2258 adc= 331.99
    {    5787,       135     }, // r=  2570 adc= 361.68
    {    6290,       130     }, // r=  2934 adc= 393.15
    {    6820,       125     }, // r=  3357 adc= 426.25
    {    7373,       120     }, // r=  3852 adc= 460.80
    {    7945,       115     }, // r=  4433 adc= 496.54
    {    8531,       110     }, // r=  5116 adc= 533.16
    {    9125,       105     }, // r=  5921 adc= 570.31
    {    9722,       100     }, // r=  6875 adc= 607.60
    {   10314,        95     }, // r=  8007 adc= 644.61
    {   10895,        90     }, // r=  9356 adc= 680.92
    {   11458,        85     }, // r= 10968 adc= 716.13
    {   11998,        80     }, // r= 12903 adc= 749.86
    {   12509,        75     }, // r= 15234 adc= 781.80
    {   12987,        70     }, // r= 18051 adc= 811.66
    {   13428,        65     }, // r= 21469 adc= 839.27
    {   13832,        60     }, // r= 25635 adc= 864.50
    {   14197,        55     }, // r= 30732 adc= 887.30
    {   14523,        50     }, // r= 36995 adc= 907.68
    {   14811,        45     }, // r= 44725 adc= 925.72
    {   15064,        40     }, // r= 54309 adc= 941.52
    {   15284,        35     }, // r= 66249 adc= 955.23
    {   15472,        30     }, // r= 81195 adc= 967.02
    {   15633,        25     }, // r=100000 adc= 977.08
    {   15769,        20     }, // r=123783 adc= 985.58
    {   15883,        15     }, // r=154025 adc= 992.71
    {   15978,        10     }, // r=192694 adc= 998.64
    {   16057,         5     }, // r=242427 adc=1003.54
    {   16121,         0     }, // r=306773 adc=1007.56


The ADC values in the table are multiplied by 16 because Marlin uses oversampling to give four more bits of precision. The old tables just multiplied the integer ADC value by 16 but I multiply it before rounding it to an integer so the table has the same precision as the oversampled ADC reading.

The script can also take ADC values as parameters instead of resistances. This allows you to calibrate a thermistor in situ. If you set the temperature to a value in an existing table and let it settle and then measure it with a thermocouple you know that the ADC value for the measured temperature is the value in the table for the set temperature. You can then produce a new more accurate table.

Two days after I put it on Github ErikZalm added a new script to the official version of Marlin to do exactly the same thing: createTemperatureLookupMarlin.py, amazing coincidence! It is different code but I think it uses exactly the same maths to find the three coefficients using simultaneous equations that I lifted from here