CN102479269B - Method for converting fluid load into solid load - Google Patents
Method for converting fluid load into solid load Download PDFInfo
- Publication number
- CN102479269B CN102479269B CN 201010565341 CN201010565341A CN102479269B CN 102479269 B CN102479269 B CN 102479269B CN 201010565341 CN201010565341 CN 201010565341 CN 201010565341 A CN201010565341 A CN 201010565341A CN 102479269 B CN102479269 B CN 102479269B
- Authority
- CN
- China
- Prior art keywords
- unit
- solid structure
- solid
- structure unit
- node
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
The invention relates to a method for converting a fluid load into a solid load. The method sequentially comprises the following steps of: 1, establishing a fluid field computing model; 2, establishing a solid computing model; 3, collecting the area of a stress surface; 4, establishing local coordinate systems; 5, establishing a transformation matrix among the coordinate systems; 6, obtaining solid structure units corresponding to pneumatic grid nodes; 7, obtaining a concentrated force; 8, establishing isoparametric units; 9, establishing a mapping relationship between the solid structure units and the corresponding isoparametric units; 10, establishing an area infinitesimal; 11, obtaining all node forces in the local coordinate systems according to a principle of minimum potential energy; 12, converting the node forces into an integral coordinate system; and 13, obtaining the node force of each solid structure unit. In the invention, according to the principle of minimum potential energy, the conversion from the fluid load to the solid load can be realized well, so that higher precision can be achieved and the method has a higher actual application value.
Description
Technical field
The present invention relates to the conversion method of a kind of stream-immobilized lotus, particularly relate to the conversion method of the open stream of a kind of high precision-immobilized lotus.
Background technology
When the structural design of aircraft, the solid coupling of stream is the problem that must consider.Difference due to fluid domain and solid domain governing equation characteristics, more to use the Euler coordinate system in fluid mechanics, the motion state of main research spatial point, find the solution with the Fluid Computation mechanics method, and custom adopts the Lagrangian coordinate system in solid mechanics, main research material motion of point adopts the computation structure dynamic method to find the solution, and its method of mainly using is finite element method.Its common characteristics are all to carry out that grid division is discrete finds the solution to finding the solution the territory.
Along with being showing improvement or progress day by day of spationautics, aerodynamic configuration and the version of aircraft become increasingly complex, and the form of aircraft bearing load also becomes increasingly complex.How solving the load of fluid aerodynamic grid to the conversion of solid structure unit load, is a very important problem.In the practical engineering calculation for aircraft, the numerical evaluation in flow field requires computing grid to be stretched over the enough general goals of space relative model characteristic length from aircraft surface, and solid structure unit calculation requirement computing grid extends to inside from aircraft surface.On the other hand, Flow Field Numerical Calculation generally wants thin at body surface slope variation general goal grid, and solid structure unit requirement body surface grid will be divided more as far as possible, in order to can obtain easily stiffness matrix.The traditional method of processing Similar Problems at industrial circles such as Aero-Space comprises that pneumatic structure conode method, choose method, spline surface method of interpolation at 3.These three kinds of popular disposal routes have certain defect, and precision is not enough.Design a kind of classic method defective that overcomes so need badly, realize the conversion method of the stream that high precision transforms between different grid systems-immobilized lotus.
Summary of the invention
The technical problem to be solved in the present invention is to provide the conversion method of a kind of high-precision open stream-immobilized lotus.
For solving the problems of the technologies described above, the conversion method of a kind of stream of the present invention-immobilized lotus comprises the following steps successively:
Exact shape and the physical size of step 1, collection solid structure are set up JASMINE; Gather the physical parameter in the residing flow field of solid structure, the stream field computation model is divided aerodynamic grid, obtains the coordinate of each aerodynamic grid node and the aerodynamic force of each aerodynamic grid node;
Exact shape and the physical size of step 2, the solid structure that gathers according to step 1 are set up solid computation model and global coordinate system; The solid computation model is divided the node of solid structure unit and solid structure unit, determined the solid structure unit under each node, the coordinate of node under global coordinate system of each solid structure unit; Described solid structure unit comprises solid element and plate shell unit;
Step 3, when the solid structure unit is solid element, differentiate the stress surface of solid structure unit, gather the area of stress surface; When the solid structure unit is the plate shell unit, gather the area of stress surface;
Step 4, set up a local coordinate system on the middle face of each solid structure unit;
Step 5, obtain the direction cosine of each local coordinate system, set up the transition matrix between global coordinate system and each local coordinate system;
Step 6, obtain the aerodynamic grid node to the distance of each solid structure unit; The aerodynamic grid node is the mean value that aerodynamic grid is pressed heart node distance of affiliated node to the solid structure unit to the distance of each solid structure unit; Obtain from the nearest solid structure unit of aerodynamic grid nodal distance; Repeat this step, find each aerodynamic grid node corresponding from its nearest solid structure unit;
Step 7, the aerodynamic force of each aerodynamic grid node be multiply by the area of the stress surface of the corresponding solid structure of this aerodynamic grid node unit, obtain concentrated force;
Step 8, set up each solid structure unit etc. the ginseng unit, and the node of the ginseng such as establishment unit, establish the solid structure unit node corresponding with this solid structure unit etc. the corresponding relation between the node of ginseng unit, establish each and wait ginseng unit displacement shape function accordingly under natural system of coordinates;
Step 9, set up each solid structure unit and its correspondence etc. the mapping relations of ginseng unit, obtain the Jacobian of mapping relations;
Step 10, set up under natural system of coordinates etc. the area element of ginseng unit, Jacobian be multiply by the area element that waits the ginseng unit under natural system of coordinates, obtain the area element of the solid structure unit under global coordinate;
Step 11, utilize minimum potential energy principal, concentrated force is converted into the nodal force of solid structure unit under local coordinate system;
Comprise: make aerodynamic force on virtual displacement arbitrarily to solid structure work and nodal force on virtual displacement, the solid structure work being equated arbitrarily, make virtual displacement and actual displacement have identical displacement shape function;
When the solid structure unit is the plate shell unit, make concentrated force act on the center that grade corresponding to plate shell unit joined the unit, with etc. the transposed matrix of displacement shape function of ginseng unit multiply by the concentrated force array, obtain waiting the nodal force of ginseng unit, according to the corresponding relation of the node of the node that waits the ginseng unit and solid structure unit, obtain the nodal force of the plate shell unit corresponding with waiting the ginseng unit under local coordinate system;
When the solid structure unit is solid element, make concentrated force act on the center that grade corresponding to solid element joined the unit, with etc. the transposed matrix of displacement shape function of ginseng unit multiply by the area element that the concentrated force array multiply by the solid structure unit under global coordinate again, and the surface area of solid element stress surface is quadratured, obtain joining with grade under local coordinate system the nodal force of solid element corresponding to unit;
Step 12, the local coordinate system of determining according to step 5 and the transition matrix of global coordinate system, the nodal force that step 11 is obtained is transformed under global coordinate system, obtains all nodal forces under global coordinate system;
Step 13, will be corresponding to the nodal force stack of same solid structure unit, obtain the nodal force of each solid structure unit.
Step 4 is set up a local coordinate system on the middle face of each solid structure unit, be to utilize the rule of setting up of finite element unit local coordinate to set up.
The present invention according to minimum potential energy principal, can realize the load transfer between fluid aerodynamic grid and solid structure grid by introducing the viewpoint of energy well, reaches higher precision, has higher actual application value.
The present invention has open characteristics, has very high engineering using value.
Embodiment
The present invention comprises the following steps successively:
Exact shape and the physical size of step 1, collection solid structure are set up JASMINE; Gather the physical parameter in the residing flow field of solid structure, the stream field computation model is divided aerodynamic grid, obtains the coordinate of each aerodynamic grid node and the aerodynamic force of each aerodynamic grid node;
Exact shape and the physical size of step 2, the solid structure that gathers according to step 1 are set up solid computation model and global coordinate system; The solid computation model is divided the node of solid structure unit and solid structure unit, determined the solid structure unit under each node, the coordinate of node under global coordinate system of each solid structure unit; The solid structure unit comprises solid element and plate shell unit;
Step 3, when the solid structure unit is solid element, differentiate the stress surface of solid structure unit, gather the area of stress surface; When the solid structure unit is the plate shell unit, gather the area of stress surface;
Step 4, according to the rule of setting up of finite element unit local coordinate, set up a local coordinate system on the middle face of each solid structure unit;
Step 5, obtain the direction cosine of each local coordinate system, set up the transition matrix between global coordinate system and each local coordinate system;
Step 6, obtain the aerodynamic grid node to the distance of each solid structure unit; The aerodynamic grid node is the mean value that aerodynamic grid is pressed heart node distance of affiliated node to the solid structure unit to the distance of each solid structure unit; Obtain from the nearest solid structure unit of aerodynamic grid nodal distance; Repeat this step, find each aerodynamic grid node corresponding from its nearest solid structure unit;
Step 7, the aerodynamic force of each aerodynamic grid node be multiply by the area of the stress surface of the corresponding solid structure of this aerodynamic grid node unit, obtain concentrated force;
Step 8, set up each solid structure unit etc. the ginseng unit, and the node of the ginseng such as establishment unit, establish the solid structure unit node corresponding with this solid structure unit etc. the corresponding relation between the node of ginseng unit, establish each and wait ginseng unit displacement shape function accordingly under natural system of coordinates;
Step 9, set up each solid structure unit and its correspondence etc. the mapping relations of ginseng unit, obtain the Jacobian of mapping relations;
Step 10, set up under natural system of coordinates etc. the area element of ginseng unit, Jacobian be multiply by the area element that waits the ginseng unit under natural system of coordinates, obtain the area element of the solid structure unit under global coordinate;
Step 11, make aerodynamic force on virtual displacement arbitrarily to solid structure work and nodal force on virtual displacement, the solid structure work being equated arbitrarily, make virtual displacement and actual displacement have identical displacement shape function;
When the solid structure unit is the plate shell unit, make concentrated force act on the center that grade corresponding to plate shell unit joined the unit, with etc. the transposed matrix of displacement shape function of ginseng unit multiply by the concentrated force array, obtain waiting the nodal force of ginseng unit, according to the corresponding relation of the node of the node that waits the ginseng unit and solid structure unit, obtain the nodal force of the plate shell unit corresponding with waiting the ginseng unit under local coordinate system;
When the solid structure unit is solid element, make concentrated force act on the center that grade corresponding to solid element joined the unit, with etc. the transposed matrix of displacement shape function of ginseng unit multiply by the area element that the concentrated force array multiply by the solid structure unit under global coordinate again, and the surface area of solid element stress surface is quadratured, obtain the nodal force of the solid element corresponding with waiting the ginseng unit under local coordinate system;
Step 12, the local coordinate system of determining according to step 5 and the transition matrix of global coordinate system, the nodal force that step 12 is obtained is transformed under global coordinate system, obtains all nodal forces under global coordinate system;
Step 13, will be corresponding to the nodal force stack of same solid structure unit, obtain the nodal force of each solid structure unit.
The present invention is applicable to flow analytical calculation and structural static strength, fatigue resistance, vibration and other analytical calculations of solid coupling.
Claims (2)
1. the conversion method of stream-immobilized lotus comprises the following steps successively:
Exact shape and the physical size of step 1, collection solid structure are set up JASMINE; Gather the physical parameter in the residing flow field of solid structure, the stream field computation model is divided aerodynamic grid, obtains the coordinate of each aerodynamic grid node and the aerodynamic force of each aerodynamic grid node;
Exact shape and the physical size of step 2, the solid structure that gathers according to step 1 are set up solid computation model and global coordinate system; The solid computation model is divided the node of solid structure unit and solid structure unit, determined the solid structure unit under each node, the coordinate of node under global coordinate system of each solid structure unit; Described solid structure unit comprises solid element and plate shell unit;
Step 3, when the solid structure unit is solid element, differentiate the stress surface of solid structure unit, gather the area of stress surface; When the solid structure unit is the plate shell unit, gather the area of stress surface;
Step 4, set up a local coordinate system on the middle face of each solid structure unit;
Step 5, obtain the direction cosine of each local coordinate system, set up the transition matrix between global coordinate system and each local coordinate system;
Step 6, obtain the aerodynamic grid node to the distance of each solid structure unit; The aerodynamic grid node is the mean value that aerodynamic grid is pressed heart node distance of affiliated node to the solid structure unit to the distance of each solid structure unit; Obtain from the nearest solid structure unit of aerodynamic grid nodal distance; Repeat this step, find each aerodynamic grid node corresponding from its nearest solid structure unit;
Step 7, the aerodynamic force of each aerodynamic grid node be multiply by the area of the stress surface of the corresponding solid structure of this aerodynamic grid node unit, obtain concentrated force;
Step 8, set up each solid structure unit etc. the ginseng unit, and the node of the ginseng such as establishment unit, establish the solid structure unit node corresponding with this solid structure unit etc. the corresponding relation between the node of ginseng unit, establish each and wait ginseng unit displacement shape function accordingly under natural system of coordinates;
Step 9, set up each solid structure unit and its correspondence etc. the mapping relations of ginseng unit, obtain the Jacobian of mapping relations;
Step 10, set up under natural system of coordinates etc. the area element of ginseng unit, Jacobian be multiply by the area element that waits the ginseng unit under natural system of coordinates, obtain the area element of the solid structure unit under global coordinate;
Step 11, utilize minimum potential energy principal, concentrated force is converted into the nodal force of solid structure unit under local coordinate system;
Comprise: make aerodynamic force on virtual displacement arbitrarily to solid structure work and nodal force on virtual displacement, the solid structure work being equated arbitrarily, make virtual displacement and actual displacement have identical displacement shape function;
When the solid structure unit is the plate shell unit, make concentrated force act on the center that grade corresponding to plate shell unit joined the unit, with etc. the transposed matrix of displacement shape function of ginseng unit multiply by the concentrated force array, obtain waiting the nodal force of ginseng unit, according to the corresponding relation of the node of the node that waits the ginseng unit and solid structure unit, obtain the nodal force of the plate shell unit corresponding with waiting the ginseng unit under local coordinate system;
When the solid structure unit is solid element, make concentrated force act on the center that grade corresponding to solid element joined the unit, with etc. the transposed matrix of displacement shape function of ginseng unit multiply by the area element that the concentrated force array multiply by the solid structure unit under global coordinate again, and the surface area of solid element stress surface is quadratured, obtain the nodal force of the solid element corresponding with waiting the ginseng unit under local coordinate system;
Step 12, the local coordinate system of determining according to step 5 and the transition matrix of global coordinate system, the nodal force that step 11 is obtained is transformed under global coordinate system, obtains all nodal forces under global coordinate system;
Step 13, will be corresponding to the nodal force stack of same solid structure unit, obtain the nodal force of each solid structure unit.
2. the conversion method of a kind of stream according to claim 1-immobilized lotus is characterized in that: described step 4 is set up a local coordinate system on the middle face of each solid structure unit, is to utilize the rule of setting up of finite element unit local coordinate to set up.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010565341 CN102479269B (en) | 2010-11-29 | 2010-11-29 | Method for converting fluid load into solid load |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010565341 CN102479269B (en) | 2010-11-29 | 2010-11-29 | Method for converting fluid load into solid load |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102479269A CN102479269A (en) | 2012-05-30 |
| CN102479269B true CN102479269B (en) | 2013-05-22 |
Family
ID=46091913
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201010565341 Active CN102479269B (en) | 2010-11-29 | 2010-11-29 | Method for converting fluid load into solid load |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102479269B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106709179B (en) * | 2016-12-21 | 2020-04-07 | 中国燃气涡轮研究院 | Blade aerodynamic force distribution data processing method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6611736B1 (en) * | 2000-07-01 | 2003-08-26 | Aemp Corporation | Equal order method for fluid flow simulation |
| CN101706832A (en) * | 2009-11-25 | 2010-05-12 | 哈尔滨工业大学 | Optimization design method of fibre enhanced composite material marine propeller blade |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007052784A1 (en) * | 2005-11-07 | 2007-05-10 | Keio University | Numerical structure analyzer by load transfer path method |
-
2010
- 2010-11-29 CN CN 201010565341 patent/CN102479269B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6611736B1 (en) * | 2000-07-01 | 2003-08-26 | Aemp Corporation | Equal order method for fluid flow simulation |
| CN101706832A (en) * | 2009-11-25 | 2010-05-12 | 哈尔滨工业大学 | Optimization design method of fibre enhanced composite material marine propeller blade |
Non-Patent Citations (4)
| Title |
|---|
| 一种新型航天器有效载荷承力结构的优化分析;王骁峰等;《导弹与航天运载技术》;20071031(第5期);全文 * |
| 王骁峰等.一种新型航天器有效载荷承力结构的优化分析.《导弹与航天运载技术》.2007,(第5期), |
| 结构流-固冲击屈曲研究进展;顾王明等;《力学进展》;19960225;第26卷(第1期);全文 * |
| 顾王明等.结构流-固冲击屈曲研究进展.《力学进展》.1996,第26卷(第1期), |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102479269A (en) | 2012-05-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106021644B (en) | The method for determining mixed dimensional modelling interface constraint equation coefficient | |
| CN103838913B (en) | The Finite Element of the curved bridge of curved box girder | |
| CN104036095A (en) | Regional-decomposition based high-precision coupling fast-calculation method for complex-shape flow field | |
| CN106446432B (en) | A kind of solution the optimal of material large deformation transports non-mesh method | |
| CN106934185A (en) | A kind of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid | |
| CN109446471B (en) | Fluid-solid coupling interface data transmission method considering load uncertainty | |
| CN102446241A (en) | Wing surface structural rigidity simulation method | |
| CN110348110A (en) | A kind of bolt connected piece rigidity adaptive recognition methods | |
| CN104281730A (en) | Great-rotating-deformation plate shell structure dynamic response finite element analysis method | |
| CN105205293A (en) | Method and system for acquiring aerodynamic load of aircraft component | |
| CN110020460A (en) | It is bolted flanged (FLGD) cylindrical shell structure frequency response function Uncertainty Analysis Method | |
| CN105912007A (en) | Differential geometry nonlinear control method of spatial mechanical arm anti-interference attitude stabilization | |
| CN114943167B (en) | Method, system, medium and equipment for calculating wall surface distance of structural grid | |
| CN102479269B (en) | Method for converting fluid load into solid load | |
| CN203203742U (en) | Box type six-component balance | |
| Pan et al. | Investigation on the choice of boundary conditions and shape functions for flexible multi-body system | |
| CN109459206A (en) | Ground experiment unsteady aerodynamic force loading method | |
| CN104807606A (en) | MATLAB-STM32 hybrid power test system and test method thereof | |
| CN105160092A (en) | Thermal environment interpolation method suitable for thermal protection system transient temperature field calculation | |
| CN104239588B (en) | A kind of structural-functional analysis method based on ultra-unit construction method | |
| CN202255706U (en) | Three-dimensional force sensor | |
| CN101813555A (en) | Soft abrasive flow field testing method based on level set | |
| CN102436514B (en) | Method for updating fluid-solid interaction grid | |
| CN111008417B (en) | Wind-induced vibration analysis method for continuous structure with large slenderness ratio | |
| CN104484498B (en) | It is determined that effectively pushing the solar sail transient dynamic analysis method that acceleration is lost |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |