CN107656136B - Solar wing output power numerical calculation method under shielding working condition - Google Patents
Solar wing output power numerical calculation method under shielding working condition Download PDFInfo
- Publication number
- CN107656136B CN107656136B CN201710652069.6A CN201710652069A CN107656136B CN 107656136 B CN107656136 B CN 107656136B CN 201710652069 A CN201710652069 A CN 201710652069A CN 107656136 B CN107656136 B CN 107656136B
- Authority
- CN
- China
- Prior art keywords
- battery
- shielding
- calculating
- voltage
- current
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R21/00—Arrangements for measuring electric power or power factor
- G01R21/06—Arrangements for measuring electric power or power factor by measuring current and voltage
Abstract
The invention provides a solar wing output power numerical calculation method under a shielding working condition, which is used for meeting the requirements of various types of spacecrafts on rapid and accurate analysis and calculation of solar wing output power at various stages of development. (1) Preprocessing a cell sheet distribution image of the solar cell array, and storing the preprocessed cell sheet distribution image as a universal text format; (2) reading the universal text format file generated in the step (1) and establishing a geometric relation model of the solar cell array; (3) reading shielding pixel information in a video or picture format, completing the correspondence of the pixel information and the geometric relation model established in the step (2), completing the setting of a shielding factor K of all the cells in the model, and recording an illumination angle theta of a corresponding working condition, namely an included angle between a solar ray and a cell normal line, which ranges from 0 degree to 90 degrees; (4) inputting output characteristic parameters of the battery piece during positive illumination under the working temperature condition; (5) and (4) calculating the total output current of the solar wing according to the output current of each battery string obtained in the step (4), and calculating the total output power of the solar wing by multiplying the total output current by the bus voltage.
Description
Technical Field
The invention relates to a solar wing output power numerical calculation method under a shielding working condition, and belongs to the technical field of spacecrafts.
Background
When the spacecraft runs in orbit, the spacecraft body can shield the solar wing along with the change of the attitude or the illumination angle, so that the output power of the solar wing is reduced. In the design stage and the in-orbit operation process of the spacecraft, the calculation of the output power of the solar wing under the shielding working condition is carried out, and the method plays an important role in confirming the energy balance of the spacecraft and evaluating the task completion capability of the spacecraft.
Under the shielding working condition, the output power of the solar wing is not only related to the area and the position of the shadow, but also related to the output characteristics of the solar wing, the series-parallel connection design, the sheet distribution mode, the illumination angle and other factors.
The traditional output power calculation work usually adopts a manual method, the calculation precision and the calculation efficiency are very low through simple conversion of shadow area or the number of affected matrixes, only static working conditions can be usually calculated, and the method cannot be suitable for the condition that the shadow shape and the shadow depth are very complicated.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the method solves the problems that the existing manual calculation method is low in calculation efficiency and precision and cannot be suitable for the working conditions with complicated shadow shapes and shadow depths, provides a solar wing output power numerical calculation method under the shielding working conditions, realizes the calculation of the solar wing output power under the complicated dynamic shielding working conditions, and meets the requirements of various stages of spacecraft development on the rapid and accurate analysis and calculation of the solar wing output power.
The technical solution of the invention is as follows:
a solar wing output power numerical value calculation method under a shielding working condition comprises the following steps:
(1) preprocessing a cell sheet distribution diagram of the solar cell array, storing the preprocessed cell sheet distribution diagram as a universal text format, and then entering the step (2);
(2) reading the universal text format file generated in the step (1), and establishing a geometric relation model of the solar cell array, wherein the geometric relation model comprises the following information: the position information of all the battery pieces, the battery pieces and the battery strings are in corresponding relation; the corresponding relation between the battery string and the array; then entering the step (3);
(3) reading shielding pixel information in a video or picture format, completing the correspondence of the pixel information and the geometric relation model established in the step (2), completing the setting of a shielding factor K of all the cells in the model, and recording an illumination angle theta of a corresponding working condition, namely an included angle between a solar ray and a cell normal line, which ranges from 0 degree to 90 degrees; in the calculation, the voltage characteristic of the battery piece is assumed to be kept unchanged when the shielding and the illumination angle are changed; then entering the step (4);
(4) inputting output characteristic parameters of the battery piece in the positive illumination under the working temperature condition: including an open circuit voltage VocShort-circuit current IscMaximum operating point voltage VmpMaximum operating point current ImpBypass diode conduction voltage drop Vd(ii) a According to the working voltage parameters of the battery strings, the output current I of each battery string is calculated in an iterative wayout(ii) a Then entering step (5);
(5) and (4) calculating the total output current of the solar wing according to the output current of each battery string obtained in the step (4), and calculating the total output power of the solar wing by multiplying the total output current by the bus voltage.
Furthermore, the position information of the battery piece adopts the coordinates (X) of the lower left corner pointL，YL) And coordinates of the upper right corner point (X)R，YR) To describe.
Furthermore, the corresponding relationship between the battery piece and the battery string is distinguished by adopting color information.
Further, the corresponding relation between the battery string and the subarray is distinguished by adopting layer information.
Furthermore, the range of the shielding factor K is 0-1, which represents the influence degree of shielding on the current parameter in the output characteristic.
Compared with the prior art, the invention has the advantages that: the method of the invention considers various factors related to output power calculation under the shielding working condition, constructs an output power calculation model by utilizing the solar wing cell piece distribution diagram and the cell piece characteristic parameters, and automatically completes the calculation of the output power of the solar wing by inputting the shadow information of the image or video format. By the method, the problem of dynamic and rapid calculation of the output power of the spacecraft with the complex configuration when the solar wing is shielded can be solved, and calculation support can be provided for task analysis, solar wing design, on-orbit working planning and the like of the spacecraft.
1) The calculation model is directly generated after the design cloth sheet diagram of the solar wing is preprocessed, so that the modeling is simple and rapid, and the combination with the spacecraft design service is tight;
2) the output power calculation under the condition of mixed occlusion of different depths of the shadow is realized by setting occlusion factors in the calculation model and using the shadow depths distinguished by different colors in the occlusion information to set the occlusion factors;
3) compared with the traditional manual calculation mode, each cell and each shadow pixel of the solar wing are analyzed by the method, the output current calculation result is accurate to be below 1mA, and the calculation result precision is high.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a schematic diagram of calculation of equivalent shading factors of the battery plate.
Detailed Description
The invention will be further described with reference to the accompanying drawings.
As shown in FIG. 1, the method for calculating the output power of the solar wing under the shading working condition comprises the following steps:
(1) preprocessing a cell sheet distribution diagram of the solar cell array, storing the cell sheet distribution diagram in a universal text format, and then entering the step (2);
(2) reading the text file generated in the step (1) by adopting a program, establishing a geometric relation model of the solar cell array, wherein the calculation model comprises the following information: position information of all cells (using coordinates of lower left corner (X)L，YL) And coordinates of the upper right corner point (X)R，YR) Description); the corresponding relation between the battery piece and the battery string (distinguished by adopting color information); and (4) corresponding relation between the battery string and the array (distinguishing by adopting layer information).Then entering the step (3);
(3) reading shielding pixel information in a video or picture format, completing the correspondence of the pixel information and a geometric model, completing the setting of shielding factors K of all the cells in the model, and recording the illumination angle theta (the included angle between the sunlight and the normal line of the cell, the range is 0-90 degrees) of the corresponding working condition; the range of the shielding factor is 0-1, representing the influence degree of shielding on current parameters in the output characteristic, and assuming that the voltage characteristic of the battery piece is kept unchanged during the shielding and illumination angle variation in the calculation; then entering the step (4);
(4) inputting output characteristic parameters of the battery piece in the positive illumination under the working temperature condition: including an open circuit voltage VocShort-circuit current IscMaximum operating point voltage VmpMaximum operating point current ImpBypass diode conduction voltage drop Vd(ii) a According to the working voltage parameters (including bus voltage and line voltage drop) of the battery strings, the output current I of each battery string is calculated in an iterative mannerout(ii) a Then entering step (5);
(5) and calculating the total output current of the solar wing according to the output currents of all the battery strings, and calculating the total output power of the solar wing by multiplying the total output current by the bus voltage.
The solar cell slice distribution drawing is a drawing which is drawn by CAD software (such as AutoCAD software) and reflects the size and the slice distribution position information of the solar cell slice; preprocessing a sheet layout, namely distinguishing an array and a battery string in the array by adopting different layers and color information in the sheet layout, and then exporting the sheet layout into a file with a text format (such as a DXF format of AutoCAD);
the read shielding information adopts different colors to represent whether shielding exists or not and the shielding shadow depth; according to the position information of the battery pieces in the calculation model, after conversion and rounding completion are carried out according to the proportional relation between the cloth piece image and the shadow image, corresponding shadow information is retrieved according to pixel positions, and the shielding factors K of all the battery pieces are calculated and set;
the occlusion information read in step (3) represents whether occlusion exists or not and the shadow depth of the occlusion by using different colors (for example, the occlusion information can be obtained by using different colorsWhite is adopted to represent no shielding, black is adopted to represent full shielding, and gray level color is adopted to represent half shielding with different transmittances); the principle of calculating the equivalent shading factor of each cell is shown in fig. 2, wherein the thick line in the figure represents the actual cell position, and the small square represents the pixel in the shading map. The first step is that according to the position information of the battery piece in the calculation model, after conversion and rounding completion are carried out according to the proportional relation and the coordinate conversion relation between the piece distribution graph and the shadow graph, a pixel set of the whole pixel in the position range of the battery piece is obtained, and the pixel set is totally m rows and n columns as shown in the figure. Reading the color of the pixel, and obtaining the shielding factor k of each pixel according to the corresponding table of the color and the shading factorij(i.e., the corresponding shading factor for the pixel in row i and column j), the range of values is 0-1. Calculating the shielding equivalent factor of the battery piece (the calculation result range is 0-1):
the calculation process of the step (4) adopts current parameters to iterate, firstly, the output voltages of all the battery pieces in the battery string are calculated through current, and the calculation formula is as follows:
c2＝(Vmp/Voc-1)/ln(1-Imp/Isc) (2)
where I is the output current for iterative calculation, VcellOther variables for the calculated cell output voltage are described above. It should be noted that the calculation of the output voltage of the cell using the current in equation (4) coexists in two cases: a) calculating a current value of I or lessscCalculating output voltage by adopting a current-voltage curve of the battery piece when the Kcos theta is reached; b) calculating a current value greater than IscWhen Kcos theta, the battery plate worksIn the bypass diode conduction mode, the output voltage is the conduction voltage drop (negative value is adopted) of the bypass diode.
And then summing the voltages of all the battery plates of the battery string to obtain a corresponding calculated voltage Vc of the battery of the string, and comparing the calculated voltage Vc with the voltage of the working point of the parallel battery string (namely the bus voltage + the line voltage drop). When Vc is smaller than the voltage of the working point, the result proves that I is larger than the real output current, and I is reduced and recalculated; when Vc is larger than the voltage of the working point, the result proves that I is smaller than the real output current, and I is adjusted to be larger for recalculation. Since the value range of I is [0, Isccosθ]Therefore, the dichotomy can be adopted to continuously iterate the approximation result in the interval (note: because the calculation formula of the battery voltage has very obvious nonlinear characteristics, the calculation formula cannot be normally converged by using the Newton iteration method, so that the calculation is not recommended) until the Ik+1-IkIf | < Δ (where Δ is a convergence threshold of the current, and may be selected according to the calculation accuracy requirement), I may be regarded ask+1I is the battery stringout。
Those skilled in the art will appreciate that those matters not described in detail in the present specification are well known in the art.
Claims (8)
1. A solar wing output power numerical value calculation method under a shielding working condition is characterized by comprising the following steps:
(1) preprocessing a cell sheet distribution diagram of the solar cell array, storing the preprocessed cell sheet distribution diagram as a universal text format, and then entering the step (2);
(2) reading the universal text format file generated in the step (1), and establishing a geometric relation model of the solar cell array, wherein the geometric relation model comprises the following information: the position information of all the battery pieces, the battery pieces and the battery strings are in corresponding relation; the corresponding relation between the battery string and the array; then entering the step (3);
(3) reading shielding pixel information in a video or picture format, completing the correspondence of the shielding pixel information and the geometric relation model established in the step (2), completing the setting of shielding factors K of all the battery pieces in the model, and recording an illumination angle theta of corresponding working conditions, namely an included angle between a solar ray and a battery piece normal line, wherein the range is 0-90 degrees; in the calculation, the voltage characteristic of the battery piece is assumed to be kept unchanged when the shielding and the illumination angle are changed; then entering the step (4);
(4) inputting output characteristic parameters of the battery piece in the positive illumination under the working temperature condition: including an open circuit voltage VocShort-circuit current IscMaximum operating point voltage VmpMaximum operating point current ImpBypass diode conduction voltage drop Vd(ii) a Iteratively calculating the output current I of each battery string according to the working voltage parameters of the battery stringsoutThe working voltage parameters of the battery string comprise bus voltage and line voltage drop; then entering step (5);
(5) and (4) calculating the total output current of the solar wing according to the output current of each battery string obtained in the step (4), and calculating the total output power of the solar wing by multiplying the total output current by the bus voltage.
2. The method according to claim 1, wherein the position information of the solar cell is the coordinates of the lower left corner (X)L，YL) And coordinates of the upper right corner point (X)R，YR) To describe.
3. The method for calculating the solar wing output power value under the shielding working condition according to claim 1, further comprising distinguishing the correspondence between the battery pieces and the battery strings by using color information.
4. The method for calculating the solar wing output power value under the shading working condition according to claim 1, further comprising the step of distinguishing the corresponding relation between the battery string and the partial array by using layer information.
5. The method for calculating the output power value of the solar wing under the shielding working condition according to any one of claims 1 to 4, wherein the cell piece distribution diagram is drawn by CAD software and reflects the size and the distribution position information of the cell piece.
6. The method for calculating the output power value of the solar wing under the shielding working condition according to any one of claims 1 to 4, wherein the read shielding pixel information adopts different colors to represent whether shielding exists and the shielding depth; and according to the position information of the battery pieces in the geometric relation model, converting according to the proportional relation between the sheet distribution graph and the shadow graph, completing the rounding, searching corresponding shadow information according to pixel positions, and calculating and setting the shielding factors K of all the battery pieces.
7. The method for calculating the output power value of the solar wing under the shielding working condition according to any one of claims 1 to 4, wherein current parameters are adopted to iteratively calculate the output current of each battery string, specifically: calculating the output voltage of all the battery pieces in the battery string through current, and summing the working voltages of the battery pieces to obtain the output voltage of the battery string; and when the difference between the output voltage calculation result and the working point voltage of the battery string, namely the bus voltage + the line voltage drop, is less than the iterative convergence threshold value delta V, the calculation is considered to be converged, and the current value at the moment is the output current of the battery string.
8. The solar wing output power numerical calculation method under the shielding working condition according to claim 7, wherein the output voltage of the cell calculated by adopting the current coexists in two conditions:
a) calculating a current value of I or lessscxKxcos θ, wherein IscCalculating the output voltage by adopting a current-voltage curve of the battery piece, wherein K is a shielding factor of the battery piece, and theta is an illumination angle of a corresponding working condition;
b) calculating a current value greater than IscAnd multiplying the multiplied by the cos theta, the battery piece works in a bypass diode conduction mode, and the output voltage is the conduction voltage drop of the bypass diode.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710652069.6A CN107656136B (en) | 2017-08-02 | 2017-08-02 | Solar wing output power numerical calculation method under shielding working condition |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710652069.6A CN107656136B (en) | 2017-08-02 | 2017-08-02 | Solar wing output power numerical calculation method under shielding working condition |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107656136A CN107656136A (en) | 2018-02-02 |
| CN107656136B true CN107656136B (en) | 2020-04-07 |
Family
ID=61128255
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710652069.6A Active CN107656136B (en) | 2017-08-02 | 2017-08-02 | Solar wing output power numerical calculation method under shielding working condition |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107656136B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3270303B2 (en) * | 1995-07-26 | 2002-04-02 | キヤノン株式会社 | Battery power supply device characteristic measuring device and measuring method |
| CN102684555A (en) * | 2011-03-17 | 2012-09-19 | 南京优阳光伏技术有限公司 | Component-level power optimization technology applied to photovoltaic power generation system |
| CN102928714A (en) * | 2012-11-02 | 2013-02-13 | 北京航空航天大学 | Moonlet sun array life forecast method based on I-V curve and energy balance |
| CN103353716A (en) * | 2013-05-16 | 2013-10-16 | 河海大学 | Topological structure capable of realizing photovoltaic array reconstruction and method for solving maximum power point |
| CN103944510A (en) * | 2014-05-06 | 2014-07-23 | 河海大学常州校区 | Photovoltaic module output characteristic judgment method |
| CN106324631A (en) * | 2016-07-28 | 2017-01-11 | 北京空间飞行器总体设计部 | Remote sensing satellite energy balance constraint analysis system and method |
| CN106575872A (en) * | 2014-05-28 | 2017-04-19 | 库珀技术公司 | Distributed low voltage power systems |
-
2017
- 2017-08-02 CN CN201710652069.6A patent/CN107656136B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3270303B2 (en) * | 1995-07-26 | 2002-04-02 | キヤノン株式会社 | Battery power supply device characteristic measuring device and measuring method |
| CN102684555A (en) * | 2011-03-17 | 2012-09-19 | 南京优阳光伏技术有限公司 | Component-level power optimization technology applied to photovoltaic power generation system |
| CN102928714A (en) * | 2012-11-02 | 2013-02-13 | 北京航空航天大学 | Moonlet sun array life forecast method based on I-V curve and energy balance |
| CN103353716A (en) * | 2013-05-16 | 2013-10-16 | 河海大学 | Topological structure capable of realizing photovoltaic array reconstruction and method for solving maximum power point |
| CN103944510A (en) * | 2014-05-06 | 2014-07-23 | 河海大学常州校区 | Photovoltaic module output characteristic judgment method |
| CN106575872A (en) * | 2014-05-28 | 2017-04-19 | 库珀技术公司 | Distributed low voltage power systems |
| CN106324631A (en) * | 2016-07-28 | 2017-01-11 | 北京空间飞行器总体设计部 | Remote sensing satellite energy balance constraint analysis system and method |
Non-Patent Citations (1)
| Title |
|---|
| 局部阴影遮挡的太阳电池组件输出特性实验研究;张臻 等;《太阳能学报》;20120131;第33卷(第1期);第5-12页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107656136A (en) | 2018-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10331805B2 (en) | Methods and systems for determining obstruction setbacks in a photovoltaic system | |
| CN107742171B (en) | Photovoltaic power station power generation power prediction method based on mobile shadow image recognition | |
| CN104657563A (en) | Intelligent design platform and intelligent design method of photovoltaic power station | |
| Lyden et al. | Modelling, parameter estimation and assessment of partial shading conditions of photovoltaic modules | |
| Zhang et al. | Modelling and estimating performance for PV module under varying operating conditions independent of reference condition | |
| Geisz et al. | Characterization of multiterminal tandem photovoltaic devices and their subcell coupling | |
| Carneiro et al. | Assessment of solar irradiance on the urban fabric for the production of renewable energy using LIDAR data and image processing techniques | |
| Hu et al. | A methodology for calculating photovoltaic field output and effect of solar tracking strategy | |
| CN111259550B (en) | Grid search and improved NM simplex algorithm-based photovoltaic model updating method | |
| Lappalainen et al. | Experimental comparison between various fitting approaches based on RMSE minimization for photovoltaic module parametric identification | |
| CN108008633A (en) | Irradiation level comprising a variety of Changes in weather and photovoltaic module coordinate incidence relation method for building up | |
| Elhammoudy et al. | Dandelion Optimizer algorithm-based method for accurate photovoltaic model parameter identification | |
| CN107656136B (en) | Solar wing output power numerical calculation method under shielding working condition | |
| CN110782024A (en) | Photo-thermal electric field solar direct normal radiation prediction method based on convolutional neural network | |
| CN113437938B (en) | Photovoltaic array output power calculation method considering field characteristics due to regional differences | |
| CN112561181A (en) | Photovoltaic power generation prediction system based on Unet network and foundation cloud picture | |
| Yang et al. | Fast simulation modeling and multiple-PS fault diagnosis of the PV array based on I–V curve conversion | |
| Li et al. | Development of spacecraft solar array electrical performance simulation system | |
| JP2015007959A (en) | Solar battery design apparatus, solar battery design method, and solar battery manufacturing method | |
| Passow et al. | PlantPredict: solar performance modeling made simple | |
| Osorio et al. | A MATLAB/GUI for Photovoltaic Modules Performance Simulations Based on Two-Diode Model | |
| Stewart | Solar Array Analysis Suite | |
| Looney et al. | Machine Learning-based Classification of Spectral Conditions for High-Throughput Indoor Testing of Photovoltaic Modules | |
| CN107340799A (en) | A kind of maximum power point of photovoltaic array tracking based on shadow Detection | |
| KR101405354B1 (en) | Method for PV cell 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 |