CN111009729B - High-density integrated active phased array T/R assembly arrangement method based on machine, electricity and heat - Google Patents

High-density integrated active phased array T/R assembly arrangement method based on machine, electricity and heat Download PDF

Info

Publication number
CN111009729B
CN111009729B CN201911096411.4A CN201911096411A CN111009729B CN 111009729 B CN111009729 B CN 111009729B CN 201911096411 A CN201911096411 A CN 201911096411A CN 111009729 B CN111009729 B CN 111009729B
Authority
CN
China
Prior art keywords
attenuation
region
sub
amplitude
components
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
Application number
CN201911096411.4A
Other languages
Chinese (zh)
Other versions
CN111009729A (en
Inventor
杨京
杨慧杰
范乃康
吕鹏
张建波
李岩
张凯
黄超
马俊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Institute of Space Radio Technology
Original Assignee
Xian Institute of Space Radio Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian Institute of Space Radio Technology filed Critical Xian Institute of Space Radio Technology
Priority to CN201911096411.4A priority Critical patent/CN111009729B/en
Publication of CN111009729A publication Critical patent/CN111009729A/en
Application granted granted Critical
Publication of CN111009729B publication Critical patent/CN111009729B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01QANTENNAS, i.e. RADIO AERIALS
    • H01Q1/00Details of, or arrangements associated with, antennas
    • H01Q1/50Structural association of antennas with earthing switches, lead-in devices or lightning protectors
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01QANTENNAS, i.e. RADIO AERIALS
    • H01Q21/00Antenna arrays or systems
    • H01Q21/0006Particular feeding systems
    • H01Q21/0025Modular arrays
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01QANTENNAS, i.e. RADIO AERIALS
    • H01Q23/00Antennas with active circuits or circuit elements integrated within them or attached to them

Landscapes

  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

The invention provides a high-density integrated active phased array T/R assembly arrangement method based on machinery, electricity and heat, which comprises the following implementation steps: 1. dividing the T/R component arrangement area into different amplitude weighting sub-areas according to the amplitude weighting value; 2. determining a 0 attenuation area, a maximum attenuation area and a interchange area, and T/R components respectively arranged in the 0 attenuation area and the maximum attenuation area and arranging the T/R components in the corresponding areas; 3. dividing the interchange area into a small attenuation interchange area, a large attenuation interchange area and a full interchange area, and arranging the rest T/R components to be arranged in the three areas according to the principle that the internal heat consumption is smaller than the external heat consumption; 4. the T/R assembly position adjustment is performed according to the total thickness of each column of T/R assemblies along the extrusion direction, so that the total thickness of each column is consistent. The invention can improve the reliability of the phased array system, improve the heat flux density distribution and improve the dynamic range of amplitude distribution.

Description

High-density integrated active phased array T/R assembly arrangement method based on machine, electricity and heat
Technical Field
The invention belongs to the technical field of antennas, and relates to a phased array T/R assembly arrangement method, in particular to a high-density integrated active phased array T/R assembly arrangement method based on machinery, electricity and heat.
Background
With the development of microwave technology and the continuous improvement of the demand, the scale of the active phased array is larger and larger, and the integration level is higher and higher. At present, a high-density integrated active phased array T/R assembly widely adopts a hard connection form, and the hard connection form has the risk of extrusion damage caused by assembly and thermal stress. And because the integration level is high, the heat flux density distribution is unbalanced, the temperature gradient is large, and the long-term stability and reliability of the system are influenced. In addition, in order to satisfy the designed array weighting value, the amplitudes of different radio frequency channels need to be trimmed, so that the dynamic range of the feed amplitude is lost, and the weighting effect of the large dynamic range of amplitude weighting is not favorable.
Disclosure of Invention
The invention aims to solve the problems and designs an arrangement method which can comprehensively consider the mechanical size, heat consumption and electrical performance of a T/R assembly. The method can effectively solve the risk of extrusion damage existing in the hard connection of the T/R components, realize the effective control of heat flow density distribution, and utilize the difference between the T/R components to furthest improve the dynamic range of the feeding amplitude under the condition of ensuring that all units meet the required feeding amplitude excitation, thereby being beneficial to the realization of array weighting. Thereby improving the safety, stability and long-term reliability of the system.
The design idea of the invention is as follows: dividing the T/R arrangement region into a 0 attenuation region, a maximum attenuation region and an interchange region; firstly, arranging the T/R components in the 0 attenuation region according to the structural size, the heat consumption and the electrical performance parameters of the T/R components, and then arranging the T/R components in the maximum attenuation region; dividing the interchange area into a small attenuation interchange area, a large attenuation interchange area and a full interchange area, and arranging the rest T/R assemblies to be arranged in the three areas according to the principle that the internal heat consumption is smaller than the external heat consumption; and finally, adjusting the position of the T/R assembly according to the total thickness of each column of the T/R assembly along the extrusion direction, so that the total thickness of each column is consistent.
According to the design concept, the technical scheme for achieving the purpose of the invention comprises the following steps:
step (1), the array scale is M multiplied by N, the radio frequency channel scale of each T/R component connected with the antenna unit is P multiplied by Q, and the array needs S multiplied by T/R components (wherein
Figure GDA0003117345910000011
Indicating rounding up X). Determining the amplitude weighted value Amp of the whole array obtained by optimization, wherein the maximum attenuation of the amplitude of the whole array obtained by optimization is AmaxdB, i.e. AmaxMax (| Amp |). Attenuator maximum gear attenuation of T/R assembly is A'maxdB, the attenuation step amount is Δ AdB.
Step (2), dividing the whole array area corresponding to the amplitude weighted value Amp into amplitude weighted sub-areas theta corresponding to S multiplied by T/R components according to the arrangement form of the T/R componentsst. The amplitude weighted sub-region ΘstCorresponding amplitude weighted value is AmpstI.e. the amplitude-weighted sub-region theta of row s and column tstWeighted value Amp ofstThe amplitude weighting value of each rf channel can be expressed as:
Figure GDA0003117345910000021
wherein P1, Q.
And (3) determining a 0 attenuation area according to the array amplitude weighted value Amp, determining and sequencing T/R components which can be arranged in the 0 attenuation area, and finally performing T/R component arrangement on the 0 attenuation area.
Step (4), counting the minimum value of the gain of the distributed T/R assembly radio frequency channel at 0dB attenuation 0 DEG phase shift corresponding to the position of 0dB attenuation with the amplitude weighted value of 0dB attenuation in the 0 attenuation region
Figure GDA0003117345910000022
And (5) determining a maximum attenuation area and an interchange area according to the array amplitude weighted value Amp, determining and sequencing T/R components which can be arranged in the maximum attenuation area, and finally arranging the T/R components in the maximum attenuation area.
And (6) dividing the interchange area into a small attenuation interchange area, a large attenuation interchange area and a full interchange area according to the Amp, the T/R assemblies arranged in the attenuation area 0, the T/R assemblies arranged in the maximum attenuation area and the gain characteristics of the rest T/R assemblies.
And (7) arranging the T/R components of the small attenuation interchange area, the large attenuation interchange area and the full interchange area according to the principle that the internal heat consumption is smaller than the external heat consumption.
And (8) counting the total thickness of each column of T/R assemblies along the extrusion direction, and adjusting the positions of the partial T/R assemblies in different columns according to the principle that the total thickness of each column is consistent.
Further, the step (1) of determining the optimized amplitude weighted value of the whole array as Amp includes the following steps:
step 1, if the phased array is a receiving phased array, the amplitude weighted value Amp is a receiving amplitude attenuation value Amp obtained by optimizationRI.e. Amp ═ AmpR
Step 2, if the phased array is a transmitting phased array, the amplitude weighted value Amp is a transmitting amplitude attenuation value Amp obtained by optimizationTI.e. Amp ═ AmpT
Step 3, if the phased array is a transmitting-receiving shared phased array, the amplitude weighted value Amp receives the dynamic range D of the amplitude weighted valueRAnd transmit amplitude weight dynamic range DTIs determined, i.e. is
Figure GDA0003117345910000031
Further, the step (3) of determining the 0 attenuation region according to the array amplitude weighting value Amp includes the following steps: setting 0 attenuation region judgment threshold
Figure GDA0003117345910000032
For the amplitude weighted sub-region ΘstThere is P e {1Satisfy } and Q ∈ { 1.,. Q }
Figure GDA0003117345910000033
Then the amplitude weighted sub-region ΘstBelonging to the 0 attenuation region. All S1, S, T1, T are judged to obtain the 0 attenuation region of the whole array.
Further, the step (3) of determining the T/R components capable of being arranged in the 0 attenuation region and sequencing the T/R components in the 0 attenuation region comprises the following steps:
step 1, setting a T/R component judgment threshold value of 0 attenuation region
Figure GDA0003117345910000034
Step 2, the amplitude weighting sub-region Θ for the 0 attenuation regionstWeighted value Amp ofstWhere all channels are 0dB attenuation, all T/R components can be used for this 0 attenuation magnitude weighted sub-region ΘstThe arrangement of (a) and (b).
Amplitude weighted sub-region Θ for the 0 attenuation regionstWeighted value Amp ofstThe condition that the channel attenuation value is not 0dB exists in the memory, when AmpstAll the 0dB attenuation positions in the radio frequency channel of the T/R component correspond to the minimum value of the gain when the 0dB attenuation is 0 DEG phase shift
Figure GDA0003117345910000035
Subtracting AmpstAll non-0 dB attenuation positions in the T/R component correspond to the gains of the radio frequency channels of the T/R component when the 0dB attenuation is 0 DEG phase shift
Figure GDA0003117345910000036
Are all less than or equal to
Figure GDA0003117345910000037
And AmpstAmplitude weighting values of all non-0 dB fading channels
Figure GDA0003117345910000038
(wherein
Figure GDA0003117345910000039
Figure GDA00031173459100000310
) Sum, i.e. when Q satisfies the following equation, P1
Figure GDA00031173459100000311
The T/R component can be used for the 0 attenuation magnitude weighting sub-region ΘstThe arrangement of (a) and (b).
When able to be arranged in the 0 attenuation amplitude weighting sub-region ΘstNumber of T/R components
Figure GDA00031173459100000312
Number of amplitude weighted sub-regions smaller than 0 attenuation region
Figure GDA00031173459100000313
Then choose in T/R component that cannot be used for 0 attenuation amplitude weighted sub-region
Figure GDA00031173459100000314
Individual weighted value AmpstMinimum value of gain at 0dB attenuation 0 degree phase shift of all 0dB attenuation positions
Figure GDA00031173459100000315
The largest T/R component is used for this 0 attenuation magnitude weighted sub-region arrangement. Statistics for the 0 attenuation amplitude weighting sub-region ΘstAll T/R components arranged at a weighted value AmpstThe minimum value of the gain of all T/R radio frequency channels at the corresponding position of 0dB when the phase is shifted by 0 DEG under the attenuation of 0dB
Figure GDA00031173459100000316
And the minimum value of the T/R module heat consumption PV according to the gain
Figure GDA00031173459100000317
Sequencing the T/R components to obtainThe sub-region ΘstCorresponding T/R arrangement sequence
Figure GDA00031173459100000318
Step 3, for phased arrays with different transceiving working states, G in step 20
Figure GDA00031173459100000319
And PV is selected according to the following method: if the phased array is a receive phased array, G0And
Figure GDA00031173459100000320
are reception state gains, and PV is the reception heat rate of the R module. If the phased array is a transmit phased array, G0And
Figure GDA0003117345910000041
both are gains at the emission operating point, and PV is the emission heat rate of the T/R module. If the phased array is a transmit-receive shared phased array, when DR≥DTWhen, G0And
Figure GDA0003117345910000042
are all receiving state gains, otherwise G0And
Figure GDA0003117345910000043
are the gains at the transmit operating point. And the heat consumption PV of the transceiving shared phased array is the emission heat consumption of the T/R assembly.
Further, the step of arranging the T/R modules in the attenuation region 0 in the step (3) is as follows:
step 1, counting each amplitude weighting sub-region theta in an attenuation region 0stWeighted value AmpstTotal number of channels attenuated by 0dB
Figure GDA0003117345910000044
When different amplitude weights the sub-region ΘstIs/are as follows
Figure GDA0003117345910000045
When not identical, the amplitude weighting sub-region ΘstIs arranged in the order of
Figure GDA0003117345910000046
The arrangement is performed in the order from large to small. When the amplitude weights the sub-region ΘstIs/are as follows
Figure GDA0003117345910000047
If, as such, the sub-regions do not have channels with non-0 dB attenuation, i.e.
Figure GDA0003117345910000048
When the sub-regions are the same and equal to P multiplied by Q, the sequence of arrangement among the sub-regions is randomly selected; if the sub-regions have channels with attenuation which is not 0dB, counting weighted values of all the channels with attenuation which is not 0dB in the sub-regions, and counting the number of the channels with the same attenuation according to the sequence of the attenuation from small to large, namely, firstly counting the number of the channels with the minimum attenuation in the attenuation which is not 0dB, and preferentially arranging the channels with the maximum attenuation if the number of the channels with the minimum attenuation is larger than the attenuation if the number of the channels with the minimum attenuation is still the same, and preferentially arranging the channels with the same number of the channels until the maximum attenuation A in the sub-regions is countedmaxdB, if the channel number of the maximum attenuation is the same, the weighting sub-region theta isstThe sequence of arrangement is randomly selected.
Step 2, determining each amplitude weighting sub-region theta in the 0dB attenuation region according to the step 1stThe arrangement order of the sub-regions is arranged in the range weighting sub-region thetastAccording to the T/R module
Figure GDA0003117345910000049
And the determined T/R components sequentially select the T/R components to be arranged.
Step 3, when the different amplitude weights the sub-region thetastIs/are as follows
Figure GDA00031173459100000410
Are all the same and equal to P × QI.e. different sub-regions ΘstWhen the weighted values of all the channels are 0dB attenuation, the position of the T/R components is adjusted according to the principle that the T/R components with high heat consumption are far away from the geometric center of the array surface and the T/R components with low heat consumption are close to the center after the T/R components are arranged according to the method. Firstly, counting the heat consumption PV of 0dB attenuation of weighted values of all channels of the T/R assembly before position interchangestAnd the corresponding sub-region ΘstIs a distance D between the geometric center of the array surface and the geometric center of the array surfacest. According to the distance DstPV (photovoltaic cell) in the T/R assembly which is arranged from large to small and needs to be subjected to position adjustment according to heat consumptionstThe T/R components are sequentially selected from large to small for position interchange. If there are multiple sub-regions ΘstD of (A)stUnder the same condition, selecting the T/R assemblies with the minimum heat consumption in the same number with the sub-regions from the rest T/R assemblies to be subjected to position adjustment, and randomly arranging the T/R assemblies in the sub-regions thetastIn (1).
Further, in the step (4), the minimum value of the gain of the distributed T/R assembly radio frequency channel at 0dB attenuation 0 DEG phase shift corresponding to the position where all amplitude weighted values in the 0dB attenuation area are 0dB attenuation is counted
Figure GDA00031173459100000411
The steps are as follows:
for all S1, 1., S, T1., T and P1., P, Q1., Q,
Figure GDA00031173459100000412
the minimum value of gain of 0dB attenuation 0 DEG phase shift of distributed T/R assembly radio frequency channels corresponding to all positions with 0dB attenuation amplitude weighted value in 0 attenuation region
Figure GDA0003117345910000051
If the phased array is a receive phased array, then
Figure GDA0003117345910000052
Is the receive gain. If the phased array is a transmit phased array, then
Figure GDA0003117345910000053
Is the gain at the transmit operating point. If the phased array is a transmit-receive shared phased array, when DR≥DTWhen the temperature of the water is higher than the set temperature,
Figure GDA0003117345910000054
gain for receiving state, otherwise
Figure GDA0003117345910000055
Is the gain at the transmit operating point.
Further, the step (5) of determining the maximum attenuation region and the interchange region is as follows:
step 1, counting all radio frequency channels of all the rest T/R components at the maximum gear A 'of the attenuator'maxMaximum in dB-attenuated 0 DEG phase-shifted gain
Figure GDA0003117345910000056
Wherein if the phased array is a receive phased array,
Figure GDA0003117345910000057
is the receive gain. If the phased array is a transmit phased array,
Figure GDA0003117345910000058
is the gain at the transmit operating point. If the phased array is a transmit-receive shared phased array, when DR≥DTWhen the temperature of the water is higher than the set temperature,
Figure GDA0003117345910000059
to receive gain, otherwise
Figure GDA00031173459100000510
Is the gain at the transmit operating point.
Step 2, setting a maximum attenuation area judgment threshold value
Figure GDA00031173459100000511
Step 3, weighting the amplitudeRegion ΘstIf P ∈ { 1.,. P } and Q ∈ { 1.,. Q } satisfy the following expression
Figure GDA00031173459100000512
The amplitude weighted sub-region ΘstBelonging to the region of maximum attenuation. And judging all S1, T1, S, T, and T to obtain the maximum attenuation area of the whole array.
Step 4, amplitude weighting sub-region theta which does not belong to the 0 attenuation region or the maximum attenuation regionstBelonging to the interchange area.
Step 5, if there is no amplitude weighting sub-region thetastBelonging to the maximum attenuation region, the remaining regions are all interchange regions.
Further, the step of determining and sorting the T/R components capable of being arranged in the maximum attenuation region in the step (5) is as follows:
step 1, setting a T/R component judgment threshold value of a maximum attenuation area
Figure GDA00031173459100000513
Step 2, the amplitude weighting sub-region theta for the maximum attenuation regionstWeighted value Amp ofstAll channels therein are AmaxdB attenuation, when a certain attenuation level alpha epsilon {0, delta Amax,...A′maxMake
Figure GDA00031173459100000514
And Amp corresponding to T/R componentstAll of A inmaxMaximum value of gain of dB attenuation channel when alpha dB attenuates 0 DEG phase shift
Figure GDA00031173459100000515
(wherein P1, Q) satisfies the following equation
Figure GDA00031173459100000516
The T/R component can be used for the maximum attenuation magnitude weighting sub-region ΘstArranging;
amplitude weighting sub-region Θ for the region of maximum attenuationstWeighted value Amp ofstThe channel attenuation value is different from AmaxdB, when there is a certain attenuation level α ∈ {0, Δ amax,...A′maxMake
Figure GDA0003117345910000061
Amp corresponding to T/R componentstAll of A inmaxMaximum value of gain of dB attenuation channel when alpha dB attenuates 0 DEG phase shift
Figure GDA0003117345910000062
All not AmaxdB attenuation channel at maximum gear A 'of T/R assembly attenuator'maxGain in dB attenuation at 0 ° phase shift
Figure GDA0003117345910000063
AmaxAnd all channel weights not subject to maximum attenuation
Figure GDA0003117345910000064
(wherein P1, Q) satisfies the following equation
Figure GDA0003117345910000065
The T/R component can be used for the maximum attenuation magnitude weighting sub-region ΘstThe arrangement of (a) and (b).
When able to be arranged in the maximum attenuation amplitude weighting sub-region ΘstNumber of T/R components
Figure GDA0003117345910000066
Less than the maximum attenuation region ΘstNumber of amplitude weighted sub-regions of
Figure GDA0003117345910000067
Then, the maximum attenuation amplitude weighting sub-region Θ cannot be arranged in the maximum attenuation amplitude weighting sub-regionstIn the T/R component of
Figure GDA0003117345910000068
Individual weighted value AmpstIn AmaxdB attenuation channel at A'maxMaximum gain in dB attenuation at 0 degree phase shift
Figure GDA0003117345910000069
The smallest T/R component is used for this maximum attenuation magnitude weighted sub-region arrangement. Statistics for the maximum attenuation amplitude weighting sub-region ΘstAll T/R components arranged at a weighted value AmpstIs AmaxRadio frequency channel at dB corresponding position of A'maxMaximum value of gain in dB attenuation 0 degree phase shift
Figure GDA00031173459100000610
And the T/R module heat loss PV, and according to the maximum value of the gain
Figure GDA00031173459100000611
The T/R components are sorted. Using the above method to weight the amplitude sub-region Θ that can be arranged within all maximum attenuation regionsstThe T/R components of (a) are sequenced to obtain the sub-region thetastCorresponding T/R arrangement sequence
Figure GDA00031173459100000612
And 3, if the phased array is a receiving phased array, the gain in the step 2 is the receiving gain, and the heat consumption is the receiving heat consumption. If the phased array is a transmitting phased array, the gains in the step 2 are the gains at the transmitting working point, and the heat consumption is the transmitting heat consumption. If the phased array is a transmit-receive shared phased array, when DR≥DTAnd if not, the gains in the step 2 are the gains at the transmitting working point. The heat consumption of the transmitting and receiving shared phased array is the transmitting heat consumption.
Further, the step of arranging the T/R module for the maximum attenuation region in step (5) is as follows:
step 1, counting each amplitude weighting sub-region theta in the maximum attenuation regionstWeighted value AmpstIs AmaxTotal number of channels in dB attenuation
Figure GDA00031173459100000613
When different amplitude weights the sub-region ΘstIs/are as follows
Figure GDA00031173459100000614
When not identical, the amplitude weighting sub-region ΘstIs arranged in the order of
Figure GDA00031173459100000615
The arrangement is performed in the order from large to small. When the amplitude weights the sub-region ΘstIs/are as follows
Figure GDA00031173459100000616
If, as such, these subregions do not have amaxChannels with dB attenuation, i.e.
Figure GDA00031173459100000617
When the sub-regions are the same and equal to P multiplied by Q, the sequence of arrangement among the sub-regions is randomly selected; if these subregions have a value other than AmaxdB-attenuated channels, then all non-A in these sub-regions are countedmaxThe weighted value of dB attenuation channels is counted according to the order of the attenuation from large to small, that is, the number of channels with the same attenuation is counted firstlymaxThe channel number with the maximum attenuation in dB attenuation is distributed preferentially if the channel number of the attenuation is more, the channel number with the maximum attenuation smaller than the attenuation is counted if the channel number is still the same, the channel number with the same channel number is distributed preferentially if the channel number is more, and the process is repeated until the minimum attenuation in the sub-regions is counted, if the channel number of the minimum attenuation is still the same, the weighted sub-region theta is obtainedstThe sequence of arrangement is randomly selected.
Step 2, determining each amplitude weighting sub-region theta in the maximum attenuation region according to the step 1stThe arrangement order of the sub-regions is arranged in the range weighting sub-region thetastAccording to the T/R module
Figure GDA0003117345910000071
And selecting T/R for arrangement according to the determined arrangement sequence of T/R.
Step 3, when the different amplitude weights the sub-region thetastIs/are as follows
Figure GDA0003117345910000072
All equal to PxQ, i.e. different sub-regions thetastThe weighted value of all channels is maximum AmaxWhen dB attenuates, the position of the T/R assembly is adjusted according to the principle that the T/R assembly with high heat consumption is far away from the geometric center of the array surface and the T/R assembly with low heat consumption is close to the center after the T/R assembly is arranged according to the method. Namely, firstly, the heat consumption PV of the T/R assembly before position interchange is countedstAnd the corresponding sub-region ΘstIs a distance D between the geometric center of the array surface and the geometric center of the array surfacest. According to the distance DstIn the T/R module which is arranged and is to be subjected to position adjustment according to heat consumption PV in the order from big to smallstThe T/R components are sequentially selected from large to small for position interchange. If a plurality of sub-regions ΘstD of (A)stIf the T/R components are the same as the sub-regions, selecting the T/R components with the minimum heat consumption in the rest T/R components to be subjected to position adjustment, and randomly arranging the T/R components in the sub-regions thetastIn (1).
Further, the step (6) of dividing the interchange area into a small attenuation interchange area, a large attenuation interchange area and a full interchange area comprises the following steps:
step 1, counting the minimum value of all radio frequency channels of all the rest T/R components in gain of 0dB attenuation and 0 DEG phase shift
Figure GDA0003117345910000073
All radio frequency channels of all remaining T/R components are at attenuator maximum gear A'maxMaximum in dB-attenuated 0 DEG phase-shifted gain
Figure GDA0003117345910000074
All radio frequency channels of all remaining T/R components are at attenuator maximum gear A'maxMinimum in dB-attenuated 0 DEG phase-shifted gain
Figure GDA0003117345910000075
If the phased array is a receiving phased array, the gains in the step 1 are all receiving gains. And if the phased array is a transmitting phased array, the gains in the step 1 are the gains at the transmitting working point. If the phased array is a transmit-receive shared phased array, when DR≥DTAnd if not, the gain in the step 1 is the gain at the transmitting working point.
Step 2, setting a judgment threshold value of a small attenuation interchange area
Figure GDA0003117345910000078
When there is P e { 1.,. P } and Q e { 1.,. Q } in the interchange area satisfy
Figure GDA0003117345910000076
Then the amplitude weighted sub-region ΘstBelonging to the exchange region of small attenuation.
Step 3, setting a judgment threshold value of a large attenuation interchange area
Figure GDA0003117345910000079
When there is P e { 1.,. P } and Q e { 1.,. Q } in the interchange area satisfy
Figure GDA0003117345910000077
Then the amplitude weighted sub-region ΘstBelonging to the exchange area of large attenuation.
Step 4, for theta in the interchange area which neither belongs to the small attenuation interchange area nor the large attenuation interchange areastIs a fully interchanged region.
Further, the step (7) of arranging the remaining T/R modules to be arranged in the small attenuation interchange area, the large attenuation interchange area and the full interchange area includes the steps of:
step 1, counting all sub-regions to be distributed thetastIs a distance D between the geometric center of the array surface and the geometric center of the array surfacest
Step 2, according to DstThe T/R components are arranged in sequence from small to large by using the following method:
setting T/R component judgment threshold of small attenuation interchange area
Figure GDA0003117345910000081
And large attenuation interchange area T/R component judgment threshold
Figure GDA0003117345910000082
② if the amplitude weighting sub-region theta to be arrangedstBelonging to the small attenuation interchange area, the amplitude weighting sub-area ΘstWeighted value of AmpstFor all P e { 1.,. P } and Q e { 1.,. Q }, there is a
Figure GDA0003117345910000083
So that
Figure GDA0003117345910000084
Subtracting the weighted value Amp corresponding to the T/R componentstIn
Figure GDA0003117345910000085
Attenuating gain at 0 ° phase shift
Figure GDA0003117345910000086
Are all less than or equal to the amplitude weighted value of the channel
Figure GDA0003117345910000087
And
Figure GDA0003117345910000088
when the difference of (a) is greater than or equal to 1, i.e., when Q satisfies the following equation for all P, i.e., P, Q, i.e., 1
Figure GDA0003117345910000089
The T/R component can be used for the small attenuation swap area ΘstThe arrangement of (a) and (b). When able to be arranged in the small attenuation interchange area thetastNumber of T/R components
Figure GDA00031173459100000810
Number of amplitude weighted sub-regions smaller than the small attenuation interchange region
Figure GDA00031173459100000811
Then choose in T/R module that can' T be used for small attenuation interchange area
Figure GDA00031173459100000812
Individual weighted value AmpstMinimum value of gain at 0dB attenuation 0 degree phase shift of all 0dB attenuation positions
Figure GDA00031173459100000813
The largest T/R-component is used for the arrangement of the small attenuation interchange area. Counting the heat consumption PV of all T/R modules capable of being used in the regionstSelective heat loss PVstThe lowest is arranged in the amplitude-weighted sub-region Θst
Third if the amplitude weighting sub-region theta to be arrangedstBelonging to the region of large attenuation interchange, the amplitude-weighted sub-region ΘstWeighted value of AmpstFor all P e { 1.,. P } and Q e { 1.,. Q }, there is a
Figure GDA00031173459100000814
So that
Figure GDA00031173459100000815
Subtracting the weighted value Amp corresponding to the T/R componentstIn
Figure GDA00031173459100000816
Attenuating gain at 0 ° phase shift
Figure GDA00031173459100000817
Are all larger than or equal to the amplitude weighted value of the channel
Figure GDA00031173459100000818
And
Figure GDA00031173459100000819
when the sum of (a) and (b) is equal to 1, Q satisfies the following equation
Figure GDA00031173459100000820
The T/R component can be used for the large attenuation interchange area ΘstThe arrangement of (a) and (b). When able to be arranged in the small attenuation interchange area thetastNumber of T/R components
Figure GDA00031173459100000821
Number of amplitude weighted sub-regions smaller than the large attenuation interchange region
Figure GDA00031173459100000822
Then choose in T/R module that can' T be used for large attenuation interchange area
Figure GDA0003117345910000091
Individual weighted value AmpstIn AmaxdB attenuation channel at A'maxMaximum gain in dB attenuation at 0 degree phase shift
Figure GDA0003117345910000092
The smallest T/R-component is used for the arrangement of this large attenuation interchange area. Counting the heat consumption PV of all T/R modules capable of being used in the regionstSelective heat loss PVstThe lowest is arranged in the amplitude-weighted sub-region Θst
Fourthly, if the amplitude weighting sub-region theta to be arrangedstBelonging to the fully-interchanged region, all the remaining T/R devices can be arranged in the region. Counting the heat consumption PV of all T/R modules capable of being used in the regionstSelectingHeat loss PVstThe lowest is arranged in the amplitude-weighted sub-region Θst
Step 3, if there are multiple subregions DstIf the sub-regions are equal, the arrangement sequence of the sub-regions is randomly determined, and the arrangement method of each sub-region is the same as that in the step 2.
If the phased array is a receiving phased array, the gains in the step 2 are all receiving gains, and the heat consumption is the receiving heat consumption. If the phased array is a transmitting phased array, the gains in the step 2 are the gains at the transmitting working point, and the heat consumption is the transmitting heat consumption. If the phased array is a transmit-receive shared phased array, when DR≥DTAnd if not, the gains in the step 2 are the gains at the transmitting working point. The heat consumption of the transmitting and receiving shared phased array is the transmitting heat consumption.
And 4, repeating the steps from the step 2 to the step 3 until all the T/R components are arranged.
Further, the adjusting of the position of the T/R assembly according to the total thickness of each column of T/R assemblies along the extrusion direction in the step (8) comprises the following steps:
and step 1, setting T columns and S rows of T/R assemblies in the extrusion direction. The average thickness of all T/R components (S X T)
Figure GDA0003117345910000093
The desired total thickness of each column of T/R assemblies is then
Figure GDA0003117345910000094
Step 2, counting the total thickness { H) of all the rows of T/R components1,H2,...,Ht,...,HT}. Setting a total thickness tolerance threshold value epsilon of each row of extrusion directionsH
Step 3, calculating the difference value { delta H) between the actual total thickness and the expected total thickness of each column of T/R1,ΔH2,...,ΔHt,...,ΔHTThat is to say
ΔHt=Ht-HExpWherein (T ═ 1, 2.. T)
In the step 4, the step of mixing,find { Δ H1,ΔH2,...,ΔHt,...,ΔHTMaximum value of
Figure GDA0003117345910000095
Namely, it is
Figure GDA0003117345910000096
Figure GDA0003117345910000097
And minimum value
Figure GDA0003117345910000098
Namely, it is
Figure GDA0003117345910000099
Wherein the maximum value is located at column t1 and the minimum value is located at column t2, such that
Figure GDA00031173459100000910
Step 5, listing all T/R components in the T1 th column
Figure GDA00031173459100000911
And all T/R modules in column T2
Figure GDA00031173459100000912
Calculating the difference deltaH 'between the actual total thickness of the T1 th column and the expected total thickness of the T2 th column and any T/R assembly of the T1 th column and any T/R assembly of the T2 th column after being interchanged't1And Δ H't2. Order to
Figure GDA00031173459100000913
Δ H 'after interchanging the e th T/R Assembly from column T1 and the f th T/R Assembly from column T2't1And Δ H't2Maximum value of absolute value, i.e.
Figure GDA0003117345910000101
Marking according to the type of the region where the T/R component is positioned
Figure GDA0003117345910000102
And
Figure GDA0003117345910000103
interchange type mark LefAs shown in Table 12-1, all e ∈ {1, 2.. The P } and f ∈ {1, 2.. The P } combinations are traversed to get information about
Figure GDA0003117345910000104
Matrix MatirxΔHAnd about LefMatrix MatirxL. In MatirxLIn finding interchange type flag minimum value Lmin=min(MatirxL)。
TABLE 12-1 interchange type flag
Figure GDA0003117345910000105
Figure GDA0003117345910000106
Figure GDA0003117345910000107
Step 6, find interchange type mark LminCorresponding matrix Matirx of the positionΔHMinimum value of Medium element Δ H'minI.e. Δ H'min=min(MatirxΔH|Lmin) Where e ∈ {1, 2., P }, and f ∈ {1, 2., P }. And recording the minimum value Δ H'minCorresponding interchangeable T/R assembly
Figure GDA0003117345910000108
And
Figure GDA0003117345910000109
step 7, when Δ Hmax>ΔH′minWhen, to
Figure GDA00031173459100001010
And
Figure GDA00031173459100001011
the components are subjected to position interchange, so that the position of the components is exchanged by delta Hmax=ΔH′minAnd sets the T/R component interchange blacklist to the whitelist (0, 0).
When Δ Hmax≤ΔH′minWhen the position is not exchanged, the T/R component exchange blacklist is set as (B)1,B2) In which B is1=t1,B2T2, i.e. B-th of T/R module1Column and B2Columns are prohibited from interchanging. At this time, { Δ H is found out of the list where the T/R device interchange blacklists1,ΔH2,...,ΔHt,...,ΔHTMaximum value of
Figure GDA0003117345910000111
Column t1, and minimum value
Figure GDA0003117345910000112
The t2 th columns, i.e., t1 and t2, satisfy the following formula
t1≠B1And t2 ≠ B2
Calculate the time of
Figure GDA0003117345910000113
Then repeating the above 5 th step to the 7 th step until the { H }1,H2,...,HTThere are no interchangeable T/R components in any two columns.
Step 8, if any two columns have no interchangeable T/R components, enabling L in step 6min=Lmin+1. Then repeating the above steps 6 to 7 until Lmin=15。
Step 9, repeating the steps from step 3 to step 8 until delta Hmax≤εHOr { H1,H2,...,HTAnd when any two columns in the T/R module are not interchanged, stopping the position adjustment of the T/R module.
Compared with the prior art, the invention has the following characteristics and advantages:
the invention comprehensively considers the factors of mechanical size, heat consumption, electrical property and the like of the T/R component. The T/R components with high gain are arranged in the area with less weight attenuation, the T/R components with low gain are arranged in the area with more weight attenuation, and the dynamic range of the feeding amplitude is improved by using the difference between the T/R components. The method of arranging the high heat consumption T/R components at the outer side and arranging the low heat consumption components at the inner side can realize effective control of heat flow density distribution. Finally, the total thickness of each column of T/R components along the extrusion direction is consistent by adjusting the individual T/R components, and the risk of extrusion damage existing in the hard connection of the T/R components can be effectively solved. Thereby improving the safety, stability and long-term reliability of the system.
Drawings
FIG. 1 is a flow chart of an implementation of the method of the present invention.
FIG. 2 is a schematic structural diagram of a T/R module according to an embodiment of the present invention.
FIG. 3 tabulates array receive amplitude weighted attenuation values (dB) for an embodiment of the present invention.
FIG. 4 tabulates the T/R component receive gain (dB) for an embodiment of the present invention.
FIG. 5 tabulates the emission heat rate of the T/R assembly of an embodiment of the invention.
FIG. 6 tabulates the thickness of the T/R assembly of an embodiment of the invention.
Figure 7 shows in tabular form the amplitude weighted sub-regions of an embodiment of the present invention.
Fig. 8 tabulates the amplitude weighted subregion types of embodiments of the present invention, where "C1" represents a 0 attenuation region, "C2" represents a maximum attenuation region, "C3" represents a small attenuation interchange region, "C4" represents a large attenuation interchange region, and "C5" represents a full attenuation interchange region.
FIG. 9 tabulates the attenuation region Θ of an embodiment of the present invention at 0stUsable T/R module and
Figure GDA0003117345910000121
FIG. 10 is a table showing the results of the arrangement of the T/R elements in the 0 attenuation region according to the embodiment of the present invention.
FIG. 11 tabulates the maximum attenuation region Θ of an embodiment of the present inventionstUsable T/R module and
Figure GDA0003117345910000122
FIG. 12 is a table showing the results of the placement of the T/R elements in the maximum attenuation region for an embodiment of the present invention.
FIG. 13 is a table showing the results of the placement of T/R components in the swap area according to an embodiment of the present invention.
FIG. 14 tabulates the T/R case thickness results (mm) after the interchange area T/R assembly arrangement of an embodiment of the present invention.
FIG. 15 is a table showing the arrangement of the T/R elements after the position adjustment according to the embodiment of the present invention.
FIG. 16 tabulates T/R assembly housing thickness results (mm) after position adjustment for embodiments of the present invention.
FIG. 17 is a table showing the T/R module emission heat rate (W) after position adjustment for an embodiment of the present invention.
Detailed Description
Referring to fig. 1, the present invention is given as an example below. The array adopts a uniform planar array distributed by a 20-element-20-element rectangular grid with a rectangular caliber. The rf channels connecting each T/R element to the antenna element are 1 x 4 in size, as shown in fig. 2. The maximum gear attenuation of the attenuator of the T/R component is 15.5dB, and the attenuation step amount is 0.5 dB. The array is a transmit-receive shared array, wherein the receive weighted attenuation values are shown in fig. 3, and the receive amplitude weighted values have a dynamic range of 15.5 dB. The transmission adopts a constant amplitude weighted value, and the dynamic range of the transmission amplitude weighted value is 0 dB. The T/R module 0dB attenuation 0 degree phase shift receiving gain, 15.5dB attenuation 0 degree phase shift receiving gain as shown in figure 4. The T/R module emission heat dissipation is shown in FIG. 5. The T/R assembly thickness is shown in FIG. 6.
Example (b): the T/R component arrangement is carried out on the transmitting-receiving shared array by referring to the method of the invention.
The method comprises the following steps:
the array is M × N (M is 20, N is 20), the radio frequency channel connected to the antenna unit is P × Q (P is 1, Q is 4) per T/R component, and the number of T/R components required by the array is S × T (S is 20, T is 5). Attenuator maximum gear attenuation of T/R assembly is A'max15.5dB, the attenuation step Δ a is 0.5 dB. The dynamic range D of the receiving amplitude weight valueR15.5dB, dynamic range of transmit amplitude weights D T0 dB. Due to DR≥DTAnd Amp is the receive array weight (as shown in figure 3). The maximum attenuation of the amplitude of the whole array obtained by optimization is Amax=15.5dB。
Step two:
dividing the whole array region corresponding to the amplitude weighted value Amp into amplitude weighted sub-regions theta corresponding to 20-5T/R components according to the arrangement form of the T/R componentsst. The amplitude weighted sub-region ΘstAs shown in fig. 7.
Step three:
determining a 0 attenuation area according to the array amplitude weighted value Amp, determining and sequencing T/R components which can be arranged in the 0 attenuation area, and finally performing T/R component arrangement on the 0 attenuation area:
(3a) setting 0 attenuation region judgment threshold
Figure GDA0003117345910000131
For the amplitude weighted sub-region ΘstIf q ∈ { 1., 4} satisfies the following equation
Figure GDA0003117345910000132
The amplitude weighted sub-region ΘstBelonging to the 0 attenuation region. A decision is made for all s 1, 20, t1, 0 attenuation regions for the entire array, as shown in fig. 8.
(3b) Setting 0 attenuation region T/R component judgment threshold
Figure GDA00031173459100001314
Due to amplitudeWeighted sub-region Θ10 3And Θ11 3All channels within the weighted value of (a) are 0dB attenuated, so the T/R component can be used for the arrangement of the 0 attenuated magnitude weighted sub-region.
Due to the amplitude weighted sub-region Θ9 3And Θ12 3When the channel attenuation value is not 0dB, the AmpstAll the 0dB attenuation positions in the radio frequency channel of the T/R component correspond to the minimum value of the gain when the 0dB attenuation is 0 DEG phase shift
Figure GDA0003117345910000133
Subtracting AmpstAll non-0 dB attenuation positions in the T/R component correspond to the gains of the radio frequency channels of the T/R component when the 0dB attenuation is 0 DEG phase shift
Figure GDA0003117345910000134
Are all less than or equal to
Figure GDA0003117345910000135
And AmpstAmplitude weighting values of all non-0 dB fading channels
Figure GDA0003117345910000136
(wherein
Figure GDA0003117345910000137
) Sum, i.e. when all q 1, 4 satisfy the following equation
Figure GDA0003117345910000138
The T/R component can be used for the 0 attenuation magnitude weighting sub-region ΘstThe arrangement of (a) and (b).
(3c) Statistics can be used for the 0 attenuation amplitude weighting sub-region ΘstAll T/R components arranged at a weighted value AmpstThe minimum value of the gain of all T/R radio frequency channels at the corresponding position of 0dB when the phase is shifted by 0 DEG under the attenuation of 0dB
Figure GDA0003117345910000139
And the minimum value of the T/R module heat consumption PV according to the gain
Figure GDA00031173459100001310
The T/R components are sorted. Using the above method to weight the amplitude sub-region Θ that can be arranged within all 0 attenuation regionsstThe T/R components of (a) are sequenced to obtain the sub-region thetastCorresponding T/R arrangement sequence
Figure GDA00031173459100001311
Can be arranged at theta9 3,Θ10 3,Θ11 3And theta12 3And their corresponding order
Figure GDA00031173459100001312
And
Figure GDA00031173459100001313
as shown in fig. 9.
(3d)Θ10 3And Θ11 3All channels within the weighted value of (a) are 0dB attenuation, so the theta is preferentially arranged10 3And Θ11 3Region T/R component, and arranging theta9 3And Θ12 3And (4) a region T/R component. The result of the arrangement of the T/R elements in the attenuation region of 0 is shown in FIG. 10.
Step four:
counting the minimum value of the receiving gain of 0dB attenuation 0 degree phase shift of the distributed T/R assembly radio frequency channel corresponding to the position with 0dB attenuation in the 0 attenuation area
Figure GDA0003117345910000141
Step five:
(5a) counting all radio frequency channels of all residual T/R components at maximum gear A 'of attenuator'maxMaximum in dB-attenuated 0 DEG phase-shifted receive gain
Figure GDA0003117345910000142
Setting a maximum attenuation region judgment threshold
Figure GDA0003117345910000143
(5b) For the amplitude weighted sub-region ΘstIf q ∈ { 1., 4} satisfies the following equation
Figure GDA0003117345910000144
The amplitude weighted sub-region ΘstBelonging to the region of maximum attenuation. A decision is made for all s 1, 20, t1, 5, resulting in the maximum attenuation region for the entire array. Amplitude-weighted sub-region Θ not belonging to the 0 attenuation region nor to the maximum attenuation regionstBelonging to the interchange area. The region of maximum attenuation is shown in figure 8.
(5c) Setting a maximum attenuation region T/R component judgment threshold
Figure GDA0003117345910000145
Amplitude weighting sub-region Θ of the region of maximum attenuation2 1,Θ2 5,Θ19 1And theta19 5Weighted value Amp ofstAll channels therein are AmaxdB attenuation, when there is some attenuation step alpha ∈ {0, Δ Amax,...A′maxMake
Figure GDA0003117345910000146
And Amp corresponding to T/R componentstAll of A inmaxMaximum value of gain of dB attenuation channel when alpha dB attenuates 0 DEG phase shift
Figure GDA0003117345910000147
(wherein q is 1.., 4.) satisfying the following equation
Figure GDA0003117345910000148
The T/R component can be used for the maximum attenuation magnitude weighting sub-region Θ2 1,Θ2 5,Θ19 1And theta19 5Arranging;
amplitude weighting sub-region Θ of the region of maximum attenuation1 1,Θ3 1,Θ4 1,Θ1 5,Θ3 5,Θ4 5,Θ20 1,Θ18 1,Θ17 1,Θ20 5,Θ18 5And theta17 5Weighted value Amp ofstThe channel attenuation value is different from AmaxdB, when there is some attenuation step α ∈ {0, Δ Amax,...A′maxMake
Figure GDA0003117345910000149
Amp corresponding to T/R componentstAll of A inmaxMaximum value of gain of dB attenuation channel when alpha dB attenuates 0 DEG phase shift
Figure GDA00031173459100001410
All not AmaxdB attenuation channel at maximum gear A 'of T/R assembly attenuator'maxGain in dB attenuation at 0 ° phase shift
Figure GDA00031173459100001411
AmaxAnd all channel weights not subject to maximum attenuation
Figure GDA00031173459100001412
(wherein q is 1, 4) each of the above-mentioned formulas
Figure GDA00031173459100001413
The T/R component can be used for the maximum attenuation magnitude weighting sub-region Θ1 1,Θ3 1,Θ4 1,Θ1 5,Θ3 5,Θ4 5,Θ20 1,Θ18 1,Θ17 1,Θ20 5,Θ18 5And theta17 5The arrangement of (a) and (b).
(5d) Using the above method for webs that can be laid out in all areas of maximum attenuationDegree weighted sub-region ΘstThe T/R components of (a) are sequenced to obtain the sub-region thetastCorresponding T/R arrangement sequence
Figure GDA0003117345910000151
Can be arranged at theta2 1,Θ2 5,Θ19 1,Θ19 5,Θ1 1,Θ3 1,Θ4 1,Θ1 5,Θ3 5,Θ4 5,Θ20 1,Θ18 1,Θ17 1,Θ20 5,Θ18 5And theta17 5The T/R modules of (a) and their corresponding sequence are shown in fig. 11(a) to 11 (d).
Θ2 1,Θ2 5,Θ19 1And Θ19 5All channels within the weighted value of (a) are 15.5dB attenuation, so the order of Θ is preferentially arranged2 1,Θ2 5,Θ19 1And Θ19 5And the other amplitude weighting subregion T/R assemblies are arranged. The result of the maximum attenuation region T/R module arrangement is shown in FIG. 12.
Step six:
(6a) counting the minimum value of the gain of all radio frequency channels of all the rest T/R components at 0dB attenuation 0 DEG phase shift
Figure GDA0003117345910000152
All radio frequency channels of all remaining T/R components are at attenuator maximum gear A'maxMaximum in dB-attenuated 0 DEG phase-shifted gain
Figure GDA0003117345910000153
All radio frequency channels of all remaining T/R components are at attenuator maximum gear A'maxMinimum in dB-attenuated 0 DEG phase-shifted gain
Figure GDA0003117345910000154
(6b) Setting judgment threshold of small attenuation interchange area
Figure GDA0003117345910000155
Satisfy for any q ∈ { 1., 4} in the swap region
Figure GDA0003117345910000156
Then the amplitude weighted sub-region ΘstBelonging to the exchange region of small attenuation.
(6c) Setting a judgment threshold value of a large attenuation interchange area
Figure GDA0003117345910000157
Satisfy for any q ∈ { 1., 4} in the swap region
Figure GDA0003117345910000158
Then the amplitude weighted sub-region ΘstBelonging to the exchange area of large attenuation.
(6d) For Θ in the exchange region that neither belongs to the small nor the large attenuation exchange regionstIs a fully interchanged region. The interchange area of the entire array is shown in fig. 8.
Step seven:
(7a) counting all sub-regions to be arranged thetastIs a distance D between the geometric center of the array surface and the geometric center of the array surfacest
(7b) According to DstThe T/R components are arranged in sequence from small to large by using the following method:
setting T/R component judgment threshold of small attenuation interchange area
Figure GDA0003117345910000159
And large attenuation interchange area T/R component judgment threshold
Figure GDA00031173459100001510
② if the amplitude weighting sub-region theta to be arrangedstBelonging to the small attenuation interchange area, the amplitude weighting sub-area ΘstWeighted value of AmpstFor any q e { 1.,. 4}, there is a
Figure GDA00031173459100001511
So that
Figure GDA00031173459100001512
Subtracting the weighted value Amp corresponding to the T/R componentstIn
Figure GDA00031173459100001513
Attenuating gain at 0 ° phase shift
Figure GDA00031173459100001514
Are all less than or equal to the amplitude weighted value of the channel
Figure GDA00031173459100001515
And
Figure GDA00031173459100001516
when the difference of (a) is equal to 1, that is, when all q's satisfy the following equation (4)
Figure GDA00031173459100001517
The T/R component can be used for the small attenuation swap area ΘstThe arrangement of (a) and (b). Counting the heat consumption PV of all T/R modules capable of being used in the regionstSelective heat loss PVstThe lowest is arranged in the amplitude-weighted sub-region Θst
Third if the amplitude weighting sub-region theta to be arrangedstBelonging to the region of large attenuation interchange, the amplitude-weighted sub-region ΘstWeighted value of AmpstFor any q e { 1.,. 4}, there is a
Figure GDA0003117345910000161
So that
Figure GDA0003117345910000162
Subtracting the weighted value Amp corresponding to the T/R componentstIn
Figure GDA0003117345910000163
Attenuating gain at 0 ° phase shift
Figure GDA0003117345910000164
Are all larger than or equal to the amplitude weighted value of the channel
Figure GDA0003117345910000165
And
Figure GDA0003117345910000166
when the sum of (a) and (b), i.e., when all q 1, 4 satisfy the following equation
Figure GDA0003117345910000167
The T/R component can be used for the large attenuation interchange area ΘstThe arrangement of (a) and (b). Counting the heat consumption PV of all T/R modules capable of being used in the regionstSelective heat loss PVstThe lowest is arranged in the amplitude-weighted sub-region Θst
Fourthly, if the amplitude weighting sub-region theta to be arrangedstBelonging to the fully-interchanged region, all the remaining T/R devices can be arranged in the region. Heat consumption PV of all T/R assembliesstSelective heat loss PVstThe lowest is arranged in the amplitude-weighted sub-region Θst
(7c) If there are a plurality of sub-regions DstAnd if the sub-regions are equal to each other, the arrangement sequence among the sub-regions is randomly determined, and the arrangement method of each sub-region is the same as that in (7 b).
(7d) And (4) repeating the steps (7b) to (7c) until all the T/R assemblies are arranged. The results of the arrangement of the T/R elements in the interchanged regions are shown in FIG. 13.
Step eight:
(8a) the T/R extrusion direction is along the X direction in the figure 2, and the T/R extrusion direction has 5 columns and 20 rows of T/R assemblies. The average thickness of all T/R components (S X T)
Figure GDA0003117345910000168
The desired total thickness of each column of T/R assemblies is then
Figure GDA0003117345910000169
(8b) Count the total thickness { H ] of all column T/R components1,H2,...,Ht,...,HTAs shown in fig. 14. Setting a total thickness tolerance threshold value epsilon of each row of extrusion directionsH=0.500mm。
(8c) Calculating the difference between the actual total thickness and the expected total thickness of each column of T/R { Delta H1,ΔH2,...,ΔHt,...,ΔHTThat is to say
ΔHt=Ht-HExpWherein (T ═ 1, 2.. T)
(8d) Find { Δ H1,ΔH2,...,ΔHt,...,ΔHTMaximum value of
Figure GDA00031173459100001610
Namely, it is
Figure GDA00031173459100001611
Figure GDA00031173459100001612
And minimum value
Figure GDA00031173459100001613
Namely, it is
Figure GDA00031173459100001614
Wherein the maximum value is located at column t1 and the minimum value is located at column t2, such that
Figure GDA00031173459100001615
(8e) Listing all T/R Components in column T1
Figure GDA00031173459100001616
And all T/R modules in column T2
Figure GDA00031173459100001617
Calculating any T/R component in the T1 th column and any T/R in the T2 th columnDifference Δ H 'of actual total thickness from desired total thickness for columns t1 and t2 after assembly interchange't1And Δ H't2. Order to
Figure GDA0003117345910000171
Δ H 'after interchanging the e th T/R Assembly from column T1 and the f th T/R Assembly from column T2't1And Δ H't2Maximum value of absolute value, i.e.
Figure GDA0003117345910000172
Marking according to the type of the region where the T/R component is positioned
Figure GDA0003117345910000173
And
Figure GDA0003117345910000174
interchange type mark LefAs shown in Table 12-1. All e {1, 2.. eta.,. P } and f e {1, 2.. eta.,. P } combinations are traversed to obtain information about
Figure GDA0003117345910000175
Matrix MatirxΔHAnd about LefMatrix MatirxL. In MatirxLIn finding interchange type flag minimum value Lmin=min(MatirxL)。
TABLE 12-1 interchange type flag
Figure GDA0003117345910000176
Figure GDA0003117345910000177
Figure GDA0003117345910000178
(8f) Finding interchange type flag LminCorresponding matrix Matirx of the positionΔHMinimum value of Medium element Δ H'minI.e. Δ H'min=min(MatirxΔH|Lmin) Where e ∈ {1, 2., P }, and f ∈ {1, 2., P }. And recording the minimum value Δ H'minCorresponding interchangeable T/R assembly
Figure GDA0003117345910000179
And
Figure GDA00031173459100001710
(8g) when Δ Hmax>ΔH′minWhen, to
Figure GDA0003117345910000181
And
Figure GDA0003117345910000182
the components are subjected to position interchange, so that the position of the components is exchanged by delta Hmax=ΔH′minAnd sets the T/R component interchange blacklist to the whitelist (0, 0).
When Δ Hmax≤ΔH′minWhen the position is not exchanged, the T/R component exchange blacklist is set as (B)1,B2) In which B is1=t1,B2T2, i.e. B-th of T/R module1Column and B2Columns are prohibited from interchanging. At this time, { Δ H is found out of the list where the T/R device interchange blacklists1,ΔH2,...,ΔHt,...,ΔHTMaximum value of
Figure GDA0003117345910000183
Column t1, and minimum value
Figure GDA0003117345910000184
The t2 th columns, i.e., t1 and t2, satisfy the following formula
t1≠B1And t2 ≠ B2
Calculate the time of
Figure GDA0003117345910000185
Subsequently repeating the above (8e) to step (8g) until { H }1,H2,...,HTThere are no interchangeable T/R components in any two columns.
(8h) If there is no interchangeable T/R component in any two columns, let L in (8f)min=Lmin+1. Subsequently repeating the above (8f) to step (8g) until Lmin=15。
(8i) Repeating the steps (8c) to (8H) until Δ Hmax≤εHOr { H1,H2,...,HTAnd when any two columns in the T/R module are not interchanged, stopping the position adjustment of the T/R module.
The position adjustment of the T/R assembly is completed after 6 times of T/R assembly interchange. The 6 times of interchanging sequence are A22 to A64, A13 to A55, A2 to A69, A60 to A5, A32 to A59 and A51 to A65. The arrangement result of the T/R assembly after the position adjustment is shown in FIG. 15, the thickness result of the T/R shell of the assembly after the position adjustment is shown in FIG. 16, and the emission heat consumption of the T/R assembly after the position adjustment is shown in FIG. 17.
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 (12)

1. A high-density integrated active phased array T/R component arrangement method based on machine, electricity and heat is characterized by comprising the following steps:
step (1), the array scale is M multiplied by N, the radio frequency channel scale of each T/R component connected with the antenna unit is P multiplied by Q, and the array needs S multiplied by T/R components in total, wherein
Figure FDA0003117345900000014
Represents rounding up to X; determining the amplitude weighted value Amp of the whole array obtained by optimization, wherein the maximum attenuation of the amplitude of the whole array obtained by optimization is AmaxdB, i.e. AmaxMax (| Amp |), attenuator maximum gear attenuation of T/R assembly is A'maxdB, the attenuation step size is Δ AdB;
step (2), dividing the whole array area corresponding to the amplitude weighted value Amp into amplitude weighted sub-areas theta corresponding to S multiplied by T/R components according to the arrangement form of the T/R componentsst(ii) a The amplitude weighted sub-region ΘstCorresponding amplitude weighted value is AmpstI.e. the amplitude-weighted sub-region theta of row s and column tstWeighted value Amp ofstThe amplitude weighting value of each rf channel can be expressed as:
Figure FDA0003117345900000011
wherein P1, Q;
Figure FDA0003117345900000012
is AmpstThe amplitude weighted value of the p-th row and the q-th column;
step (3), determining a 0 attenuation area according to the array amplitude weighted value Amp, determining and sequencing T/R components which can be arranged in the 0 attenuation area, and finally performing T/R component arrangement on the 0 attenuation area;
step (4), counting the minimum value of the gain of the distributed T/R assembly radio frequency channel at 0dB attenuation 0 DEG phase shift corresponding to the position of 0dB attenuation with the amplitude weighted value of 0dB attenuation in the 0 attenuation region
Figure FDA0003117345900000013
Step (5), determining a maximum attenuation area and an interchange area according to the array amplitude weighted value Amp, determining and sequencing T/R components which can be arranged in the maximum attenuation area, and finally arranging the T/R components in the maximum attenuation area;
step (6), dividing the interchange area into a small attenuation interchange area, a large attenuation interchange area and a full interchange area according to Amp, T/R assemblies arranged in the 0 attenuation area, T/R assemblies arranged in the maximum attenuation area and the gain characteristics of the rest T/R assemblies;
step (7), arranging the T/R components of the small attenuation interchange area, the large attenuation interchange area and the full interchange area according to the principle that the internal heat consumption is smaller than the external heat consumption;
and (8) counting the total thickness of each column of T/R assemblies along the extrusion direction, and adjusting the positions of the partial T/R assemblies in different columns according to the principle that the total thickness of each column is consistent.
2. The method for arranging the high-density integrated active phased array T/R component based on the mechanical, electrical and thermal methods as claimed in claim 1, wherein the optimized amplitude weighted value of the whole array determined in the step (1) is Amp, comprising the following steps:
(2a) if the phased array is a receiving phased array, the amplitude weighted value Amp is a receiving amplitude attenuation value Amp obtained by optimizationRI.e. Amp ═ AmpR
(2b) If the phased array is a transmitting phased array, the amplitude weighted value Amp is an optimized transmitting amplitude attenuation value AmpTI.e. Amp ═ AmpT
(2c) If the phased array is a transmitting-receiving shared phased array, the amplitude weighted value Amp is in a dynamic range D according to the received amplitude weighted valueRAnd transmit amplitude weight dynamic range DTIs determined, i.e. is
Figure FDA0003117345900000021
3. The method for arranging the high-density integrated active phased array T/R component based on the mechanical, electrical and thermal characteristics as claimed in claim 1, wherein the step (3) of determining the 0 attenuation region according to the array amplitude weighting value Amp comprises the following steps:
(3a) setting 0 attenuation region judgment threshold
Figure FDA0003117345900000022
(3b) For the amplitude weighted sub-region ΘstIf P ∈ { 1.,. P } and Q ∈ { 1.,. Q } satisfy the following expression
Figure FDA0003117345900000023
The amplitude weighted sub-region ΘstBelong to the 0 attenuation region; all S1, S, T1, T are judged to obtain the 0 attenuation region of the whole array.
4. The method for arranging the T/R components of the high-density integrated active phased array based on the mechanical, electrical and thermal characteristics as claimed in claim 1, wherein the step (3) of determining the T/R components capable of being arranged in the 0 attenuation region and sequencing the T/R components in the 0 attenuation region comprises the following steps:
(4a) setting 0 attenuation region T/R component judgment threshold
Figure FDA0003117345900000024
(4b) Amplitude weighted sub-region Θ for the 0 attenuation regionstWeighted value Amp ofstWhere all channels are 0dB attenuation, all T/R components can be used for this 0 attenuation magnitude weighted sub-region ΘstArranging;
amplitude weighted sub-region Θ for the 0 attenuation regionstWeighted value Amp ofstThe condition that the channel attenuation value is not 0dB exists in the memory, when AmpstAll the 0dB attenuation positions in the radio frequency channel of the T/R component correspond to the minimum value of the gain when the 0dB attenuation is 0 DEG phase shift
Figure FDA0003117345900000025
Subtracting AmpstAll non-0 dB attenuation positions in the T/R component correspond to the gains of the radio frequency channels of the T/R component when the 0dB attenuation is 0 DEG phase shift
Figure FDA0003117345900000026
Are all less than or equal to
Figure FDA0003117345900000027
And AmpstAmplitude weighting values of all non-0 dB fading channels
Figure FDA0003117345900000028
(wherein
Figure FDA0003117345900000029
Figure FDA0003117345900000031
Sum, i.e. when Q satisfies the following equation, P1
Figure FDA0003117345900000032
The T/R component can be used for the 0 attenuation magnitude weighting sub-region ΘstArranging;
when able to be arranged in the 0 attenuation amplitude weighting sub-region ΘstNumber of T/R components
Figure FDA0003117345900000033
Number of amplitude weighted sub-regions smaller than 0 attenuation region
Figure FDA0003117345900000034
Then choose in T/R component that cannot be used for 0 attenuation amplitude weighted sub-region
Figure FDA0003117345900000035
Individual weighted value AmpstMinimum value of gain at 0dB attenuation 0 degree phase shift of all 0dB attenuation positions
Figure FDA0003117345900000036
The largest T/R component is used for arranging the 0 attenuation amplitude weighting subareas; statistics for the 0 attenuation amplitude weighting sub-region ΘstAll T/R components arranged at a weighted value AmpstThe minimum value of the gain of all T/R radio frequency channels at the corresponding position of 0dB when the phase is shifted by 0 DEG under the attenuation of 0dB
Figure FDA0003117345900000037
And the T/R module heat loss PV and according to the increaseMinimum value of benefit
Figure FDA0003117345900000038
Sequencing the T/R components to obtain the sub-region thetastCorresponding T/R arrangement sequence
Figure FDA0003117345900000039
(4c) For phased arrays of different transmit-receive operating states, G in (4b)0
Figure FDA00031173459000000310
And PV is selected according to the following method: if the phased array is a receive phased array, G0And
Figure FDA00031173459000000311
the gain is the receiving state, and PV is the receiving heat consumption of the R component; if the phased array is a transmit phased array, G0And
Figure FDA00031173459000000312
the gain is the gain at the emission working point, and PV is the emission heat consumption of the T/R assembly; if the phased array is a transmit-receive shared phased array, when DR≥DTWhen, G0And
Figure FDA00031173459000000313
are all receiving state gains, otherwise G0And
Figure FDA00031173459000000314
the gains are all at the transmit operating point; and the heat consumption PV of the transceiving shared phased array is the emission heat consumption of the T/R assembly.
5. The method for arranging high-density integrated active phased array T/R components based on mechanical, electrical and thermal according to claim 1, wherein the step of arranging the T/R components in the attenuation 0 region in the step (3) is as follows:
(5a) counting each amplitude weighting sub-region theta in 0 attenuation regionstWeighted value AmpstTotal number of channels attenuated by 0dB
Figure FDA00031173459000000315
When different amplitude weights the sub-region ΘstIs/are as follows
Figure FDA00031173459000000316
When not identical, the amplitude weighting sub-region ΘstIs arranged in the order of
Figure FDA00031173459000000317
Arranging in the order from big to small; when the amplitude weights the sub-region ΘstIs/are as follows
Figure FDA00031173459000000318
If, as such, the sub-regions do not have channels with non-0 dB attenuation, i.e.
Figure FDA00031173459000000319
When the sub-regions are the same and equal to P multiplied by Q, the sequence of arrangement among the sub-regions is randomly selected; if the sub-regions have channels with attenuation which is not 0dB, counting weighted values of all the channels with attenuation which is not 0dB in the sub-regions, and counting the number of the channels with the same attenuation according to the sequence of the attenuation from small to large, namely, firstly counting the number of the channels with the minimum attenuation in the attenuation which is not 0dB, and preferentially arranging the channels with the maximum attenuation if the number of the channels with the minimum attenuation is larger than the attenuation if the number of the channels with the minimum attenuation is still the same, and preferentially arranging the channels with the same number of the channels until the maximum attenuation A in the sub-regions is countedmaxdB, if the channel number of the maximum attenuation is the same, the weighting sub-region theta isstThe sequence of arrangement is randomly selected;
(5b) respective amplitude weighted sub-regions Θ within the 0dB attenuation region determined in accordance with (5a) abovestThe arrangement order of the/R components is arranged, and the amplitude weighted sub-region theta is arrangedstAccording to the T/R module
Figure FDA0003117345900000041
The determined T/R components are sequentially selected to be arranged;
(5c) when different amplitude weights the sub-region ΘstIs/are as follows
Figure FDA0003117345900000042
All equal to PxQ, i.e. different sub-regions thetastWhen the weighted values of all the channels are attenuated by 0dB, the position of the T/R components is adjusted according to the principle that the T/R components with high heat consumption are far away from the geometric center of the array surface and the T/R components with low heat consumption are close to the center after the T/R components are arranged according to the method; firstly, counting the heat consumption PV of 0dB attenuation of weighted values of all channels of the T/R assembly before position interchangestAnd the corresponding sub-region ΘstIs a distance D between the geometric center of the array surface and the geometric center of the array surfacest(ii) a According to the distance DstPV (photovoltaic cell) in the T/R assembly which is arranged from large to small and needs to be subjected to position adjustment according to heat consumptionstSequentially selecting T/R components from large to small for position interchange; if there are multiple sub-regions ΘstD of (A)stUnder the same condition, selecting the T/R assemblies with the minimum heat consumption in the same number with the sub-regions from the rest T/R assemblies to be subjected to position adjustment, and randomly arranging the T/R assemblies in the sub-regions thetastIn (1).
6. The method as claimed in claim 1, wherein the minimum value of the gain of 0dB attenuation 0 ° phase shift of the rf channel of the arranged T/R module corresponding to the position with 0dB attenuation of all amplitude weighting values in the attenuation 0 region counted in step (4) is calculated
Figure FDA0003117345900000043
The steps are as follows:
for all S1, 1., S, T1., T and P1., P, Q1., Q,
Figure FDA0003117345900000044
the minimum value of gain of 0dB attenuation 0 DEG phase shift of distributed T/R assembly radio frequency channels corresponding to all positions with 0dB attenuation amplitude weighted value in 0 attenuation region
Figure FDA0003117345900000045
If the phased array is a receive phased array, then
Figure FDA0003117345900000046
Is the receive gain; if the phased array is a transmit phased array, then
Figure FDA0003117345900000047
Is the gain at the transmit operating point; if the phased array is a transmit-receive shared phased array, when DR≥DTWhen the temperature of the water is higher than the set temperature,
Figure FDA0003117345900000048
gain for receiving state, otherwise
Figure FDA0003117345900000049
Is the gain at the transmit operating point.
7. The method for arranging the high-density integrated active phased array T/R assembly based on the mechanical, electrical and thermal characteristics as claimed in claim 1, wherein the step (5) of determining the maximum attenuation region and the interchange region comprises the following steps:
(7a) counting all radio frequency channels of all residual T/R components at maximum gear A 'of attenuator'maxMaximum in dB-attenuated 0 DEG phase-shifted gain
Figure FDA00031173459000000410
Wherein, if the phased array is a receive phased array,
Figure FDA00031173459000000411
is the receive gain; if the phased array is a transmit phased array,
Figure FDA00031173459000000412
gain at the transmit operating point; if the phased array is a transmit-receive shared phased array, when DR≥DTWhen the temperature of the water is higher than the set temperature,
Figure FDA00031173459000000413
to receive gain, otherwise
Figure FDA00031173459000000414
Is the gain at the transmit operating point;
(7b) setting a maximum attenuation region judgment threshold
Figure FDA00031173459000000415
(7c) For the amplitude weighted sub-region ΘstWhen P ∈ {1, …, P } and Q ∈ {1, …, Q } satisfy the following expression
Figure FDA00031173459000000416
The amplitude weighted sub-region ΘstBelongs to the maximum attenuation region; judging all S1, T1, S, T, and T to obtain the maximum attenuation area of the whole array;
(7d) amplitude-weighted sub-region Θ not belonging to the 0 attenuation region nor to the maximum attenuation regionstBelonging to a interchange area;
(7e) without the amplitude weighting sub-region ΘstBelonging to the maximum attenuation region, the remaining regions are all interchange regions.
8. The method for arranging high-density integrated active phased array T/R components based on mechanical, electrical and thermal according to claim 6, wherein the step of determining and ordering T/R components capable of being arranged in the maximum attenuation region in the step (5) is as follows:
(8a) setting a maximum attenuation region T/R component judgment threshold
Figure FDA0003117345900000051
(8b) Amplitude weighting sub-region Θ for the region of maximum attenuationstWeighted value Amp ofstAll channels therein are AmaxdB attenuation, when a certain attenuation level alpha epsilon {0, delta Amax,...A′maxMake
Figure FDA0003117345900000052
And Amp corresponding to T/R componentstAll of A inmaxMaximum value of gain of dB attenuation channel when alpha dB attenuates 0 DEG phase shift
Figure FDA0003117345900000053
(wherein P1, Q) satisfies the following equation
Figure FDA0003117345900000054
The T/R component can be used for the maximum attenuation magnitude weighting sub-region ΘstArranging;
amplitude weighting sub-region Θ for the region of maximum attenuationstWeighted value Amp ofstThe channel attenuation value is different from AmaxdB, when there is a certain attenuation level α ∈ {0, Δ amax,...A′maxMake
Figure FDA0003117345900000055
Amp corresponding to T/R componentstAll of A inmaxMaximum value of gain of dB attenuation channel when alpha dB attenuates 0 DEG phase shift
Figure FDA0003117345900000056
All not AmaxdB attenuation channel at maximum gear A 'of T/R assembly attenuator'maxGain in dB attenuation at 0 ° phase shift
Figure FDA0003117345900000057
AmaxAnd all channel weights not subject to maximum attenuation
Figure FDA0003117345900000058
(wherein P1, Q) satisfies the following equation
Figure FDA0003117345900000059
The T/R component can be used for the maximum attenuation magnitude weighting sub-region ΘstArranging;
when able to be arranged in the maximum attenuation amplitude weighting sub-region ΘstNumber of T/R components
Figure FDA00031173459000000510
Less than the maximum attenuation region ΘstNumber of amplitude weighted sub-regions of
Figure FDA00031173459000000511
Then, the maximum attenuation amplitude weighting sub-region Θ cannot be arranged in the maximum attenuation amplitude weighting sub-regionstIn the T/R component of
Figure FDA00031173459000000512
Individual weighted value AmpstIn AmaxdB attenuation channel at A'maxMaximum gain in dB attenuation at 0 degree phase shift
Figure FDA00031173459000000513
The smallest T/R component is used for the arrangement of the maximum attenuation amplitude weighting subarea; statistics for the maximum attenuation amplitude weighting sub-region ΘstAll T/R components arranged at a weighted value AmpstIs AmaxRadio frequency channel at dB corresponding position of A'maxdB attenuation 0 degree shiftMaximum value of gain of phase time
Figure FDA0003117345900000061
And the T/R module heat loss PV, and according to the maximum value of the gain
Figure FDA0003117345900000062
Sequencing the T/R components to obtain the sub-region thetastCorresponding T/R arrangement sequence
Figure FDA0003117345900000063
(8c) If the phased array is a receiving phased array, the gains in (8b) are all receiving gains, and the heat consumption is the receiving heat consumption; if the phased array is a transmitting phased array, the gains in the step (8b) are the gains at the transmitting working point, and the heat consumption is the transmitting heat consumption; if the phased array is a transmit-receive shared phased array, when DR≥DTIf so, the gains in (8b) are all receiving gains, otherwise, the gains in (8b) are all gains at the transmitting working point; the heat consumption of the transmitting and receiving shared phased array is the transmitting heat consumption.
9. The method for arranging the T/R components of the high-density integrated active phased array based on the mechanical, electrical and thermal aspects of claim 8, wherein the step of arranging the T/R components in the maximum attenuation region in the step (5) comprises the following steps:
(9a) counting each amplitude weighting sub-region theta in the maximum attenuation regionstWeighted value AmpstIs AmaxTotal number of channels in dB attenuation
Figure FDA0003117345900000064
When different amplitude weights the sub-region ΘstIs/are as follows
Figure FDA0003117345900000065
When not identical, the amplitude weighting sub-region ΘstIs arranged in the order of
Figure FDA0003117345900000066
Arranging in the order from big to small; when the amplitude weights the sub-region ΘstIs/are as follows
Figure FDA0003117345900000067
If, as such, these subregions do not have amaxChannels with dB attenuation, i.e.
Figure FDA0003117345900000068
When the sub-regions are the same and equal to P multiplied by Q, the sequence of arrangement among the sub-regions is randomly selected; if these subregions have a value other than AmaxdB-attenuated channels, then all non-A in these sub-regions are countedmaxThe weighted value of dB attenuation channels is counted according to the order of the attenuation from large to small, that is, the number of channels with the same attenuation is counted firstlymaxThe channel number with the maximum attenuation in dB attenuation is distributed preferentially if the channel number of the attenuation is more, the channel number with the maximum attenuation smaller than the attenuation is counted if the channel number is still the same, the channel number with the same channel number is distributed preferentially if the channel number is more, and the process is repeated until the minimum attenuation in the sub-regions is counted, if the channel number of the minimum attenuation is still the same, the weighted sub-region theta is obtainedstThe sequence of arrangement is randomly selected;
(9b) each amplitude weighted sub-region Θ within the region of maximum attenuation determined in accordance with (9a) abovestThe arrangement order of the sub-regions is arranged in the range weighting sub-region thetastAccording to the T/R module
Figure FDA0003117345900000069
Selecting T/R for arrangement according to the determined arrangement sequence of T/R;
(9c) when different amplitude weights the sub-region ΘstIs/are as follows
Figure FDA00031173459000000610
All equal to PxQ, i.e. different sub-regions thetastThe weighted value of all channels is maximum AmaxIn dB attenuation, according to the methodAfter the T/R components are arranged, position adjustment is carried out according to the principle that the T/R components with high heat consumption are far away from the geometric center of the array surface and the T/R components with low heat consumption are close to the center; namely, firstly, the heat consumption PV of the T/R assembly before position interchange is countedstAnd the corresponding sub-region ΘstIs a distance D between the geometric center of the array surface and the geometric center of the array surfacest(ii) a According to the distance DstIn the T/R module which is arranged and is to be subjected to position adjustment according to heat consumption PV in the order from big to smallstSequentially selecting T/R components from large to small for position interchange; if a plurality of sub-regions ΘstD of (A)stIf the T/R components are the same as the sub-regions, selecting the T/R components with the minimum heat consumption in the rest T/R components to be subjected to position adjustment, and randomly arranging the T/R components in the sub-regions thetastIn (1).
10. The method of claim 1, wherein the step (6) of dividing the interchange area into the small attenuation interchange area, the large attenuation interchange area and the full interchange area comprises the steps of:
(10a) counting the minimum value of the gain of all radio frequency channels of all the rest T/R components at 0dB attenuation 0 DEG phase shift
Figure FDA0003117345900000071
All radio frequency channels of all remaining T/R components are at attenuator maximum gear A'maxMaximum in dB-attenuated 0 DEG phase-shifted gain
Figure FDA0003117345900000072
All radio frequency channels of all remaining T/R components are at attenuator maximum gear A'maxMinimum in dB-attenuated 0 DEG phase-shifted gain
Figure FDA0003117345900000073
Wherein if the phased array is a receive phased array, the gains in (10a) are all receive gains; if the phased array is a transmit phased array, then the gains in (10a) are all transmitGain at the operating point; if the phased array is a transmit-receive shared phased array, when DR≥DTIf so, the gains in (10a) are all receiving gains, otherwise, the gains in (10a) are all gains at a transmitting working point;
(10b) setting judgment threshold of small attenuation interchange area
Figure FDA0003117345900000074
When there is P e { 1.,. P } and Q e { 1.,. Q } in the interchange area satisfy
Figure FDA0003117345900000075
Then the amplitude weighted sub-region ΘstBelonging to a small attenuation interchange area;
(10c) setting a judgment threshold value of a large attenuation interchange area
Figure FDA0003117345900000076
When there is P e { 1.,. P } and Q e { 1.,. Q } in the interchange area satisfy
Figure FDA0003117345900000077
Then the amplitude weighted sub-region ΘstBelonging to a large attenuation interchange area;
(10d) for Θ in the exchange region that neither belongs to the small nor the large attenuation exchange regionstIs a fully interchanged region.
11. The method of claim 6, wherein the step (7) of arranging the remaining T/R components to be arranged in the small attenuation interchange area, the large attenuation interchange area and the full interchange area comprises the steps of:
(11a) counting all sub-regions to be arranged thetastIs a distance D between the geometric center of the array surface and the geometric center of the array surfacest
(11b) According to DstThe T/R components are arranged in sequence from small to large by using the following method:
setting up a small attenuationInterchangeable region T/R component judgment threshold
Figure FDA0003117345900000078
And large attenuation interchange area T/R component judgment threshold
Figure FDA0003117345900000079
② if the amplitude weighting sub-region theta to be arrangedstBelonging to the small attenuation interchange area, the amplitude weighting sub-area ΘstWeighted value of AmpstFor all P e { 1.,. P } and Q e { 1.,. Q }, there is a
Figure FDA0003117345900000081
So that
Figure FDA0003117345900000082
Subtracting the weighted value Amp corresponding to the T/R componentstIn
Figure FDA0003117345900000083
Gain in dB attenuation at 0 ° phase shift
Figure FDA0003117345900000084
Are all less than or equal to the amplitude weighted value of the channel
Figure FDA0003117345900000085
And
Figure FDA0003117345900000086
when the difference of (a) is greater than or equal to 1, i.e., when Q satisfies the following equation for all P, i.e., P, Q, i.e., 1
Figure FDA0003117345900000087
The T/R component can be used for the small attenuation swap area ΘstArranging; when it can be arranged in the small attenuation interchange areaDomain thetastNumber of T/R components
Figure FDA0003117345900000088
Number of amplitude weighted sub-regions smaller than the small attenuation interchange region
Figure FDA0003117345900000089
Then choose in T/R module that can' T be used for small attenuation interchange area
Figure FDA00031173459000000810
Individual weighted value AmpstMinimum value of gain at 0dB attenuation 0 degree phase shift of all 0dB attenuation positions
Figure FDA00031173459000000811
The largest T/R component is used for the arrangement of the small attenuation interchange area; counting the heat consumption PV of all T/R modules capable of being used in the regionstSelective heat loss PVstThe lowest is arranged in the amplitude-weighted sub-region Θst
Third if the amplitude weighting sub-region theta to be arrangedstBelonging to the region of large attenuation interchange, the amplitude-weighted sub-region ΘstWeighted value of AmpstFor all P e { 1.,. P } and Q e { 1.,. Q }, there is a
Figure FDA00031173459000000812
So that
Figure FDA00031173459000000813
Subtracting the weighted value Amp corresponding to the T/R componentstIn
Figure FDA00031173459000000814
Attenuating gain at 0 ° phase shift
Figure FDA00031173459000000815
Are all larger than or equal to the amplitude weighted value of the channel
Figure FDA00031173459000000816
And
Figure FDA00031173459000000817
when the sum of (a) and (b) is equal to 1, Q satisfies the following equation
Figure FDA00031173459000000818
The T/R component can be used for the large attenuation interchange area ΘstArranging; when able to be arranged in the small attenuation interchange area thetastNumber of T/R components
Figure FDA00031173459000000819
Number of amplitude weighted sub-regions smaller than the large attenuation interchange region
Figure FDA00031173459000000820
Then choose in T/R module that can' T be used for large attenuation interchange area
Figure FDA00031173459000000821
Individual weighted value AmpstIn AmaxdB attenuation channel at A'maxMaximum gain in dB attenuation at 0 degree phase shift
Figure FDA00031173459000000822
The smallest T/R-module is used for the arrangement of the large attenuation interchange area; counting the heat consumption PV of all T/R modules capable of being used in the regionstSelective heat loss PVstThe lowest is arranged in the amplitude-weighted sub-region Θst
Fourthly, if the amplitude weighting sub-region theta to be arrangedstBelonging to a full interchange area, all the rest T/R components can be arranged in the area; counting the heat consumption PV of all T/R modules capable of being used in the regionstSelective heat loss PVstThe lowest is arranged in the amplitude-weighted sub-region Θst
(11c) If there are a plurality of sub-regions DstIf the sub-regions are equal to each other, the arrangement sequence among the sub-regions is randomly determined, and the arrangement method of each sub-region is the same as that in the step (11 b);
if the phased array is a receiving phased array, the gains in (11b) are all receiving gains, and the heat consumption is the receiving heat consumption; if the phased array is a transmitting phased array, the gains in the step (11b) are the gains at the transmitting working point, and the heat consumption is the transmitting heat consumption; if the phased array is a transmit-receive shared phased array, when DR≥DTIf so, the gains in (11b) are all receiving gains, otherwise, the gains in (11b) are all gains at the transmitting working point; the heat consumption of the transmitting and receiving common phased array is the transmitting heat consumption;
(11d) and (4) repeating the steps (11b) to (11c) until all the T/R assemblies are arranged.
12. The method for arranging the T/R components of the high-density integrated active phased array based on the mechanical, electrical and thermal characteristics as claimed in claim 1, wherein the step (8) of adjusting the position of the T/R components according to the total thickness of each column of the T/R components along the extrusion direction comprises the following steps:
(12a) let the total number of T columns and the total number of S rows of T/R modules in the extrusion direction, and the average thickness of all T/R modules (S multiplied by T) be
Figure FDA0003117345900000091
The desired total thickness of each column of T/R assemblies is then
Figure FDA0003117345900000092
(12b) Count the total thickness { H ] of all column T/R components1,H2,...,Ht,...,HTAnd setting a total thickness tolerance threshold value epsilon of each column in the extrusion directionH
(12c) Calculating the difference between the actual total thickness and the expected total thickness of each column of T/R { Delta H1,ΔH2,...,ΔHt,...,ΔHTThat is to say
ΔHt=Ht-HExpWherein (T ═ 1, 2.. T)
(12d) Find { Δ H1,ΔH2,...,ΔHt,...,ΔHTMaximum value of
Figure FDA0003117345900000093
Namely, it is
Figure FDA0003117345900000094
Figure FDA0003117345900000095
And minimum value
Figure FDA0003117345900000096
Namely, it is
Figure FDA0003117345900000097
Wherein the maximum value is located at column t1 and the minimum value is located at column t2, such that
Figure FDA0003117345900000098
(12e) Listing all T/R Components in column T1
Figure FDA0003117345900000099
And all T/R modules in column T2
Figure FDA00031173459000000910
Calculating the difference deltaH 'between the actual total thickness of the T1 th column and the expected total thickness of the T2 th column and any T/R assembly of the T1 th column and any T/R assembly of the T2 th column after being interchanged't1And Δ H't2Let us order
Figure FDA00031173459000000911
Δ H 'after interchanging the e th T/R Assembly from column T1 and the f th T/R Assembly from column T2't1And Δ H't2Maximum value of absolute value, i.e.
Figure FDA00031173459000000912
Marking according to the type of the region where the T/R component is positioned
Figure FDA00031173459000000913
And
Figure FDA00031173459000000914
interchange type mark LefTraversing all e ∈ {1, 2.,. P } and f ∈ {1, 2.,. P } combinations to get about
Figure FDA00031173459000000915
Matrix MatirxΔHAnd about LefMatrix MatirxL(ii) a In MatirxLIn finding interchange type flag minimum value Lmin=min(MatirxL),
Figure FDA00031173459000000916
Figure FDA0003117345900000101
(12f) Finding interchange type flag LminCorresponding matrix Matirx of the positionΔHMinimum value of Medium element Δ H'minI.e. Δ H'min=min(MatirxΔH|Lmin) Where e ∈ {1,2, …, P }, f ∈ {1,2, …, P }, and the minimum value Δ H 'is recorded'minCorresponding interchangeable T/R assembly
Figure FDA0003117345900000102
(12g) When Δ Hmax>ΔH′minWhen, to
Figure FDA0003117345900000103
And
Figure FDA0003117345900000104
the components are subjected to position interchange, so that the position of the components is exchanged by delta Hmax=ΔH′minSetting a T/R component interchange blacklist as a blank list (0, 0);
when Δ Hmax≤ΔH′minWhen the position is not exchanged, the T/R component exchange blacklist is set as (B)1,B2) In which B is1=t1,B2T2, i.e. B-th of T/R module1Column and B2Column inhibit swapping; at this time, { Δ H is found out of the list where the T/R device interchange blacklists1,ΔH2,…,ΔHt,…,ΔHTMaximum value of
Figure FDA0003117345900000105
Column t1, and minimum value
Figure FDA0003117345900000106
The t2 th columns, i.e., t1 and t2, satisfy the following formula
t1≠B1And t2 ≠ B2
Calculate the time of
Figure FDA0003117345900000107
Subsequently repeating the above (12e) to step (12g) until { H }1,H2,...,HTThere are no interchangeable T/R components in any two columns;
(12h) if there is no interchangeable T/R component in any two columns, let L in (12f)min=Lmin+1, followed by repeating steps (12f) to (12g) above until Lmin=15;
(12i) Repeating the steps (12c) to (12H) until Δ Hmax≤εHOr { H1,H2,…,HTAnd when any two columns in the T/R module are not interchanged, stopping the position adjustment of the T/R module.
CN201911096411.4A 2019-11-11 2019-11-11 High-density integrated active phased array T/R assembly arrangement method based on machine, electricity and heat Active CN111009729B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911096411.4A CN111009729B (en) 2019-11-11 2019-11-11 High-density integrated active phased array T/R assembly arrangement method based on machine, electricity and heat

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911096411.4A CN111009729B (en) 2019-11-11 2019-11-11 High-density integrated active phased array T/R assembly arrangement method based on machine, electricity and heat

Publications (2)

Publication Number Publication Date
CN111009729A CN111009729A (en) 2020-04-14
CN111009729B true CN111009729B (en) 2021-08-06

Family

ID=70112282

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911096411.4A Active CN111009729B (en) 2019-11-11 2019-11-11 High-density integrated active phased array T/R assembly arrangement method based on machine, electricity and heat

Country Status (1)

Country Link
CN (1) CN111009729B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104052515A (en) * 2014-05-13 2014-09-17 成都雷电微力科技有限公司 High-integrity TR radio frequency module
WO2016053501A1 (en) * 2014-10-03 2016-04-07 Raytheon Company Transmit/receive daughter card wth integral circulator
CN106772256A (en) * 2016-12-20 2017-05-31 中国航空工业集团公司雷华电子技术研究所 A kind of Connectors for Active Phased Array Radar antenna Antenna Subarray Division
BR112017015658A2 (en) * 2015-01-23 2018-03-20 Huawei Tech Co Ltd phase control for antenna network
CN108387878A (en) * 2018-06-01 2018-08-10 中国人民解放军陆军工程大学石家庄校区 A kind of phased-array radar TR components automatic testing equipment and method

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7671696B1 (en) * 2006-09-21 2010-03-02 Raytheon Company Radio frequency interconnect circuits and techniques
US10419062B2 (en) * 2016-10-12 2019-09-17 Massachusetts Institute Of Technology Simultaneous transmit and receive with digital phased arrays
CN107275806B (en) * 2017-05-19 2019-11-12 北京空间飞行器总体设计部 A kind of phased array antenna front method of weighting
CN109980367B (en) * 2019-03-28 2020-12-29 中国人民解放军陆军工程大学 Array antenna with rapid self-repairing capability and self-repairing method thereof

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104052515A (en) * 2014-05-13 2014-09-17 成都雷电微力科技有限公司 High-integrity TR radio frequency module
WO2016053501A1 (en) * 2014-10-03 2016-04-07 Raytheon Company Transmit/receive daughter card wth integral circulator
BR112017015658A2 (en) * 2015-01-23 2018-03-20 Huawei Tech Co Ltd phase control for antenna network
CN106772256A (en) * 2016-12-20 2017-05-31 中国航空工业集团公司雷华电子技术研究所 A kind of Connectors for Active Phased Array Radar antenna Antenna Subarray Division
CN108387878A (en) * 2018-06-01 2018-08-10 中国人民解放军陆军工程大学石家庄校区 A kind of phased-array radar TR components automatic testing equipment and method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Superconducting Sub-array Module as T/R Module for X-band Active Phased Array Antenna;Kenta Iijima;《2015 IEEE Radar Conference》;20150625;全文 *
一种Ka频段瓦片式TR组件子阵集成方案;赵青;《电讯技术》;20120720;全文 *
基于随机阵列的相控阵T/R组件排布方法;王罗胜斌;《雷达科学与技术》;20151015;全文 *

Also Published As

Publication number Publication date
CN111009729A (en) 2020-04-14

Similar Documents

Publication Publication Date Title
CN1748340B (en) Feed network for simultaneous generation of narrow and wide beams with a rotational-symmetric antenna
CN104270820B (en) United vertical beam control and power distribution method in the extensive mimo systems of 3D
CN107403991A (en) Super ultra wide band AESA system and method
CN101359947A (en) Broadcast wave bean shaping method and apparatus for multi-antenna array system
CN103167512B (en) Antenna for base station angle of declination defining method, device and base station equipment
CN1129237C (en) Digital wave beam forming method and module in radio communication system and its array receiver
CN101938837A (en) Inter-cell interference suppression method and device
US6781560B2 (en) Phased array antenna including archimedean spiral element array and related methods
CN112733334B (en) Multi-band common-aperture array antenna sparse optimization method and system
CN205863410U (en) A kind of multifrequency Shared aperture is combined phased array antenna structure
CN105842666B (en) Radar Subarray partition optimization method based on difference algorithm
CN111009729B (en) High-density integrated active phased array T/R assembly arrangement method based on machine, electricity and heat
CN109800474A (en) A kind of array antenna energy concocting method based on genetic algorithm
CN103715503B (en) A kind of multiple sector antenna and communication system
CN103812547A (en) Communication device and method and device for beam forming
CN105896102A (en) Subarray antenna structure suitable for millimeter wave LOS MIMO and design method
US10581501B2 (en) Flexible analog architecture for sectorization
CN112736486A (en) Broadband RCS (radar cross section) reduced diffuse reflection structure and design method thereof
CN109379119B (en) Transmission method for carrying out SWIPT by utilizing vertical dimension for 3D large-scale MIMO system
CN103956585A (en) Phased-array antenna unit with flat-top directional diagram characteristics
CN115719015A (en) Design method of target radiation array structure
CN206148623U (en) Wide band omnidirectional antenna
CN113851854A (en) Wide-angle scanning lens terminal antenna and scanning mode thereof
CN110687527B (en) Method for improving directivity of wide-beam transmitting transducer and transducer
CN1160890C (en) Digital wave beam forming method and module in radio communication system and its array receiver

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