The present disclosure relates to the technical field of photo-curing three dimensional printing simulation, in particular to a photo-curing three dimensional printing preview method.
Three Dimensional Printing (3DP) is a novel rapid prototyping technology developed by relying on electronic information, mechanical processing, sensors and new material development, has the incomparable advantages of design and manufacturing integration, short production cycle, low cost and the like compared with a forced forming technology like traditional planing and milling, is known as “manufacturing technology with industrial revolution significance”, and has been widely used in tissue engineering, industrial manufacturing, energy materials, aerospace and other fields. Compared with other 3DP technologies, photo-curing three dimensional printing has the advantages of high forming accuracy, low cost, and diversified design of functional photosensitive resins, thereby making it one of the most popular 3DP technologies at present.
Photo-curing three dimensional printing uses UV light to irradiate thin layers of liquid photosensitive resins for stacking molding, and light-induced thermosetting liquid resins undergo cross-linking reaction between layers under UV light irradiation. Therefore, printing settings that affect the reaction kinetics, such as a slicing layer thickness, the printing time of each layer and the UV light source light intensity will affect the printing effect and quality of a finished product, the debugging process wastes more time and labor than 2D printing, and materials are prone to be wasted. Compared with 2D printing editing software with a built-in print preview function, such as OFFICE series, 3DP, especially photo-curing three dimensional printing control software involving cross-linking reaction, needs to introduce a printing preview function in addition to a slicing function to help users guide printing settings and even automatically design and compare printing setting solutions based on imported graphics.
The technical problem to be solved by the present disclosure is to provide a photo-curing three dimensional printing preview method, which simulates a resin layer conversion rate of a printed part under printing settings by combining a 3DP layer-by-layer stacking mechanism and a photo-curing cross-linking kinetic principle so as to characterize a preview effect.
The technical solution adopted by the present disclosure to solve the technical problem is as follows: a photo-curing three dimensional printing preview method is provided, and includes the following steps:
(1) obtaining a relationship curve between a conversion rate and time of a photosensitive resin system under different curing conditions by measuring a relationship between an absorbance and time of the photosensitive resin system under different curing conditions;
(2) fitting an experimental curve by using an iterative solution method based on a photo-curing kinetic equation, obtaining a photo-curing kinetic constant of the photosensitive resin system, and establishing a photo-curing reaction kinetic differential equation set;
(3) solving the photo-curing reaction kinetic differential equation set by using a fixed-step Euler method, and obtaining a spatial distribution and time relationship curve of a content of each component of the photosensitive resin system in a photo-curing reaction process;
(4) measuring a relationship curve between an incident light transmissivity and a thickness of a cured product before and after a photo-curing reaction of the photosensitive resin system through an ultraviolet spectrum; and
(5) performing simulated printing according to a structure of a three dimensional printing product, a light source intensity, a layer thickness and single-layer printing time, obtaining spatial distribution of the content of each component in each layer structure in the three dimensional printing product after printing is completed by combining with the spatial distribution and time relationship curve of the content of each component of the photosensitive resin system in the photo-curing reaction process, and completing photo-curing three dimensional printing preview.
The step (1) is specifically as follows: measuring a change situation of an absorbance of functional groups with time in 3-7 groups of photosensitive resins under the different curing conditions by using an in-situ Fourier transform infrared spectrometer, and converting the change situation into a change curve of the conversion rate with time through an internal standard method, wherein, the different curing conditions include different UV light source intensities and different photoinitiator contents.
The photo-curing kinetic constant obtained in the step (2) includes a chain initiation constant, a chain growth constant and a termination constant.
The photo-curing reaction kinetic differential equation set in the step (2) is obtained by simultaneously obtaining reaction rate equations including initiator decomposition, chain initiation, chain growth, and chain termination.
Due to the adoption of the above-mentioned technical solution, the present disclosure has the following advantages and positive effects compared with the prior art: the present disclosure can obtain a kinetic parameter input value of an unknown or arbitrary formula photosensitive resin system by a method for combining an infrared spectrum and the UV spectrum to characterize a photosensitive resin stacking curing property, thereby making a simulation result more realistic. Combining a classical free radical polymerization kinetic model introduced with oxygen inhibition and a light attenuation principle, a method for realizing the simulation of a photo-curing three dimensional printing preview effect provided by the present disclosure makes up for the blank of three dimensional printing control software on the aspect of the preview function.
The present disclosure will be further described below with reference to specific embodiments. It should be understood that these embodiments are merely illustrative of the present disclosure and are not intended to limit the scope of the present disclosure. In addition, it should be understood that after reading the content taught by the present disclosure, those skilled in the art can make various changes or modifications to the present disclosure, and these equivalent forms also fall within the scope defined by the appended claims of the present application.
An implementation of the present disclosure relates to a photo-curing three dimensional printing preview method. The method is based on a free radical cross-linking reaction model, since a three dimensional printing process is usually carried out in an air atmosphere, and oxygen captures active free radicals generated in the reaction process, an oxygen inhibition model is introduced into classical chain initiation, chain growth and chain termination mechanisms; and a fixed-step Euler method is used to solve a kinetic differential equation set, a time-varying value of each variable is derived, and according, a simulated three dimensional printing preview effect is achieved. The method specifically includes the following steps.
(1) A relationship curve between a conversion rate and time of a photosensitive resin system under different curing conditions is obtained by measuring a relationship between an absorbance and time of the photosensitive resin system under different curing conditions, specifically: a change situation of an absorbance of functional groups with time in 3-7 groups of photosensitive resins under the different curing conditions (different UV light source intensities, different photoinitiator contents and the like) is measured by using an in-situ Fourier transform infrared spectrometer, and converted into a change curve of the conversion rate with time through an internal standard method.
(2) An experimental curve is fitted by using an iterative solution method based on a photo-curing kinetic equation according to the relationship curve between the conversion rate and time of the photosensitive resin system under the different curing conditions, and a photo-curing kinetic constant of the photosensitive resin system is obtained by calculation. For the same curing system, the kinetic constant (a chain initiation constant, a chain growth constant, and a termination constant) that conform to each experimental curve may be calculated according to the change curve of the conversion rate with time under several input sets of different settings. The parameters of the free radical photopolymerization of photosensitive resins involved in the calculation may be set according to experimental measurement results and the actual needs of photo-curing printing, specifically including calculation parameters such as a calculation step size, a grid size, and an area size, curing parameters such as a light source intensity, a layer thickness and an exposure time of single layer, and physical properties such as resin monomer and initiator densities, a molar absorption coefficient, a temperature, and an initial molar concentration.
(3) The photo-curing reaction kinetic differential equation set is solved by using a fixed-step Euler method, and a spatial distribution and time relationship curve of a content of each component of the photosensitive resin system in a photo-curing reaction process is obtained.
(4) A relationship curve between an incident light transmissivity and a thickness of a cured product before and after a photo-curing reaction of the photosensitive resin system is measured through a UV spectrum.
(5) Simulated printing is performed according to a structure of a three dimensional printing product, the light source intensity, the layer thickness, and single-layer printing time; spatial distribution of a content of each component in each layer structure in the three dimensional printing product after printing is completed is obtained by combining with the spatial distribution and time relationship curve of the content of each component of the photosensitive resin system in the photo-curing reaction process, and photo-curing three dimensional printing preview is completed.
In the process of simulating layer-by-layer curing and superposition, the present disclosure collects the change situation of the absorbance of the functional groups with time through an infrared spectrum, converts it into the change curve of the conversion rate with time through the internal standard method, and then obtains the kinetic constant, thereby obtaining the spatial distribution of the content of each component in the reaction process, the value may be represented by data or by a gradient color, and the spatial distribution is reflected by displaying a spatial content diagram of any component according to a selected section.
The present disclosure will be further described below through a specific embodiment.
Step 1, a photosensitive resin curing infrared spectrum is collected.
Epoxy acrylic resins containing 1%, 2% and 3% of 2,4,6-trimethylbenzoyl-diphenylphosphine oxide are prepared separately and divided into three groups for experiments. For each group of resins, an in-situ photo-curing Fourier transform infrared spectrometer is adopted to achieve in-situ infrared spectrum detection by UV light curing through an optical path system: UV light is projected onto the surface of a sample through an optical fiber, and an energy density of a light source is measured by a UV light power meter to realize closed-loop control. The liquid resin is placed between two potassium bromide tablets with a thickness of 1 mm and is surrounded by a ring member to fix a thickness and prevent leakage of a liquid monomer. An infrared light source enters a detector after passing through the sample to be tested, and a computer measures the infrared spectrum of the sample at regular intervals. Installation of an in-situ photo-curing Fourier transform infrared spectrometer optical path system and a liquid resin sample is shown in
Step 2, photosensitive resin kinetic parameters are calculated.
According to infrared spectrum data of the monomer, a characteristic peak is selected, a monomer conversion curve is calculated according to the decrease of an absorption peak relative to an initial state, and a monomer conversion curve with different curing parameters is obtained according to experimental results of different initiator concentrations. Calculation parameters such as a calculation step size, a grid size, and an area size, curing parameters such as a light source intensity, a layer thickness and an exposure time of single layer, physical properties such as resin monomer and initiator densities, a molar absorption coefficient, a temperature, and an initial molar concentration, and estimated values of reaction kinetic constants such as a chain initiation constant, a chain growth constant, and a termination constant are set. A calculation program automatically calculates experimental values of the reaction kinetic constants according to input parameters. A monomer infrared spectrum absorption characteristic peak-monomer concentration curve is shown in
Step 3, a photo-curing three dimensional printing curing process is calculated.
A three dimensional model file is input, the calculation parameters and the curing parameters are set, and the program simulates the curing process according to the settings, finally outputs time and space distribution data of various parameters including the monomer conversion rate, the light intensity, and the initiator concentration, and generates a three dimensional printing molding effect preview according to the calculation data. The comparison of the curing thickness of a three dimensional printing layer under different additive concentrations refers to
Number | Date | Country | Kind |
---|---|---|---|
201911012217.3 | Oct 2019 | CN | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CN2019/118315 | 11/14/2019 | WO |