Experimenteel ontwerp voor de ontwikkeling van 3D-printprocessen: van orthogonale screening en responsvlakmodellering tot kalibratie van het procesvenster

Inleiding: waarom parameteroptimalisatie steeds in cirkels blijft draaien

Bij de praktische afstelling van metaal-laserpoederbedfusie (LPBF) en polymeersinterprocessen is het meest voorkomende scenario als volgt: een engineer verhoogt het laservermogen van 250 W naar 280 W, de dichtheid stijgt met 0,3%, en het vermogen wordt verder opgevoerd. Bij 300 W ontstaat een keyhole-defect, waarna de instelling weer wordt verlaagd. Na een ronde met meer dan tien proefstukken resteert slechts een gefragmenteerde ervaringscurve; bij een andere machine of een nieuwe poederbatch moet alles opnieuw worden gedaan. Het probleem is niet dat engineers zich onvoldoende inspannen, maar dat de methode waarbij één factor per keer wordt gewijzigd, interacties niet kan verwerken. Vermogen en scansnelheid bepalen samen de volumetrische energiedichtheid, terwijl de laagdikte van het uitgespreide poeder de warmteafvoer van het smeltbad verandert. De drie factoren zijn onderling gekoppeld, waardoor het volledige beeld nooit zichtbaar wordt wanneer telkens slechts één variabele verandert.

Design of Experiments (DoE) biedt precies de technische aanpak voor dergelijke problemen met gekoppelde factoren: met een minimaal aantal proefstukken kunnen hoofdeffecten, interactie-effecten en responsvlakken gelijktijdig worden geschat, inclusief een vrij te geven procesvenster met betrouwbaarheidsinterval. Dit artikel biedt voor 3D-printprocesontwikkeling een volledig uitvoerbare methode, van orthogonale screening en responsvlakmodellering tot vensterverificatie.

1. Maak eerst onderscheid tussen drie soorten proeven: screening, modellering en verificatie

Procesontwikkeling kan niet met één enkele proef worden afgerond. Gewoonlijk verloopt zij in drie niveaus, die verschillen in doel, aantal factoren en aantal proefstukken.

Het eerste niveau bestaat uit screeningproeven, met als doel de werkelijk significante 2 tot 4 factoren te selecteren uit 6 tot 10 kandidaatfactoren. Typische kandidaatfactoren zijn laservermogen, scansnelheid, scanafstand, laagdikte, voorverwarmingstemperatuur, scanstrategie (strookbreedte), rotatiehoek tussen lagen en beschermgasstroomsnelheid. In de screeningsfase volstaan een fractioneel factoriëel ontwerp of een Plackett-Burman-ontwerp met 8 tot 12 proeven; als responsen worden alleen de globale indicatoren dichtheid en oppervlakteruwheid bekeken.

Het tweede niveau bestaat uit modelleringsproeven, met als doel een responsvlakmodel op te bouwen voor de geselecteerde 2 tot 3 kritische factoren. Veelgebruikte methoden zijn een centraal samengesteld ontwerp (CCD) of een Box-Behnken-ontwerp. Een CCD met drie factoren omvat 6 axiale punten, 8 kubuspunten en 6 herhalingen in het centrum, in totaal 20 proeven. De responsen worden uitgebreid met dichtheid, treksterkte, maatafwijking en oppervlakteruwheid Ra.

Het derde niveau bestaat uit verificatieproeven. Binnen het door het model voorspelde optimale gebied worden 3 tot 5 bevestigingspunten opnieuw geprint, waarna wordt gecontroleerd of het verschil tussen voorspelde en gemeten waarden binnen ±5% blijft. Na succesvolle verificatie worden de grenzen van het venster met 10% aangescherpt als parameterbereik voor de productieproceskaart.

2. Orthogonale proeven: met zo weinig mogelijk proefstukken een zo groot mogelijke ruimte afdekken

Wanneer het aantal factoren beperkt is en alle factoren hetzelfde aantal niveaus hebben, is een orthogonale tabel de meest economische keuze. Neem LPBF met roestvast staal 316L als voorbeeld, met de drie factoren laservermogen P, scansnelheid v en scanafstand h, elk op drie niveaus. Een volledig factoriëel ontwerp vereist 27 proeven, terwijl een L9(3^4)-orthogonale tabel slechts 9 proeven vereist.

Het concrete ontwerp is: P op 200/240/280 W, v op 800/1000/1200 mm/s en h op 0,09/0,11/0,13 mm. Na het plannen van de proeven volgens L9 wordt voor elke niveaucombinatie een kubisch proefstuk van 10×10×10 mm geprint. De dichtheid wordt gemeten met de Archimedes-methode en de oppervlakteruwheid Ra wordt tegelijkertijd geregistreerd.

Een bereikanalyse is voldoende om de factoren snel te rangschikken: bereken voor elke factor de gemiddelde respons op elk niveau; het bereik R = max(gemiddelde respons) − min(gemiddelde respons). Hoe groter R, hoe significanter de invloed van de factor. Indien R(P)=2,1%, R(v)=1,4% en R(h)=0,6%, is de volgorde van invloed vermogen > snelheid > afstand. Voer vervolgens variantieanalyse (ANOVA) uit en gebruik de F-toets om significantie te beoordelen: factoren met een F-waarde boven de kritische waarde en een p-waarde kleiner dan 0,05 gaan door naar de volgende modelleringsronde.

De beperking van orthogonale proeven is dat zij veronderstellen dat er geen interacties tussen factoren bestaan. Wanneer een sterke koppeling tussen vermogen en snelheid wordt vermoed — de volumetrische energiedichtheid E = P/(v·h·t) is zelf al een productrelatie — moet in de modelleringsfase een responsvlakontwerp worden gebruikt in plaats van orthogonale tabellen geforceerd toe te passen.

3. Responsvlakmethode: van discrete punten naar een continu oppervlak

De kern van de responsvlakmethode (RSM) is het fitten van de relatie tussen respons en factoren met een kwadratische veelterm:

Y = β0 + Σβi·xi + Σβii·xi² + Σβij·xi·xj + ε

De kruisterm βij is de gekwantificeerde uitdrukking van het interactie-effect. In een CCD-ontwerp met drie factoren is de axiale afstand α een belangrijke parameter. Voor een roteerbaar ontwerp geldt α = 2^(k/4) = 2^(3/4) ≈ 1,682, waarbij k het aantal factoren is. De axiale punten liggen in de gecodeerde ruimte op ±1,682, de kubuspunten op ±1, en het centrumpunt wordt 6 keer herhaald om de zuivere fout te schatten.

Na het fitten moeten drie controles worden uitgevoerd: ten eerste de lack-of-fit-toets; een p-waarde groter dan 0,05 betekent dat er geen significante afwijking bestaat tussen de modelvorm en de gegevens. Ten tweede moeten de determinatiecoëfficiënt R² en de aangepaste R² worden beoordeeld; het verschil tussen beide moet kleiner zijn dan 0,2, anders kan sprake zijn van overfitting. Ten derde moeten de normaliteit van de residuen en de homogeniteit van de variantie worden gecontroleerd. In de praktijk komt het vaak voor dat R² 0,95 bereikt, maar de aangepaste R² slechts 0,86 bedraagt; in dat geval moeten niet-significante termen van de derde orde worden verwijderd en moet het model opnieuw worden gefit.

Optimalisatie met meerdere responsen is een uitdaging bij 3D-printtoepassingen: de dichtheid moet worden gemaximaliseerd, de ruwheid geminimaliseerd en de maatafwijking moet zo dicht mogelijk bij nul liggen. Deze drie doelen zijn vaak onderling strijdig. In de techniek wordt een desirability-functie gebruikt om elke respons naar een bereik van 0 tot 1 om te zetten. Het geometrische gemiddelde levert vervolgens de totale desirability D op, en de parametercombinatie wordt bepaald waar D maximaal is. De weging wordt toegekend volgens de prioriteiten in de klantacceptatie; voor luchtvaartonderdelen kan bijvoorbeeld dichtheid een gewicht van 0,5 krijgen, ruwheid 0,3 en maatafwijking 0,2.

4. Kalibratie van het procesvenster: maak van het optimale punt een vrij te geven bereik

Het vinden van het optimale punt is slechts de eerste stap; voor serieproductie is een bereik nodig. Een bruikbare methode is een overlappende contourgrafiek (overlay contour plot): teken voor elke respons een contour op basis van de acceptatiedrempel — dichtheid ≥ 99,5%, Ra ≤ 12 μm en maatafwijking ≤ ±0,1 mm — en gebruik het overlappende gebied van alle contouren als procesvenster.

Nadat het venster is vastgesteld, moet een veiligheidsmarge worden aangehouden. Het advies is om de grenzen van het venster 10% naar binnen te verschuiven. Als het model bijvoorbeeld een haalbaar vermogensbereik van 235–265 W geeft, wordt op de proceskaart 238–262 W vermeld. De reden hiervoor is dat variatie tussen poederbatches, fluctuaties in de zuiverheid van het beschermgas en vermogensverval van het optische systeem het werkelijke procespunt geleidelijk kunnen verschuiven. Door marge aan te houden kan de procescapabiliteitsindex Cpk van 1,0 naar meer dan 1,33 worden verhoogd.

Voor vensterverificatie worden drie batches gebruikt, elk met 5 stuks. Per batch wordt een groep geprint in het centrum van het venster, net binnen de bovengrens en net binnen de ondergrens. Wanneer drie opeenvolgende batches volledig voldoen en de verschuiving van het gemiddelde tussen batches kleiner is dan een derde van de standaardafwijking van de respons, kan het venster als stabiel worden beoordeeld. De verificatiegegevens worden tevens in het procesregister opgenomen als basislijn voor latere traceerbaarheid.

5. Veelvoorkomende fouten en controlelijst

Fout één: het aantal proefstukken beschouwen als de volledige kostenpost. In werkelijkheid zijn de totale kosten van 9 orthogonale proeven plus 20 CCD-proeven veel lager dan die van 30 blinde rondes met één factor tegelijk, terwijl de uitkomst een herbruikbaar model is.

Fout twee: blokeffecten negeren. Als 20 CCD-proeven over drie dagen worden uitgevoerd, kunnen veranderingen in machinestatus, poederbatch en temperatuur en vochtigheid van de omgeving in de residuen terechtkomen. De juiste aanpak is om de datum als blokfactor in het ontwerp op te nemen, of herhalingen van centrumpunten binnen elk blok te plannen om drift te bewaken.

Fout drie: alleen naar dichtheid kijken. Een procespunt met een dichtheid van 99,5% voldoet niet noodzakelijk aan de mechanische eisen; gebrek-aan-versmelting-defecten kunnen bij dichtheidsmetingen worden onderschat. Het wordt aanbevolen om ten minste treksterkte en rek gelijktijdig te testen en, waar nodig, metallografisch onderzoek en steekproefcontroles met industriële CT toe te voegen.

Uitvoeringscontrolelijst: dekken de factorniveaus het bereik waarin de apparatuur stabiel kan werken; zijn herhalingen van centrumpunten gepland om de zuivere fout te schatten; is een lack-of-fit-toets uitgevoerd; sluiten de gewichten van meerdere responsen aan op de acceptatiecriteria van de klant; bevat het venster 10% marge; dekken de verificatiebatches de venstergrenzen af?

Samenvatting

Experimenteel ontwerp verandert de ontwikkeling van 3D-printprocessen van empirisch trial-and-error in een kwantificeerbare en reproduceerbare technische activiteit: orthogonale proeven voltooien de factorselectie in 9 tot 12 tests, de responsvlakmethode bouwt met ongeveer 20 tests een kwadratisch model met interactie-effecten op, de desirability-functie lost conflicten tussen meerdere doelen op en overlappende contouren leveren een vrij te geven procesvenster met veiligheidsmarge. Voor lifecyclemanagementplatforms zoals lantu3D, die het traject van ontwerptekening tot fysieke oplevering ondersteunen, beperkt de waarde van DoE zich niet tot parameterafstelling: het legt proceskennis vast in doorzoekbare, herbruikbare en verifieerbare parametermodellen, zodat elke procesontwikkeling het vertrekpunt voor de volgende wordt in plaats van telkens opnieuw vanaf nul te beginnen.