CN102682168A - Method for determining dimension of rectangular beam section of low-speed flutter main beam model of airplane - Google Patents

Method for determining dimension of rectangular beam section of low-speed flutter main beam model of airplane Download PDF

Info

Publication number
CN102682168A
CN102682168A CN2012101462225A CN201210146222A CN102682168A CN 102682168 A CN102682168 A CN 102682168A CN 2012101462225 A CN2012101462225 A CN 2012101462225A CN 201210146222 A CN201210146222 A CN 201210146222A CN 102682168 A CN102682168 A CN 102682168A
Authority
CN
China
Prior art keywords
rectangular beam
section
cross
width
inertia
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
Application number
CN2012101462225A
Other languages
Chinese (zh)
Other versions
CN102682168B (en
Inventor
罗务揆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Aircraft Design and Research Institute of AVIC
Original Assignee
Xian Aircraft Design and Research Institute of AVIC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Aircraft Design and Research Institute of AVIC filed Critical Xian Aircraft Design and Research Institute of AVIC
Priority to CN201210146222.5A priority Critical patent/CN102682168B/en
Publication of CN102682168A publication Critical patent/CN102682168A/en
Application granted granted Critical
Publication of CN102682168B publication Critical patent/CN102682168B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention belongs to the field of structural mechanics and relates to a method for determining the dimension of a rectangular beam section of a low-speed flutter main beam model of an airplane. The method is characterized by adapting to the situation that the half width a of a rectangular section is smaller than the half height b of the rectangular section or the situation that J/Ix<1.69, wherein Ix refers to vertical inertia moment and J refers to polar inertia moment. The method for determining the dimension of the rectangular beam section with a lug piece comprises the following steps of: calculating the half width ar and the half height br of the rectangular section when the rectangular beam without the lug piece has the vertical inertia moment Ix of a preset value and the polar inertia moment J of a preset value; adjusting the polar inertia moment J with the preset value; calculating the half width a and the half height b; and calculating the half width 1 of the rectangular beam section with the lug piece when the rectangular beam with the lug piece has lateral inertia moment Iy of a preset value and the thickness of the lug piece is t. The method disclosed by the invention has the advantages of increasing the rigidity precision of the section of the model, reducing the uncertainty of model design, shortening the time for determining the dimension of the section, and increasing the design efficiency of the flutter model.

Description

A kind of aircraft low speed flutter girder model rectangular beam sectional dimension is confirmed method
Technical field
The invention belongs to the structural mechanics field, relate to a kind of aircraft low speed flutter girder model rectangular beam sectional dimension and confirm method.
Background technology
For the airplane flutter modelling of high aspect ratio wing, single-beam structure of models form is widely used in structures such as wing, fuselage, empennage.In order to satisfy the designing requirement of three-way rigidity, the square-section of band auricle is an ideal designs form.
Usually need obtain the sectional dimension that meets the demands through the method that examination is gathered in the past.The method that examination is gathered has following shortcoming: the first, the adjustment sectional dimension data needs designer's experience; The dimensional data that provides by rule of thumb often error is very big; Even the situation of the size that meets design requirement can appear being difficult to adjust; Influence the precision of model section rigidity, increased the uncertainty of modelling.The second, to confirm sectional dimension through method of trial and error, need carry out the manual work adjustment of many rounds and judge that the time is long that efficient is low, has a strong impact on the modelling cycle.
Referring to Chinese patent " band auricle rectangular beam sectional dimension confirm method " (application number 201110232655.8), need not examination and gather, only need provide and comprise the auricle thickness t, vertically to moment of inertia I x, side direction moment of inertia I yWith polar moment of inertia J in interior input data, can obtain sectional dimension data, comprise width 2a, height 2b and the cross-sectional width 2l of rectangle, but the applicable elements of this method is J/I x>=1.69.Work as J/I x<1.69, design error can increase.
Summary of the invention
The objective of the invention is: propose a kind of aircraft low speed flutter girder model rectangular beam sectional dimension and confirm method; To improve the precision of model section rigidity, reduce the uncertainty of modelling, shorten the time of confirming sectional dimension; Improve the design efficiency of flutter model, and be applicable to cross section J/I x<1.69 o'clock requirement.
Technical solution of the present invention is: aircraft low speed flutter girder model rectangular beam cross section is the rectangular beam of band auricle, guarantee band auricle rectangular beam have predetermined value vertically to moment of inertia I x, side direction moment of inertia I yWith polar moment of inertia J and auricle thickness be under the situation of t, confirm width 2a, height 2b and the cross-sectional width 2l of its rectangle, it is characterized in that, confirm that the step of band auricle rectangular beam sectional dimension is following:
1, calculate not with the auricle rectangular beam have predetermined value vertically to moment of inertia I xSquare-section half width a during with polar moment of inertia J rWith cross section half height b r:
1.1, set the initial value α of iteration variable α 0=0.93 ~ 0.99:
1.2, according to computes intermediate variable β and iteration variable α 1:
24 ( 2 &pi; ) 5 &alpha; 0 &beta; 3 - 4 &beta; 2 + J I x = 0 &alpha; 1 = tanh &pi; 2 &beta; + 1 243 tanh 3 &pi; 2 &beta; . . . [ 1 ]
Attention: when finding the solution β, relate to finding the solution of simple cubic equation, the flat-sawn method is adopted in suggestion, and β ∈ (0,1);
1.3, calculate iterative value α 1With initial value α 0Scale error e 1:
e 1 = &alpha; 1 - &alpha; 0 &alpha; 0 . . . [ 2 ]
If | e 1|<0.000001, square-section half width a then rWith cross section half height b rFor:
a r = 3 I x &beta; 3 4 4 b r = 3 I x 4 &beta; 4 . . . [ 3 ]
Step 1.1 finishes; Otherwise, carry out step 1.4;
1.4, make α 01, the method for repeating step 1.2~step 1.3, through iterative computation, up to | e 1|<0.000001, and obtain square-section half width a at this moment rWith cross section half height b r, step 1 finishes;
2, predetermined polar moment of inertia J is adjusted: go out the adjusted value J of polar moment of inertia J according to computes 1,
Figure BDA00001628974700024
Note answering guaranteed conditions: t<a r;
3, calculate half width a and half height b: change the J in the formula [1] into J 1, calculate half width a and half height b according to the described method of step 1 then;
4, calculate band auricle rectangular beam cross section half-breadth l:
4.1, according to computes intermediate variable δ:
&delta; = 3 I y - 4 a 3 b 2 t . . . [ 5 ]
4.2, calculate band auricle rectangular beam cross section half-breadth l:
l = &delta; + a 3 3 . . . [ 6 ]
So far, obtain half width a, half height b and cross section half-breadth l, and then obtain width 2a, height 2b and the cross-sectional width 2l of corresponding rectangular beam with the auricle rectangular beam.
Advantage of the present invention is: improved the precision of model section rigidity, reduced the uncertainty of modelling, shortened the time of definite sectional dimension, improved the design efficiency of flutter model, and satisfied J/I x<1.69 requirement.One embodiment of the present of invention are compared with the present examination method of gathering, and the present invention confirms that the time of sectional dimension is merely 3 minutes, and the present examination method of gathering needs 4 hours, and the time spent of the present invention is merely examination and gathers 1/80th of method.
Description of drawings
Fig. 1 is a schematic cross section of not being with the rectangular beam of auricle.A among the figure rBe the half width of the rectangle that calculates of step 1 of the present invention, b rIt is half height of the rectangle that calculates of step 1 of the present invention.The initial point 0 of the two-dimensional coordinate system among Fig. 1 is the center of rectangle, and the x axle is parallel to the Width of rectangle, and positive dirction is towards the right side, and the positive dirction of y axle up.
Fig. 2 is the schematic cross section of the rectangular beam of band auricle.Two-dimensional coordinate system among Fig. 2 is identical with Fig. 1.
Embodiment
Explain further details in the face of the present invention down.Referring to Fig. 1,2, band auricle rectangular beam sectional dimension is confirmed method, guarantee band auricle rectangular beam have predetermined value vertically to moment of inertia I x, side direction moment of inertia I yWith polar moment of inertia J and auricle thickness be under the situation of t, confirm width 2a, height 2b and the cross-sectional width 2l of its rectangle, it is characterized in that, confirm that the step of band auricle rectangular beam sectional dimension is following:
1, calculate not with the auricle rectangular beam have predetermined value vertically to moment of inertia I xSquare-section half width a during with polar moment of inertia J rWith cross section half height b r:
1.1, set the initial value α of iteration variable α 0=0.93 ~ 0.99:
1.2, according to computes intermediate variable β and iteration variable α 1:
24 ( 2 &pi; ) 5 &alpha; 0 &beta; 3 - 4 &beta; 2 + J I x = 0 &alpha; 1 = tanh &pi; 2 &beta; + 1 243 tanh 3 &pi; 2 &beta; . . . [ 1 ]
Attention: when finding the solution β, relate to finding the solution of simple cubic equation, the flat-sawn method is adopted in suggestion, and β ∈ (0,1);
1.3, calculate iterative value α 1With initial value α 0Scale error e 1:
e 1 = &alpha; 1 - &alpha; 0 &alpha; 0 . . . [ 2 ]
If | e 1|<0.000001, square-section half width a then rWith cross section half height b rFor:
a r = 3 I x &beta; 3 4 4 b r = 3 I x 4 &beta; 4 . . . [ 3 ]
Step 1 finishes; Otherwise, carry out step 1.4;
1.4, make α 01, the method for repeating step 1.2~step 1.3, through iterative computation, up to | e 1|<0.000001, and obtain square-section half width a at this moment rWith cross section half height b r, step 1 finishes;
2, predetermined polar moment of inertia J is adjusted: go out the adjusted value J of polar moment of inertia J according to computes 1,
Figure BDA00001628974700034
Note answering guaranteed conditions: t<a r;
3, calculate half width a and half height b: change the J in the formula [1] into J 1, calculate half width a and half height b according to the described method of step 1 then;
4, calculate band auricle rectangular beam cross section half-breadth l:
4.1, according to computes intermediate variable δ:
&delta; = 3 I y - 4 a 3 b 2 t . . . [ 5 ]
4.2, calculate band auricle rectangular beam cross section half-breadth l:
l = &delta; + a 3 3 . . . [ 6 ]
So far, obtain half width a, half height b and cross section half-breadth l, and then obtain width 2a, height 2b and the cross-sectional width 2l of corresponding rectangular beam with the auricle rectangular beam.
Principle of work of the present invention is: derive and the parameter correction through the mechanics of materials and Elasticity fundamental formular; (seeing the logical Higher Education Publishing House 1997 of " Elasticity " Yang Gui and " mechanics of materials " single brightness ancestral Higher Education Publishing House 1999) obtained a kind ofly directly to obtain the method with auricle rectangular beam sectional dimension from cross section property; For definite sectional dimension, this is a kind of reverse mentality of designing and method, therefore; Than need try to gather the method that obtains sectional dimension through manual work in the past; Efficient and precision have obtained improving greatly, and than the reverse engineer method that proposes before, this method more is applicable to J/I x<1.69 conditions.
Embodiment
To the method for the invention, calculate checking.
Provide three groups of aircraft low speed flutter girder model band auricle rectangular beam cross sections, P1~P3 cross section is respectively embodiment 1~embodiment 3.Table 1 has provided the geometrical property predetermined value of three embodiment, i.e. desired value, and desired value J/I x<1.69.For P1, make t=2.0mm, for P2 and P3, make t=1.0mm, carry out the sectional dimension design.Table 1 has also provided the geometrical property of using the design section that the inventive method obtains, the i.e. error of design load, and design load.Table 2 has provided the size design value of three embodiment.
The design load of the cross section geometric characteristic of table 1 is the sectional dimension corresponding according to table 2, is calculated by FEMAP v9.31.The desired value of comparing, the error of design section characteristic value all is not more than 2%, says from engineering viewpoint, and this is a kind of high-precision result of calculation.
Table 1 cross section geometric characteristic, unit: mm 4
I x I y J
The P1 desired value 227635.2 213104.0 18696.0
The P1 design load 227761.0 213117.0 18535.4
The P1 error amount 0.06% 0.01% -0.86%
The P2 desired value 3528.3 52734.5 5835.9
The P2 design load 3531.3 52745.6 5892.3
The P2 error amount 0.09% 0.02% 0.97%
The P3 desired value 2276.35 2131.04 236.96
The P3 design load 2282.3 2132.4 232.7
The P3 error amount 0.26% 0.06% -1.82%
Table 2 cross section geometry design load, unit: mm.
t 2a 2b 2l
P1 2.00 9.67 65.62 107.72
P2 1.00 14.08 14.43 84.12
P3 1.00 3.30 20.24 29.19

Claims (1)

1. an aircraft low speed flutter girder model rectangular beam sectional dimension is confirmed method, and aircraft low speed flutter girder model rectangular beam cross section is the rectangular beam of band auricle, guarantee band auricle rectangular beam have predetermined value vertically to moment of inertia I x, side direction moment of inertia I yWith polar moment of inertia J and auricle thickness be under the situation of t, and satisfy J/I x<1.69, confirm to it is characterized in that width 2a, height 2b and the cross-sectional width 2l of its rectangle, confirm that the step of band auricle rectangular beam sectional dimension is following:
1.1, calculate not with the auricle rectangular beam have predetermined value vertically to moment of inertia I xSquare-section half width a during with polar moment of inertia J rWith cross section half height b r:
1.1.1, set the initial value α of iteration variable α 0=0.93 ~ 0.99:
1.1.2, according to computes intermediate variable β and iteration variable α 1:
24 ( 2 &pi; ) 5 &alpha; 0 &beta; 3 - 4 &beta; 2 + J I x = 0 &alpha; 1 = tanh &pi; 2 &beta; + 1 243 tanh 3 &pi; 2 &beta; . . . [ 1 ]
Attention: when finding the solution β, relate to finding the solution of simple cubic equation, the flat-sawn method is adopted in suggestion, and β ∈ (0,1);
1.1.3, calculate iterative value α 1With initial value α 0Scale error e 1:
e 1 = &alpha; 1 - &alpha; 0 &alpha; 0 . . . [ 2 ]
If | e 1|<0.000001, square-section half width a then rWith cross section half height b rFor:
a r = 3 I x &beta; 3 4 4 b r = 3 I x 4 &beta; 4 . . . [ 3 ]
Step 1.1 finishes; Otherwise, carry out step 1.1.4;
1.1.4, make α 01, the method for repeating step 1.1.2~step 1.1.3, through iterative computation, up to | e 1|<0.000001, and obtain square-section half width a at this moment rWith cross section half height b r, step 1.1 finishes;
1.2, predetermined polar moment of inertia J is adjusted: go out the adjusted value J of polar moment of inertia J according to computes 1,
Figure FDA00001628974600014
Note answering guaranteed conditions: t<a r;
1.3, calculate half width a and half height b: change the J in the formula [1] into J 1, calculate half width a and half height b according to the described method of step 1.1 then;
1.4, calculate band auricle rectangular beam cross section half-breadth l:
1.4.1, according to computes intermediate variable δ:
&delta; = 3 I y - 4 a 3 b 2 t . . . [ 5 ]
1.4.2, calculate band auricle rectangular beam cross section half-breadth l:
l = &delta; + a 3 3 . . . [ 6 ]
So far, obtain half width a, half height b and cross section half-breadth l, and then obtain width 2a, height 2b and the cross-sectional width 2l of corresponding rectangular beam with the auricle rectangular beam.
CN201210146222.5A 2012-05-11 2012-05-11 Method for determining dimension of rectangular beam section of low-speed flutter main beam model of airplane Active CN102682168B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210146222.5A CN102682168B (en) 2012-05-11 2012-05-11 Method for determining dimension of rectangular beam section of low-speed flutter main beam model of airplane

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210146222.5A CN102682168B (en) 2012-05-11 2012-05-11 Method for determining dimension of rectangular beam section of low-speed flutter main beam model of airplane

Publications (2)

Publication Number Publication Date
CN102682168A true CN102682168A (en) 2012-09-19
CN102682168B CN102682168B (en) 2014-03-19

Family

ID=46814088

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210146222.5A Active CN102682168B (en) 2012-05-11 2012-05-11 Method for determining dimension of rectangular beam section of low-speed flutter main beam model of airplane

Country Status (1)

Country Link
CN (1) CN102682168B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102968525A (en) * 2012-11-12 2013-03-13 中国航空工业集团公司西安飞机设计研究所 Method for determining overweight ratio of plane flutter model

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102411652A (en) * 2011-08-15 2012-04-11 中国航空工业集团公司西安飞机设计研究所 Method for determining sectional dimension of rectangular beam with lugs
CN102419163A (en) * 2011-08-15 2012-04-18 中国航空工业集团公司西安飞机设计研究所 Method for determining sectional dimension of thin-walled rectangular hollow beam with lugs

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102411652A (en) * 2011-08-15 2012-04-11 中国航空工业集团公司西安飞机设计研究所 Method for determining sectional dimension of rectangular beam with lugs
CN102419163A (en) * 2011-08-15 2012-04-18 中国航空工业集团公司西安飞机设计研究所 Method for determining sectional dimension of thin-walled rectangular hollow beam with lugs

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
胡志勇,等。: "一种新的颤振计算方法", 《振动工程学报》 *
谢怀强: "xx飞机垂尾低速颤振模型的设计、制造与试验验证", 《第八届全国空气弹性学术交流会论文》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102968525A (en) * 2012-11-12 2013-03-13 中国航空工业集团公司西安飞机设计研究所 Method for determining overweight ratio of plane flutter model

Also Published As

Publication number Publication date
CN102682168B (en) 2014-03-19

Similar Documents

Publication Publication Date Title
CN102567596B (en) Parametric modeling method of face gear teeth
CN102968524B (en) A kind of modeling method of two-dimentional variable curvature section bar part process model
CN102411652B (en) Method for determining sectional dimension of rectangular beam with lugs
CN101644915B (en) Computer-control polishing method based on removal function prediction model
CN104697462B (en) A kind of aerial blade surface-type feature parameter extracting method based on axis
CN103902785A (en) Structure finite element model correcting method based on multi-element uncertainty
CN102419163A (en) Method for determining sectional dimension of thin-walled rectangular hollow beam with lugs
Ito et al. Unstructured Mesh Generation for Viscous Flow Computations.
CN103324778A (en) Ground load determination method of multi-fulcrum airplane
CN105109705A (en) Method for calculating rigidity of aircraft airfoil surface structure
CN109726437B (en) Cabin door pneumatic load equivalent node force processing method
CN102867097A (en) Method for designing photo-cure quickly formed wind tunnel model in consideration of influence of static elastic deformation
CN110210130B (en) Shape optimization method for I-beam two-dimensional model
CN104899345A (en) Method for determining laser peening forming process parameter of complex curved-surface-shaped workpiece
CN103399996A (en) Multi-target topological optimization design method for flexible mechanism for fast servo tool rest
CN105697625A (en) Design method of few-leaf parabolic isostress steel plate spring provided with ends of different structures
CN102661730B (en) Method for determining hollow girder cross-section size of airplane high-speed vibration girder model
CN103577635A (en) Fitting method of blade surface data
CN107169171A (en) A kind of EMUs head dummy design method
CN111551343B (en) Design method for full-speed domain aerodynamic characteristic wind tunnel test of rocket with grid rudder
CN102682168B (en) Method for determining dimension of rectangular beam section of low-speed flutter main beam model of airplane
CN105718619A (en) Method for determining fuel quality characteristics of aircraft based on finite element method
CN104102774A (en) Method for designing linearity of mounting surface for rolling linear guide of precision machine tool
CN103902786A (en) Root extending segment appearance parameterization design method of turbine air-cooling moving blade
CN102661729B (en) Method for confirming sectional dimension of I-shaped hollow beam of high-speed fluttering model of airplane

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