CN106156386A - A kind of slow test for housing reinforced structure and predicting method - Google Patents
A kind of slow test for housing reinforced structure and predicting method Download PDFInfo
- Publication number
- CN106156386A CN106156386A CN201510169408.6A CN201510169408A CN106156386A CN 106156386 A CN106156386 A CN 106156386A CN 201510169408 A CN201510169408 A CN 201510169408A CN 106156386 A CN106156386 A CN 106156386A
- Authority
- CN
- China
- Prior art keywords
- measuring point
- test
- finite element
- file
- result
- 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
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
This technology belongs to slow test and indication field, is specifically related to a kind of slow test for housing reinforced structure and predicting method.By structure static(al) indication simulation result under different loads effect, different boundary processing mode is contrasted with test data.On the one hand, solve slow test and the systematicness contrast of static(al) indication emulation, comprehensive assessment correlated error, solve individual characteristics point small data in pole to be analyzed control methods in the past;On the other hand, by relative analysis under these different stages, accumulation different boundary analog form, difference solve mode, different tests method etc. to actual experimental or simulation calculation factor of influence data.Another further aspect, by contrast cloud atlas, curve, similarity can reaction test intuitively and indication between relation, by the assessment of above-mentioned data be may decide that test whether proceed, with or without testpieces or personal safety as well as the property safety problem.
Description
Technical field
This technology belongs to slow test and indication field, is specifically related to a kind of static(al) for housing reinforced structure
Test and predicting method.
Background technology
Carrier rocket structural strength and rigidity in development process are important examination objects, and we mainly lead at present
Cross 3 kinds of means to analyze: 1, Engineering Algorithm analysis;2, finite element indication simulation analysis;3, entity static(al)
Test examination.Engineering Algorithm analysis is mainly used in the design initial stage and assesses structure macro strength and rigidity, its
Feature is that application is simple, low cost, calculating time are few, almost without visualization display, low precision, analysis
Result general trend is correct, but, product this for carrier rocket, it is desirable to structural units's quality high-mechanic
Ratio, has increasing need for designing fine structureization, and Engineering Algorithm cannot meet subsequent request.And for limited
Unit's simulation calculation is in addition to having low cost, time short feature, it is also possible to visual Simulation structure is at ring
The variation tendency of the numerical value such as full structural stress, strain, displacement, acceleration and magnitude under the load effect of border.But
It is, during simulation calculation, due to the simplification of model, the constraint on border, the equivalence of material, analysis
Theoretical difference causes the feasibility of simulation result to be worth assessment.Can be true for true slow test
The deformation such as model configuration malformation under bearing environmental load, destruction, unstability.But, owing to it is actual
Operating process can run into following situation: 1, cost is high;2, measurement data is limited, it is impossible to so closing
Heart point measures;3, measurement space is limited, to measure sensor cannot layout area (surface curvature is too big,
Little closing space) it is to obtain measurement data;4, when load applying or boundary constraint, inevitably
Need other auxiliary mould, the usual rigidity of these frocks and true environment lower boundary rigidity by part variation, but
It is that intensity, the rigidity regional area of structure are affected greatly, can not ignore by border factor;5, other environment is former
Because of inevitable, such as temperature Change, the factor of wind load change.
Summary of the invention
It is an object of the invention to: a kind of slow test for housing reinforced structure and predicting method are provided,
For the design of aerospace craft with analyze and provide stronger support.
Technical scheme is as follows: a kind of slow test for housing reinforced structure and predicting method,
Comprise the following steps:
Step one: set up FEM (finite element) model, according to the physical dimension of model, material, boundary condition, load
Set up corresponding FEM (finite element) model, and computational analysis, obtain finite element result file;
Step 2: according to measuring point space coordinates XYZ, obtains correspondence position finite element node serial number, thus obtains
Measuring point information, measuring point information is made up of two files: finite element node serial number file, finite element result file,
Step 3: by measuring point space coordinates XYZ, obtain the structure determination numbering of correspondence position;Thus obtain
Test measuring point information, test measuring point information is made up of two files: test structure numbering file, result of the test
File,
Step 4: search and test measuring point geometric position distance in the FEM (finite element) model set up in step one
Near node, then completes one-to-one relationship by the degree of freedom of the degree of freedom of measuring point Yu finite element node;
Step 5: by the finite element result file in calculation procedure two and the result of the test file in step 3
Modal assurance criterion MAC, checking finite element indication is the most consistent, wherein with test deformation tendency:
In formula, ΦeRepresent test data matrix, ΦaRepresent indication emulation FEM calculation matrix, wherein, MAC value
In the range of 0-1, as MAC=0, represent that finite element indication is entirely different with test deformation tendency, work as MAC=1
Time, represent that both deformation tendency are identical;
Step 6: draw comparison diagram
Step 7: FEM updating
Finite element result in step one is carried out the MAC value in step 5 with the result of the test in step 3
Contrast, then contrasted by the comparison diagram drawn in step 6, revised by comparison diagram version
FEM (finite element) model, thus reach MAC as close possible to 1, makes result of calculation and result of the test similarity
Greatly.
In described step 2, the sequence that measuring point is carried out by finite element node serial number file according to measuring point order, for
NX1 matrix, exports ess-strain displacement by finite element node serial number file, and finite element result file is by nX3
Matrix description, wherein n is measure-point amount, and every string is stated by three data, and first is measuring point x coordinate,
Second is y-coordinate, and the 3rd is measuring point stress, strain or displacement, for less than 3 row or the survey of 3 row
Point, if a curve can be plotted as under space, then measuring point information is by 2Xn matrix description, and first is classified as
Measuring point is along the distance of curve, and second is classified as measuring point stress, strain or displacement size;If can not paint under space
It is made as a curve, then terminates.
In described step 3, the sequence that measuring point is carried out by test structure numbering file according to measuring point order, for nX1
Matrix, exports ess-strain displacement by test structure numbering file, and result of the test file is retouched by nX3 matrix
Stating, wherein n is measure-point amount, and every string is stated by three data, and first is measuring point x coordinate, second
Individual for y-coordinate, the 3rd is measuring point stress, strain or displacement, for less than 3 row or the measuring points of 3 row,
If a curve can be plotted as under space, then test measuring point information is by 2Xn matrix description, and first is classified as
Measuring point is along the distance of curve, and second is classified as measuring point stress, strain or displacement size;If can not paint under space
It is made as a curve, then terminates.
In described step 6, when being depicted as 3 row, 3 row and above mesh shape on measuring point direction in space,
Selecting cloud atlas to draw comparison diagram, wherein, x coordinate is spatially transverse size;Y-coordinate is axial dimension;Cloud atlas
Color table is shown as stress value size, when only 1,2 row or 1,2 row on measuring point direction in space, selects road
Footpath correlation curve draws comparison diagram, and wherein, x-axis is path relative distance, and y-axis is the stress of measuring point, answers
Become or displacement.When contrasting a certain characteristic point test data, select and load rank or load as x-axis, y-axis
For measuring point stress, strain or displacement.
The remarkable result of the present invention is: by structure in different loads effect, different boundary processing mode
Lower static(al) indication simulation result contrasts with test data.On the one hand, slow test and static(al) indication emulation are solved
Systematicness contrast, comprehensive assessment correlated error, solve in the past can only be to pole individual characteristics point small data minute
Analysis control methods;On the other hand, by relative analysis under these different stages, accumulate different boundary simulation side
Formula, difference solve mode, different tests method etc. to actual experimental or simulation calculation factor of influence data.Again
On the one hand, by contrast cloud atlas, curve, similarity can reaction test intuitively and indication between relation,
By the assessment of above-mentioned data being may decide that whether test proceeds, with or without testpieces or person property peace
Full problem.
Detailed description of the invention
A kind of slow test for housing reinforced structure and predicting method, comprise the following steps
Step one: set up FEM (finite element) model, according to the physical dimension of model, material, boundary condition, load
Set up corresponding FEM (finite element) model, and computational analysis, obtain finite element result file;
Step 2: according to measuring point space coordinates XYZ, obtains correspondence position finite element node serial number, thus obtains
Measuring point information, measuring point information is made up of two files: finite element node serial number file, finite element result file,
The sequence that measuring point is carried out by finite element node serial number file according to measuring point order, for nX1 matrix, by limited
Unit node serial number file output ess-strain displacement, finite element result file by nX3 matrix description, wherein n
For measure-point amount, every string is stated by three data, and first is measuring point x coordinate, and second is y-coordinate,
3rd is measuring point stress, strain or displacement, for less than 3 row or the measuring point of 3 row, if can be in space
Under be plotted as a curve, then measuring point information is by 2Xn matrix description, and first is classified as the measuring point distance along curve,
Second is classified as measuring point stress, strain or displacement size;If a curve can not be plotted as under space, then
Terminate;
Step 3: by measuring point space coordinates XYZ, obtain the structure determination numbering of correspondence position;Thus obtain
Test measuring point information, test measuring point information is made up of two files: test structure numbering file, result of the test
File.The sequence that measuring point is carried out by test structure numbering file according to measuring point order, for nX1 matrix, passes through
Test structure numbering file output ess-strain displacement, result of the test file by nX3 matrix description, wherein n
For measure-point amount, every string is stated by three data, and first is measuring point x coordinate, and second is y-coordinate,
3rd is measuring point stress, strain or displacement.For less than 3 row or the measuring point of 3 row, if can be in space
Under be plotted as a curve, then test measuring point information is by 2Xn matrix description, and first is classified as measuring point along curve
Distance, second is classified as measuring point stress, strain or displacement size;If a curve can not be plotted as under space,
Then terminate;
Step 4: search and test measuring point geometric position distance in the FEM (finite element) model set up in step one
Near node, then completes one-to-one relationship by the degree of freedom of the degree of freedom of measuring point Yu finite element node;
Step 5: by the finite element result file in calculation procedure two and the result of the test file in step 3
Modal assurance criterion MAC, checking finite element indication is the most consistent, wherein with test deformation tendency:
In formula, ΦeRepresent test data matrix, ΦaRepresent indication emulation FEM calculation matrix.Wherein, MAC value
In the range of 0-1, as MAC=0, represent that finite element indication is entirely different with test deformation tendency, work as MAC=1
Time, represent that both deformation tendency are identical;
Step 6: draw comparison diagram
When 3 row, 3 row and above mesh shape can be depicted as on measuring point direction in space, cloud atlas is selected to draw
Comparison diagram.Wherein, x coordinate is spatially transverse size;Y-coordinate is axial dimension;Cloud atlas color table is shown as should
Force value size.When only 1,2 row or 1,2 row on measuring point direction in space, path correlation curve is selected
Drawing comparison diagram, wherein, x-axis is path relative distance, and y-axis is the stress of measuring point, strain or displacement.?
When contrasting a certain characteristic point test data, select and load rank or load as x-axis, y-axis be measuring point stress,
Strain or displacement;
Step 7: FEM updating
Finite element result in step one is carried out the MAC value in step 5 with the result of the test in step 3
Contrast, then contrasted by the comparison diagram drawn in step 6, revised by comparison diagram version
FEM (finite element) model, so that MAC value is close to 1, makes result of calculation close to result of the test.
Claims (4)
1. the slow test for housing reinforced structure and predicting method, it is characterised in that: include following
Step:
Step one: set up FEM (finite element) model, according to the physical dimension of model, material, boundary condition, load
Set up corresponding FEM (finite element) model, and computational analysis, obtain finite element result file;
Step 2: according to measuring point space coordinates XYZ, obtains correspondence position finite element node serial number, thus obtains
Measuring point information, measuring point information is made up of two files: finite element node serial number file, finite element result file,
Step 3: by measuring point space coordinates XYZ, obtain the structure determination numbering of correspondence position;Thus obtain
Test measuring point information, test measuring point information is made up of two files: test structure numbering file, result of the test
File,
Step 4: search and test measuring point geometric position distance in the FEM (finite element) model set up in step one
Near node, then completes one-to-one relationship by the degree of freedom of the degree of freedom of measuring point Yu finite element node;
Step 5: by the finite element result file in calculation procedure two and the result of the test file in step 3
Modal assurance criterion MAC, checking finite element indication is the most consistent, wherein with test deformation tendency:
In formula, ΦeRepresent test data matrix, ΦaRepresent indication emulation FEM calculation matrix, wherein, MAC value
In the range of 0-1, as MAC=0, represent that finite element indication is entirely different with test deformation tendency, work as MAC=1
Time, represent that both deformation tendency are identical;
Step 6: draw comparison diagram
Step 7: FEM updating
Finite element result in step one is carried out the MAC value in step 5 with the result of the test in step 3
Contrast, then contrasted by the comparison diagram drawn in step 6, revised by comparison diagram version
FEM (finite element) model, so that MAC value is close to 1, makes result of calculation close to result of the test.
A kind of slow test for housing reinforced structure the most according to claim 1 and predicting method,
It is characterized in that: in described step 2, measuring point is carried out by finite element node serial number file according to measuring point order
Sequence, for nX1 matrix, exports ess-strain displacement, finite element result by finite element node serial number file
File is by nX3 matrix description, and wherein n is measure-point amount, and every string is stated by three data, and first is
Measuring point x coordinate, second is y-coordinate, and the 3rd is measuring point stress, strain or displacement, for less than 3
Row or the measuring point of 3 row, if a curve can be plotted as under space, then measuring point information is retouched by 2Xn matrix
Stating, first is classified as the measuring point distance along curve, and second is classified as measuring point stress, strain or displacement size;If no
A curve can be plotted as under space, then terminate.
A kind of slow test for housing reinforced structure the most according to claim 1 and predicting method,
It is characterized in that: in described step 3, the row that measuring point is carried out by test structure numbering file according to measuring point order
Sequence, for nX1 matrix, exports ess-strain displacement by test structure numbering file, result of the test file by
NX3 matrix description, wherein n is measure-point amount, and every string is stated by three data, sits for measuring point x for first
Mark, second is y-coordinate, and the 3rd is measuring point stress, strain or displacement, for less than 3 row or 3 row
Measuring point, if a curve can be plotted as under space, then test measuring point information by 2Xn matrix description,
First is classified as the measuring point distance along curve, and second is classified as measuring point stress, strain or displacement size;If can not be
It is plotted as a curve under space, then terminates.
A kind of slow test for housing reinforced structure the most according to claim 1 and predicting method,
It is characterized in that: in described step 6, when can be depicted as on measuring point direction in space 3 row, 3 row and above
Mesh shape, selects cloud atlas to draw comparison diagram, and wherein, x coordinate is spatially transverse size;Y-coordinate is axial
Size;Cloud atlas color table is shown as stress value size, when only 1,2 row or 1,2 row on measuring point direction in space
Time, select path correlation curve to draw comparison diagram, wherein, x-axis is path relative distance, and y-axis is measuring point
Stress, strain or displacement.When contrasting a certain characteristic point test data, select and load rank or load work
For x-axis, y-axis is measuring point stress, strain or displacement.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510169408.6A CN106156386B (en) | 2015-04-10 | 2015-04-10 | A kind of slow test and predicting method for shell reinforced structure |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510169408.6A CN106156386B (en) | 2015-04-10 | 2015-04-10 | A kind of slow test and predicting method for shell reinforced structure |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106156386A true CN106156386A (en) | 2016-11-23 |
CN106156386B CN106156386B (en) | 2019-09-13 |
Family
ID=57336664
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510169408.6A Active CN106156386B (en) | 2015-04-10 | 2015-04-10 | A kind of slow test and predicting method for shell reinforced structure |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106156386B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108256214A (en) * | 2018-01-16 | 2018-07-06 | 滨州学院 | A kind of orthogonal stiffened panel calculating method of stiffness of aviation and device |
CN108334653A (en) * | 2017-10-20 | 2018-07-27 | 北京空天技术研究所 | A kind of Static Analysis Model of Micro-machined modification method, update the system and Static Strength Analysis method |
CN109359360A (en) * | 2018-09-30 | 2019-02-19 | 国家超级计算天津中心 | A kind of structural stress processing method based on local feature |
CN110717273A (en) * | 2019-10-11 | 2020-01-21 | 内蒙古第一机械集团股份有限公司 | Technological process simulation boundary condition construction method |
CN111639457A (en) * | 2020-06-23 | 2020-09-08 | 广州电力机车有限公司 | Design method and test method of electric wheel test tool of mining dump truck |
CN112067265A (en) * | 2020-08-19 | 2020-12-11 | 上海航天精密机械研究所 | Testing device suitable for static force of boundary capable of sliding in one direction |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020132068A1 (en) * | 2001-03-16 | 2002-09-19 | Onanov Herman George | Superstability of thin-walled structures |
CN104239639A (en) * | 2014-09-18 | 2014-12-24 | 上海理工大学 | Method for simplifying parameter modification of topological optimization model of transmission shell |
CN104484526A (en) * | 2014-12-16 | 2015-04-01 | 中国第一汽车股份有限公司 | Method for improving finite element analysis accuracy of transmission case |
-
2015
- 2015-04-10 CN CN201510169408.6A patent/CN106156386B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020132068A1 (en) * | 2001-03-16 | 2002-09-19 | Onanov Herman George | Superstability of thin-walled structures |
CN104239639A (en) * | 2014-09-18 | 2014-12-24 | 上海理工大学 | Method for simplifying parameter modification of topological optimization model of transmission shell |
CN104484526A (en) * | 2014-12-16 | 2015-04-01 | 中国第一汽车股份有限公司 | Method for improving finite element analysis accuracy of transmission case |
Non-Patent Citations (3)
Title |
---|
李佰灵等: "基于响应面方法的多目标有限元模型修正技术研究", 《强度与环境》 * |
李双等: "某飞行器三角翼动力学有限元模型与试验数据相关分析原理及实现", 《强度与环境》 * |
郝鹏: "面向新一代运载火箭的网格加筋柱壳结构优化研究", 《中国博士学位论文全文数据库工程科技II辑》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108334653A (en) * | 2017-10-20 | 2018-07-27 | 北京空天技术研究所 | A kind of Static Analysis Model of Micro-machined modification method, update the system and Static Strength Analysis method |
CN108256214A (en) * | 2018-01-16 | 2018-07-06 | 滨州学院 | A kind of orthogonal stiffened panel calculating method of stiffness of aviation and device |
CN109359360A (en) * | 2018-09-30 | 2019-02-19 | 国家超级计算天津中心 | A kind of structural stress processing method based on local feature |
CN109359360B (en) * | 2018-09-30 | 2022-11-11 | 国家超级计算天津中心 | Structural stress processing method based on local features |
CN110717273A (en) * | 2019-10-11 | 2020-01-21 | 内蒙古第一机械集团股份有限公司 | Technological process simulation boundary condition construction method |
CN110717273B (en) * | 2019-10-11 | 2023-03-17 | 内蒙古第一机械集团股份有限公司 | Technological process simulation boundary condition construction method |
CN111639457A (en) * | 2020-06-23 | 2020-09-08 | 广州电力机车有限公司 | Design method and test method of electric wheel test tool of mining dump truck |
CN111639457B (en) * | 2020-06-23 | 2023-06-09 | 广州电力机车有限公司 | Design method and test method of electric wheel test fixture of mining dump truck |
CN112067265A (en) * | 2020-08-19 | 2020-12-11 | 上海航天精密机械研究所 | Testing device suitable for static force of boundary capable of sliding in one direction |
Also Published As
Publication number | Publication date |
---|---|
CN106156386B (en) | 2019-09-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106156386A (en) | A kind of slow test for housing reinforced structure and predicting method | |
Miller et al. | Understanding local deformation in metallic polycrystals using high energy X-rays and finite elements | |
CN103267683B (en) | Method for determining remaining life of heat-resisting metal material | |
Niesłony et al. | A study of compatibility between two classical fatigue curve models based on some selected structural materials | |
JP6655105B2 (en) | Manufacturing auxiliary equipment | |
Kowalczyk | Identification of mechanical parameters of composites in tensile tests using mixed numerical-experimental method | |
CN110807285B (en) | Engineering transport vehicle frame fatigue analysis method, system and device | |
US20080294397A1 (en) | Analytical method, recording medium, and analyzing apparatus | |
Buljak et al. | Assessment of residual stresses and mechanical characterization of materials by “hole drilling” and indentation tests combined and by inverse analysis | |
Spagnuolo et al. | Mesoscale modeling and experimental analyses for pantographic cells: effect of hinge deformation | |
Lyngdoh et al. | Bayesian updating of structural model with a conditionally heteroscedastic error distribution | |
Calvin et al. | A new effective stress intensity factor approach to determine thickness-independent fatigue crack growth rate curves | |
Firat et al. | Improving the accuracy of stamping analyses including springback deformations | |
Pavlov et al. | Development of mass and cost efficient grid-stiffened and lattice structures for space applications | |
Carrere et al. | Efficient structural computations with parameters uncertainty for composite applications | |
Morse et al. | Practical conformance evaluation in the measurement of flexible parts | |
Yu et al. | Measurement of the M-integral for a hole in an aluminum plate or strip | |
CN113962134A (en) | Strain monitoring method for linear superposition of condition generation type countermeasure network and load strain | |
Lipina et al. | Implementation of polycarbonate material mechanical properties of rapid prototyping into system CREO, laboratory verification of the results | |
Cosseron et al. | Optimized gauging for tire–rim loading identification | |
US6935186B2 (en) | Fatigue safety factor testing method and fatigue safety factor testing apparatus | |
CN105426660B (en) | A kind of Forecasting Methodology of machine components table planar residual stress field | |
JP4631207B2 (en) | Deformed shape calculation device and program | |
dos Reis | Constitutive parameters identification of metal sheets using circular bulge tests | |
Wankhede et al. | SPSA: An image processing based software for single point strain analysis |
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 |