Disclosure of Invention
It is an object of the present invention to address the technical deficiencies of the prior art by providing a method for calculating the average power capacity of a dielectric integrated suspension line for guiding SISL a transmission line design based on power considerations.
The technical scheme adopted for realizing the purpose of the invention is as follows:
a method of calculating an average power capacity of a media integrated suspension line, comprising:
s1, obtaining an attenuation constant of a SISL transmission line through simulation, and calculating a first expression of heat flow phi gen generated in the SISL transmission line, wherein phi gen = -delta P, delta P represents the power variation of the length delta z of the SISL transmission line, delta P= -2 alpha P delta z, P=P 0e-2αz, alpha is a power flow function in the transmission direction (z direction) of the consumable transmission line through simulation of a full-wave simulator, P 0 represents input power, and the corresponding coordinate z=0;
s2, calculating and SISL a relation between the outer surface temperature T s of the transmission line and the heat flow phi gen, and obtaining a second expression of the heat flow phi gen:
Wherein n and m are the number of surface layers associated with convection and radiation boundaries, respectively, h c,i is the heat transfer coefficient, A i is the surface area associated with convection, A j is the surface area associated with radiation, the surface area A i and A j over the length of Δz on the SISL transmission line are calculated as the product of the width or height of the SISL transmission line and Δz, wherein the area of the upper and lower horizontal surfaces is the product of the width and Δz, the area of the left and right vertical surfaces is the product of the height and Δz, σ is the blackbody radiation constant, ζ is the emissivity of the surface layer associated with radiation, T ∞ is the ambient temperature, where AndThe upper scale of (2) is the power of 4;
S3, calculating the relation between the highest temperature T w in the SISL transmission line and the heat flow phi gen to obtain a third expression of the heat flow phi gen
Wherein Φ cond denotes that K air is the thermal conductivity of air, W e Δz is the area of thermal conduction, W e denotes that the width of the SISL transmission line equivalent to the width of the parallel plate thermal model after being widened by the influence of the substrate supporting the signal line, h=b/2, b is the total height of the upper and lower cavities,Representing total radiation heat exchange quantity in the upper cavity and the lower cavity, wherein ζ c is emissivity of a signal wire, ζ d is emissivity of a medium, ζ in is emissivity of materials on the upper surface of the bottom PCB and the lower surface of the top PCB, F 12 is view angle factor from a circuit layer to a corresponding cover plate, w air represents cavity width in the transmission wire, w is width of the signal wire, whereinAndThe upper scale of (2) is the power of 4;
s4, calculating SISL the average power capacity of the transmission line
Two equations related to the input power P 0 and the external surface temperature T s of the SISL transmission line are obtained through simultaneous equations of the first expression, the second expression and the third expression of the heat flow phi gen, the two obtained equations form an equation set, the corresponding input power P 0 when the highest temperature inside SISL reaches the glass-state softening temperature T g of the medium is obtained through solving the equation set, and the input power P 0 is the average power capacity of the SISL transmission line.
The method of the invention realizes the calculation of the corresponding relation between SISL surface temperature and input power through heat transfer analysis, establishes SISL parallel plate thermal model for internal heat transfer analysis through analogy of heat flow field and electrostatic field, realizes the calculation of the corresponding relation between the internal highest temperature and input power of SISL under the condition of composite medium heat transfer through numerical modeling and heat transfer analysis, and can realize the rapid calculation of the average power capacity of SISL transmission lines.
The invention provides a method for rapidly calculating SISL average power capacity of a transmission line, which can rapidly estimate the maximum continuous wave power which can be born when the SISL transmission line is adopted for circuit design and has important reference significance when a SISL platform is utilized for designing a high-power circuit.
Detailed Description
The invention is described in further detail below with reference to the drawings and the specific examples. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
As shown in fig. 1, in the method for calculating SISL the average power capacity of the transmission line according to the embodiment of the present invention, the average power capacity of the SISL transmission line can be quickly calculated. The method comprises the following specific steps:
Step one, obtaining the attenuation constant of SISL transmission lines through simulation, and calculating the heat flow generated in SISL transmission lines.
The attenuation constant alpha of the uniform SISL transmission line is obtained through simulation of a full-wave simulator, and when the input power is P 0, the power of the SISL transmission line at the position with the length z is reduced as follows:
P=P0e-2αz (1)
The amount of power change over the length Δz of SISL transmission line is:
ΔP=-2αPΔz (2)
the heat flow Φ gen generated in SISL transmission lines is therefore:
Φgen=-ΔP (3)
And step two, calculating SISL the relation between the outer surface temperature T s of the transmission line and the heat flow phi gen.
Heat dissipation of SISL transmission lines includes convection and radiation, and the self-packaging of SISL transmission lines is approximately isothermal due to the presence of copper cladding and a large number of metallized vias.
The heat exchange coefficient h c of the convective heat exchange has the following calculation formula:
1) Coefficient of heat transfer with hot side facing horizontally up
2) Heat transfer coefficient with heat face down
3) Heat transfer coefficient with vertical hot face
Wherein, T ∞ is the ambient temperature, w box is the width of the SISL self-packaging, and h box is the height of the SISL self-packaging. The temperature T s of the outer surface of SISL transmission line is related to the heat flow Φ gen as follows:
wherein n and m are the number of surface layers related to convection and radiation boundaries, h c,i is a heat exchange coefficient, a i is a surface layer area related to convection, a j is a surface layer area related to radiation, the surface layer areas a i and a j over the Δz length on the SISL transmission line are calculated as the product of the width or height of the SISL transmission line and Δz, wherein the area of the upper and lower horizontal surfaces is the product of the width and Δz, the area of the left and right vertical surfaces is the product of the height and Δz, σ is a blackbody radiation constant, ζ is the emissivity of the surface layer related to radiation, T ∞ is the ambient temperature, and the formulae T s and T ∞ are marked to the power of 4;
And thirdly, calculating SISL the relation between the internal highest temperature T w of the transmission line and the heat flow phi gen.
SISL are mainly located in the signal wire and the supporting medium around the signal wire, and the heat converted by the partial loss is transferred to the outer surface through the composite medium formed by the upper cavity, the lower cavity and the supporting medium of the signal wire. Typically SISL the cavity height is insufficient to produce convective heat transfer, and the internal heat transfer is only by conduction and radiation, both paths being in parallel.
When SISL is filled with air, the thermal flow field can be analogically compared with the electrostatic field, and in order to accurately calculate the area of heat conduction outwards from the signal line and the supporting medium around the signal line, SISL is equivalent to a parallel plate thermal model according to the principle that the thermal flow field is analogically compared with the electrostatic field, and the thermal flow field width calculation formula is as follows:
Wherein h=b/2, b is the total height of the upper and lower cavities, Is the characteristic impedance of the air-filled SISL. When SISL is configured to meet the stripline range,Can be approximated by the following equation
1) When w/(b-t) is not more than 0.35, t/b is not more than 0.25, and d ' '/d ' is not more than 0.11:
2) When w/(b-t) is more than or equal to 0.35 and t/b is less than or equal to 0.25:
where t is the signal line thickness, d″ is the small one of w and t, and d' is the large one of w and t.
Since the thermal conductivity of the substrate under the signal line is much higher than that of air, the calculated width W e of the parallel plate thermal model is widened to two sides under the condition of only air filling, and the widened width is related to the cavity height h air, the cavity width W air, the thermal conductivity K sub of the substrate supporting the signal line, the line width W of the signal line and the thickness h sub of the substrate supporting the signal line, the invention provides fitting the multiple lambda of the widened width W e of the parallel plate thermal model of the air filling SISL through numerical modeling, and the width of the parallel plate model influenced by the substrate supporting the signal line is W e=λ·we.
A large amount of data of the corresponding relation between the input power of SISL and the internal highest temperature is simulated as a sample through thermal simulation software, and a mathematical model is adopted to fit a widened and widened multiple lambda of the width w e of the parallel plate thermal model of the air filling SISL transmission line, wherein lambda is related to the cavity height h air, the cavity width w air, the substrate thermal conductivity K sub for supporting the signal line, the line width w of the signal line and the substrate thickness h sub for supporting the signal line.
Wherein the heat transfer heat can be calculated by the following formula:
Where K air is the thermal conductivity of air and W e Δz is the area of thermal conduction.
The total radiant heat exchange in the upper and lower cavities can be approximated by the following equation:
Xi c is the emissivity of the signal wire, xi d is the emissivity of the medium, xi in is the emissivity of the material on the upper surface of the bottom PCB and the material on the lower surface of the top PCB, F 12 is the visual angle factor from the circuit layer to the corresponding cover plate, and the calculation formula is that
Finally, the energy conservation includes:
And step four, calculating SISL average power capacity.
Substituting equations (3) - (6) into equation (7) yields an equation for input power and surface temperature:
Substituting the formulas (3) and (8) - (15) into the formula (16) to obtain a second equation about the input power and the surface temperature, wherein the obtained two equations form an equation set, and the corresponding input power when the highest temperature in SISL reaches the glass-state softening temperature T g of the medium can be obtained by solving the equation set, and the input power is the average power capacity of SISL.
The average power calculation method of the transmission line can be realized by a computer program realizing processing model, corresponding parameters are input into the model, then corresponding results are output by the model, the change of the transmitted temperature is realized by changing the input power, and when the highest temperature in SISL reaches the glass-state softening temperature T g of the medium, the corresponding input power, namely the average power capacity of SISL. Therefore, when the SISL transmission line is adopted for circuit design, the maximum continuous wave power which can be born by the circuit can be rapidly estimated, so that the instruction on circuit design is realized, the designed circuit meets the corresponding power requirement, and the circuit has important reference significance when the SISL platform is used for designing a high-power circuit.
In the following, a section of SISL transmission line operating at 3GHz is taken as an example, and its maximum transmission power is calculated. The S parameters are shown in fig. 5. The plates constituting SISL transmission lines all use FR4, the relative dielectric constant epsilon r =4.4, and tan delta=0.02. The glass transition temperature of conventional FR4 is 140 ℃, and thus the ultimate operating temperature of SISL transmission lines is 140 ℃. The thickness of the five layers of plates is 0.3mm, 0.6mm, 0.127mm, 0.6mm and 0.3mm, and the copper thickness is 0.035mm. The working environment is room temperature T ∞ =22℃.
SISL the transmission line was packaged from width W box =20 mm, height h box = 1.997mm, cavity width W air =10 mm, cavity height h air =0.6 mm, b=1.397 mm, signal line width w=1.75 mm, thermal conductivities of fr4 and copper are 0.294W/(m· ℃ and 400W/(m· ℃ respectively), emissivity is 0.9 and 0.1 respectively.
Step one, obtaining the attenuation constant of SISL transmission lines through simulation.
As shown in FIG. 6, the 3GHz attenuation constant is 0.186N p/m, and the heat flow rate generated by loss in the transmission line is calculated by substituting the values (2) and (3)
Φgen=2αPΔz=2×0.186×P0×Δz (17)
And step two, calculating SISL the relation between the outer surface temperature T s of the transmission line and the heat flow phi gen.
The ambient temperature T ∞ is 295K (22 ℃), and substituting data of formulae (4) - (6) and dimensions, materials, etc. into (7) yields:
And thirdly, calculating SISL the relation between the internal highest temperature T w of the transmission line and the heat flow phi gen.
Because SISL in the embodiment is sized to match the stripline, the characteristic impedance can be calculated as follows from equations (11), (12)
Thus, the width of the parallel plate thermal model is
Using a single formula as a model function of lambda, one expression of lambda fitted by using a plurality of groups of data of the corresponding relation between input power and maximum temperature as a sample is
The application range of the expression is that the line width W of the signal line is 0.1mm-50Ω, the cavity height h air is 0.3mm-2mm, the cavity width W air:max{4w,w+4hair } -30mm, the thickness h sub of the signal line supporting medium is 0.127mm, the thermal conductivity of the signal line supporting medium is K sub, and the range is 0.2W/(m·DEG C) -0.3W/(m·DEG C)
From this, W can be calculated e
Substituting the glass-state softening temperature T g of the medium into T w, wherein the heat flow of heat conduction is
The viewing angle factor of the heat exchange of radiation in SISL transmission lines can be calculated as
The total radiation heat exchange quantity in the upper cavity and the lower cavity is
Finally, the formulas (23) and (25) are substituted into the formula (16) to obtain
And step four, calculating SISL the average power capacity of the transmission line.
Substituting formula (17) into formula (18) to obtain
Substituting formula (17) into formula (26) to obtain
By solving the system of equations consisting of equations (27) and (28), P 0=58.96W,Ts = 76.48 ℃,58.96W is the average power capacity of the SISL transmission line in the example.
In order to verify the accuracy of the calculation result, thermal simulation software is adopted to verify the calculation result.
Modeling simulation is carried out in a full-wave simulator, the simulation model and the result are imported into thermal simulation software, heat dissipation boundary conditions are set to be convection and radiation, and the convection heat transfer coefficients of the upper surface, the lower surface and the side surface are respectively obtained by formulas (1) - (7).
The maximum temperature of SISL is 140 ℃ by varying the input power P0. As shown in fig. 7, when SISL reaches the glass softening temperature, the input power is 64.10W, the corresponding surface temperature is 77.78 ℃, and the loss data transfer accuracy and meshing of the full wave simulator to the thermal simulation software are part of the reasons for temperature and power deviation, which are very close to the calculation results in the examples, thus proving the accuracy of the calculation method of the invention.
While the fundamental and principal features of the invention and advantages of the invention have been shown and described, it will be apparent to those skilled in the art that the invention is not limited to the details of the foregoing exemplary embodiments, but may be embodied in other specific forms without departing from the spirit or essential characteristics thereof;
The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein.
Furthermore, it should be understood that although the present disclosure describes embodiments, not every embodiment is provided with a separate embodiment, and that this description is provided for clarity only, and that the disclosure is not limited to the embodiments described in detail below, and that the embodiments described in the examples may be combined as appropriate to form other embodiments that will be apparent to those skilled in the art.