CN113419120A - Method and system for measuring thermal resistance of dielectric film and metal interface - Google Patents
Method and system for measuring thermal resistance of dielectric film and metal interface Download PDFInfo
- Publication number
- CN113419120A CN113419120A CN202110500544.4A CN202110500544A CN113419120A CN 113419120 A CN113419120 A CN 113419120A CN 202110500544 A CN202110500544 A CN 202110500544A CN 113419120 A CN113419120 A CN 113419120A
- Authority
- CN
- China
- Prior art keywords
- dielectric film
- thermal
- response current
- thermal response
- interface
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/003—Environmental or reliability tests
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
The invention relates to a method and a system for measuring thermal resistance of a dielectric film and a metal interface, wherein a thermal pulse experiment is carried out after the dielectric film and a metal foil are tightly attached to obtain an actually measured thermal response current signal curve; and establishing a heat transfer simulation model of the dielectric film and the metal foil, and when the fitting degree of a theoretical thermal response signal current curve and an actually measured thermal response current signal curve of the simulation model meets a preset convergence condition, taking the interface thermal resistance value of the simulation model as the interface thermal resistance of the dielectric film and the metal foil. Compared with the prior art, the method has the advantages that the implementation process is simple, the heat transfer simulation model is established after the actually measured thermal response current is obtained, the interface thermal resistance value of the simulation model is continuously adjusted until the fitting degree of the actually measured thermal response current and the theoretical thermal response current of the simulation model meets the convergence condition, the measurement speed is high, the accuracy is high, the operation is convenient, the dielectric film is not damaged, and the interface thermal resistance between the independent dielectric film material with the thickness of several micrometers and metal can be measured at high precision.
Description
Technical Field
The invention relates to the field of measurement of interface thermal resistance, in particular to a method and a system for measuring the interface thermal resistance of a dielectric film and metal.
Background
In an electronic device, the safety and stability of the device operation are generally ensured by mounting a heat sink, a heat dissipation plate, and the like to transfer internal heat. The dielectric film is widely used as a functional material for microelectronics, electrical appliance insulation and energy storage capacitors, so that it is very important to accurately master thermal resistance information between the dielectric film and the metal heat sink to effectively manage the heat of an electronic device. The classical theoretical calculation model of the interface thermal resistance is as follows:
R12 k/σ=c(H/p)n
wherein R is12The contact thermal resistance of two contact materials is shown, k is the harmonic mean value of the thermal conductivities of the material 1 and the material 2, sigma is the surface roughness, H is the microhardness of the material, p is the contact pressure, and c and n are empirical constants. According to the formula, the theoretical calculation of the thermal resistance model is complex, a plurality of variables which are difficult to determine exist, and the interface thermal resistance is difficult to determine through the theoretical calculation, so that various experimental measurements for representing the interface thermal resistance are provided at present.
The experimental measurement is roughly classified into a steady-state method and a transient method. The steady state method comprises the following steps: maintaining a certain temperature difference on the two contact samples, measuring the axial temperature values of the two samples, and extrapolating the temperature values to the contact interface by the Fourier law so as to obtain the temperature difference on the interface; steady State method D5470-06 is often used for bulk materials according to ASTM Standard with a measurable Range of 10-7~10-6m2K/W, the measurement accuracy is good, but the time consumption is long, and the method is not suitable for medium films with the thickness of only a few micrometers. The transient measurement method mainly comprises a laser flash method, a heat reflection method, a TPS method and the like, most of the transient methods directly detect the temperature distribution of a sample to reversely deduce the interface thermal resistance, but for a dielectric film with the thickness of only a few micrometers, the accurate detection of the temperature in the thickness direction has certain difficulty, and the problem of insufficient precision of a detected temperature signal is difficult to avoid.
Therefore, how to measure the interface thermal resistance between the dielectric film and the metal heat sink and meet the requirement of measurement accuracy becomes an urgent problem to be solved.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provides a method and a system for measuring the interface thermal resistance between a medium film and a metal, the method and the system are simple in implementation process, the metal foil is selected as a heat sink, a heat transfer simulation model is established after the actually measured thermal response current of the medium film and the metal foil in a thermal pulse experiment is measured, the interface thermal resistance value in the heat transfer simulation model is continuously adjusted until the fitting degree of the theoretical thermal response current and the actually measured thermal response current of the heat transfer simulation model meets the convergence condition, the interface thermal resistance value of the heat transfer simulation model is used as the interface thermal resistance between the medium film and the metal foil, the measurement speed is high, the accuracy is high, the operation is convenient, the medium film is not damaged, and the interface thermal resistance between an independent medium film material with the thickness of several micrometers and the metal can be measured at high accuracy.
The purpose of the invention can be realized by the following technical scheme:
a method for measuring thermal resistance of a dielectric film and a metal interface comprises tightly attaching the dielectric film and a metal foil, performing a heat pulse experiment, and measuring to obtain an actually measured thermal response current signal curve; establishing a heat transfer simulation model of the dielectric film and the metal foil; and adjusting the interface thermal resistance value of the heat transfer simulation model, calculating to obtain a theoretical thermal response signal current curve, and taking the interface thermal resistance value of the heat transfer simulation model as the actual interface thermal resistance value of the dielectric film and the metal foil when the fitting degree of the actually measured thermal response current signal curve and the theoretical thermal response signal current curve meets the preset convergence condition.
Further, obtaining the measured thermal response current signal curve comprises the following steps:
a1: after the medium film is subjected to metallization treatment, the medium film is tightly attached to the metal foil, and direct-current voltages are externally connected to two sides of the medium film, so that an electric field which is uniformly distributed exists in the medium film;
a2: installing the dielectric film and the metal foil in a shielding box, applying pulse laser to heat the dielectric film, and collecting the actually measured thermal response current of the dielectric film;
a3: and converting the actually measured thermal response current into a frequency domain signal through fast Fourier transform to obtain an actually measured thermal response current signal curve.
Furthermore, the thermal response current signal is affected by the external dc voltage circuit and the pre-current amplifier, in step a3, distortion compensation is performed on the measured thermal response current to obtain a calibrated measured thermal response current, and then the calibrated measured thermal response current is converted into a frequency domain signal to obtain a measured thermal response current signal curve, where the formula is as follows:
Im~(f)=I0~(f)*H~(f)
wherein, Im~(f) Representing the collected measured thermally responsive current, I0~(f) Represents the thermal response current in the ideal case, H~(f) Is a transfer function determined from the effect of the external circuit and amplifier bandwidth limitations on the thermally responsive current.
Further, establishing a heat transfer simulation model of the dielectric film and the metal foil specifically comprises: the method comprises the steps of obtaining initial temperature distribution, thickness, thermal diffusion coefficient, thermal conductivity coefficient, relative dielectric coefficient, temperature coefficient and thermal expansion coefficient of the relative dielectric coefficient of a dielectric film, obtaining initial temperature distribution, thickness, thermal diffusion coefficient and thermal conductivity coefficient of a metal foil, obtaining parameters of laser pulse and electric field in a thermal pulse experiment process, establishing a double-layer structure of the dielectric film and the metal foil in simulation software to serve as a heat transfer simulation model, initializing interface thermal resistance values of the heat transfer simulation model, and simulating disturbance of an externally applied direct current electric field and laser pulse in the double-layer structure in the simulation software.
Further, the step of obtaining the theoretical thermal response current signal curve comprises the following steps:
b1: establishing a one-dimensional heat conduction equation of a double-layer structure of a heat transfer simulation model:
the boundary conditions are as follows:
the initial conditions were:
θ1=T1(x)(0<x<x1)(t=0)
θ2=T2(x)(x1<x<x2)(t=0)
wherein, theta1Is the temperature distribution of the dielectric film, a1Is the thermal diffusion coefficient, θ, of the dielectric film2Is the temperature distribution of the metal foil, a2Is the thermal diffusion coefficient of the metal foil, x is the position from the dielectric film to the thickness direction of the metal foil, x1Equal to the thickness, x, of the dielectric film2Is equal to the sum of the thicknesses of the dielectric film and the metal foil, t is time, k1Is the thermal conductivity, k, of the dielectric film2Is the thermal conductivity of the metal foil, h12Is the heat exchange coefficient, T, of the dielectric film and the metal foil1(x) Is the initial temperature distribution, T, of the dielectric film2(x) Is the initial temperature distribution of the metal foil, eta is the absorption coefficient of the pulse laser, q is the energy of the pulse laser, and delta (t) is the Gaussian distribution function of the pulse laser;
b2: solving the one-dimensional heat conduction equation to obtain the theoretical calculation of the temperature distribution of the dielectric filmTheoretical calculation of temperature distribution based on dielectric filmObtaining the temperature change inside the dielectric film according to the following formulaCalculating to obtain theoretical thermal response current:
g(x)=ε0εr(αε-αx)E(x)
wherein I (t) represents a thermal response current, A represents an area of the dielectric film irradiated by the laser, d represents a thickness of the dielectric film, Δ θ represents a temperature change inside the dielectric film, g (x) represents a distribution function, E (x) represents an electric field distribution inside the dielectric film, ε0Denotes the vacuum dielectric constant,. epsilonrDenotes the relative dielectric coefficient, alpha, of the dielectric filmεTemperature coefficient, alpha, representing the relative dielectric coefficientxRepresents the thermal expansion coefficient of the dielectric film;
b3: and performing fast Fourier transform on the theoretical thermal response current and converting the theoretical thermal response current into a frequency domain signal to obtain a theoretical thermal response current signal curve.
Further, adjusting the interface thermal resistance value of the heat transfer simulation model and determining the actual interface thermal resistance value of the dielectric thin film and the metal foil specifically comprises:
c1: obtaining a theoretical thermal response current signal curve based on the current interface thermal resistance value of the heat transfer simulation model;
c2: fitting the theoretical thermal response current signal curve with the actually measured thermal response current signal curve, if the fitting degree does not meet the preset convergence condition, adjusting the interface thermal resistance value in the heat transfer simulation model, and repeating the step C1, otherwise, executing the step C3;
c3: and taking the interface thermal resistance value of the heat transfer simulation model as the actual interface thermal resistance value of the dielectric film and the metal foil.
Further, a Nelder-Mead simplex algorithm is used to fit the theoretical thermal response current signal curve to the measured thermal response current signal curve in step C2.
Further, the preset convergence condition is a default convergence condition in the Nelder-Mead simplex algorithm.
A system for measuring thermal interface resistance of a dielectric film and a metal comprises:
the actual measurement module is used for tightly attaching the dielectric film to the metal foil, performing a heat pulse experiment and measuring to obtain an actual measurement heat response current signal curve;
the modeling module is used for establishing a double-layer structure of the dielectric film and the metal foil according to the parameters of the dielectric film and the metal foil and the parameters of the heat pulse experiment;
the calculation module is used for calculating a theoretical thermal response current signal curve of the heat transfer simulation model;
and the fitting module is used for fitting the actually measured thermal response current signal curve and the theoretical thermal response current signal curve and judging whether the fitting degree meets the preset convergence condition or not.
Furthermore, in the actual measurement module, the actual measurement thermal response current of the thermal pulse experiment is obtained, distortion compensation is performed on the actual measurement thermal response current to obtain a calibrated actual measurement thermal response current, and then the calibrated actual measurement thermal response current is converted into a frequency domain signal to obtain an actual measurement thermal response current signal curve.
Compared with the prior art, the invention has the following beneficial effects:
(1) the method has the advantages that the realization process is simple, after the actually measured thermal response current of the dielectric film and the metal foil in the thermal pulse experiment is measured, the heat transfer simulation model is established, the interface thermal resistance value in the heat transfer simulation model is continuously adjusted until the fitting degree of the theoretical thermal response current and the actually measured thermal response current of the heat transfer simulation model meets the convergence condition, the interface thermal resistance value of the heat transfer simulation model is used as the interface thermal resistance of the dielectric film and the metal foil, the measurement speed is high, the accuracy is high, the operation is convenient, the dielectric film is not damaged, and the interface thermal resistance between the independent dielectric film material with the thickness of several micrometers and the metal can be measured with high accuracy.
(2) The influence of the pre-current amplifier on the actually measured thermal response current when the external direct current voltage and the thermal response current are collected is considered, distortion compensation is carried out on the actually measured thermal response current measured in the thermal pulse experiment, and more accurate measurement is achieved.
(3) Establishing a heat transfer simulation model, simulating a heat pulse experiment in simulation software, calculating the temperature distribution of the medium film through the established heat conduction equation, calculating by using a first class of Fredholm integral equation according to the temperature distribution to obtain theoretical thermal response current, and calculating the theoretical thermal response current of the heat transfer simulation model under different interface thermal resistances.
Drawings
FIG. 1 is a flow chart of a method for measuring thermal resistance at the interface between a dielectric film and a metal;
fig. 2 is a fitting result of a measured thermal response current signal curve and a theoretical thermal response signal current curve.
Detailed Description
The invention is described in detail below with reference to the figures and specific embodiments. The present embodiment is implemented on the premise of the technical solution of the present invention, and a detailed implementation manner and a specific operation process are given, but the scope of the present invention is not limited to the following embodiments.
Example 1:
a method for measuring thermal resistance of a dielectric film and a metal interface belongs to a transient measurement method, is convenient to operate, has high measurement speed, does not damage a sample, and can effectively measure the thermal resistance of the interface between the dielectric film and a metal foil, as shown in figure 1, the method mainly comprises the following steps:
tightly attaching the dielectric film and the metal foil, performing a heat pulse experiment and measuring to obtain an actually measured thermal response current signal curve;
establishing a heat transfer simulation model of the dielectric film and the metal foil;
and adjusting the interface thermal resistance value of the heat transfer simulation model, calculating to obtain a theoretical thermal response signal current curve, and taking the interface thermal resistance value of the heat transfer simulation model as the actual interface thermal resistance value of the dielectric film and the metal foil when the fitting degree of the actually measured thermal response current signal curve and the theoretical thermal response signal current curve meets the preset convergence condition.
A system for measuring thermal interface resistance of a dielectric film and a metal comprises:
the actual measurement module is used for tightly attaching the dielectric film to the metal foil, performing a heat pulse experiment to obtain an actual measurement thermal response current of the heat pulse experiment, performing distortion compensation on the actual measurement thermal response current to obtain a calibrated actual measurement thermal response current, and then converting the calibrated actual measurement thermal response current into a frequency domain signal to obtain an actual measurement thermal response current signal curve;
the modeling module is used for establishing a double-layer structure of the dielectric film and the metal foil according to the parameters of the dielectric film and the metal foil and the parameters of the heat pulse experiment;
the calculation module is used for calculating a theoretical thermal response current signal curve of the heat transfer simulation model;
and the fitting module is used for fitting the actually measured thermal response current signal curve and the theoretical thermal response current signal curve and judging whether the fitting degree meets the preset convergence condition or not.
In this embodiment, the dielectric film is a BOPP dielectric film with a thickness of 5.8 μm, the metal foil is a copper foil with a thickness of 2mm, related parameters (thermal conductivity, thermal diffusivity, etc.) of the dielectric film can be measured in advance, for example, the multilayer dielectric films are stacked to measure the related parameters, the metal foil is a copper foil or a copper sheet, the related parameters of which can be directly obtained by looking up the parameters of the metal copper, and obtaining the actually measured thermal response current signal curve includes the following steps:
a1: after the medium film is subjected to metallization treatment, the medium film is tightly attached to the metal foil, and direct-current voltages are externally connected to two sides of the medium film, so that an electric field which is uniformly distributed exists in the medium film;
a2: installing the dielectric film and the metal foil in a shielding box, applying short pulse laser (25ns,1Hz) to heat the dielectric film, and collecting the actually measured thermal response current of the dielectric film through a signal collection device of a thermal pulse experiment;
the thermal pulse method generates thermal disturbance by applying thermal pulse to the dielectric film, the temperature distribution of the dielectric film is changed, the thermal pulse drives charges in the internal propagation process of the dielectric film to trigger a measurable thermal response current, and the characteristics of the thermal response current can effectively represent the information of the temperature distribution in the dielectric film.
The heat flow flows into the metal foil after passing through the medium film, and the air gap between the medium film and the metal foil causes the air with low heat conductivity to be reserved, so that the heat transfer of the heat flow between the two material interfaces is blocked, and the heat flow is formedAn interfacial thermal resistance R is provided. The air gap not only hinders heat conduction, but also affects electrical contact between an external power source and the dielectric film, resulting in poor signal-to-noise ratio in response to current. Therefore, in experimental measurement, thermal interface materials such as a thermal coupling agent and the like need to be added between the metal foil and the dielectric film, thermal resistance is reduced, and signal to noise ratio is improved, so as to obtain effective measurement dataintThe model of (a) is expressed as:
Rint=Rcontact1+Rcond+Rcontact2
Rcond=dg/kTims
thermal interface resistance RintMainly comprises an upper interface contact thermal resistance Rcontact1Lower interface thermal contact resistance Rcontact2And self thermal resistance R of thermal interface materialcond,dgIs the thickness of the thermal interface material, kTimsIs the thermal conductivity of the thermal interface material.
And (3) while acquiring a measured thermal response current signal curve, recording parameters of laser pulse and an electric field in a thermal pulse experiment, if the laser is pulse laser, the frequency is 1Hz, the pulse width is 25ns, and the electric field intensity E in the dielectric film is 34.5k V/mm in the thermal pulse experiment, so that preparation is made for subsequently establishing a heat transfer simulation model.
A3: and converting the actually measured thermal response current into a frequency domain signal through fast Fourier transform to obtain an actually measured thermal response current signal curve.
The signal of the thermal response current can be influenced by an external direct-current voltage circuit and a front-end current amplifier, in order to realize more accurate determination, step a3 further comprises distortion compensation of the measured thermal response current to obtain the calibrated measured thermal response current, and then the calibrated measured thermal response current is converted into a frequency domain signal to obtain a measured thermal response current signal curve, wherein the formula is as follows:
Im~(f)=I0~(f)*H~(f)
wherein, Im~(f) Representing the collected measured thermally responsive current, I0~(f) Represents the thermal response current in the ideal case, H~(f) Is a transfer function determined according to the influence of the external circuit and the amplifier bandwidth limitation on the thermal response current, the transfer function H~(f) Distortion compensation transfer functions commonly used by those skilled in the art may be used.
The establishment of the heat transfer simulation model of the dielectric film and the metal foil comprises the following specific steps: the method comprises the steps of obtaining parameters such as initial temperature distribution, thickness, thermal diffusion coefficient, thermal conductivity, relative dielectric coefficient, temperature coefficient and thermal expansion coefficient of a dielectric film, obtaining parameters such as initial temperature distribution, thickness, thermal diffusion coefficient and thermal conductivity coefficient of a metal foil, obtaining parameters of laser pulse and electric field in a thermal pulse experiment process, establishing a double-layer structure of the dielectric film and the metal foil in simulation software to serve as a heat transfer simulation model, initializing interface thermal resistance of the heat transfer simulation model, and simulating disturbance of an externally applied direct current electric field and the laser pulse in the double-layer structure in the simulation software.
In the heat transfer simulation model, the interface thermal resistance value is first set empirically, and in this embodiment, the interface thermal resistance value is initially set to 1 × 10 empirically-6K·m2And W, continuously adjusting the interface thermal resistance value according to the fitting degree of the theoretical thermal response current signal curve and the actually measured thermal response current signal curve.
In this example, the dielectric film had a thickness of 5.8 μm and a density of 910kg/m3The specific heat capacity is 1900J/(kg.K), the thermal conductivity, i.e., the thermal conductivity is 0.165W/(m.K), and the thermal diffusivity is 0.95X 10-7m2(s) the relative dielectric constant (60Hz) of the dielectric thin film was 2.2, and the temperature coefficient of the relative dielectric constant of the dielectric thin film was 3.31X 10-41/K, coefficient of thermal expansion of the dielectric film is 1.35 multiplied by 10-4 1/K;
The metal foil has a thickness of 2mm and a density of 8960kg/m3The specific heat capacity is 384.6J/(kg.K), the thermal conductivity, i.e., the thermal conductivity is 401W/(m.K), and the thermal diffusivity is 1.16X 10-4m2/s。
The method for acquiring the theoretical thermal response current signal curve comprises the following steps:
b1: establishing a one-dimensional heat conduction equation of a double-layer structure of a heat transfer simulation model:
the boundary conditions are as follows:
the initial conditions were:
θ1=T1(x)(0<x<x1)(t=0)
θ2=T2(x)(x1<x<x2)(t=0)
wherein, theta1Is the temperature distribution of the dielectric film, a1Is the thermal diffusion coefficient, θ, of the dielectric film2Is the temperature distribution of the metal foil, a2Is the thermal diffusion coefficient of the metal foil, x is the position from the dielectric film to the thickness direction of the metal foil, x1Equal to the thickness, x, of the dielectric film2Equal to the sum of the thicknesses of the dielectric film and the metal foil, wherein x is 0, namely the upper surface of the dielectric film, and x is x1On the lower surface of the dielectric film, or on the upper surface of the metal foilFace, x ═ x2Where is the lower surface of the metal foil, t is time, k1Is the thermal conductivity, k, of the dielectric film2Is the thermal conductivity of the metal foil, h12Is the heat exchange coefficient between the dielectric film and the metal foil, which is determined based on the interface thermal resistance value, T1(x) Is the initial temperature distribution, T, of the dielectric film2(x) Is the initial temperature distribution of the metal foil, eta is the absorption coefficient of the pulse laser, q is the energy of the pulse laser, and delta (t) is the Gaussian distribution function of the pulse laser;
the temperature distribution refers to the temperature distribution of different thicknesses and different time, the initial temperature distribution can be measured in advance, the temperature of each thickness of the dielectric film and the metal foil is room temperature, and after the laser pulse heating is applied, the temperature of the dielectric film at different thicknesses can change until the thermal equilibrium is reached and then the temperature is kept unchanged.
B2: the principal methods for solving the one-dimensional heat conduction equation include separation variable method, Green function method, Laplace transform, finite element method, etc., and the theoretical calculation of the temperature distribution of the obtained dielectric filmTheoretical calculation of temperature distribution based on dielectric filmObtaining the temperature change in the medium film, wherein the heat diffusion depth of the laser in the medium of the medium film is far less than the radius of the laser light target, the heat conduction problem can be approximated to a one-dimensional heat conduction problem, the medium film is a non-polar medium, and the calculation of the thermal response current relates to a first Fredholm integral equation as follows:
g(x)=ε0εr(αε-αx)E(x)
wherein I (t) represents a thermal response current, A represents an area of the dielectric film irradiated by the laser, d represents a thickness of the dielectric film, and Delta theta represents the dielectric filmInternal temperature change, g (x) represents a distribution function, E (x) represents an electric field distribution within the dielectric film,. epsilon0Representing a vacuum dielectric constant of 8.8537X 10-12F/m,εrDenotes the relative dielectric coefficient, alpha, of the dielectric filmεTemperature coefficient, alpha, representing the relative dielectric coefficientxRepresents the thermal expansion coefficient of the dielectric film;
b3: and performing fast Fourier transform on the theoretical thermal response current and converting the theoretical thermal response current into a frequency domain signal to obtain a theoretical thermal response current signal curve.
Adjusting the interface thermal resistance value of the heat transfer simulation model and determining the actual interface thermal resistance value of the dielectric film and the metal foil specifically as follows:
c1: obtaining a theoretical thermal response current signal curve based on the current interface thermal resistance value of the heat transfer simulation model;
c2: fitting the theoretical thermal response current signal curve and the actually measured thermal response current signal curve by using a Nelder-Mead simplex algorithm, wherein the preset convergence condition is a default convergence condition in the Nelder-Mead simplex algorithm, if the fitting degree does not meet the preset convergence condition, the interface thermal resistance value in the heat transfer simulation model is adjusted, and the step C1 is repeated, otherwise, the step C3 is executed;
c3: and taking the interface thermal resistance value of the heat transfer simulation model as the actual interface thermal resistance value of the dielectric film and the metal foil.
The Nelder-Mead simplex algorithm can process the minimum solving problem of a multivariable function, fitting a theoretical thermal response current signal curve and an actually measured thermal response current signal curve and adjusting interface thermal resistance are carried out, the theoretical thermal response current signal curve closest to the actually measured thermal response current signal curve and the interface thermal resistance value of the theoretical thermal response current signal curve can be quickly and accurately found, MATLAB is provided with an fmisearch () function, the fmisearch () function is directly used for solving to obtain the interface thermal resistance value, and the default convergence condition in the algorithm is used as the convergence condition.
In the whole heat transfer process, the thermal resistance has the function of preventing heat flow from being rapidly transferred to the metal foil and slowing down the heat dissipation of the medium film, so that the interface thermal resistance RintThe main contribution of (1) is mainly reflected in the process of cooling a dielectric film sample, the action on time domain response current is mainly reflected in the characteristic that the thermal response current is attenuated to a negative value region from a peak value, the signal characteristic of a relatively low frequency band is mainly concerned in a frequency domain, the fitting result is shown in figure 2, the fitting of a real part and an imaginary part of an actually measured thermal response current signal curve and a theoretical thermal response current signal curve is considered, finally, a theoretical thermal response current signal curve with the highest fitting degree with the actually measured thermal response current signal curve is found, and the interface thermal resistance value corresponding to the theoretical thermal response current signal curve is 7.39 multiplied by 10-6K·m2W, the actual interface thermal resistance value of the dielectric film and the metal foil is 7.39 multiplied by 10-6K·m2/W。
The foregoing detailed description of the preferred embodiments of the invention has been presented. It should be understood that numerous modifications and variations could be devised by those skilled in the art in light of the present teachings without departing from the inventive concepts. Therefore, the technical solutions available to those skilled in the art through logic analysis, reasoning and limited experiments based on the prior art according to the concept of the present invention should be within the scope of protection defined by the claims.
Claims (10)
1. A method for measuring thermal resistance of a dielectric film and a metal interface is characterized in that the dielectric film is tightly attached to a metal foil, a heat pulse experiment is carried out, and an actually measured thermal response current signal curve is obtained through measurement; establishing a heat transfer simulation model of the dielectric film and the metal foil; and adjusting the interface thermal resistance value of the heat transfer simulation model, calculating to obtain a theoretical thermal response signal current curve, and taking the interface thermal resistance value of the heat transfer simulation model as the actual interface thermal resistance value of the dielectric film and the metal foil when the fitting degree of the actually measured thermal response current signal curve and the theoretical thermal response signal current curve meets the preset convergence condition.
2. The method for measuring thermal interface resistance between a dielectric film and a metal as claimed in claim 1, wherein obtaining a measured thermal response current signal curve comprises the steps of:
a1: after the medium film is subjected to metallization treatment, the medium film is tightly attached to the metal foil, and direct-current voltages are externally connected to two sides of the medium film, so that an electric field which is uniformly distributed exists in the medium film;
a2: applying pulse laser to heat the medium film, and collecting the actually measured thermal response current of the medium film;
a3: and converting the actually measured thermal response current into a frequency domain signal through fast Fourier transform to obtain an actually measured thermal response current signal curve.
3. The method for measuring thermal resistance of a dielectric film-metal interface according to claim 2, wherein in step a3, the method further comprises performing distortion compensation on the measured thermal response current to obtain a calibrated measured thermal response current, and then converting the calibrated measured thermal response current into a frequency domain signal to obtain a measured thermal response current signal curve, wherein the formula is as follows:
Im~(f)=I0~(f)*H~(f)
wherein, Im~(f) Representing the collected measured thermally responsive current, I0~(f) Represents the thermal response current in the ideal case, H~(f) Is a transfer function determined from the effect of the external circuit and amplifier bandwidth limitations on the thermally responsive current.
4. The method for measuring the thermal interface resistance between the dielectric film and the metal as claimed in claim 1, wherein the establishing of the heat transfer simulation model of the dielectric film and the metal foil specifically comprises: the method comprises the steps of obtaining initial temperature distribution, thickness, thermal diffusion coefficient, thermal conductivity coefficient, relative dielectric coefficient, temperature coefficient and thermal expansion coefficient of the relative dielectric coefficient of a dielectric film, obtaining initial temperature distribution, thickness, thermal diffusion coefficient and thermal conductivity coefficient of a metal foil, obtaining parameters of laser pulse and electric field in a thermal pulse experiment process, establishing a double-layer structure of the dielectric film and the metal foil in simulation software to serve as a heat transfer simulation model, initializing interface thermal resistance values of the heat transfer simulation model, and simulating disturbance of an externally applied direct current electric field and laser pulse in the double-layer structure in the simulation software.
5. The method for measuring thermal interface resistance between a dielectric film and a metal as claimed in claim 4, wherein the step of obtaining a theoretical thermal response current signal curve comprises the steps of:
b1: establishing a one-dimensional heat conduction equation of a double-layer structure of a heat transfer simulation model:
the boundary conditions are as follows:
the initial conditions were:
θ1=T1(x) (0<x<x1) (t=0)
θ2=T2(x) (x1<x<x2) (t=0)
wherein, theta1Is the temperature distribution of the dielectric film, a1Is the thermal diffusion coefficient, θ, of the dielectric film2Is the temperature distribution of the metal foil, a2Is the thermal diffusion coefficient of the metal foil, x is the position from the dielectric film to the thickness direction of the metal foil, x1Equal to the thickness, x, of the dielectric film2Is equal to the sum of the thicknesses of the dielectric film and the metal foil, t is time, k1Is the thermal conductivity, k, of the dielectric film2Is the thermal conductivity of the metal foil, h12Is the heat exchange coefficient, T, of the dielectric film and the metal foil1(x) Is the initial temperature distribution, T, of the dielectric film2(x) Is the initial temperature distribution of the metal foil, eta is the absorption coefficient of the pulse laser, q is the energy of the pulse laser, and delta (t) is the Gaussian distribution function of the pulse laser;
b2: solving the one-dimensional heat conduction equation to obtain the theoretical calculation of the temperature distribution of the dielectric filmTheoretical calculation of temperature distribution based on dielectric filmObtaining the temperature change inside the dielectric film, and calculating according to the following formula to obtain the theoretical thermal response current:
g(x)=ε0εr(αε-αx)E(x)
wherein I (t) represents a thermal response current, A represents an area of the dielectric film irradiated by the laser, d represents a thickness of the dielectric film, Δ θ represents a temperature change inside the dielectric film, g (x) represents a distribution function, E (x) represents an electric field distribution inside the dielectric film, ε0Denotes the vacuum dielectric constant,. epsilonrDenotes the relative dielectric coefficient, alpha, of the dielectric filmεTemperature coefficient, alpha, representing the relative dielectric coefficientxRepresents the thermal expansion coefficient of the dielectric film;
b3: and performing fast Fourier transform on the theoretical thermal response current and converting the theoretical thermal response current into a frequency domain signal to obtain a theoretical thermal response current signal curve.
6. The method for measuring the thermal resistance of the interface between the dielectric film and the metal as claimed in claim 1, wherein the step of adjusting the interface thermal resistance value of the heat transfer simulation model and determining the actual interface thermal resistance value of the dielectric film and the metal foil is specifically as follows:
c1: obtaining a theoretical thermal response current signal curve based on the current interface thermal resistance value of the heat transfer simulation model;
c2: fitting the theoretical thermal response current signal curve with the actually measured thermal response current signal curve, if the fitting degree does not meet the preset convergence condition, adjusting the interface thermal resistance value in the heat transfer simulation model, and repeating the step C1, otherwise, executing the step C3;
c3: and taking the interface thermal resistance value of the heat transfer simulation model as the actual interface thermal resistance value of the dielectric film and the metal foil.
7. The method as claimed in claim 6, wherein the Nelder-Mead simplex algorithm is used to fit the theoretical thermal response current signal curve to the measured thermal response current signal curve in step C2.
8. The method as claimed in claim 7, wherein the predetermined convergence condition is a default convergence condition in Nelder-Mead simplex algorithm.
9. A system for measuring thermal interface resistance between a dielectric film and a metal, which is based on the method for measuring thermal interface resistance between the dielectric film and the metal as claimed in any one of claims 1 to 8, and comprises:
the actual measurement module is used for tightly attaching the dielectric film to the metal foil, performing a heat pulse experiment and measuring to obtain an actual measurement heat response current signal curve;
the modeling module is used for establishing a double-layer structure of the dielectric film and the metal foil according to the parameters of the dielectric film and the metal foil and the parameters of the heat pulse experiment;
the calculation module is used for calculating a theoretical thermal response current signal curve of the heat transfer simulation model;
and the fitting module is used for fitting the actually measured thermal response current signal curve and the theoretical thermal response current signal curve and judging whether the fitting degree meets the preset convergence condition or not.
10. The system for measuring thermal resistance of a dielectric film and a metal interface of claim 9, wherein the actual measurement module obtains an actual measurement thermal response current of a thermal pulse experiment, performs distortion compensation on the actual measurement thermal response current to obtain a calibrated actual measurement thermal response current, and then converts the calibrated actual measurement thermal response current into a frequency domain signal to obtain an actual measurement thermal response current signal curve.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110500544.4A CN113419120B (en) | 2021-05-08 | 2021-05-08 | Method and system for measuring thermal resistance of dielectric film and metal interface |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110500544.4A CN113419120B (en) | 2021-05-08 | 2021-05-08 | Method and system for measuring thermal resistance of dielectric film and metal interface |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113419120A true CN113419120A (en) | 2021-09-21 |
CN113419120B CN113419120B (en) | 2022-10-25 |
Family
ID=77712168
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110500544.4A Active CN113419120B (en) | 2021-05-08 | 2021-05-08 | Method and system for measuring thermal resistance of dielectric film and metal interface |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113419120B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114264695A (en) * | 2022-01-28 | 2022-04-01 | 同济大学 | Method and system for measuring heat conductivity coefficient of trace liquid |
CN114460131A (en) * | 2022-01-28 | 2022-05-10 | 同济大学 | Method and device for measuring cross-scale solid heat conductivity coefficient |
Citations (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5015950A (en) * | 1989-07-21 | 1991-05-14 | Iowa State University Research Foundation, Inc. | Method and apparatus for determining thermal resistance and structural integrity of coatings on conducting materials by monitoring electrical conductance of the underlying material upon localized heating of the overlying coating |
US20020080850A1 (en) * | 2000-10-17 | 2002-06-27 | Nat'l Inst. Of Advanced Industrial Sci. And Tech | Method for measuring thermal diffusivity and interface thermal resistance |
US20030165179A1 (en) * | 2002-03-01 | 2003-09-04 | Danley Robert L. | System and method for calibrating contact thermal resistances in differential scanning calorimeters |
JP2007019327A (en) * | 2005-07-08 | 2007-01-25 | Kansai Electric Power Co Inc:The | High heat-resistant film capacitor |
JP2010243482A (en) * | 2009-03-19 | 2010-10-28 | National Institute For Materials Science | Thin film thermophysical property measuring instrument and method of measuring thermal conductivity and interface thermal resistance using this measuring instrument |
CN102590724A (en) * | 2012-01-08 | 2012-07-18 | 中国石油大学(华东) | Method for accurately measuring interface thermal resistance of semiconductor thin film |
CN103245694A (en) * | 2013-05-13 | 2013-08-14 | 北京工业大学 | Method for measuring thermal contact resistance between semiconductor device and contact material |
CN103759851A (en) * | 2014-02-18 | 2014-04-30 | 悉雅特万科思自动化(杭州)有限公司 | Temperature compensation method used for thermal resistance collecting module |
CN107490595A (en) * | 2017-09-21 | 2017-12-19 | 北京工业大学 | A kind of normal load influences the assay method and device of relation on faying face thermal contact resistance |
US20180017511A1 (en) * | 2015-02-06 | 2018-01-18 | Saint-Gobain Isover | Determination of the thermal resistance of a wall |
CN107688039A (en) * | 2017-07-14 | 2018-02-13 | 昆明理工大学 | The test system and its method of testing of light sheet material thermal conductivity factor and interface resistance |
CN108226218A (en) * | 2016-12-09 | 2018-06-29 | 上海大学 | A kind of thermal resistance measurement method and system of electronic device |
CN108363849A (en) * | 2018-01-31 | 2018-08-03 | 电子科技大学 | A kind of method for extracting thermal resistance and system |
CN108535313A (en) * | 2018-02-11 | 2018-09-14 | 中国矿业大学 | A method of with interface resistance between heat-pole method two solids of measurement |
CN108593707A (en) * | 2018-05-29 | 2018-09-28 | 哈尔滨工业大学深圳研究生院 | A kind of measurement method and device of GaN epitaxy wafer interface thermal resistance |
CN208297403U (en) * | 2018-05-29 | 2018-12-28 | 哈尔滨工业大学深圳研究生院 | A kind of measuring device of GaN epitaxy wafer interface thermal resistance |
CN109283216A (en) * | 2018-10-12 | 2019-01-29 | 广州特种承压设备检测研究院 | A kind of measurement method and device of grapheme material interface resistance |
CN109557129A (en) * | 2018-10-29 | 2019-04-02 | 同济大学 | A kind of measurement method of film thermal diffusion coefficient |
CN109959676A (en) * | 2017-12-14 | 2019-07-02 | 核工业西南物理研究院 | A kind of graphite and graphite film material thermal contact resistance test method |
CN110057862A (en) * | 2019-04-24 | 2019-07-26 | 苏州浪潮智能科技有限公司 | A kind of method and system using unary linear regression equation test thermal contact resistance |
CN209821110U (en) * | 2019-03-05 | 2019-12-20 | 佛山市联动科技实业有限公司 | Thermal resistance testing device |
CN111505399A (en) * | 2020-04-28 | 2020-08-07 | 同济大学 | Method for measuring space charge distribution of polymer dielectric film |
CN111537561A (en) * | 2020-06-17 | 2020-08-14 | 清华大学 | Method and system for measuring interface thermal resistance |
US20210048400A1 (en) * | 2019-08-16 | 2021-02-18 | Illumina, Inc. | Method for measuring thermal resistance between a thermal component of an instrument and a consumable |
CN112415046A (en) * | 2020-10-30 | 2021-02-26 | 同济大学 | System and method for measuring longitudinal thermal diffusion coefficient of film based on medium detector |
-
2021
- 2021-05-08 CN CN202110500544.4A patent/CN113419120B/en active Active
Patent Citations (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5015950A (en) * | 1989-07-21 | 1991-05-14 | Iowa State University Research Foundation, Inc. | Method and apparatus for determining thermal resistance and structural integrity of coatings on conducting materials by monitoring electrical conductance of the underlying material upon localized heating of the overlying coating |
US20020080850A1 (en) * | 2000-10-17 | 2002-06-27 | Nat'l Inst. Of Advanced Industrial Sci. And Tech | Method for measuring thermal diffusivity and interface thermal resistance |
US20030165179A1 (en) * | 2002-03-01 | 2003-09-04 | Danley Robert L. | System and method for calibrating contact thermal resistances in differential scanning calorimeters |
JP2007019327A (en) * | 2005-07-08 | 2007-01-25 | Kansai Electric Power Co Inc:The | High heat-resistant film capacitor |
JP2010243482A (en) * | 2009-03-19 | 2010-10-28 | National Institute For Materials Science | Thin film thermophysical property measuring instrument and method of measuring thermal conductivity and interface thermal resistance using this measuring instrument |
CN102590724A (en) * | 2012-01-08 | 2012-07-18 | 中国石油大学(华东) | Method for accurately measuring interface thermal resistance of semiconductor thin film |
CN103245694A (en) * | 2013-05-13 | 2013-08-14 | 北京工业大学 | Method for measuring thermal contact resistance between semiconductor device and contact material |
CN103759851A (en) * | 2014-02-18 | 2014-04-30 | 悉雅特万科思自动化(杭州)有限公司 | Temperature compensation method used for thermal resistance collecting module |
US20180017511A1 (en) * | 2015-02-06 | 2018-01-18 | Saint-Gobain Isover | Determination of the thermal resistance of a wall |
CN108226218A (en) * | 2016-12-09 | 2018-06-29 | 上海大学 | A kind of thermal resistance measurement method and system of electronic device |
CN107688039A (en) * | 2017-07-14 | 2018-02-13 | 昆明理工大学 | The test system and its method of testing of light sheet material thermal conductivity factor and interface resistance |
CN107490595A (en) * | 2017-09-21 | 2017-12-19 | 北京工业大学 | A kind of normal load influences the assay method and device of relation on faying face thermal contact resistance |
CN109959676A (en) * | 2017-12-14 | 2019-07-02 | 核工业西南物理研究院 | A kind of graphite and graphite film material thermal contact resistance test method |
CN108363849A (en) * | 2018-01-31 | 2018-08-03 | 电子科技大学 | A kind of method for extracting thermal resistance and system |
CN108535313A (en) * | 2018-02-11 | 2018-09-14 | 中国矿业大学 | A method of with interface resistance between heat-pole method two solids of measurement |
CN108593707A (en) * | 2018-05-29 | 2018-09-28 | 哈尔滨工业大学深圳研究生院 | A kind of measurement method and device of GaN epitaxy wafer interface thermal resistance |
CN208297403U (en) * | 2018-05-29 | 2018-12-28 | 哈尔滨工业大学深圳研究生院 | A kind of measuring device of GaN epitaxy wafer interface thermal resistance |
CN109283216A (en) * | 2018-10-12 | 2019-01-29 | 广州特种承压设备检测研究院 | A kind of measurement method and device of grapheme material interface resistance |
CN109557129A (en) * | 2018-10-29 | 2019-04-02 | 同济大学 | A kind of measurement method of film thermal diffusion coefficient |
CN209821110U (en) * | 2019-03-05 | 2019-12-20 | 佛山市联动科技实业有限公司 | Thermal resistance testing device |
CN110057862A (en) * | 2019-04-24 | 2019-07-26 | 苏州浪潮智能科技有限公司 | A kind of method and system using unary linear regression equation test thermal contact resistance |
US20210048400A1 (en) * | 2019-08-16 | 2021-02-18 | Illumina, Inc. | Method for measuring thermal resistance between a thermal component of an instrument and a consumable |
CN113287007A (en) * | 2019-08-16 | 2021-08-20 | 伊鲁米纳公司 | Method for measuring thermal resistance between thermal component and consumable of instrument |
EP3881060A1 (en) * | 2019-08-16 | 2021-09-22 | Illumina Inc | Method for measuring thermal resistance between a thermal component of an instrument and a consumable |
CN111505399A (en) * | 2020-04-28 | 2020-08-07 | 同济大学 | Method for measuring space charge distribution of polymer dielectric film |
CN111537561A (en) * | 2020-06-17 | 2020-08-14 | 清华大学 | Method and system for measuring interface thermal resistance |
CN112415046A (en) * | 2020-10-30 | 2021-02-26 | 同济大学 | System and method for measuring longitudinal thermal diffusion coefficient of film based on medium detector |
Non-Patent Citations (7)
Title |
---|
D. P. H. HASSELMAN, KIMBERLY Y. DONALDSON, FRED D. BARLOW, AICHA: "Interfacial Thermal Resistance and Temperature Dependence of Three Adhesives for Electronic Packaging", 《IEEE TRANSACTIONS ON COMPONENTS AND PACKAGING TECHNOLOGIES》 * |
JIAN CHEN, WEI ZHANG, ZHILI FENG, WAYNE CAI: "Determination of thermal contact conductance between thin metal sheets of battery tabs", 《INTERNATIONAL JOURNAL OF HEAT AND MASS TRANSFER》 * |
傅德舟,郑飞虎,安振连,张冶文: "基于瞬态放电电流法的聚丙烯薄膜中空间电荷密度的动态监测", 《内蒙古科技大学学报》 * |
朱长纯等: "碳纳米管薄膜电阻的热敏特性测试", 《西安交通大学学报》 * |
李想,吴钢,周刚,毕柯,汤智胤,张青枝,马计: "基于Labview的激光光热法界面热阻测量系统", 《低温工程》 * |
王镇,莫德锋.汪洋,徐红艳,李雪,余利泉, 李俊,王小坤: "界面热阻测量方法及影响因素研究进展", 《工程热物理学报》 * |
蔡浩原等: "高阻值纳米薄膜材料的热电特性测量", 《光学精密工程》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114264695A (en) * | 2022-01-28 | 2022-04-01 | 同济大学 | Method and system for measuring heat conductivity coefficient of trace liquid |
CN114460131A (en) * | 2022-01-28 | 2022-05-10 | 同济大学 | Method and device for measuring cross-scale solid heat conductivity coefficient |
CN114460131B (en) * | 2022-01-28 | 2023-08-29 | 同济大学 | Method and device for measuring trans-scale solid heat conductivity coefficient |
CN114264695B (en) * | 2022-01-28 | 2023-08-29 | 同济大学 | Method and system for measuring micro-liquid heat conductivity coefficient |
Also Published As
Publication number | Publication date |
---|---|
CN113419120B (en) | 2022-10-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113419120B (en) | Method and system for measuring thermal resistance of dielectric film and metal interface | |
Braun et al. | On the steady-state temperature rise during laser heating of multilayer thin films in optical pump–probe techniques | |
CN108614005B (en) | Method and system for testing thermophysical property parameters based on multilayer composite material | |
CN104034749B (en) | Based on the method for testing of thermal contact resistance between the layer material of 3 ω methods | |
CN109557129B (en) | Method for measuring film thermal diffusion coefficient | |
CN108169279A (en) | One kind is based on VO2The thermal conductivity of thin film measuring device and method of film | |
CN111537561B (en) | Method and system for measuring interface thermal resistance | |
CN107037079B (en) | A kind of support beam type MEMS compound sensor and its preparation and test method | |
CN108226219B (en) | Method for detecting heat generation uniformity of film resistor | |
CN112415046A (en) | System and method for measuring longitudinal thermal diffusion coefficient of film based on medium detector | |
CN113092523B (en) | Device and method for testing heat conduction performance of film material | |
CN114264695B (en) | Method and system for measuring micro-liquid heat conductivity coefficient | |
Siroka et al. | Comparison of thin film heat flux gauge technologies emphasizing continuous-duration operation | |
CN106053529A (en) | Device and method for measuring porous metal material heat conductivity coefficient through comparison plate | |
CN103698357A (en) | Thermal conductivity and thermal diffusivity sensor based on MEMS double heater | |
CN114460131B (en) | Method and device for measuring trans-scale solid heat conductivity coefficient | |
CN109269682A (en) | A kind of caliberating device and scaling method of heat flow transducer | |
Aberle et al. | Implementation of an in-situ infrared calibration method for precise heat transfer measurements on a linear cascade | |
Jha et al. | Visualization and measurement of natural convection boundary layer by particle image velocimetry | |
CN114384118B (en) | Thermal diffusion coefficient measuring method and device for dielectric film with substrate | |
CN114371379A (en) | Method and system for measuring space charge injection threshold electric field | |
CN114660127A (en) | Material identification sensor and method for identifying material attribute by using same | |
CN111157573B (en) | Measuring device and measuring method for film thermal conductivity | |
CN117368255A (en) | System and method for testing heat conductivity coefficient of filiform or thin-film material | |
Parler et al. | Predicting operating temperature and expected lifetime of aluminum-electrolytic bus capacitors with thermal modeling |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |