CN105160135A - Method for calculating consistency of test data and analytical data of test piece - Google Patents

Abstract

The invention discloses a method for calculating the consistency of the test data and the analytical data of a test piece. The method comprises the following steps: step 1, selecting a test point on the test piece, performing a static test on the test piece under a preset load condition and obtain the test data of the test point; step 2, obtaining the calculated analytical data of the test piece under the preset load condition; step 3, obtaining a position relation of the test point in the step 1 to the test piece in the step 2, and establishing a corresponding relation table of the test point in the step 1 and the test piece in the step 2; and step 4, drawing a correlation curve chart of calculation and testing by calculating the correlation coefficient, mean-square deviation, relative deviation value and the like of the test data and the analytical data. A structural static analysis method provided by the invention realizes comparative calculation of the test measurement result and the structural analysis result of the test piece, contributes to the improvement of the quality of structural analysis and testing, and provides reference data for structural design.

Description

A kind of method calculating testpieces test figure and analyze data consistency

Technical field

The present invention relates to static construction analysis technical field, particularly relate to a kind of method calculating testpieces test figure and analyze data consistency.

Background technology

Static structural analysis experiment can produce a large amount of test measurement data, and structural design needs to carry out a large amount of structural analysis and computations, generates computation structure data.

Such as, aircaft configuration slow test produces a large amount of test measurement data, and structural design needs to carry out a large amount of structural analysis and computations, generates calculation result data.Aircaft configuration slow test data and analysis data consistency computing method, with structural analysis and computation result, calculating is compared to the test measurement result (mainly displacement and stress) at flight test vehicle examination position, and to result of calculation analysis, reason is looked in test and analysis two, modified test method and modeling method respectively, make both more reach unanimity, contribute to improving structure analysis and test mass, and provide reference data for structural design.

In prior art, still there is no a kind of method to carry out above-mentioned work.

Therefore, wish a kind of technical scheme to overcome or at least alleviate at least one above-mentioned defect of prior art.

Summary of the invention

The object of the present invention is to provide a kind of static structural analysis method to overcome or at least alleviate at least one the above-mentioned defect in prior art.

For achieving the above object, the invention provides a kind of method calculating testpieces test figure and analyze data consistency, for calculating aircraft or aircraft components, described method comprises the steps: step 1: the testing site on optional test part, under predetermined load condition, carry out the slow test of testpieces, and obtain the test figure of testing site; Step 2: obtain the computational analysis data of described testpieces under described predetermined load condition with calculation and analysis methods; Step 3: obtain the position relationship of the described testing site in described step 1 relative to the described testpieces in described step 2, and mapping table is set up in the described testing site in described step 1 and the testpieces in described step 2; Step 4: by calculating test figure and analyzing the relative coefficient of data, mean square deviation and relative deviation amount etc., the correlation curve figure of draw calculation and test.

Preferably, the testing site on described optional test part specifically comprises: displacement measurement point and stress measurement point.

Preferably, the method for selecting of institute's displacement measurement point is: the point of intersection selecting described each position of testpieces; The method for selecting of described stress measurement point is: be limited to and be selected in unidirectional tension and compression and plane stressed zone and/or be selected in STRESS VARIATION gradient shoulder.

Preferably, described predetermined load condition is: the load capacity choosing 67% of the design load of described testpieces carries out slow test.

Preferably, described calculation and analysis methods comprises analytical method and finite element method.

Preferably, described finite element method comprises: the finite element model setting up testpieces, and obtains the computational analysis data of described testpieces under described predetermined load condition; Described testing site in the finite element model of described testpieces is obtained by described step 3.

Preferably, the mapping table in described step 3 at least comprises following data:

The corresponding test period in described finite element model in testing site in described step 1, element numerals, element type, testing site coordinate, material behavior corresponding to testing site.

Preferably, described testing site is multiple.

Preferably, described test figure functional arrangement and described computational analysis data function figure comprise σ _j~ σ _scurve, Δ σ-σ _scurve, Russian curve ρ ~ | σ _s|, σ-X launches the multiple curve such as distribution, P-Δ/ε probability density distribution curve, and displacement DJ ~ DS curve displacement curve, the figures such as D-X displacement measurement curve.

Present invention also offers a kind of computing system, described computing system is used for the method calculating testpieces test figure and analyze data consistency as above, and described computing system comprises: Database Unit, and described Database Unit is used for storage function; Call unit, described call unit is for calling the function in described Database Unit; Function generation unit, described function generation unit is for generating test figure functional arrangement and computational analysis data function figure; Function contrast unit, described function contrast unit is used for described test figure functional arrangement and described analysis data function figure to contrast; Model architecture unit, described model architecture unit is used for the model of framework testpieces.

Static structural analysis method provided by the present invention, achieve to testpieces test measurement result (mainly displacement, stress or strain) and results of structural analysis comparison calculate, and to result of calculation analysis, reason is looked in test and analysis two, modified test method and modeling method respectively, makes both more reach unanimity.Contribute to improving structure analysis and test mass, and provide reference data for structural design.

Accompanying drawing explanation

Fig. 1 is according to the schematic flow sheet of calculating testpieces test figure of the present invention with the method for analysis data consistency.

Embodiment

For making object of the invention process, technical scheme and advantage clearly, below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is further described in more detail.In the accompanying drawings, same or similar label represents same or similar element or has element that is identical or similar functions from start to finish.Described embodiment is the present invention's part embodiment, instead of whole embodiments.Be exemplary type below by the embodiment be described with reference to the drawings, be intended to for explaining the present invention, and can not limitation of the present invention be interpreted as.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creation type work prerequisite, belong to the scope of protection of the invention.Below in conjunction with accompanying drawing, embodiments of the invention are described in detail.

In describing the invention; it will be appreciated that; term " " center ", " longitudinal direction ", " transverse direction ", "front", "rear", "left", "right", " vertically ", " level ", " top ", " end " " interior ", " outward " etc. instruction orientation or position relationship be based on orientation shown in the drawings or position relationship; be only the present invention for convenience of description and simplified characterization; instead of instruction or imply indication device or element must have specific orientation, with specific azimuth configuration and operation, therefore can not be interpreted as limiting the scope of the invention.

Fig. 1 is according to the schematic flow sheet of calculating testpieces test figure of the present invention with the method for analysis data consistency.

The method of calculating testpieces test figure as shown in Figure 1 and analysis data consistency is for calculating aircraft or aircraft components.

Calculating testpieces test figure as shown in Figure 1 comprises the steps: step 1 with the method analyzing data consistency: the testing site on optional test part, carries out the slow test of testpieces, and obtain the test figure of testing site under predetermined load condition; Step 2: obtain the computational analysis data of testpieces under predetermined load condition with calculation and analysis methods; Step 3: obtain the position relationship of testing site in step 1 relative to the described testpieces in step 2, and mapping table is set up in the testing site in step 1 and the testpieces in step 2; Step 4: by calculating test figure and analyzing the relative coefficient of data, mean square deviation and relative deviation amount etc., the correlation curve figure of draw calculation and test.

Be understandable that, departure (specifically comprising relative coefficient, mean square deviation and relative deviation amount) in above-mentioned step 4 is consistance value (trend), when this departure, more hour (this relative coefficient is close to 1, relative deviation amount more hour), its trend is convergent, thus consistance is better, then poorer on the contrary.

In the present embodiment, the testing site on optional test part specifically comprises: displacement measurement point and stress measurement point.Namely in the present embodiment, measurement point and stress measurement point is selected.

In the present embodiment, the method for selecting of displacement measurement point is: the point of intersection at each position of Selection experiment part; The method for selecting of stress measurement point is: be limited to and be selected in unidirectional tension and compression and plane stressed zone and/or be selected in STRESS VARIATION gradient shoulder.

Be understandable that, above-mentioned displacement measurement point may be selected in beam, rib, frame crossing point of axes place.Be understandable that, above-mentioned stress measurement point can also be chosen according to other conditions, or chooses condition in conjunction with each and comprehensively choose.Such as, stress measurement point can also be selected, by unidirectional tension and compression and plane stressed zone, should not be selected in triaxiality district; Select, in the less region of STRESS VARIATION gradient, should not area of stress concentration be selected in; Collapsing load should not be selected in and be less than maneuvering load part; Be selected in the preferred point of respective element as far as possible.

In the present embodiment, predetermined load condition is: the load capacity choosing 67% of the design load of testpieces carries out slow test.Why select above-mentioned predetermined load condition, because consider the impact of testpieces construction weight and loading lever weight, measure inaccurate when being less than 30% load, the load of participating in calculating is that the load at the zero point to 67% calculated according to the slope of the linearity range between 30% ~ 67% is poor.

In the present embodiment, calculation and analysis methods comprises analytical method and finite element method.

In the present embodiment, adopt finite element method, finite element method comprises: the finite element model setting up testpieces, and obtains the computational analysis data of testpieces under predetermined load condition; Described testing site in the finite element model of testpieces is obtained by described step 3.

In the present embodiment, the mapping table in step 3 at least comprises following data:

The corresponding test period in finite element model in testing site in step 1, element numerals (finite element element numerals, finite element element type), element type, testing site coordinate, material behavior corresponding to testing site.

Advantageously, in order to make test figure accurate, testing site is set to multiple.

In the present embodiment, test figure functional arrangement and computational analysis data function figure comprise calculated stress with test measure stress relation line, σ _j~ σ _scurve, Δ σ-σ _scurve, Russian curve ρ ~ | σ _s|, σ-X launches the multiple curve such as distribution, P-Δ/ε probability density distribution curve, and displacement DJ ~ DS curve displacement curve, the figures such as D-X displacement measurement curve.

Be understandable that, above-mentioned various function curves according to the difference of the physical quantity of the difference of tested testpieces and institute's preanalysis, and separately or can combine and are suitable for each above-mentioned function curve.

Particularly, convenient in order to describe, provide the concrete meaning of each curve below:

Symbol definition: σ _jcalculated stress value, σ _sstress value is measured in test, D _jdisplacement calculating value, D _sshift value is measured in test, and n statistics is counted.

Calculated stress measures stress relation line with test

If Stress calculation value σ _jwith test measurements σ _spass be:

σ _s＝kσ _j+Δσ………………………………………………………(1)

Wherein k is defined as consistency coefficient, and best expectation value is 1, and Δ σ is error, and establishes itself and stress intensity to have nothing to do, and Normal Distribution, then according to the principle of least square, the k making the mean square value of Δ σ minimum can be released as follows,

Stress consistency coefficient k:

(to displacement by σ _j, σ _sbe changed to D _j, D _s, lower same) ... (2)

Stress root-mean-square error ε:

$\mathrm{\ϵ}=\sqrt{\frac{\underset{i=1}{\overset{n}{\Σ}}{({\mathrm{\σ}}_{si}-{\mathrm{k\σ}}_{ji})}^2}{(n-1)}}\mathrm{...}\left(3\right)$

Stress related coefficient γ:

$\mathrm{\γ}=\frac{\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\σ}}_{ji}{\mathrm{\σ}}_{si}}{\sqrt{\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\σ}}_{ji}^2\·\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\σ}}_{si}^2}}\mathrm{...}\left(4\right)$

σ _j-σ _scurve (95% envelope):

σ _s＝kσ _j+2ε…………………………………………………………(5)

σ _s＝kσ _j-2ε…………………………………………………………(6)

Russia's curve:

$\mathrm{\ρ}~\left|{\mathrm{\σ}}_S\right|,\mathrm{\ρ}=\left|\frac{{\mathrm{\σ}}_s}{{\mathrm{\σ}}_j}\right|-\mathrm{1...}\left(7\right)$

Probability density curve (unimodal bell):

Δσ _i＝σ _si-kσ _ji………………………………………………………(8)

With for spacing, divide 14 sections: (-∞ ,-3), (-3 ,-2.5) ... (2,2.5), (2.5,3), (3, ∞), adds up the probability dropped in every interval ε: р=2m/n

N: statistics number, m: drop on the statistics in district is counted.

Fuselage tangent plane along the circumferential direction launches stress distribution

In a certain tangent plane of fuselage, the stress distribution along the circumferential direction launched, adopts and pays vertical leaf curve:

$\mathrm{\σ}=a_0+a_{1}cos\frac{2\mathrm{\π}x}{L}+b_1\sin \frac{2\mathrm{\π}x}{L}+a_{2}cos\frac{2\mathrm{\π}x}{L}+b_2\sin \frac{2\mathrm{\π}x}{L}+\mathrm{......}\left(9\right)$

Wherein: $a_0=\frac{\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\σ}}_i}{n}$

$a_k=\frac{2\underset{i=1}{\overset{n}{\Σ}}\[{\mathrm{\σ}}_{i}cos\frac{2{\mathrm{k\πx}}_i}{L}\]}{n}$

$b_k=\frac{2\underset{i=1}{\overset{n}{\Σ}}\[{\mathrm{\σ}}_{i}sin\frac{2{\mathrm{k\πx}}_i}{n}\]}{n},k=1...m$

The each measurement point of x launches coordinate

L circumference overall length

Matching square error is:

$\mathrm{\ϵ}=\sqrt{\frac{{(\mathrm{\σ}-{\mathrm{\σ}}_i)}^2}{n-2m-1}}\mathrm{...}\left(10\right)$

Fuselage is along the displacement curve approximating method in course (Y-X) or the span (Y-Z) direction

If displacement calculated value D _jicoordinate X _i, test measurements D _si, with the calculated curve of this matching be:

D _J(X)＝D _j0+D _J1X+D _J2X ²+…+D _JmX ^m…………………………(11)

D _s(X)＝D _js+D _s1X+D _s2X ²+…+D _smX ^m…………………………(12)

If function: D _s(X)=K _dd _j(X) ... (13)

Depending on D _j(X) for known, be D to the test measurements of i point _si, can obtain according to least square method:

$K_D=\frac{\underset{i=1}{\overset{n}{\Σ}}{D}_j\left(X_i\right)D_{si}}{\underset{i=1}{\overset{n}{\Σ}}{\[D_j\left(X_i\right)\]}^2}\mathrm{...}\left(14\right)$

K _dclose to 1 time, illustrate that consistance is better.

Be understandable that, above-mentioned concrete equation and the computation process schematically illustrating wherein several curve, other formula do not elaborate herein.Be understandable that, other formula are prior art, therefore repeat no more.

Be understandable that, data required in the computational analysis data in above-mentioned step 2 i.e. each curve, such as, calculated stress data etc.

Additionally provide a kind of computing system in the present invention, this computing system is used for realizing above-described embodiment and preferred implementation, has carried out namely repeating no more of explanation in above-mentioned.As used below, term " unit or module " is for can realize the software of predetermined function and/or the combination of hardware.Although the device described by following examples realizes better with software, but, by hardware, or the implementation of the combination of software and hardware also may and conceived.

Present invention also offers a kind of computing system, described computing system realizes the method calculating testpieces test figure and analyze data consistency as above, and described computing system comprises:

Database Unit, described Database Unit is used for storage function;

Call unit, described call unit is for calling the function in described Database Unit;

Function generation unit, described function generation unit is for generating test figure functional arrangement and computational analysis data function figure;

Function contrast unit, described function contrast unit is used for described test figure functional arrangement and described analysis data function figure to contrast;

Model architecture unit, described model architecture unit is used for the model of framework testpieces.

In embodiment provided by the present invention, should be understood that disclosed relevant apparatus and method can realize in other way.Such as, device described above is only schematic, such as, the division of described module and unit, be only a kind of logic function to divide, actual can have other dividing mode when realizing, such as multiple unit or assembly can in conjunction with or another system can be integrated into, or some features can be ignored, or do not perform.Another point, shown or discussed coupling each other or direct-coupling or communication connection can be by some interfaces, the indirect coupling of device or unit or communicate overlooking, can be electrically, machinery or other form.

The described unit illustrated as separating component or can may not be and physically separates, and the parts as unit display can yes or no physical location, namely can be positioned at a place, or also can be distributed in multiple network element.Some or all of unit wherein can be selected according to the actual needs to realize the object of the present embodiment scheme.

In addition, each functional unit in an embodiment of the present invention can be integrated in a processing unit, also can be that the independent physics of unit exists, also can two or more unit in a unit integrated.Above-mentioned integrated unit both can adopt the form of hardware to realize, and the form of SFU software functional unit also can be adopted to realize.

If described integrated unit using the form of SFU software functional unit realize and as independently production marketing or use time, can be stored in a computer read/write memory medium, based on such understanding, the part that technical scheme of the present invention contributes to prior art in essence in other words or all or part of of this technical scheme can embody with the form of software product, this computer software product is stored in a storage medium, comprises all or part of step of some instructions in order to make computer processor perform method described in each embodiment of the present invention.And aforesaid storage medium comprise USB flash disk, portable hard drive, ROM (read-only memory), random access memory, magnetic disc or CD etc. various can be program code stored medium.

Finally it is to be noted: above embodiment only in order to technical scheme of the present invention to be described, is not intended to limit.Although with reference to previous embodiment to invention has been detailed description, those of ordinary skill in the art is to be understood that: it still can be modified to the technical scheme described in foregoing embodiments, or carries out equivalent replacement to wherein portion of techniques feature; And these amendments or replacement, do not make the essence of appropriate technical solution depart from the spirit and scope of various embodiments of the present invention technical scheme.

Claims (10)

1. calculate testpieces test figure and the method analyzing data consistency, for calculating aircraft or aircraft components, it is characterized in that, described method comprises the steps:

Step 1: the testing site on optional test part, carries out the slow test of testpieces under predetermined load condition, and obtains the test figure of testing site;

Step 2: obtain the computational analysis data of described testpieces under described predetermined load condition with calculation and analysis methods;

Step 3: obtain the position relationship of the described testing site in described step 1 relative to the described testpieces in described step 2, and mapping table is set up in the described testing site in described step 1 and the testpieces in described step 2;

Step 4: by calculating test figure and analyzing the relative coefficient of data, mean square deviation and relative deviation amount etc., the correlation curve figure of draw calculation and test.

2. the method calculating testpieces test figure and analyze data consistency as claimed in claim 1, it is characterized in that, the testing site on described optional test part specifically comprises: displacement measurement point and stress measurement point.

3. the method calculating testpieces test figure and analyze data consistency as claimed in claim 2, it is characterized in that, the method for selecting of institute's displacement measurement point is: the point of intersection selecting described each position of testpieces;

The method for selecting of described stress measurement point is: be limited to and be selected in unidirectional tension and compression and plane stressed zone and/or be selected in STRESS VARIATION gradient shoulder.

4. the method calculating testpieces test figure and analyze data consistency as claimed in claim 1, it is characterized in that, described predetermined load condition is: the load capacity choosing 67% of the design load of described testpieces carries out slow test.

5. the method calculating testpieces test figure and analyze data consistency as claimed in claim 1, it is characterized in that, described calculation and analysis methods comprises analytical method and finite element method.

6. the method calculating testpieces test figure and analyze data consistency as claimed in claim 5, it is characterized in that, described finite element method comprises: the finite element model setting up testpieces, and obtains the computational analysis data of described testpieces under described predetermined load condition; Described testing site in the finite element model of described testpieces is obtained by described step 3.

7. the method calculating testpieces test figure and analyze data consistency as claimed in claim 6, it is characterized in that, the mapping table in described step 3 at least comprises following data:

The corresponding test period in described finite element model in testing site in described step 1, element numerals, element type, testing site coordinate, material behavior corresponding to testing site.

8. the method calculating testpieces test figure and analyze data consistency as claimed in claim 1, it is characterized in that, described testing site is multiple.

9. the method calculating testpieces test figure and analyze data consistency as claimed in claim 1, it is characterized in that, described test figure functional arrangement and described computational analysis data function figure comprise σ _j~ σ _scurve, Δ σ-σ _scurve, Russian curve ρ ~ | σ _s|, σ-X launches the multiple curve such as distribution, P-Δ/ε probability density distribution curve, and displacement DJ ~ DS curve displacement curve, the figures such as D-X displacement measurement curve.

10. a computing system, described computing system is used for calculating testpieces test figure and the method analyzing data consistency as in one of claimed in any of claims 1 to 9, and it is characterized in that, described computing system comprises:

Database Unit, described Database Unit is used for storage function;

Call unit, described call unit is for calling the function in described Database Unit;

Function generation unit, described function generation unit is for generating test figure functional arrangement and computational analysis data function figure;

Function contrast unit, described function contrast unit is used for described test figure functional arrangement and described analysis data function figure to contrast;

Model architecture unit, described model architecture unit is used for the model of framework testpieces.