The present invention relates to a method for determining the behaviour of an electrical or electronic power component, more particularly of an inductor, operating under typical operating conditions of power generators, wherein that inductor is used as an energy periodic temporary storage element, to carry out a static energy conversion process.
The present invention also relates to a measuring station and a system configured to implement such a method.
Before going into the merits of the present invention it is noted that explicit reference will be made hereinafter to an inductor as an element of periodic temporary storage of energy. However, it must be kept in mind that the present invention can be applied in a completely analogous way also to capacitors as will become clear in the following.
The operating conditions of an inductor which are of particular interest for the present invention are those imposed by a periodic charge-discharge process. In said process, a zero mean square wave time-varying voltage vL(t) is applied to the inductor terminals, for example by means of active and passive semiconductor devices, which operate as electronic switches that alternately connect the terminals of the inductor to electronic components configured to impose a positive or negative voltage value, thus obtaining in a first step, an energy transfer from an energy source towards the inductor and, in the second step, an energy transfer from the inductor to a load.
The electrical behaviour of the inductor, resulting from this periodic charge-discharge process, can be determined based on the trend of the current iL(t) flowing therethrough, which current varies over time. The current flowing through such an inductor during the periodic charge-discharge process is comparable to a periodic triangular wave, with a positive slope in the charging step of the inductor and a negative slope in the discharge step of the inductor. The triangular wave of the current flowing through the inductor can be subjected to more or less accentuated slope variations, in each of the two steps (charge and/or discharge), due to the more or less abrupt occurrence of the saturation phenomenon of the inductor electronic core, as the variable current flowing therethrough varies.
As a rule, the behaviour of a commercially available inductor is experimentally determined by its manufacturer, using different types of techniques, with the aim of quantifying one of its characteristic parameters, specifically the self-inductance coefficient L. The self-inductance coefficient is usually represented in graphical form as a function of the variable average current flowing through the inductor and parametrized as a function of the environmental temperature.
Traditionally, a manufacturer of a commercially available inductor also obtains, from the results of experimental tests that he performs on the component, useful information for the formulation of heuristic loss equations (for example through Steinmetz formulas), depending on the conditions in which such component can operate, for example when the frequency varies.
Over time, techniques and apparatuses have been developed for determining the behaviour of inductors, in relation to the varying conditions in which they operate.
These determination techniques with their respective apparatuses, however, suffer from many drawbacks. For example, in some cases they are not informative, as they provide parameters that do not assume unique values, as they are subject to variations according to the operational parameters of the components. In particular, reference is made to the inductor’s self-induction coefficient L, expressed at most as a function of the average current flowing therethrough and the environmental temperature. The self-induction coefficient L is traditionally provided as it is usually used to calculate the performance of inductors, such as peak-to-peak current ripple and losses, although, in fact, it does not assume a unique value as it varies according to, for example, the working frequency, the voltage and the duty-cycle, as well as the current and the temperature.
Grids of parameter values are also traditionally used to generate signals to be sent to the inductors during the execution of experimental tests thereon. Those grids cover multidimensional domains, the exploration of which is carried out through nested scanning cycles of each quantity, which nested scanning cycles may include — if not appropriately restricted — n-tuple of operational quantities corresponding to unsustainable operational conditions of the inductor being evaluated and, therefore, such that they compromise the test results due to inductor failures that may occur and the consequent possible damage to the hardware used for the test. This type of approach is also very demanding from a computational point of view.
The determination of the behaviour of the inductors is also traditionally performed by carrying out tests on these components in small-signal sinusoidal regime, which tests however do not reflect the real working conditions of the inductors in power generators.
Complex architectures of test systems are also adopted to determine the behaviour of these inductors, which architectures require a programmable electronic load for the imposition of the desired current on the inductor under test and for the automation of the test cycles.
Some examples of the above traditional apparatuses and methods are, for example, taught in EP2807589B1, US7307412B1 and US20120316817A1. Further traditional methods for determining a characteristic parameter of a device are described in US 5793640 and in US 2016/0132625.
In any case, all the methods traditionally used provide information that does not express in a direct and easy to understand way, the relationship between the operational conditions to which an inductor is normally subjected (for example, the frequency of oscillation of the voltage, average current and the environmental temperature) and its response parameters such as, for example, current ripple, dissipated power and surface temperature.
The need is therefore felt to improve the state of the art in the field of electrical or electronic power components and, thus, the main object of the present invention is to provide a method for determining the behaviour of an electrical or electronic power component which is safe, in the sense that its implementation does not cause faults in the electrical or electronic power component under test or in the used hardware, that is less demanding from a computational point of view, that is reliable and provides for more easily interpretable and useful results for the purpose, compared to traditional methods.
A further object of the present invention is to provide for a measuring station, configured for the implementation of the aforementioned method.
Not the least object of the present invention is to provide a system for the implementation of the method for determining the behaviour of an electrical or electronic power component, which is easy to implement and use.
It is a specific object of the present invention a method for determining the behaviour of one electrical or electronic power component, with respect to a working limit condition, the method comprising the following operational steps:
According to another aspect of the invention, said operational parameters of interest for said electrical or electronic power component can be selected among:
According to a further aspect of the invention, said three-dimensional mathematical space can be defined by assigning to each one of said operational parameters of interest of said electrical or electronic power component, one corresponding mathematical axis of an Euclidean three-dimensional mathematical space.
According to an additional aspect of the invention, said field of exploration can be defined by selecting:
According to another aspect of the invention, the exploration of said three-dimensional mathematical space can be carried out, starting from one starting point Pn of said field of exploration, along said directions of exploration, in an orderly way, and according to one exploration rule such that the coordinates of an explored n-th point Pn of said three-dimensional mathematical space differ from those of a previous point Pn-1, for the value of one or more of its components.
According to a further aspect of the invention, the coordinates of said points Pn of said field of exploration can have values comprised between one minimum value and one maximum value of respective operational parameters and the finite number of said points Pn of said field of exploration can be a function, for each one of said directions of exploration of said three-dimensional mathematical space, of
According to an additional aspect of the invention, for each direction of said directions for exploration:
According to another aspect of the invention, the compliance or not of said working limit condition in any point Pn of said three-dimensional mathematical space can depend on the value of at least one response parameter of said set R of response parameters of said electrical or electronic power component, the response parameters being determined based on the response detected from said electrical or electronic power component, after the application of said stimulus to its terminals, and on a pre-set search logic.
According to a further aspect of the invention, said response parameter of said set R of response parameters of said electrical or electronic power component can be one among:
According to an additional aspect of the invention, said working limit condition of said electrical or electronic power component can be a function of a prefixed threshold value for each response parameter.
According to another aspect of the invention, said at least one response parameter can be compliant with said at least one working limit condition if its value is lower than or equal to a respective threshold value.
According to another aspect of the invention, said step C can comprise:
According to an additional aspect of the invention, the stored coordinates of said point P*n of said mathematical space, which point is comprised in the neighbourhood of said point Pn of said field of exploration, can be calculated through an interpolation formula, optionally a linear one.
According to another aspect of the invention, said step D can comprises applying to said points thereby stored, through said data control and processing unit, at least one Genetic Programming or Grammatical Evolution algorithm.
It is further a specific object of the invention, a measuring station, configured for implementing the method for determining the behaviour of an electrical or electronic power component as described above, comprising:
According to another aspect of the invention, said one stimulus generating device and said detecting device can be obtained through three power converter stages operatively connected in cascade according to an Opposition Method, so as to:
According to a further aspect of the invention, said three power converter stages can comprise an Input Stage, a Test Stage, and an Output Stage, wherein the Test Stage is connected between the Input Stage and the Output stage and is further configured to be connected to said electrical or electronic power component.
According to an additional aspect of the invention, said measuring station can comprise downstream of said Output Stage one switching element for the current output from said Output Stage, toward the input of the Test Stage or an external load.
According to a further aspect of the invention, said Input Stage can be configured to work in closed loop and provide in input to said Test Stage one direct voltage, through the adjustment of the output voltage thereof, and said Test stage is configured to work in open loop and provide for a switching frequency of said zero mean and square wave voltage, through the adjustment of its own frequency and duty cycle.
According to another aspect of the invention, said Output stage can be configured to operate in closed loop and impose at the output of said Test Stage a direct current, through the adjustment of its own input average current.
According to an additional aspect of the invention, said measuring station can comprise:
According to a further aspect of the invention, said positioning system can comprise:
According to a further aspect of the invention, said elastic anchoring group can comprise one couple of elastically charged screws, configured for being screwed on said printed circuit, passing through suitable openings obtained in said positioning plates.
According to an additional aspect of the invention, said conductive paths can have a polygonal configuration, optionally a trapezoidal configuration.
According to another aspect of the invention, said detecting device can comprise at least one temperature sensor, operatively connected to said data control and processing unit, and configured for detecting said surface temperature of said electrical or electronic power component and for transmitting it to said data control and processing unit.
It is also a specific object of the invention a system for implementing the method for determining the behaviour of an electrical or electronic power component as described above, comprising at least one measuring station as described above and at least one remote processing unit, wherein said measuring station and said remote processing unit are operatively connected to each other, optionally through one network, optionally through one cable or wirelessly, and wherein said data control and processing unit of said measuring station is configured to send to said at least one remote processing unit said coordinates of said points P*n thereby stored, and said at least one remote processing unit is configured to apply at least one mathematical algorithm, optionally a Genetic Programming or Grammatical Evolution algorithm, providing for in output a description, in analytical form, of the locus of said points P*n, and therefore the locus of the corresponding operational parameters that determine a response from said electrical or electronic power component that meets at least one working limit condition.
The present invention will be now described, for illustrative but not limiting purposes, according to its preferred embodiments, with particular reference to the Figures of the accompanying drawings, wherein:
In the enclosed Figures the same reference numbers will be used for similar elements.
Before going into the merits of the description, it is specified that in this description the terms “inductor”, “electrical or electronic power component”, “electrical or electronic component” or “component” may be used in a completely equivalent way, that is in any case clear from the context.
It is also specified that, according to the method of the present invention, an inductor operational in an environment having temperature Ta, can be subjected to a stimulus comprising:
With particular reference to
The response of an inductor operational at an environmental temperature Ta, to which the stimulus described above has been applied, i.e. a variable voltage vL(t) as described above and a certain average current IL,av, can be described through the following parameters, again shown in
This said, the method according to the present invention, indicated by reference number 1 is provided for determining the behaviour of an electrical or electronic power component (represented in
More specifically, according to a preferred embodiment of the method of the invention, the operational parameters of interest for the inductor, the behaviour of which is to be determined, are three and are selected from:
On the basis of three of the above operational parameters of interest for the inductor 2, the three-dimensional mathematical space 3 is defined, according to method 1 of the invention, by assigning to each of these operational parameters of interest, a corresponding mathematical axis of a Euclidean three-dimensional mathematical space 3.
Then, method 1 according to the invention comprises a second step B, which is always performed for example by means of a data control and processing unit, wherein an exploration field 4 of the three-dimensional mathematical space 3 is defined. At step B of the invention method 1, a working limit condition of interest for the inductor and a set R of response parameters for that inductor are also defined, as will be seen below.
In this regard, according to a preferred embodiment of the invention, the exploration field 4 is defined by selecting:
According to a particularly advantageous aspect of the invention, the coordinates of the points Pn of the exploration field 4 assume values between a minimum value and a maximum value of the respective operational parameters (Veq, fsw, IL). Specifically, therefore:
Within the finite set of points Pn that make up the exploration field 4 there will be, therefore, a point with minimum coordinates, for example given by (Veq,min, fsw,min, IL,min) and a point of coordinates maximum, for example (Veq,max, fsw,max, IL,max) at a given environmental temperature Ta.
The finite number of points Pn included in the exploration field 4 is a function, for each of the exploration directions of the three-dimensional mathematical space 3, of:
In other words, the number of samples (NVeq, Nfsw, NIL)and the offset (ΔVeq, Δfsw, ΔIL) for each one of the three directions of exploration, determine the coordinates, within the mathematical space 3, of the points Pn included in the exploration field 4.
Specifically, according to method 1 of the present invention, for each one of the three directions of exploration of the exploration field 4, the number of samples (NVeq, Nfsw, NIL) can be predetermined and the offset (ΔVeq, Δfsw, ΔIL) can be also prefixed and constant, or the number of samples (NVeq, Nfsw, NIL) can be prefixed while the offset (ΔVeq,i, Δfsw,j, ΔIL,k) with i = 1,..,NVeq-1, j = 1,..,Nfsw-1, k = 1,..,NIL-1, between two consecutive samples (corresponding to two consecutive points Pn-1 and Pn of the exploration field 4) can vary, optionally according to a selection between a logarithmic, power, trigonometric, transcendent or predetermined numerical series function. According to a further variant of method 1 of the invention, the number of samples (NVeq, Nfsw, NIL) may not be predetermined and may depend on an offset (ΔVeq,i, Δfsw,j, ΔIL,k,) with i = 1, ..NVeq-1, j = 1,..,Nfsw-1, k = 1,..,NIL-1, that varies between one sample and the next one, which is calculated each time during the exploration of the exploration field 4, for each point Pn, optionally based on a value assumed by a set R of response parameters (which will be better described below), calculated at previously explored points Pn-1 and Pn-2 of the exploration field 4.
At the end of step B, the method 1 of the invention comprises at step C to explore the three-dimensional mathematical space 3, starting from a starting point included in the exploration field 4 defined as above, in search of points P of the three-dimensional mathematical space 3 that meet at least one working limit condition for inductor 2. As will be seen in detail below, the exploration of the three-dimensional mathematical space 3 takes by means of:
More particularly, the exploration of the exploration field 4 is carried out in an orderly way (as will be better disclosed below), starting from one starting point Pn of said exploration field 4, along the directions of exploration of the ordered set of directions for exploration and according to one exploration rule such that the coordinates of an explored n-th point Pn of said three-dimensional mathematical space 3 differ from those of a previous point Pn-1, for the value of one or more of its components. Thus, for example:
According to a variation of the method 1 of the present invention, the exploration of the exploration field 4 can also be performed according to a different exploration rule for which the coordinates of an n-th point Pn differ from those of the previous point Pn-1, for the value assumed by two or all three components. For example, according to a variation of the method of the present invention, the exploration of the exploration field 4 can be performed over parallel planes wherein, then, the points Pn of the exploration field 4, which lie on a plane parallel to the one on which two axes of three-dimensional mathematical space 3 lie, are first explored starting from a starting point of minimum coordinates (for example starting from P1 having coordinates (Veq,min, fsw,min, IL,min) until a point Pn of maximum coordinates in that plane is reached (thus arriving at a point of coordinates (Veq,max, fsw,max, IL,min) and then passing to the adjacent parallel plane and exploring this adjacent parallel plane starting from a point with minimum coordinates in that plane Pn+1= (Veq,min, fsw,min, IL,min+ΔIL). It is understood that the one provided is a purely indicative and non-limiting example for the present invention and the person skilled in the art will have no difficulty in understanding how it is possible to explore the points of the exploration field 4 also in other ways, for example starting from a point of maximum coordinates for the field of exploration or in a plane thereof as described above.
With regard to the working limit condition, it is specified that whether or not that working limit condition is met in any point Pn of the three-dimensional mathematical space 3 (including those of the exploration field 4) depends on the value assumed by at least one response parameter of the set R of response parameters of the inductor 2, following the application to its terminals of at least one stimulus as described above (i.e. comprising a square wave voltage vL(t) and an average current IL), and by a pre-set search logic, as explained below.
The response parameters that are comprised in the above said set R are one or more between:
It should be noted that, according to a preferred embodiment of the invention, the aforementioned response parameters can be detected by a detecting device included in a measuring station which also forms part of the present invention, and which will be described below. More particularly, the peak-to-peak variation of the time-varying current(iL(t)) flowing through the inductor and the surface temperature Ts can be directly detected by respective analog components of the detecting device, while the dissipated power Pd can be calculated as a function of the stimulus variable voltage vL(t) applied to the inductor 2 and the variable current iL(t) detected by it. Alternatively, according to a variant of the present invention, the dissipated power Pd can also be determined in an alternative way. For example, according to a variant of the present invention, the dissipated power Pd can be determined in a more approximate but still reliable way, as a function of the surface temperature Ts of the inductor 2, in a way known to the man skilled in the art.
According to method 1 of the invention, the aforementioned working limit condition for inductor 2 is a function of a threshold value (ΔiLpp,ref, Ts,ref, Pd,ref) for a corresponding response parameter. More particularly, the working limit condition is met if one or more response parameters assume a value lower than or equal to the respective threshold value ΔiLpp,ref, Ts,ref, Pd,ref).
As far as the search logic is concerned, the working limit condition for the points P of the three-dimensional mathematical space 3 can concern, as mentioned above:
In other words, the exploration of three-dimensional mathematical space 3 can be carried out by considering three distinct working limit conditions according to the following different search logics:
According to an advantageous embodiment of the invention, the search logic is selected at step B or in any case before the execution of step C of invention method 1. Alternatively, the search logic can be pre-set.
Going into the merits of step C of method 1 of the present invention, it comprises:
With particular reference to the above-mentioned step C.7.2.1, according to a preferred embodiment of method 1 of the present invention, it is envisaged that the stored coordinates of each point P*n of the mathematical space 3, which is included in the neighbourhood of said point Pn of the exploration field 4, can be calculated, by means of an interpolation formula, optionally a linear one. Specifically in the aforementioned case wherein:
Of course, if the exploration of the exploration field 4 was along the second or third direction of exploration, instead of Veq and ΔVeq, the above formula would involve, respectively, fsw with Δfsw and IL with ΔIL
Once the coordinates of a point P*n have been determined, the above method 1 foresees, at step C.7.2.2, to select the previous point Pn-1 and, starting from this, a new point Pn in the mathematical space 3 along the second exploration direction, where the coordinates of the new point are given by (Veq[n-1], fsw[n-1]+ Δfsw, IL[n-1]) or ((Veq[n-1], fsw[n-1]- Δfsw, IL[n-1]) depending on the case. In fact, depending on which starting point is chosen, for the exploration of the exploration field 4, the values of the components of the coordinates of the points Pn included therein will be increased or reduced, to allow the method of the present invention to explore and select the points of three-dimensional mathematical space 3, starting from those of the exploration field 4, which comply with the working limit condition.
At step C.7.3.1 of method 1 of the invention, on the other hand, a new point Pn is selected in mathematical space 3 along the first direction of exploration, where the coordinates of the new point are determined starting from Pn and are given by (Veq[n] ]+ΔVeq, fsw[n], IL[n]) or (Veq[n] - ΔVeq, fsw[n], IL[n]) depending on the case. In this case also, in fact, depending on which starting point is chosen, for the exploration of the exploration field 4, the values of the coordinates of the points Pn included therein will be increased or reduced, to allow the method of the present invention to explore and select the points of three-dimensional mathematical space 3, starting from those of the exploration field 4, which comply with the working limit condition.
In other words, and as will also be seen in the following, the three-dimensional mathematical space 3 is explored, starting from the points of the exploration field 4, around its boundary surface which divides the three-dimensional mathematical space 3 into two half-spaces, one having points wherein the working limit condition is met and another having points wherein the working limit condition is not met.
The aforementioned step C.8 is about the case wherein a new selected point Pn of the mathematical creation space 3 is not included in the exploration field 4 has, that is, wherein a component of its coordinates is outside the range of minimum and maximum values mentioned above, for example in the case wherein the exploration of the exploration field 4 is carried out along the first exploration direction and the new point Pn has coordinates (Veq,max[n] +ΔVeq, fsw[n], IL[n]). Method 1 of the present invention comprises selecting a new point Pn of the three-dimensional mathematician, starting from this previously selected point Pn, along the second direction of exploration (step C.8.1) and, if this point Pn is included in the field exploration, the method returns to step C.2; otherwise the method 1 of the invention comprises selecting the new point Pn of the mathematical space, starting from the previously selected point Pn-1, along the third direction of exploration (step C.8.2) and, if the point direction of exploration is included in the exploration field 4, the method returns to step C.2. According to a further alternative, when the points Pn of the exploration field 4, that are located around the boundary surface have all been explored, the method can continue to the next step (step D).
As anticipated above at step C.7.1 of the invention method, if during the execution of method 1 it is found, only for the initial point Pn of the exploration field 4, that the working limit condition is not met, the method comprises searching for the first point Pn of the exploration field 4 wherein the working limit condition is met, along the first direction of exploration. In other words, at the start of the invention method 1, the exploration field 4 is explored in search of a boundary point of the three-dimensional mathematical space 3, at which the inductor 2 has a behaviour that complies with the aforementioned working limit condition.
It should be noted that thanks to the guided exploration of the exploration field 4 at step C, the method 1 of the present invention allows selecting a finite number of points (or combinations of operational parameters) within the mathematical space 3, for which the inductor 2 complies with at least one limit condition (which typically corresponds to a threshold within which inductor 2 works in compliance with the working limit condition). The number of points of the exploration field 4 at which the exploration is actually carried out can be, depending on the case, much lower than the total number of points of the same exploration field 4, so that with the method 1 of the present invention the computational complexity required to determine the behaviour of inductor 2 is considerably reduced.
Purely by way of a non-limiting example, reference is made to
In the provided example, therefore, the three-dimensional mathematical space 3 is explored in search of the points P*n (i.e., triples of operational parameters Veq, fsw, iL) which constitute the boundary of the region (half-space) of the three-dimensional mathematical space 3, wherein the aforementioned predetermined working limit condition is valid.
The operational parameters of interest are Veq, fsw and IL and vary in the ranges Veq,range = [ 1 V, 10 V ], fsw,range = [ 100 kHz, 1 MHz ] and IL,range = [ 1 A, 6 A ], with uniform offsets of the operational quantities having amplitude equal to ΔVeq = 1 V, Δfsw = 100 kHz, ΔIL = 1 A. The exploration field 4 is explored, in
In the graphs illustrated at
The result of the exploration is a list of coordinates of points P*n of mathematical space 3 (corresponding to triplets of operational parameters for the inductor 2 of interest), wherein the working limit condition ΔiLpp=ΔiLpp,ref and/or Ts=Ts,ref e/o Pd=Pd,ref in met. With reference to the example of
The method 1 of the present invention after step C, comprise the execution of step D wherein a mathematical model is determined, for example by means of the data control and processing unit, in the three-dimensional mathematical space 3, the mathematical model describing in analytic form the locus of the points P*n thus stored, thereby obtaining a locus 5 of the operational parameters which determine a response in the inductor 2 which compiles with the working limit condition.
What is obtained, therefore, is an analytical description of the surface 5 of the three-dimensional mathematical space 3 which surface divides this space into two half-spaces: a first half-space, which defines the operational conditions at which inductor 2 operates in compliance with the working limit condition, and another half-space opposite to the first half-space with respect to surface 5, which defines the operational conditions at which the inductor 2 does not operate in compliance with the working limit condition.
According to a particularly advantageous embodiment of the invention, step D of method 1 includes the application of at least one algorithm, for example of Genetic programming algorithm (J. R. Koza, Genetic Programming: On the Programming of Computers by Means of Natural Selection, MIT Press, Cambridge, 1992), to above mentioned points P*n stored during the previous step. It is quite clear, in this regard, that more than one algorithm belonging to the Genetic programming family of algorithms can be used for the execution of step D of the method according to the present invention. For example, instead of a preferred algorithm that will be disclosed by way of non-limiting example in the following, a Cartesian Genetic Programming algorithm (J.F. Miller, Cartesian Genetic Programming, Springer, ISBN 978-3-642-17310-3) or, alternatively to this algorithm family, a Grammatical Evolution algorithm (C. Ryan, M. O′Neill, Michael, J.J. Collins, Handbook of Grammatical Evolution, Springer, ISBN 978-3-319-78717-6), can be used, provided that such an algorithm is able to obtain a mathematical model that analytically describes the locus of the stored P*n points and, therefore, a locus 5 of the operational parameters which determine a response in the inductor 2 compiling with the working limit condition.
With reference to
According to a preferred example of method 1 of the present invention, the Genetic Programming algorithm allows to obtain a function that analytically describes the locus of the operational parameters that determine an inductor response characterized by a certain value of the peak-to-peak current ripple ΔiLpp=ΔiLpp,ref and/or its surface temperature Ts=Ts,ref and/or its dissipated power Pd=Pd,ref, starting from the varying over time response of inductor 2 (i.e. from the current iL(t)) experimentally acquired through a hardware system which will be described below, which is also the subject of the invention.
Therefore, as already mentioned above, the method 1 of the invention does not require the exhaustive execution of experimental tests on an entire mathematical space or three-dimensional domain of the values of the operational quantities {Veq, fsw, IL}, but only on a very reduce subset of their combinations. To provide an example, if the exploration field 4 included 10 samples for each of the operational quantities {Veq, fsw, IL}, compared to 1000 possible combinations (points Pn), the method of the present invention would test about 200 points.
A Genetic Programming algorithm that can be advantageously used in the invention method 1 can be developed ad hoc or provided by software widely used in the field, for example MatLAB. Such an algorithm, as is known, adopts for the representation of a mathematical model a tree structure like the one shown, by way of example, in
The tree that represents the mathematical model that is automatically generated by the Genetic Programming algorithm, providing as input a given set of elementary functions that operate on constant coefficients (c1, c2 e c3) and independent variables (x1 and x2).
In the case of the present invention, the pair of input variables (x1, x2) can be associated, through the data control and processing unit referred to above, to any pair of the above operational parameters ((Veq, fsw), (fsw, IL) or (Veq, IL), while the output can be associated with the remaining third operational parameter. The coefficients c1, c2, ... depend on the working limit condition of the inductor 2 under investigation, ΔiLpp=ΔiLpp,ref, or Ts=Ts,ref, or Pd=Pd,ref, for which the operational parameters triplet list has been generated, which describe a surface in the three-dimensional mathematical 3 space {Veq, fsw, IL} identified by the combinations of operational parameters at which a limit condition is compiled with. It follows that the Genetic programming algorithm can generate nine different mathematical models, three for each inductor response parameter, as follows:
The coefficients and independent variables identify the terminal nodes of the tree. The elementary functions identify the non-terminal nodes of the tree. The independent variables are the inputs to the model. The model output is the result of the sequence of operations defined by the elementary functions. Table 1 reports an exemplary, non-limiting list of possible elementary functions. A complexity index is associated with each elementary function, as indicated in Table 1 by way of a non-limiting example. The input variables are assigned a complexity index of 1. The complexity of the model is obtained as follows:
The Genetic programming algorithm (J. R. Koza, Genetic Programming: On the Programming of Computers by Means of Natural Selection, MIT Press, Cambridge, 1992) which can be advantageously implemented in the method of the present invention, once the model to be generated (M.1, ..., M.9) has been determined, operates for example on a population of at least 100 models, optionally 500 models, which evolve over time over at least 100 generations, optionally 300 generations. Each model is a combination of elementary functions and coefficients, represented by a tree with a number of nodes not exceeding 100 nodes, optionally not exceeding 50 nodes. Starting from an initial population of models generated randomly, at each generation the precision with which each model reproduces the reference data is evaluated, and the population is then replaced by a new generation by applying a selection operator (in particular the binary tournament operator), of a cross-over operator (in particular the subtree cross-over operator) and of a mutation operator (in particular the subtre¬e mutation operator), with equal cross-over probability of 80%, probability of subtree mutation equal to 18% and probability of node mutation equal to 2%.
The coefficients cx of the functions relating to the selected model can be determined, for the limit condition value {ΔiLpp,ref, Ts,ref, Pd,ref}, for example through the Levenberg - Marquardt method, which identifies the best values of the coefficients, on the basis of the error minimization criterion χ2 applied to the yn values of the dependent variable of the selected model obtained in the dataset of the independent variables on n samples of the independent variables (x1,i,x2,i) i = 1, ..., n corresponding to the list of triplets obtained through the experimental tests. For example, in the case of model M.1 the error is expressed by the following formula:
where Veq,i is the experimental i-th value of the operative variable Veq corresponding to the i-th experimental values ƒsw,i e lL,i of the operational quantities ƒsw and lL. The same formula applies to the remaining models M.2, ..., M.9.
It is evident that the advantage of models M.1, ..., M.9 consists in the fact that they allow putting into relation the operational parameters of an inductor in a direct and simpler way than traditional methods and verifying whether or not they meet certain response conditions. For example, in the case of model M.1, given Veq, fsw and lL, and given the function fveq,pd(fsw,iL) corresponding to a certain Pd,ref, graphically represented by surface 5, if Veq<fveq,pd(fsw,IL) then Pd<Pd,ref.
The model obtained according to the method of the present invention (for example fVeq,Pd(fsw,lL)) is defined on the samples of the surface it represents, for the identification of which few samples are sufficient, and which can be identified starting from the acquired experimental data limited to a neighbourhood thereof, in other words, without necessarily having to subject the inductor 2 to operational conditions which would determine a response very far from a limit operational condition of interest, of the type described above. It follows that, for the model obtained according to the method of the present invention, it is easier to obtain a good precision, since the model is valid on a portion restricted to a limited surface of the three-dimensional domain {[Veq,min, Veq,max],[fsw,min, fsw,max],[lL,min, lL,max]}.
The method 1 of the present invention can be advantageously implemented by a measuring station, indicated in
Such a measuring station 100, configured for the implementation of method 1 of the invention, includes:
The remote processing unit 200 can in turn be configured to apply at least one algorithm, optionally a Genetic Programming algorithm, to the finite set of such points P*n received by the measuring station 100 and output a description in analytical form of the locus of such points P*n, thus obtaining the locus 5 of the operational parameters which determine a response in the inductor 2 under investigation that complies with at least one working limit condition.
If desired, according to a preferred variant of the present invention, the data control and processing unit 101 can be physically distributed among several units, for example it can comprise a first data acquisition and control unit, optionally digital, indicated with the abbreviation UDAC in
According to a further variant of the measuring station 100 of the present invention, the data control and processing unit 101 can be configured to directly apply the above algorithm to the finite set of points stored in the storage 104, providing in output a description in analytical form of the locus of said points P*n thus stored, and thus obtaining the locus 5 of the operational parameters which determine a response in the inductor 2 under investigation which comply with at least one working limit condition.
Not only that, the data control and processing unit 101 is configured to handle (through a firmware stored therein) steps A to C of method 1 of the present invention, namely:
According to a preferred embodiment of the measuring station 100 of the present invention, the stimulus generation device 102 is advantageously obtained through three power converting stages operatively connected in cascade according to the so-called Opposition Method (for example as taught in F. Forest et al., “Use of opposition method in the test of high-power electronic converters,” in IEEE Transactions on Industrial Electronics, vol. 53, no. 2, pp. 530-541, April 2006), in order to subject inductor 2 under investigation to a stimulus corresponding to a triplet of the operational conditions described above.
The stimulus, as already mentioned in the introduction, comprises a zero mean square wave voltage, typical of hard-switching dc-dc PWM power converters, and a direct current, the characteristics of which, according to the present invention, are set by means of such converter stages.
More particularly, the three power converter stages comprise an Input Stage Sdl, a Test Stage SdT and an Output Stage SdU, wherein the Test Stage SdT is connected between the Input Stage Sdl and the Output stage SdU and is further configured to be connected to the inductor 2 to be investigated.
The Output Stage SdU is configured to operate in a closed loop and impose a direct current at the output of the Test Stage SdT, by adjusting its own average input current Io: the SdU, therefore, operates as a variable direct current electrical load for SdT. The output of the Output Stage SdU is therefore connected to the input of the Test Stage SdT, so that the output current of SdU is returned to the input of SdT, so as not to make it strictly necessary, but optional, the inclusion of a power dissipating element (load EL in
As can be seen from
The Input Stage Sdl is configured to work in closed loop and provide in input to the Test Stage SdT one direct voltage, through adjustment of its output voltage Vi, and the Test stage SdT is configured to work in open loop and provide for the switching frequency of the varying over time zero mean square wave voltage vL(t), through adjustment of its own frequency ft and duty cycle Dt. In other words, the Test Stage modulates the direct voltage which is imposed by the Input Stage Sdl with the direct current lL which is imposed on it by the Output stage SdU and therefore imposes to the terminals of the inductor 2 connected thereto the stimulus above described, comprising a varying over time square wave voltage vL(t) as described above and a direct current lL.
The above adjustments are made based on specific control signals emitted by the data control and processing unit 101, according to a temporal sequence determined by an algorithm (firmware) which is stored in the data control and processing unit 101 or is supplied to it, by means of another processing unit, even a remote one (for example 200 in
The detecting device 103 of the measuring station 100, according to a preferred variant of the invention represented in
Not only that, those sensors are configured to send to the data control and processing unit 101 also a measured signal of the output voltage Vi,mis of Sdl (which is one of the elements that help define the voltage vL(t) of method 1 - in particular the component VLr of vL(t) corresponds to the difference between Vi (output voltage of the Sdl) and Vo (output voltage of SdT), while VLf of vL(t) corresponds to the opposite of Vo (reference is always made to
Those signals are produced by the analog conditioning circuits incorporated in the three stages Sdl, SdT and SdU.
In this regard, the attached
As can be seen (
The Test Stage (
In the measuring station 100 of the invention, the positioning system 107 comprises:
The elastic anchoring group 1072 comprises one couple of elastically charged screws 10721, configured for being screwed on said printed circuit 105, passing through suitable openings obtained in said positioning plates 1071.
In this way, the inductor 2 is configured to be placed between the printed circuit 105 and the positioning plate 1071 and to be subjected to a pressure towards the conductive tracks 106 of the printed circuit 105, by means of the elastic anchoring group 1072 thus realizing with them an electrical contact.
On the positioning plate 1071 one of the aforementioned sensors (not shown in the drawings) is also provided, configured to detect and transmit the surface temperature Ts,mis of the inductor 2 under investigation to the data control and processing unit 101.
As can be seen from the attached figures, the surface portion 1061 of the conductive tracks 106 has a polygonal configuration, optionally a trapezoidal one, which allows electrically and easily connecting inductors 2 of various sizes to the conductive tracks 106 of the measuring station 100, without the need for welding.
The advantages of this architecture, for the measuring station 100 of the present invention, are many. First of all, it does not require a dissipative load at the output of the SdU output stage and this therefore allows its dimensions to be reduced to a minimum. It also reduces the power that the Input Stage must deliver, equal to the sum of the losses of the SdT and SdU and allows the use of a low power AC / DC mains power supply (see
Again, purely by way of non-limiting example, the aforementioned Input Stage Sdl can be designed in such a way as to have:
The Test Stage SdT can be designed to have:
The Output Stage SdU can be designed to have:
The measuring station 100 described above, according to a variant of the invention and as anticipated above, can be used in conjunction with a remote processing unit 200, in a system 1000 for the implementation of method 1 of the present invention.
In system 1000, some embodiments of which are represented purely by way of non-limiting example in
Such system 1000 is configured for carrying out the method 1 of the present invention on the power inductors 2 which operate in typical conditions of hard-switching dc-dc PWM power converters, according to a suitable sequence of operational parameters given by the invention method 1. System 1000, according to a preferred embodiment of the invention, comprises:
The architecture of system 1000 of the present invention is conceived as an open laboratory, which can be used by a single user or by a community of users and managed by an administrator. For this purpose, the software architecture can envisage three distinct types of remote processing units 200: one that operates as a System Administration (IAS), one that manages the actual execution of the measurement tests (AGE - the algorithm for determining the locus of the points of the three-dimensional mathematical space 3 that comply with the limit condition is physically executed on this unit) and a System Brokerage or Routing Unit (UBS), as illustrated in
According to a preferred embodiment of the invention, two operational modes are provided in system 1000:
According to a further variant of the system 1000 of the present invention, this can work according to a:
In view of the aforementioned description and the examples provided, the advantages offered by the method 1, the measuring station 100 and system 1000 of the present invention are clear.
In fact, thanks to method 1, the measuring station 100 and the system 1000 described above, it is possible to determine, quickly, with less computational effort than traditional methods, in a safe and reliable way, a function that describes in analytical form the locus 5 of the operational parameters (points P*n) of an electrical or electronic power component that respect the working limit condition, for example a value of the peak-to-peak current ripple ΔiLpp=ΔiLpp,ref, and/or its surface temperature Ts=Ts,ref, and/or its dissipated power Pd=Pd,ref.
The locus 5 of those operational parameters can be used by an operator, through the measuring station 100 or the system 1000 described above, to determine (for example visually by means of a display which can be operatively connected to such measuring station 100 or to a component of such system 1000) if starting from a triplet of operational parameters of the inductor under investigation, such inductor 2 complies or not with a limit condition. For example, according to a variant of the present invention, one or more mobile devices can be connected to the measuring station or system 1000, for example through the wifi, for displaying the results obtained by the method described above.
In the foregoing the preferred embodiments were described, and some modifications of the present invention were suggested, but it should be understood that those skilled in the art can make modifications and changes without departing from the relative scope of protection, as defined by the appended claims.
Thus, for example, as anticipated in the introduction, it is specified that although the present invention has been described with particular reference to power inductors, it is quite clear that it can also be applied to other electronic power components such as power capacitors, taking into account that where in this text reference is made to the voltage one should instead refer to the current and vice versa, and that where in this text reference is made to the inductance L one should instead refer to the capacitance C. In this case, the stimulus applied to the capacitor during step C.2 of method 1 would comprise a given square wave current ic(t), varying over time, and an average voltage VC,av and the response measured at step C.3 would comprise a varying over time voltage vc(t), detected at the terminals of the same capacitor. Similarly, based on the measured response, corresponding response parameters at step C.4 would be determined comprised between the voltage ripple ΔvCpp, the surface temperature Ts of the capacitor, and the average electrical power Pd dissipated by the capacitor.
Number | Date | Country | Kind |
---|---|---|---|
102020000026599 | Nov 2020 | IT | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/IB2021/060270 | 11/5/2021 | WO |