CN104750930B - A kind of determination method of ribbon parachute minimum closing in ratio - Google Patents

A kind of determination method of ribbon parachute minimum closing in ratio Download PDF

Info

Publication number
CN104750930B
CN104750930B CN201510145150.6A CN201510145150A CN104750930B CN 104750930 B CN104750930 B CN 104750930B CN 201510145150 A CN201510145150 A CN 201510145150A CN 104750930 B CN104750930 B CN 104750930B
Authority
CN
China
Prior art keywords
mrow
msub
mfrac
parachute
time
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
CN201510145150.6A
Other languages
Chinese (zh)
Other versions
CN104750930A (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.)
Beijing Institute of Space Research Mechanical and Electricity
Original Assignee
Beijing Institute of Space Research Mechanical and Electricity
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 Beijing Institute of Space Research Mechanical and Electricity filed Critical Beijing Institute of Space Research Mechanical and Electricity
Priority to CN201510145150.6A priority Critical patent/CN104750930B/en
Publication of CN104750930A publication Critical patent/CN104750930A/en
Application granted granted Critical
Publication of CN104750930B publication Critical patent/CN104750930B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

A kind of determination method of ribbon parachute minimum closing in ratio, first with the data point data of tape canopy apical pore position, is calculated the relative distance between multigroup measurement point on parachute formed symmetrical, and it is defined as to the apparent diameter of canopy apical pore；Then the Time frame variation coefficient and spatial variability coefficient and the determination of stability factor for being combined into canopy inflation of ribbon parachute are calculated respectively according to apparent diameter；Finally the stability of ribbon parachute and closing in are determined than situation according to the value of the determination of stability factor, determine the minimum closing in ratio of ribbon parachute.The present invention obtains the data point data of apical pore position according to the dynamic geometry configuration of canopy, and then judge and calculate minimum closing in ratio, the blank at present in terms of the minimum closing in of parachute is than theoretical research is filled up, and operability is stronger, brief and practical, it can be applied in engineering practice.

Description

A kind of determination method of ribbon parachute minimum closing in ratio
Technical field
The present invention relates to spacecraft parachute technology field, particularly a kind of determination side of ribbon parachute minimum closing in ratio Method.
Background technology
Parachute closing in technology is to reduce the opening shock most efficient method of big load body parachute-opening under high-speed condition.It is right For the large area ribbon parachute of the parachute-opening under high-speed condition, the multistage closing in technology for increasing ribbon parachute drag area step by step is reducing Load body opening shock, the Structural strength calls for reducing umbrella, the stability etc. that increases umbrella have remarkable result, and minimum receive Mouth is than being then the primary important evidence for closing up and designing.Close up than too small, canopy air inflow very little, is hardly formed the pneumatic outer of stabilization Shape, stable aerodynamic force can not be formed to load body；Close up than excessive, subsequent stages close up to narrow than optional scope, it is difficult to reaches To effective off-load and the purpose of optimization design.The determination of current minimum closing in ratio still lacks effective test data and theory analysis Achievement is as support, and to closing up, ratio is optimized so as to realize that the optimization to parachute global design is difficult to carry out.
The content of the invention
Present invention solves the technical problem that it is：Overcome the deficiencies in the prior art, there is provided one kind determines that ribbon parachute is minimum and received The method of mouth ratio, to determine that closing in ratios at different levels provide reference frame during engineer applied multistage closing in technology.
The present invention technical solution be：A kind of determination method of ribbon parachute minimum closing in ratio, comprises the following steps：
(1) according to ribbon parachute physical dimension, material, it is current close up than and given dynamic pressure establish ribbon parachute model after imitated Very, the data that ribbon parachute canopy different parts change over time are obtained, by two initial position messages on ribbon parachute apical pore circumference On the symmetrical position in ribbon parachute axis as one group of data point, A group data points are chosen altogether, and calculate in i-th group of data point The distance of two data points is Di(t), wherein, i=1,2,3 ... A, t be calculate 2 points apart from the time of and be limited；It is described Initial position message is that ribbon parachute canopy different parts are emulating positional information when starting；
(2) D is intercepted respectivelyi(t) the period postscript of stabilization process is Δ t on curvei, calculate A Δ tiArithmetic average It is worth and is designated as Δ t, calculates any instant air inlet equivalent diameter D in Δ tequ(t) it is equal with the air inlet equivalent diameter of stabilization process Value Dave：
Or
Wherein,For n in stabilization process period average Δ ttIndividual moment Dequ(t) value sum；
(3) the time meansquaredeviationσ of stabilization process is calculatedtWith Time frame variation coefficient Ct：
Wherein, ntNumber at the time of correspondence for data point in stabilization process period average Δ t,For n in stabilization process period average Δ tt(D corresponding to the individual momentave-Dequ(t))2Sum；
(4) the space meansquaredeviationσ that A group canopy different parts change over time data is calculatedp(t) it is
Then n in stabilization process period average Δ t is calculatedtThe number that individual moment A groups canopy different parts change over time According to mean space mean square deviationFor
Finally obtain spatial variability coefficient CpFor
(5) calculating canopy dimensional stability sex determination factor S is
S=C1·Ct+C2·Cp
And judged, wherein, C1、C2For weight factor；
(6) ifIt is then current to close up than the minimum closing in ratio for experiment ribbon parachute under current working And export, ifThen increase is current closes up than simultaneously repeat step (1)-step (5) until minimum receipts are calculated Mouth ratio, ifThen reduce current close up than simultaneously repeat step (1)-step (5) until minimum receipts are calculated Mouth ratio.
Δ t in described step (2)iComputational methods be：For the D of acquisitioni(t) it is, square along time shaft since t=0 To first cyclic curve section is searched, it is t to remember this section of curve initial time estimatestd=tstd0, time span estimate is Tf =Tf0, take the time average of this section of curve and be designated asTake second period curved section [tstd+Tf,tstd+2Tf] time it is equal Value, is designated asIf
Then think current tstdBe ribbon parachute enter steady-state period section initial time and section cycle steady-state period be Tf, And take current tstdTo Di(t) when a length of Δ t of curve end, otherwise takes Tf=1.05TfSaid process is repeated until meeting above-mentioned Rule of judgment, if T in calculating processfValue is more than 1.5 times of Tf0, then T is refetchedf=Tf0, and make tstd=tstd+Tf, continue to calculate directly To meeting above-mentioned Rule of judgment.
A >=4 in described step (1).
The present invention compared with prior art the advantages of be：
(1) present invention has filled up the blank in terms of the minimum closing in of parachute is than theoretical research at present, can be in tape More accurately and effectively controlled in umbrella design process to closing up, the closing in other umbellate forms also provides reference than design；
(2) stability and minimum closing in ratio of the dynamic geometry change of configuration of the invention according to canopy to Parachute Suit Inflating Judged and calculated, relevance is strong, the degree of accuracy is high；
(3) present invention is using the relative position of measurement point near canopy apical pore as analysis foundation, and operability is stronger, succinctly Practicality, it can be applied in engineering practice.
Brief description of the drawings
Fig. 1 is the inventive method flow chart.
Embodiment
The embodiment of the present invention is further described in detail below in conjunction with the accompanying drawings.The present invention first with Measurement point (data point) data of tape canopy apical pore position, are calculated between multigroup measurement point on parachute formed symmetrical Relative distance, and it is defined as to the apparent diameter D of canopy apical porei(t), i=1,2,3 ... A, t are to calculate two point distances Moment and to be limited；On this basis, the Time frame variation coefficient C of ribbon parachute is calculated respectivelytWith spatial variability coefficient CpAnd combine Into the determination of stability factor S of canopy inflation；Finally the stability of ribbon parachute and closing in are made than situation and sentenced according to S value It is fixed, the final minimum closing in ratio for determining ribbon parachute, after the minimum ratio that closes up is calculated, it can be closed up according to minimum and compare ribbon parachute It is designed, and obtains the band parachute jumping parameter such as priming system placement information.As shown in figure 1, a kind of ribbon parachute is minimum to close up ratio really Determine method, including step is as follows：
(1) according to the physical dimension of ribbon parachute, material parameter, close up ratio and given dynamic pressure utilize LS-DYNA softwares foundation Simulation model, calculate and obtain the data that ribbon parachute geometric shape changes over time；The space for choosing A group canopy different parts is sat Mark the data changed over time and (contain two data points in every group, it is desirable to which the initial position of the two data points is in ribbon parachute Axisymmetrical；All data point requirements are substantially distributed on canopy apical pore circumference)；Remember the distance between i-th group of symmetry number strong point For Di(t), wherein t is the measurement moment.And then obtain a series of Di(t) curve, D is rememberedi(t) it is the apparent diameter of canopy apical pore.
(2) D is interceptedi(t) stabilization time section on curve, and remember that this time is Δ ti；Calculate { Δ tiArithmetic mean of instantaneous value And Δ t is designated as, wherein Δ tiComputational methods be：For the D of acquisitioni(t) time history curve, along time shaft since t=0 Positive direction searches first cyclic curve section, and this section of curve initial time estimate is tstd=tstd0, time span estimate is Tf=Tf0, the time average of this section of curve is taken, is designated asTake [tstd+Tf,tstd+2Tf] curved section time average, be designated asIf
Then think current tstdBe ribbon parachute enter steady-state period section initial time and section cycle steady-state period be Tf, And take current tstdTo Di(t) when a length of Δ t of curve endi, otherwise take Tf=1.05TfAbove formula is computed repeatedly until in satisfaction Condition is stated, if T in calculating processfValue is more than 1.5 times of Tf0, then T is refetchedf=Tf0, and make tstd=tstd+Tf, continue calculate until Meet above-mentioned condition, obtain Δ tiValue.Calculate the air inlet equivalent diameter D at a certain moment in whole stabilization time sectionequ(t) and Air inlet equivalent diameter average D in whole stabilization time sectionave；
Consider that test data is usually time uniform discrete point, DaveIt can be estimated with following formula：
Wherein, ntFor Dequ(t) curve data point number of the curve in whole stabilization time section,For ntIt is individual D on curve data pointequ(t) value sum.
(3) D is directly utilizedequAnd D (t)aveCalculate the time meansquaredeviationσ of whole stabilization time sectiontAnd Time frame variation coefficient Ct；
Wherein, ntNumber at the time of correspondence for data point in stabilization process period average Δ t,For n in stabilization process period average Δ tt(D corresponding to the individual momentave-Dequ(t))2Sum.
(4) first with DequAnd D (t)i(t) the space meansquaredeviationσ of a certain moment A group data points is calculatedp(t), then calculate N in stabilization process period average Δ ttThe mean space for the data that individual moment A groups canopy different parts change over time is square DifferenceFinally obtain spatial variability coefficient Cp；
(5) canopy dimensional stability sex determination factor S is calculated：
S=C1·Ct+C2·Cp
Wherein, C1=0.6, C2=0.4 is weight factor.
(6) judged according to following table：
The content not being described in detail in description of the invention belongs to the known technology of those skilled in the art.

Claims (3)

  1. A kind of 1. determination method of ribbon parachute minimum closing in ratio, it is characterised in that comprise the following steps：
    (1) according to ribbon parachute physical dimension, material, it is current close up than and given dynamic pressure establish ribbon parachute model after emulated, Obtain the data that change over time of ribbon parachute canopy different parts, by two initial position messages on ribbon parachute apical pore circumference on A group data points are chosen altogether, and calculate in i-th group of data point two in the symmetrical position in ribbon parachute axis as one group of data point The distance of data point is Di(t), wherein, i=1,2,3 ... A, t be calculate 2 points apart from the time of and be limited；It is described initial Positional information is that ribbon parachute canopy different parts are emulating positional information when starting；
    (2) D is intercepted respectivelyi(t) the period postscript of stabilization process is Δ t on curvei, calculate A Δ tiArithmetic mean of instantaneous value simultaneously Δ t is designated as, calculates any instant air inlet equivalent diameter D in Δ tequ(t) and stabilization process air inlet equivalent diameter average Dave：
    <mrow> <msub> <mi>D</mi> <mrow> <mi>e</mi> <mi>q</mi> <mi>u</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>A</mi> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>A</mi> </munderover> <msub> <mi>D</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow>
    Or
    Wherein,For n in stabilization process period average Δ ttIndividual moment Dequ(t) value sum；
    (3) the time meansquaredeviationσ of stabilization process is calculatedtWith Time frame variation coefficient Ct：
    <mrow> <msub> <mi>&amp;sigma;</mi> <mi>t</mi> </msub> <mo>=</mo> <msqrt> <mfrac> <mrow> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>a</mi> <mi>l</mi> <mi>l</mi> </mrow> </munder> <msup> <mrow> <mo>(</mo> <msub> <mi>D</mi> <mrow> <mi>a</mi> <mi>v</mi> <mi>e</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>D</mi> <mrow> <mi>e</mi> <mi>q</mi> <mi>u</mi> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <msub> <mi>n</mi> <mi>t</mi> </msub> </mfrac> </msqrt> </mrow>
    <mrow> <msub> <mi>C</mi> <mi>t</mi> </msub> <mo>=</mo> <mfrac> <msub> <mi>&amp;sigma;</mi> <mi>t</mi> </msub> <msub> <mi>D</mi> <mrow> <mi>a</mi> <mi>v</mi> <mi>e</mi> </mrow> </msub> </mfrac> </mrow>
    Wherein, ntNumber at the time of correspondence for data point in stabilization process period average Δ t,For N in stabilization process period average Δ tt(D corresponding to the individual momentave-Dequ(t))2Sum；
    (4) the space meansquaredeviationσ that A group canopy different parts change over time data is calculatedp(t) it is
    <mrow> <msub> <mi>&amp;sigma;</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>A</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>D</mi> <mi>i</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>-</mo> <msub> <mi>D</mi> <mrow> <mi>e</mi> <mi>q</mi> <mi>u</mi> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <mi>A</mi> </mfrac> </msqrt> <mo>,</mo> </mrow>
    Then n in stabilization process period average Δ t is calculatedtThe data that individual moment A groups canopy different parts change over time are put down Equal space mean square deviationFor
    <mrow> <mover> <msub> <mi>&amp;sigma;</mi> <mi>p</mi> </msub> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mrow> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>a</mi> <mi>l</mi> <mi>l</mi> </mrow> </munder> <msub> <mi>&amp;sigma;</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <msub> <mi>n</mi> <mi>t</mi> </msub> </mfrac> <mo>,</mo> </mrow>
    Finally obtain spatial variability coefficient CpFor
    <mrow> <msub> <mi>C</mi> <mi>P</mi> </msub> <mo>=</mo> <mfrac> <mover> <msub> <mi>&amp;sigma;</mi> <mi>p</mi> </msub> <mo>&amp;OverBar;</mo> </mover> <msub> <mi>D</mi> <mrow> <mi>a</mi> <mi>v</mi> <mi>e</mi> </mrow> </msub> </mfrac> <mo>;</mo> </mrow>
    (5) calculating canopy dimensional stability sex determination factor S is
    S=C1·Ct+C2·Cp
    And judged, wherein, C1、C2For weight factor；
    (6) ifThen it is current close up than for tested under current working the minimum closing in of ribbon parachute than and it is defeated Go out, ifThen increase is current closes up than simultaneously repeat step (1)-step (5) up to minimum closing in ratio is calculated, IfThen reduce current close up than simultaneously repeat step (1)-step (5) until minimum closing in ratio is calculated.
  2. A kind of 2. determination method of ribbon parachute minimum closing in ratio according to claim 1, it is characterised in that：Described step (2) Δ t iniComputational methods be：For the D of acquisitioni(t), a cycle song is searched along time shaft positive direction since t=0 Line segment, it is t to remember this section of curve initial time estimatestd=tstd0, time span estimate is Tf=Tf0, take this section of curve Time average is simultaneously designated asTake second period curved section [tstd+Tf,tstd+2Tf] time average, be designated asIf
    <mrow> <mfrac> <mrow> <mo>|</mo> <mrow> <mover> <msub> <mi>F</mi> <mi>b</mi> </msub> <mo>&amp;OverBar;</mo> </mover> <mo>-</mo> <mover> <msub> <mi>F</mi> <mi>a</mi> </msub> <mo>&amp;OverBar;</mo> </mover> </mrow> <mo>|</mo> </mrow> <mover> <msub> <mi>F</mi> <mi>a</mi> </msub> <mo>&amp;OverBar;</mo> </mover> </mfrac> <mo>&amp;le;</mo> <mn>5</mn> <mi>%</mi> </mrow>
    Then think current tstdBe ribbon parachute enter steady-state period section initial time and section cycle steady-state period be Tf, and take and work as Preceding tstdTo Di(t) when a length of Δ t of curve endi, otherwise take Tf=1.05TfSaid process is repeated until meeting above-mentioned judgement Condition, if T in calculating processfValue is more than 1.5 times of Tf0, then T is refetchedf=Tf0, and make tstd=tstd+Tf, continue to calculate until full The above-mentioned Rule of judgment of foot.
  3. A kind of 3. determination method of ribbon parachute minimum closing in ratio according to claim 1 or 2, it is characterised in that：Described A >=4 in step (1).
CN201510145150.6A 2015-03-30 2015-03-30 A kind of determination method of ribbon parachute minimum closing in ratio Active CN104750930B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510145150.6A CN104750930B (en) 2015-03-30 2015-03-30 A kind of determination method of ribbon parachute minimum closing in ratio

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510145150.6A CN104750930B (en) 2015-03-30 2015-03-30 A kind of determination method of ribbon parachute minimum closing in ratio

Publications (2)

Publication Number Publication Date
CN104750930A CN104750930A (en) 2015-07-01
CN104750930B true CN104750930B (en) 2018-03-09

Family

ID=53590609

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510145150.6A Active CN104750930B (en) 2015-03-30 2015-03-30 A kind of determination method of ribbon parachute minimum closing in ratio

Country Status (1)

Country Link
CN (1) CN104750930B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111090908B (en) * 2019-11-13 2023-05-09 北京空间机电研究所 Group umbrella system design parameter selection method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103544053A (en) * 2013-10-17 2014-01-29 南京航空航天大学 Parachute aerodynamic performance analyzing method considering fabric breathability
CN104102766A (en) * 2014-06-24 2014-10-15 南京航空航天大学 Ergun theory-based prediction method for breathability of parachute cloth fabric

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103544053A (en) * 2013-10-17 2014-01-29 南京航空航天大学 Parachute aerodynamic performance analyzing method considering fabric breathability
CN104102766A (en) * 2014-06-24 2014-10-15 南京航空航天大学 Ergun theory-based prediction method for breathability of parachute cloth fabric

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"基于LS-DYNA的降落伞伞衣织物透气性参数仿真验证";贾贺等;《航天返回与遥感》;20090331;第30卷(第1期);15-20 *
"带牵顶伞的大面积环帆伞充气性能分析";尚小娟等;《航天返回与遥感》;20100831;第31卷(第4期);21-26 *
"环帆伞可变参数对设计偏差的灵敏度分析";甘和麟等;《航天返回与遥感》;20140430;第35卷(第2期);14-24 *

Also Published As

Publication number Publication date
CN104750930A (en) 2015-07-01

Similar Documents

Publication Publication Date Title
CN108108562A (en) A kind of Analytic modeling method of the Wind turbines wake flow based on Gaussian Profile
CN104070525B (en) For the method for space manipulator continuous trajectory tracking
Holden Establishment time of laminar separated flows
CN105527069A (en) Transonic wind tunnel test supporting device, installation method and application thereof
CN104836278B (en) A kind of many electricity aircraft electric loads start management method
CN104750930B (en) A kind of determination method of ribbon parachute minimum closing in ratio
CN105022919B (en) A kind of method for predicting circular metal thin plate low velocity impact dimple size
CN107944137A (en) The thermographic curve computing technique of hypersonic aircraft trajectory state multi- scenarios method
CN104344967A (en) Tread pattern groove pumping noise source recognition method
CN101864978A (en) Method for determining ventilation on-way resistance of curved tunnel
CN105353650B (en) Foundation temporarily rushes formula Asia across the super preset aperture model method of across the flow field pressure regulator valve in wind-tunnel Asia
CN205779470U (en) A kind of blunt trailing edge wind mill airfoil circulation control device
CN106347059A (en) Hub-driven electric car active suspension two-loop PID control method based on particle swarm optimization
CN106508020B (en) A kind of complex aircraft calculation of aerodynamic heating method that can be used for engineering design
CN107976908A (en) A kind of aircraft coupled dynamic stability characteristic analysis method
CN108038292A (en) A kind of efficient self-adapted method of sampling based on dual-proxy technology
CN109492234A (en) A kind of improved fluid structurecoupling interpolation method
CN105046048A (en) Solving method of takeoff performance of ground effect vehicle
CN104573265B (en) A kind of method that ringsail parachute aeroperformance is optimized
CN103939283A (en) Efficient blade special for vertical axis wind turbine
CN104166348A (en) Dynamic stability judgment method of single-slider variable-mass-center control aircraft
CN103676652A (en) Dynamic modeling method of flash evaporator
CN112926132B (en) Fixed wing airfoil aerodynamic shape design method considering influence of three-dimensional effect
CN111090908A (en) Method for selecting design parameters of group umbrella system
DRELA Viscous and inviscid inverse schemes using Newton's method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant