CN106021770A - Stability analysis method and device for airplane structure rectangular flat plate - Google Patents
Stability analysis method and device for airplane structure rectangular flat plate Download PDFInfo
- Publication number
- CN106021770A CN106021770A CN201610371614.XA CN201610371614A CN106021770A CN 106021770 A CN106021770 A CN 106021770A CN 201610371614 A CN201610371614 A CN 201610371614A CN 106021770 A CN106021770 A CN 106021770A
- Authority
- CN
- China
- Prior art keywords
- temprx
- rectangular flat
- ratio
- curve data
- value
- 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.)
- Granted
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 60
- 230000006835 compression Effects 0.000 claims abstract description 33
- 238000007906 compression Methods 0.000 claims abstract description 33
- 238000010008 shearing Methods 0.000 claims abstract description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 230000008901 benefit Effects 0.000 claims description 7
- 238000010845 search algorithm Methods 0.000 abstract 1
- 230000008569 process Effects 0.000 description 16
- 238000010586 diagram Methods 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 6
- 238000004590 computer program Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 241000283690 Bos taurus Species 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000012804 iterative process Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000006399 behavior Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The embodiment of the invention provides a stability analysis method and device for an airplane structure rectangular flat plate. The method comprises the steps that electronization is carried out on a curve chart of dual-pressure shearing and buckling related curves of the rectangular flat plate, and curve data is obtained, wherein the x-coordinates, the y-coordinates and values on the curves of the curve chart are marked as CurveRx, CurveRy and CurveRs in the curve data respectively; Rx, Ry and Rs are calculated according to the compression allowable stress sigma ccrx of the rectangular flat plate in the X direction, the compression allowable stress sigma ccry of the rectangular flat plate in the Y direction, the shearing allowable stress tau cr, the compression working stress sigma cx in the X direction, the compression working stress sigma cy in the Y direction and the shearing working stress tau; the CurveRx value, the CurveRy value and the CurveRs value meeting formula (I) are sought in the curve data based on a binary search algorithm; the safety margin MS=K-1 of the rectangular flat plate is calculated; the stability of the rectangular flat plate is analyzed according to the safety margin MS. By means of the method and device, the worked out safety margin is high in scientificity and reliability, and therefore the stability analysis result of the airplane structure rectangular flat plate is more reasonable and reliable.
Description
Technical field
The present invention relates to aircaft configuration technical field, particularly relate to the stability analysis of a kind of aircaft configuration rectangular flat
Method and apparatus.
Background technology
The design of aircaft configuration bearing stability and explication de texte work, for the load-bearing safety of structure, the life-span and
The aspect such as Weight-optimised all has important engineering significance.Wherein, " flat board double pressure shearing instability analysis " is by FEM calculation
Result combines with trial curve and empirical equation and calculates flat board stability.
And in " the structural stability design manual " that Mr. Cui Degang writes, by means of this flat board double pressure shearing instability
When the curve analyzed carries out the stability analysis of aircraft rectangular flat structure, mainly by calculated curve abscissa Rx, curve
On value Rs, the value of vertical coordinate Ry, and by use two values therein, such as Rx, Rs as input parameter, then the most permissible
Determining that on curve value Ry is corresponding with them uniquely, Ry CurveRy here represents, then margin of safety is
But, the method for this calculating margin of safety, it is assumed that be that the compression stress of X-direction is constant, shear stress is constant,
The compression stress allowing Y-direction increases, the margin of safety drawn.
Obviously, this hypothesis do not meet the actual loading situation of aircaft configuration rectangular flat, then calculating with this
To margin of safety carry out stability analysis time, it is clear that the problem that there is poor reliability.
Summary of the invention
Embodiment of the present invention technical problem to be solved is to provide the stability analysis of a kind of aircaft configuration rectangular flat
Method and apparatus, to solve in prior art when aircaft configuration rectangular flat is carried out stability analysis, does not considers that aircraft is tied
The problem of the analysis result poor reliability that the actual loading situation of structure rectangular flat is caused.
In order to solve the problems referred to above, according to an aspect of the present invention, the invention discloses a kind of aircaft configuration rectangle to put down
The method for analyzing stability of plate, including:
The curve chart of rectangular flat double pressure shear buckling correlation curve is carried out electronic disposal, obtains curve data, its
In, the value on the abscissa of described curve chart, vertical coordinate and curve be designated as respectively in described curve data CurveRx,
CurveRy and CurveRs；
The compression allowable stress σ of the X-direction according to rectangular flatccrx, the compression work stress σ of X-directioncx, the pressure of Y-direction
Contracting allowable stress σccry, the compression work stress σ of Y-directioncy, shear allowable stress τcr, shearing work stress τ calculate:
In described curve data, CurveRx, the CurveRy meeting below equation (1) is searched based on binary chop algorithm
And the value of CurveRs:
Calculate margin of safety MS=K-1 of described rectangular flat；
The stability of described rectangular flat is analyzed according to described margin of safety MS.
According to a further aspect in the invention, the invention also discloses the stability analysis dress of a kind of aircaft configuration rectangular flat
Put, including:
Electronic disposal module, for carrying out electronization to the curve chart of double pressure shear buckling correlation curves of rectangular flat
Processing, obtain curve data, wherein, the value on the abscissa of described curve chart, vertical coordinate and curve is in described curve data
It is designated as CurveRx, CurveRy and CurveRs respectively；
First computing module, for the compression allowable stress σ of the X-direction according to rectangular flatccrx, the compression work of X-direction
Make stress σcx, the compression allowable stress σ of Y-directionccry, the compression work stress σ of Y-directioncy, shear allowable stress τcr, shear work
Make stress τ to calculate:
Search module, meet below equation (1) for searching in described curve data based on binary chop algorithm
The value of CurveRx, CurveRy and CurveRs:
Second computing module, for calculating margin of safety MS=K-1 of described rectangular flat；
Analyze module, for analyzing the stability of described rectangular flat according to described margin of safety MS.
Compared with prior art, the embodiment of the present invention includes advantages below:
The embodiment of the present invention uses and figure is carried out electronic disposal, utilizes follow-up data search analysis；And utilize two
Point-score is searched technology and is searched the point meeting preset formula in the graphic, it is achieved respectively synchronize to amplify, so that calculate to load
Margin of safety there is stronger science and reliability, so that the stability analysis result to aircaft configuration rectangular flat
More rationally reliable.
Accompanying drawing explanation
Fig. 1 is the flow chart of steps of the method for analyzing stability embodiment of a kind of aircaft configuration rectangular flat of the present invention；
Fig. 2 is the schematic diagram of the rectangular flat that a kind of length-width ratio is 2 double pressure shear buckling correlation curve of the present invention；
Fig. 3 is the schematic diagram of the rectangular flat that a kind of length-width ratio is 3 double pressure shear buckling correlation curve of the present invention；
After Fig. 4 is the rectangular flat that a kind of length-width ratio is 2 double pressure shear buckling correlation curve electronic disposal of the present invention
Schematic diagram；
After Fig. 5 is the rectangular flat that a kind of length-width ratio is 3 double pressure shear buckling correlation curve electronic disposal of the present invention
Schematic diagram；
Fig. 6 is the flow process signal of the method for analyzing stability embodiment of the another kind of aircaft configuration rectangular flat of the present invention
Figure；
Fig. 7 is the steps flow chart of the method for analyzing stability embodiment of the another kind of aircaft configuration rectangular flat of the present invention
Figure；
Fig. 8 is the structured flowchart of the analysis of stability analysis apparatus embodiment of a kind of aircaft configuration rectangular flat of the present invention.
Detailed description of the invention
Understandable for enabling the above-mentioned purpose of the present invention, feature and advantage to become apparent from, real with concrete below in conjunction with the accompanying drawings
The present invention is further detailed explanation to execute mode.
With reference to Fig. 1, it is shown that the step of the method for analyzing stability embodiment of a kind of aircaft configuration rectangular flat of the present invention
Rapid flow chart, specifically may include steps of:
Step 101, carries out electronic disposal to the curve chart of rectangular flat double pressure shear buckling correlation curve, obtains curve
Data, wherein, the value on the abscissa of described curve chart, vertical coordinate and curve is designated as in described curve data respectively
CurveRx, CurveRy and CurveRs；
Wherein, utilize in the present embodiment " flat board double pressure shear stability analysis " shears about the double pressure of rectangular flat and bend
Bent correlation curve carries out stability analysis to aircaft configuration rectangular flat (being also a kind of rectangular flat structure).
Step 103, according to the compression allowable stress σ of the X-direction of rectangular flatccrx, the compression work stress σ of X-directioncx、Y
The compression allowable stress σ in directionccry, the compression work stress σ of Y-directioncy, shear allowable stress τcr, shearing work stress τ meter
Calculate:
Step 105, search in described curve data based on binary chop algorithm meet below equation (1) CurveRx,
The value of CurveRy and CurveRs:
Step 107, calculates margin of safety MS=K-1 of described rectangular flat；
Step 109, analyzes the stability of described rectangular flat according to described margin of safety MS.
Wherein, calculated margin of safety is the margin of safety of analyzed aircaft configuration rectangular flat, thus enters
One step analyzes the stability of this aircaft configuration rectangular flat.
By means of the technical scheme of the above embodiment of the present invention, the embodiment of the present invention uses and carries out figure at electronization
Reason, utilizes follow-up data search analysis；And utilize binary search technology to search the point meeting preset formula in the graphic, real
Now respectively synchronize to amplify to load, so that the margin of safety calculated has stronger science and reliability, so that
More rationally reliable to the stability analysis result of aircaft configuration rectangular flat.
In order to be better understood from the technique scheme of the present invention, it is combined unstability below in conjunction with a concrete rectangular flat and divides
The technique scheme of the present invention is described in detail by the embodiment that analysis calculates.
First, rectangular flat double pressure shear buckling correlation curve is carried out electronic disposal, obtains curve data, wherein,
Value on the abscissa of described curve chart, vertical coordinate and curve is designated as CurveRx, CurveRy in described curve data respectively
And CurveRs；In Fig. 1, Fig. 2, abscissa is Rx, and vertical coordinate is Ry, and the numeral on graph curve is Rs, i.e. with Rs for certain
The equal pitch contour of numeral, then solution procedure is as follows:
First, calculating rectangular flat stability is combined Instability Analysis when, length-width ratio (the i.e. a/ of computing board is typically wanted
B), this example gives double pressures of the a/b=3 shown in double pressure shear buckling correlograms of the a/b=2 shown in Fig. 2 and Fig. 3 and cuts
Cut flexing correlogram；
Here, the curve chart shown in Fig. 2 and Fig. 3 is carried out electronic disposal, i.e. by corresponding horizontal stroke each in curve chart
Value on coordinate, vertical coordinate and curve is all processed into the curve data of electronization, and the curve data obtained is respectively at Fig. 4 and Tu
Shown in 5.Wherein, in figures 4 and 5, the abscissa Rx (being designated as CurveRx below) in curve chart is shown in X list, from figure
2, Fig. 3 is it can be seen that Rx value is in the range of [0,1], the Y1 in Fig. 4 and Fig. 5, Y2 ... Y10 is expressed as figure in curve chart
The value (being designated as CurveRy below) of the vertical coordinate Ry on each bar curve in shape, 10 altogether；Numeral on every curve
The value (being designated as CurveRs below) of Rs defines by field in the curve data (such as data base) of electronization.
It should be noted that the curve chart the most only schematically illustrating rectangular flat that length-width ratio is 2,3
Electronic disposal process, substantially arbitrarily length-width ratio are that the electronic disposal process of the curve chart of the rectangular flat of integer is homogeneous
With, therefore, repeat no more.
Wherein, in actual applications, on the one hand, for the ease of analyzing the steady of the aircaft configuration rectangular flat of different length-width ratio
Qualitative, double pressure shear buckling correlograms of the aircaft configuration rectangular flat that multiple length-width ratios are integer can be made electricity respectively
Sub-ization processes, and obtains the multiple curve data of corresponding different length-width ratio；On the other hand, it is also possible to computational stability as required
The concrete ratio of length and width of aircaft configuration rectangular flat, select double pressure shear buckling correlograms of corresponding ratio
Carry out electronic disposal, obtain the curve data that correspondence is currently needed for the length-width ratio of the aircaft configuration rectangular flat analyzed, subtract
Few system-computed amount.
Wherein, for same rectangular flat unit, at material behavior (elastic modelling quantity, Poisson's ratio, yield limit stress)
In the case of determining with geometrical property (length and width, thickness), its allowable stress value is constant, i.e. the compression allowable stress of X-direction
σccrx, the compression allowable stress σ of Y-directionccry, the shearing allowable stress τ of rectangular flatcrFor the value determined, the compression work of X-direction
Make stress σcx, the compression work stress σ of Y-directioncy, shearing work stress τ is that FEM calculation is on this rectangular flat unit
Stress result, therefore, it can actual stress based on rectangular flat unit and calculates Rx, Ry, Rs:
Then, then based on binary chop algorithm search in response curve data meet below equation (1) CurveRx,
The value of CurveRy and CurveRs:
Specifically, when the ratio of the length and width of this aircaft configuration rectangular flat is integer, such as ratio is 2,
The curve data (that is, the data shown in Fig. 4) the most then selecting ratio to be 2 in the multiple curve data of corresponding different ratios；
And lookup meets the point (i.e. the value of CurveRx, CurveRy and CurveRs) of formula (1) from the data shown in Fig. 4.Its
In, from above-mentioned electronic disposal step it can be seen that point any one curve, the abscissa Rx of this point is in curve data
Representing with CurveRx, vertical coordinate Ry represents with CurveRy in curve data, and value Rs on curve is used in curve data
CurveRs represents.
Therefore it may only be necessary to search the point on the curve shown in a Fig. 2 in the curve data shown in Fig. 4 so that this point
The value of CurveRx, CurveRy and CurveRs respectively with Rx, Ry, the Rs's hereinbefore utilizing various Stress calculation to obtain
Corresponding ratio is identical, and ratio k now represents.
And if the value of length-width ratio a/b of this aircaft configuration rectangular flat is decimal 2.5, i.e. between 2 and 3, then need
The curve of two integer ratios (i.e. ratio be 2,3) adjacent with 2.5 is selected in the multiple curve data of corresponding different ratios
Data (that is, the curve data corresponding to Fig. 3 and Fig. 4)；And between the curve data of Fig. 3 and Fig. 4, carry out the meter of Newton interpolating method
Calculate the value searching CurveRx, CurveRy and CurveRs of meeting formula (1).
So, when corresponding with virgin curve figure in the length-width ratio facing rectangular flat length-width ratio is not mated, it is possible to use cattle
Interpolation method of pausing calculates, and has the most sufficient practicality.
It should be noted that for the value on the abscissa in the curve chart shown in Fig. 2, Fig. 3, vertical coordinate and curve
The meaning that they represent and the same meaning of actual stress based on rectangular flat unit calculated Rx, Ry, Rs, distinguish
In being theoretical value in the graph, calculated, it is actual value.
Then, it is possible to calculate margin of safety MS of rectangular flat, i.e. MS=K-1.
Finally, so that it may analyze the stability of described rectangular flat according to described margin of safety MS.Specifically, in basis
When described margin of safety MS analyzes the stability of described rectangular flat, it can be determined that whether described margin of safety MS is more than zero；As
Fruit judges more than zero, it is determined that having good stability of described rectangular flat；If it is determined that less than zero, it is determined that described rectangular flat
Poor stability.
Wherein, in one embodiment, described public affairs are met searching in described curve data based on binary chop algorithm
During the value of CurveRx, CurveRy and CurveRs of formula (1), it is first determined meeting condition is: after adjacent twice iteration,
When the absolute value of the difference of CurveRyi/Ry Yu tempRx/Rx is less than 0.001, it is believed that meet formula (1).
It is with the ratio of i & lt Rx that two way classification calculates the basic ideas of margin of safetyThe ratio searched with i+1 time
ValueDo subtraction, if the absolute value of difference is less than or equal to 0.001, then K=CurveRyi+1/ Ry, required
Margin of safety
Although it should be noted that be in the present embodiment byMake comparisons with 0.001, but
It is that the value that it will be understood by those of skill in the art that comparison other is not limited to 0.001, it is also possible to be 0.002,0.01 etc.,
It is to say, the value of comparison other can set flexibly according to the required precision of actual application scenarios herein, if required precision
Height, then the value of this comparison other is relatively small, is closer to 0, and if required precision is relatively low, then the value of this comparison other
It is less than 0.05.
Implementation step, with reference to Fig. 7, specifically includes:
The Rx that A, acquisition obtain according to various above-mentioned Stress calculation, Rs, Ry, CurveRy are to perform the S001 shown in Fig. 6
Result；
Definition logical judgment variable criterionLoop, for judging that the sign of MS is changed the most in an iterative process
Becoming, initial value is false；
Initial value definition start, end, for represent adjacent twice iteration sign change time, initial tempRx and
The tempRx of ending, represents that iteration result will produce between start and end；
Definition temporary variable tempRx, tempRs, for storing the input of the S001 file in iterative process；
Definition numberK, for calculating interim ratio；
Define Dynamic Array ListMS, for recording the difference of K/Ry Yu tempRx/Rx；
Define Dynamic Array ListRx, ListRs, the tempRx during recording each loop iteration,
tempRy。
B, make tempRx=Rx, tempRs=Rs, i.e. input calculated Rx, Rs, look into curve data and obtain first
CurveRy1, because, tempRx/Rx=1, tempRs/Rs=1.Therefore, at this moment Rx, Rs are equivalent to amplify 1 times, calculate difference
CurveRy1/Ry-tempRx/Rx, i.e. CurveRy1/Ry-1, be stored in ListMS, ifThen table
Show that the value currently found meets margin of safety computing formula (1), then margin of safety MS=CurveRy1/Ry-1.
If C looks into for the first time value CurveRy1 that curve obtains and is unsatisfactory for margin of safety formula (1), i.e. CurveRy1/
Ry-1 > 0.001, then be multiplied by the multiple of Rx, Rs same ratio value (CurveRy1/Ry+1)/2 with two way classification, and tied by multiplication
Fruit represents with tempRx, tempRs respectively, and records the difference of this CurveRy/Ry Yu tempRx/Rx performed with ListMS
(for example, number less than 0), ListRx records this tempRx, ListRs performed and records this tempRs performed.
D, it is that input continues executing with S001 and looks into curve with this tempRx, tempRs；
E is if meeting the computing formula (1) of margin of safety, then margin of safety MS=CurveRy2/Ry-(CurveRy1/
Ry+1)/2；
If F is unsatisfactory for, the most first judge last ListMS [i] (i.e. the difference of CurveRy/Ry Yu tempRx/Rx)
Symbol；
If the symbol of the last time that the symbol of G F step is adjacent with F step is identical, then decision logic variable
CriterionLoop, if false, then amplifies (CurveRy/Ry+ with two way classification by the multiple of tempRx, tempRs
TempRx/Rx)/2 times, record ListMS, ListRx, ListRs, continue executing with step D；
If the last symbol that the symbol of H F step is adjacent from F step is different, ListMS [i-1] and ListMs [i] is described
Opposite sign, if logical variable criterionLoop is false, then start=ListRx [i-1]/Rx, end=ListRx [i]/
Rx, the value of iteration will produce between satrt and end, utilize two way classification that tempRx, tempRs are amplified (start+end)/2
Times, logical variable criterionLoop becomes true, performs D；
If I is by ListMS [i] and ListMS [i-1] jack per line after H step execution, the ratio of tempRx is described the most relatively
Greatly, at this moment the value of start is constant, and end=ListRx [i]/Rx, tempRx, tempRy continue to amplify (start+end)/2
Times, perform D；
If J is by ListMS [i] and ListMS [i-1] opposite sign after H step execution, illustrate that tempRx ratio is less, at this moment
The value of start=ListRx [i]/Rx, end is constant, and tempRx, tempRy continue to amplify (start+end)/2 times, perform D；
K, until meet condition end loop.
The embodiment of the present invention will be looked in rectangular flat stable composition general formula in " structural stability design manual "
The curve looked for is looked into value method and is improved, and allows the compressive load of X-direction, the compressive load of Y-direction and shear-type load simultaneous growth,
The result of calculation making margin of safety is more reasonable, so that the stability analysis result more adduction of aircaft configuration rectangular flat
Reason, and more reference significance.
It should be noted that for embodiment of the method, in order to be briefly described, therefore it is all expressed as a series of action group
Closing, but those skilled in the art should know, the embodiment of the present invention is not limited by described sequence of movement, because depending on
According to the embodiment of the present invention, some step can use other orders or carry out simultaneously.Secondly, those skilled in the art also should
Knowing, embodiment described in this description belongs to preferred embodiment, and the involved action not necessarily present invention implements
Necessary to example.
Think corresponding with the method that the invention described above embodiment is provided, with reference to Fig. 8, it is shown that one aircaft configuration of the present invention
The structured flowchart of the analysis of stability analysis apparatus embodiment of rectangular flat, specifically can include such as lower module:
Electronic disposal module 81, for carrying out at electronization double pressure shear buckling correlograms of rectangular flat
Reason, obtains curve data, and wherein, the value on the abscissa of described curve chart, vertical coordinate and curve is divided in described curve data
It is not designated as CurveRx, CurveRy and CurveRs；
First computing module 82, for the compression allowable stress σ of the X-direction according to rectangular flatccrx, the compression of Y-direction
Allowable stress σccry, shear allowable stress τcr, the compression work stress σ of X-directioncx, the compression work stress σ of Y-directioncx, shear
Working stress τ calculates:
Search module 83, meet below equation (1) for searching in described curve data based on binary chop algorithm
The value of CurveRx, CurveRy and CurveRs:
Second computing module 84, for calculating margin of safety MS=K-1 of described rectangular flat；
Analyze module 85, for analyzing the stability of described rectangular flat according to described margin of safety MS.
Wherein, described electronic disposal module 81, for the rectangular flat that ratio the is integer difference to length and width
Make the electronic disposal of double pressure shear buckling correlogram, obtain the multiple curve data of corresponding different ratio.
Described lookup module 83, is additionally operable to when the ratio of the length and width of described rectangular flat is integer, described
The multiple curve data of corresponding different ratio selects the curve data of the ratio of the length and width of corresponding described rectangular flat
Lookup meets the value of CurveRx, CurveRy and CurveRs of described formula (1).
Optionally, described device also includes:
Select module (not shown), be used for when the ratio of the length and width of described rectangular flat is decimal, described
The multiple curve data of corresponding different ratio selects the curve data of two integer ratios adjacent with described decimal；Accordingly
Described lookup module 83, is additionally operable to carry out cattle between the curve data of two integer ratios adjacent with described decimal
The value of CurveRx, CurveRy and CurveRs of meeting described formula (1) is searched in the calculating of differential technique.
Additionally, described analysis module 85 includes:
First judges submodule (not shown), is used for judging that whether described margin of safety MS is more than zero；
First determines submodule (not shown), for if it is determined that be more than zero, it is determined that the stability of described rectangular flat
Well；
Second determines submodule (not shown), for if it is determined that be less than zero, it is determined that the stability of described rectangular flat
Difference.
Described lookup module 83 comprises the steps that
Represent submodule (not shown), for iteration variable tempRx of Rx is represented, the iteration variable of Rs is used
TempRs represents, utilizes the ith iteration result of the vertical coordinate that tempRx and tempRs check in described curve data to use
CurveRyi represents；
Two points of calculating sub module (not shown), for allowing calculated described Rx and described Rs take advantage of respectively by two way classification
It is worth in the same scale, obtains described tempRx and described tempRs；
Wherein,Wherein,For described Rx
The described ratio value taken advantage of during with described Rs loop iteration, and, when iteration first, i=1, make tempRx=Rx, tempRs
=Rs；
Search submodule (not shown), corresponding with described tempRx, tempRs for searching in described curve data
CurveRyi；
First calculating sub module (not shown), for calculating the ratio of described CurveRyi and calculated described Ry,
Obtain CurveRyi/Ry；
Second calculating sub module (not shown), for calculating the ratio of tempRx and calculated described Rx, obtains
tempRx/Rx；
Record sub module (not shown), for recording the difference of CurveRyi/Ry Yu tempRx/Rx, obtains CurveRyi/
Ry-tempRx/Rx and the absolute value of described difference, obtain | CurveRyi/Ry-tempRx/Rx |；
Second judges submodule (not shown), is used for judging whether | CurveRyi/Ry-tempRx/Rx | is less than or equal to
0.001；
3rd determines submodule (not shown), for working asTime, determine K=
CurveRyi/Ry。
Variable assignments module (not shown), is used for when | CurveRyi/Ry-tempRx/Rx | is more than 0.001, and is sentencing
After disconnected ith iteration ratio C urveRyi/Ry of CurveRyi Yu Ry and the ratio tempRx/Rx of tempRx Yu Rx difference and
When the described difference symbol of the 1st calculating is contrary, makes the ratio of the variable start tempRx Yu Rx equal to the i-th-1 time, make variable
End is equal to the ratio of i & lt tempRx with Rx；
Described two points of calculating sub module (not shown), are additionally operable to utilize binary search method allow described tempRs and described
TempRs is multiplied by the same ratio value between start and end respectively, carries out the iterative step of described two way classification, until full
FootThen exit iteration.
The embodiment of the present invention will be looked in rectangular flat stable composition general formula in " structural stability design manual "
The curve looked for is looked into value method and is improved, and allows the compressive load of X-direction, the compressive load of Y-direction and shear-type load simultaneous growth,
The result of calculation making margin of safety is more reasonable, so that the stability analysis result more adduction of aircaft configuration rectangular flat
Reason, and more reference significance.
For device embodiment, due to itself and embodiment of the method basic simlarity, so describe is fairly simple, relevant
Part sees the part of embodiment of the method and illustrates.
Each embodiment in this specification all uses the mode gone forward one by one to describe, what each embodiment stressed is with
The difference of other embodiments, between each embodiment, identical similar part sees mutually.
Those skilled in the art are it should be appreciated that the embodiment of the embodiment of the present invention can be provided as method, device or calculate
Machine program product.Therefore, the embodiment of the present invention can use complete hardware embodiment, complete software implementation or combine software and
The form of the embodiment of hardware aspect.And, the embodiment of the present invention can use one or more wherein include computer can
With in the computer-usable storage medium (including but not limited to disk memory, CD-ROM, optical memory etc.) of program code
The form of the computer program implemented.
The embodiment of the present invention is with reference to method, terminal unit (system) and computer program according to embodiments of the present invention
The flow chart of product and/or block diagram describe.It should be understood that can be by computer program instructions flowchart and/or block diagram
In each flow process and/or the flow process in square frame and flow chart and/or block diagram and/or the combination of square frame.These can be provided
Computer program instructions sets to general purpose computer, special-purpose computer, Embedded Processor or other programmable data processing terminals
Standby processor is to produce a machine so that held by the processor of computer or other programmable data processing terminal equipment
The instruction of row produces for realizing in one flow process of flow chart or multiple flow process and/or one square frame of block diagram or multiple square frame
The device of the function specified.
These computer program instructions may be alternatively stored in and can guide computer or other programmable data processing terminal equipment
In the computer-readable memory worked in a specific way so that the instruction being stored in this computer-readable memory produces bag
Including the manufacture of command device, this command device realizes in one flow process of flow chart or multiple flow process and/or one side of block diagram
The function specified in frame or multiple square frame.
These computer program instructions also can be loaded on computer or other programmable data processing terminal equipment so that
On computer or other programmable terminal equipment, execution sequence of operations step is to produce computer implemented process, thus
The instruction performed on computer or other programmable terminal equipment provides for realizing in one flow process of flow chart or multiple flow process
And/or the step of the function specified in one square frame of block diagram or multiple square frame.
Although having been described for the preferred embodiment of the embodiment of the present invention, but those skilled in the art once knowing base
This creativeness concept, then can make other change and amendment to these embodiments.So, claims are intended to be construed to
The all changes including preferred embodiment and falling into range of embodiment of the invention and amendment.
Finally, in addition it is also necessary to explanation, in this article, the relational terms of such as first and second or the like be used merely to by
One entity or operation separate with another entity or operating space, and not necessarily require or imply these entities or operation
Between exist any this reality relation or order.And, term " includes ", " comprising " or its any other variant meaning
Containing comprising of nonexcludability, so that include that the process of a series of key element, method, article or terminal unit not only wrap
Include those key elements, but also include other key elements being not expressly set out, or also include for this process, method, article
Or the key element that terminal unit is intrinsic.In the case of there is no more restriction, by wanting that statement " including ... " limits
Element, it is not excluded that there is also other identical element in including the process of described key element, method, article or terminal unit.
Method for analyzing stability and a kind of aircraft knot to a kind of aircaft configuration rectangular flat provided by the present invention above
The analysis of stability analysis apparatus of structure rectangular flat, is described in detail, the specific case principle to the present invention used herein
And embodiment is set forth, the explanation of above example is only intended to help to understand that the method for the present invention and core thereof are thought
Think；Simultaneously for one of ordinary skill in the art, according to the thought of the present invention, in specific embodiments and applications
All will change, in sum, this specification content should not be construed as limitation of the present invention.
Claims (14)
1. the method for analyzing stability of an aircaft configuration rectangular flat, it is characterised in that including:
The curve chart of rectangular flat double pressure shear buckling correlation curve is carried out electronic disposal, obtains curve data, wherein, institute
State the value on the abscissa of curve chart, vertical coordinate and curve be designated as respectively in described curve data CurveRx, CurveRy and
CurveRs；
The compression allowable stress σ of the X-direction according to rectangular flatccrx, the compression allowable stress σ of Y-directionccry, shear allowable stress
τcr, the compression work stress σ of X-directioncx, the compression work stress σ of Y-directioncy, shearing work stress τ calculate:
Search in described curve data based on binary chop algorithm meet below equation (1) CurveRx, CurveRy and
The value of CurveRs:
Calculate margin of safety MS=K-1 of described rectangular flat；
The stability of described rectangular flat is analyzed according to described margin of safety MS.
Method the most according to claim 1, it is characterised in that described double pressure shear buckling correlation curves to rectangular flat
Figure carries out electronic disposal, and the step obtaining curve data includes:
The rectangular flat that ratio is integer of length and width is made at the electronization of double pressure shear buckling correlogram respectively
Reason, obtains the multiple curve data of corresponding different ratio.
Method the most according to claim 2, it is characterised in that when the ratio of the length and width of described rectangular flat is whole
During number, the multiple curve data of described corresponding different ratios selects the ratio of the length and width of corresponding described rectangular flat
Curve data search the value of CurveRx, CurveRy and CurveRs meeting described formula (1).
Method the most according to claim 2, it is characterised in that described method also includes:
When the ratio of the length and width of described rectangular flat is decimal, at the multiple curve data of described corresponding different ratios
The curve data of two integer ratios that middle selection is adjacent with described decimal；
The calculating carrying out newton differential technique between the curve data of two integer ratios adjacent with described decimal is satisfied to search
The value of CurveRx, CurveRy and CurveRs of described formula (1).
Method the most according to claim 1, it is characterised in that described flat according to the described margin of safety MS described rectangle of analysis
The step of the stability of plate includes:
Judge that whether described margin of safety MS is more than zero；
If it is determined that more than zero, it is determined that having good stability of described rectangular flat；
If it is determined that less than zero, it is determined that the poor stability of described rectangular flat.
Method the most according to claim 1, it is characterised in that described based on binary chop algorithm in described curve data
The step of the value that lookup meets CurveRx, CurveRy and CurveRs of described formula (1) includes:
Iteration variable tempRx of Rx is represented, iteration variable tempRs of Rs is represented, utilizes tempRx and tempRs
Ith iteration result CurveRyi of the vertical coordinate checked in described curve data represents；
Allow calculated described Rx and described Rs be multiplied by identical ratio value respectively by two way classification, obtain described tempRx and
Described tempRs；
Wherein,Wherein,For described Rx and institute
The described ratio value taken advantage of when stating Rs loop iteration, and, when iteration first, i=1, make tempRx=Rx, tempRs=
Rs；
The CurveRyi corresponding with described tempRx, tempRs is searched in described curve data；
Calculate the ratio of described CurveRyi and calculated described Ry, obtain CurveRyi/Ry；
Calculate the ratio of tempRx and calculated described Rx, obtain tempRx/Rx；
The difference of record CurveRyi/Ry Yu tempRx/Rx, obtains CurveRyi/Ry-tempRx/Rx and described difference
Absolute value, obtains | CurveRyi/Ry-tempRx/Rx |；
Judge that whether | CurveRyi/Ry-tempRx/Rx | is less than or equal to 0.001；
WhenTime, determine K=CurveRyi/Ry.
Method the most according to claim 6, it is characterised in that when | CurveRyi/Ry-tempRx/Rx | is more than 0.001
Time, described search in described curve data based on binary chop algorithm CurveRx, CurveRy of meeting described formula (1) with
And the step of the value of CurveRs also includes:
Ratio C urveRyi/Ry of CurveRyi Yu Ry and the ratio tempRx/Rx of tempRx Yu Rx after judging ith iteration
Difference contrary with the described difference symbol that the 1st time calculates time, make the variable start ratio equal to tempRx and Rx of the i-th-1 time
Value, makes the variable end ratio equal to i & lt tempRx with Rx；
Utilize the same ratio value that binary search method allows described Rx and described Rs be multiplied by respectively between start and end,
Go out tempRx and tempRs, carry out the iterative step of described two way classification, until meetingThen
Exit iteration.
8. the analysis of stability analysis apparatus of an aircaft configuration rectangular flat, it is characterised in that including:
Electronic disposal module, for the curve chart of rectangular flat double pressure shear buckling correlation curve is carried out electronic disposal,
Obtaining curve data, wherein, the value on the abscissa of described curve chart, vertical coordinate and curve is remembered in described curve data respectively
For CurveRx, CurveRy and CurveRs；
First computing module, for the compression allowable stress σ of the X-direction according to rectangular flatccrx, the compression allowable stress of Y-direction
σccry, shear allowable stress τcr, the compression work stress σ of X-directioncx, the compression work stress σ of Y-directioncy, shearing work stress
τ calculates:
Search module, meet below equation (1) for searching in described curve data based on binary chop algorithm
The value of CurveRx, CurveRy and CurveRs:
Second computing module, for calculating margin of safety MS=K-1 of described rectangular flat；
Analyze module, for analyzing the stability of described rectangular flat according to described margin of safety MS.
Device the most according to claim 8, it is characterised in that described electronic disposal module, for length and width
The rectangular flat that ratio is integer make the electronic disposal of double pressure shear buckling correlogram respectively, obtain correspondence the most on year-on-year basis
The multiple curve data of value.
Device the most according to claim 9, it is characterised in that described lookup module, is additionally operable to when described rectangular flat
When the ratio of length and width is integer, corresponding described rectangle is selected to put down in the multiple curve data of described corresponding different ratios
The curve data of the ratio of the length and width of plate searches CurveRx, CurveRy and the CurveRs meeting described formula (1)
Value.
11. devices according to claim 9, it is characterised in that described device also includes:
Select module, be used for when the ratio of the length and width of described rectangular flat is decimal, at described corresponding different ratios
Multiple curve data in select the curve data of two integer ratios adjacent with described decimal；Accordingly
Described lookup module, is additionally operable to carry out newton difference between the curve data of two integer ratios adjacent with described decimal
The value of CurveRx, CurveRy and CurveRs of meeting described formula (1) is searched in the calculating of method.
12. devices according to claim 8, it is characterised in that described analysis module includes:
First judges submodule, is used for judging that whether described margin of safety MS is more than zero；
First determines submodule, for if it is determined that be more than zero, it is determined that having good stability of described rectangular flat；
Second determines submodule, for if it is determined that be less than zero, it is determined that the poor stability of described rectangular flat.
13. devices according to claim 8, it is characterised in that described lookup module includes:
Represent submodule, for iteration variable tempRx of Rx is represented, iteration variable tempRs of Rs is represented, utilize
Ith iteration result CurveRyi of the vertical coordinate that tempRx and tempRs checks in described curve data represents；
Two points of calculating sub module, for allowing calculated described Rx and described Rs take advantage of respectively in the same scale by two way classification
Value, obtains described tempRx and described tempRs；
Wherein,Wherein,For described Rx and institute
The described ratio value taken advantage of when stating Rs loop iteration, and, when iteration first, i=1, make tempRx=Rx, tempRs=
Rs；
Search submodule, for searching the CurveRyi corresponding with described tempRx, tempRs in described curve data；
First calculating sub module, for calculating the ratio of described CurveRyi and calculated described Ry, obtains CurveRyi/
Ry；
Second calculating sub module, for calculating the ratio of tempRx and calculated described Rx, obtains tempRx/Rx；
Record sub module, for recording the difference of CurveRyi/Ry Yu tempRx/Rx, obtains CurveRyi/Ry-tempRx/Rx
And the absolute value of described difference, obtain | CurveRyi/Ry-tempRx/Rx |；
Second judges submodule, is used for judging that whether | CurveRyi/Ry-tempRx/Rx | is less than or equal to 0.001；
3rd determines submodule, for working asTime, determine K=CurveRyi/Ry.
14. methods according to claim 13, it is characterised in that described lookup module also includes:
Variable assignments module, is used for when | CurveRyi/Ry-tempRx/Rx | is more than 0.001, and after judging ith iteration
The difference of the ratio tempRx/Rx of ratio C urveRyi/Ry of CurveRyi Yu Ry and tempRx Yu Rx and the institute calculated for the 1st time
State difference symbol contrary time, make the variable start ratio equal to tempRx Yu Rx of the i-th-1 time, make variable end be equal to i & lt
The ratio of tempRx Yu Rx；
Described two points of calculating sub module, are additionally operable to utilize binary search method to allow described Rx and described Rs be multiplied by respectively and are positioned at
Same ratio value between start and end, draws tempRx and tempRs, carries out the iterative step of described two way classification, until full
FootThen exit iteration.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610371614.XA CN106021770B (en) | 2016-05-30 | 2016-05-30 | The method for analyzing stability and device of aircaft configuration rectangular plate |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610371614.XA CN106021770B (en) | 2016-05-30 | 2016-05-30 | The method for analyzing stability and device of aircaft configuration rectangular plate |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106021770A true CN106021770A (en) | 2016-10-12 |
| CN106021770B CN106021770B (en) | 2019-08-09 |
Family
ID=57091858
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610371614.XA Withdrawn - After Issue CN106021770B (en) | 2016-05-30 | 2016-05-30 | The method for analyzing stability and device of aircaft configuration rectangular plate |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106021770B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106564615A (en) * | 2016-10-26 | 2017-04-19 | 北京数设科技有限公司 | Automatic plane structure sectional material partitioning method and device |
| CN107977488A (en) * | 2017-11-07 | 2018-05-01 | 西安合科软件有限公司 | Static strength check method and device |
| CN108733864A (en) * | 2017-04-25 | 2018-11-02 | 南京航空航天大学 | A kind of aircraft wing structure Global sensitivity analysis method based on support vector machines |
| CN108763624A (en) * | 2018-04-11 | 2018-11-06 | 北京航空航天大学 | A method of calculating uncertain but bounded buckling structure safety margin |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2379298A (en) * | 2001-09-04 | 2003-03-05 | Jahid Zeghiche | Analysis software for concrete-filled steel tubes |
| US20080164372A1 (en) * | 2007-01-05 | 2008-07-10 | Airbus France | Method of optimizing stiffened panels under stress |
| CN104778309A (en) * | 2015-03-19 | 2015-07-15 | 合科软件（北京）有限责任公司 | Airplane structure strength checking method and device |
| CN105197253A (en) * | 2015-08-14 | 2015-12-30 | 中国航空工业集团公司西安飞机设计研究所 | Safety margin calculation method of wing reinforced wallboard |
-
2016
- 2016-05-30 CN CN201610371614.XA patent/CN106021770B/en not_active Withdrawn - After Issue
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2379298A (en) * | 2001-09-04 | 2003-03-05 | Jahid Zeghiche | Analysis software for concrete-filled steel tubes |
| US20080164372A1 (en) * | 2007-01-05 | 2008-07-10 | Airbus France | Method of optimizing stiffened panels under stress |
| CN104778309A (en) * | 2015-03-19 | 2015-07-15 | 合科软件（北京）有限责任公司 | Airplane structure strength checking method and device |
| CN105197253A (en) * | 2015-08-14 | 2015-12-30 | 中国航空工业集团公司西安飞机设计研究所 | Safety margin calculation method of wing reinforced wallboard |
Non-Patent Citations (2)
| Title |
|---|
| 牛春匀: "《实用飞机结构应力分析及尺寸设计》", 31 December 2009 * |
| 空面导弹系统设计: "《空面导弹系统设计》", 31 December 2013 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106564615A (en) * | 2016-10-26 | 2017-04-19 | 北京数设科技有限公司 | Automatic plane structure sectional material partitioning method and device |
| CN106564615B (en) * | 2016-10-26 | 2019-03-22 | 北京数设科技有限公司 | A kind of method and device of aircaft configuration profile automatic Partitioning |
| CN108733864A (en) * | 2017-04-25 | 2018-11-02 | 南京航空航天大学 | A kind of aircraft wing structure Global sensitivity analysis method based on support vector machines |
| CN108733864B (en) * | 2017-04-25 | 2021-07-27 | 南京航空航天大学 | Aircraft wing structure global sensitivity analysis method based on support vector machine |
| CN107977488A (en) * | 2017-11-07 | 2018-05-01 | 西安合科软件有限公司 | Static strength check method and device |
| CN107977488B (en) * | 2017-11-07 | 2021-09-28 | 上海数设科技有限公司 | Static strength checking method and device |
| CN108763624A (en) * | 2018-04-11 | 2018-11-06 | 北京航空航天大学 | A method of calculating uncertain but bounded buckling structure safety margin |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106021770B (en) | 2019-08-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106021770A (en) | Stability analysis method and device for airplane structure rectangular flat plate | |
| Yang et al. | A stable node-based smoothed finite element method for metal forming analysis | |
| Keum et al. | Compatible description of tool surfaces and FEM meshes for analysing sheet forming operations | |
| EP2396736A2 (en) | Displaying multiple row and column header areas in a summary table | |
| Firat | U-channel forming analysis with an emphasis on springback deformation | |
| US11250185B2 (en) | Method and apparatus for calculating equivalent mechanical parameters of film layer etching region | |
| CN109784365A (en) | A kind of feature selection approach, terminal, readable medium and computer program | |
| Balzani et al. | Construction of statistically similar representative volume elements | |
| CN104834795B (en) | Band connection structure nonlinear contact with friction simulated behavior method and system | |
| CN104778309A (en) | Airplane structure strength checking method and device | |
| Valente et al. | Parameter identification and shape optimization: An integrated methodology in metal forming and structural applications | |
| CN117030448A (en) | Method and related device for analyzing mechanical test data based on digital image method | |
| Chow et al. | A generalized mixed isotropic-kinematic hardening plastic model coupled with anisotropic damage for sheet metal forming | |
| CN107609110B (en) | Mining method and device for maximum multiple frequent patterns based on classification tree | |
| Rangel et al. | Ftool 5.0: Nonlinear, stability and natural vibration analyses | |
| CN109635346A (en) | The analysis method for reliability of mechanical connecting structure | |
| US11164662B2 (en) | Simulation method, simulation program, and simulation device | |
| Olukunle et al. | Bernstein induced one step hybrid scheme for general solution of second order initial value problems | |
| Landkammer et al. | On inverse form finding for orthotropic plasticity | |
| de la Hoz et al. | A Sylvester-based IMEX method via differentiation matrices for solving nonlinear parabolic equations | |
| Valencia et al. | Bernstein polynomials in element-free Galerkin method | |
| Hu et al. | Comparison of explicit and implicit finite element methods in the quasistatic simulation of uniaxial tension | |
| CN106407620A (en) | ABAQUS-based engineering structure response surface stochastic finite element analysis processing method | |
| EP2325755B1 (en) | Method of searching a set of real numbers for a nearest neighbor | |
| US20130091160A1 (en) | High-Speed Statistical Processing in a Database |
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 | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 610A, block D, Zhongguancun International Incubation Park, No. 2, information road, Haidian District, Beijing 100085 Patentee after: Beijing Mowei Technology Co.,Ltd. Address before: 610A, block D, Zhongguancun International Incubation Park, No. 2, information road, Haidian District, Beijing 100085 Patentee before: BEIJING SHUSHE TECHNOLOGY Co.,Ltd. |
|
| CP01 | Change in the name or title of a patent holder | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210811 Address after: Exhibition hall A153, 1st floor, 11 Guotai Road, Yangpu District, Shanghai Patentee after: SHANGHAI DIGITAL DESIGN TECHNOLOGY Co.,Ltd. Address before: 610A, block D, Zhongguancun International Incubation Park, No. 2, information road, Haidian District, Beijing 100085 Patentee before: Beijing Mowei Technology Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| AV01 | Patent right actively abandoned |
Granted publication date: 20190809 Effective date of abandoning: 20240319 |
|
| AV01 | Patent right actively abandoned |
Granted publication date: 20190809 Effective date of abandoning: 20240319 |
|
| AV01 | Patent right actively abandoned | ||
| AV01 | Patent right actively abandoned |