CN102479269A - Method for converting fluid load into solid load - Google Patents
Method for converting fluid load into solid load Download PDFInfo
- Publication number
- CN102479269A CN102479269A CN2010105653415A CN201010565341A CN102479269A CN 102479269 A CN102479269 A CN 102479269A CN 2010105653415 A CN2010105653415 A CN 2010105653415A CN 201010565341 A CN201010565341 A CN 201010565341A CN 102479269 A CN102479269 A CN 102479269A
- 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.)
- Granted
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.Because fluid territory and solid territory governing equation characteristics is different; In fluid mechanics, more be to use the Euler coordinate system, mainly study the motion state of spatial point, find the solution with the Fluid Computation mechanics method; And custom adopts the Lagrangian coordinate system in solid mechanics; The motion of main research material point adopts the computation structure dynamic method to find the solution, and its method of mainly using is a finite element method.Its common characteristic all is to divide that grid is discrete finds the solution to finding the solution the territory.
Along with being showing improvement or progress day by day of spationautics, the aerodynamic configuration and the version of aircraft become increasingly complex, and the form of aircraft bearing load also becomes increasingly complex.How solving the conversion of the load of the pneumatic grid of fluid to the solid structure unit load, is a very important problem.In being directed against the practical engineering calculation of aircraft, numerical Calculation of 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 the solid structure unit requires the body surface grid will divide more as far as possible, so that can obtain stiffness matrix easily.The traditional method of handling similar problem 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 all have certain defective, and precision is not enough.Design a kind of classic method defective that overcomes so need badly, realize the conversion method of stream-immobilized lotus that high precision transforms between the different grid systems.
Summary of the invention
The technical matters that the present invention will solve provides 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 may further comprise the steps successively:
The 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 pneumatic grid, obtains the coordinate of each pneumatic grid node and the aerodynamic force of each pneumatic grid node;
The exact shape and the physical size of step 2, the solid structure gathered according to step 1 are set up solid computation model and global coordinate system; To the node of solid computation model division solid structure unit and solid structure unit, confirm the solid structure unit that each node is affiliated, the coordinate of the node of each solid structure unit under global coordinate system; Said solid structure unit comprises solid element and plate shell unit;
Step 3, when the solid structure unit is a 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, on the middle face of each solid structure unit, set up a local coordinate system;
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 pneumatic grid node to the distance of each solid structure unit; Pneumatic grid node is the mean value that pneumatic 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 pneumatic grid node; Repeat this step, find each pneumatic grid node pairing from its nearest solid structure unit;
Step 7, the aerodynamic force of each pneumatic grid node multiply by the area of the stress surface of the pairing solid structure of this pneumatic grid node unit, obtain concentrated force;
Step 8, set up each solid structure unit etc. the ginseng unit; And the node of ginseng such as establishment unit, establish the node of solid structure unit corresponding with this solid structure unit etc. the corresponding relation between the node of ginseng unit, establish each etc. join unit corresponding displacement shape function under natural system of coordinates;
Step 9, set up each solid structure unit and its corresponding etc. the mapping relations of ginseng unit, obtain the Jacobi determinant of mapping relations;
Step 10, set up under the natural system of coordinates etc. the area element of ginseng unit, the Jacobi determinant multiply by the area element of wait joining the unit natural system of coordinates under, obtain the area element of the solid structure unit under the global coordinate;
Step 11, utilize minimum potential energy principal, concentrated force is converted into the nodal force of solid structure unit under the local coordinate system;
Step 12, the local coordinate system of confirming according to step 5 and the transition matrix of global coordinate system, the nodal force that step 11 is obtained is transformed under the global coordinate system, obtains all nodal forces under the global coordinate system;
Step 13, will obtain the nodal force of each solid structure unit corresponding to the stack of the nodal force of same solid structure unit.
Step 11 comprises: make aerodynamic force on virtual displacement arbitrarily to solid structure work and nodal force on the virtual displacement 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 corresponding center of waiting the ginseng unit of plate shell 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 node that waits the ginseng unit and solid structure unit, obtain local coordinate system under with etc. join the nodal force of the corresponding plate shell unit in unit;
When the solid structure unit is a solid element; Make concentrated force act on the corresponding center of waiting the ginseng unit of solid element; 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 the global coordinate again; And the surface area of solid element stress surface quadratured, obtain the nodal force of solid element corresponding under the local coordinate system with waiting ginseng 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 conversion between pneumatic grid of fluid and the solid structure grid through introducing viewpoint of energy well, reaches higher precision, has higher actual application value.
The present invention has open characteristics, has very high practical applications and is worth.
Embodiment
The present invention may further comprise the steps successively:
The 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 pneumatic grid, obtains the coordinate of each pneumatic grid node and the aerodynamic force of each pneumatic grid node;
The exact shape and the physical size of step 2, the solid structure gathered according to step 1 are set up solid computation model and global coordinate system; To the node of solid computation model division solid structure unit and solid structure unit, confirm the solid structure unit that each node is affiliated, the coordinate of the node of each solid structure unit under global coordinate system; The solid structure unit comprises solid element and plate shell unit;
Step 3, when the solid structure unit is a 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, on the middle face of each solid structure unit, set up a local coordinate system;
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 pneumatic grid node to the distance of each solid structure unit; Pneumatic grid node is the mean value that pneumatic 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 pneumatic grid node; Repeat this step, find each pneumatic grid node pairing from its nearest solid structure unit;
Step 7, the aerodynamic force of each pneumatic grid node multiply by the area of the stress surface of the pairing solid structure of this pneumatic grid node unit, obtain concentrated force;
Step 8, set up each solid structure unit etc. the ginseng unit; And the node of ginseng such as establishment unit, establish the node of solid structure unit corresponding with this solid structure unit etc. the corresponding relation between the node of ginseng unit, establish each etc. join unit corresponding displacement shape function under natural system of coordinates;
Step 9, set up each solid structure unit and its corresponding etc. the mapping relations of ginseng unit, obtain the Jacobi determinant of mapping relations;
Step 10, set up under the natural system of coordinates etc. the area element of ginseng unit, the Jacobi determinant multiply by the area element of wait joining the unit natural system of coordinates under, obtain the area element of the solid structure unit under the global coordinate;
Step 11, make aerodynamic force on virtual displacement arbitrarily to solid structure work and nodal force on the virtual displacement 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 corresponding center of waiting the ginseng unit of plate shell 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 node that waits the ginseng unit and solid structure unit, obtain local coordinate system under with etc. join the nodal force of the corresponding plate shell unit in unit;
When the solid structure unit is a solid element; Make concentrated force act on the corresponding center of waiting the ginseng unit of solid element; 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 the global coordinate again; And the surface area of solid element stress surface quadratured, obtain the nodal force of the solid element corresponding under the local coordinate system with waiting ginseng unit;
Step 12, the local coordinate system of confirming according to step 5 and the transition matrix of global coordinate system, the nodal force that step 12 is obtained is transformed under the global coordinate system, obtains all nodal forces under the global coordinate system;
Step 13, will obtain the nodal force of each solid structure unit corresponding to the stack of the nodal force of same solid structure unit.
The present invention is applicable to analytical calculation and structure static strength, fatigue resistance, vibration and other analytical calculations of the solid coupling of stream.
Claims (3)
1. the conversion method of stream-immobilized lotus may further comprise the steps successively:
The 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 pneumatic grid, obtains the coordinate of each pneumatic grid node and the aerodynamic force of each pneumatic grid node;
The exact shape and the physical size of step 2, the solid structure gathered according to step 1 are set up solid computation model and global coordinate system; To the node of solid computation model division solid structure unit and solid structure unit, confirm the solid structure unit that each node is affiliated, the coordinate of the node of each solid structure unit under global coordinate system; Said solid structure unit comprises solid element and plate shell unit;
Step 3, when the solid structure unit is a 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, on the middle face of each solid structure unit, set up a local coordinate system;
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 pneumatic grid node to the distance of each solid structure unit; Pneumatic grid node is the mean value that pneumatic 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 pneumatic grid node; Repeat this step, find each pneumatic grid node pairing from its nearest solid structure unit;
Step 7, the aerodynamic force of each pneumatic grid node multiply by the area of the stress surface of the pairing solid structure of this pneumatic grid node unit, obtain concentrated force;
Step 8, set up each solid structure unit etc. the ginseng unit; And the node of ginseng such as establishment unit, establish the node of solid structure unit corresponding with this solid structure unit etc. the corresponding relation between the node of ginseng unit, establish each etc. join unit corresponding displacement shape function under natural system of coordinates;
Step 9, set up each solid structure unit and its corresponding etc. the mapping relations of ginseng unit, obtain the Jacobi determinant of mapping relations;
Step 10, set up under the natural system of coordinates etc. the area element of ginseng unit, the Jacobi determinant multiply by the area element of wait joining the unit natural system of coordinates under, obtain the area element of the solid structure unit under the global coordinate;
Step 11, utilize minimum potential energy principal, concentrated force is converted into the nodal force of solid structure unit under the local coordinate system;
Step 12, the local coordinate system of confirming according to step 5 and the transition matrix of global coordinate system, the nodal force that step 11 is obtained is transformed under the global coordinate system, obtains all nodal forces under the global coordinate system;
Step 13, will obtain the nodal force of each solid structure unit corresponding to the stack of the nodal force of same solid structure unit.
2. the conversion method of a kind of stream according to claim 1-immobilized lotus; It is characterized in that: said step 11 comprises: make aerodynamic force on virtual displacement arbitrarily to solid structure work and nodal force on the virtual displacement 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 corresponding center of waiting the ginseng unit of plate shell 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 node that waits the ginseng unit and solid structure unit, obtain local coordinate system under with etc. join the nodal force of the corresponding plate shell unit in unit;
When the solid structure unit is a solid element; Make concentrated force act on the corresponding center of waiting the ginseng unit of solid element; 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 the global coordinate again; And the surface area of solid element stress surface quadratured, obtain the nodal force of the solid element corresponding under the local coordinate system with waiting ginseng unit;
3. the conversion method of a kind of stream according to claim 1-immobilized lotus is characterized in that: said 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 true CN102479269A (en) | 2012-05-30 |
CN102479269B 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) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106709179A (en) * | 2016-12-21 | 2017-05-24 | 中国燃气涡轮研究院 | Blade aerodynamic distribution data processing method suitable for finite element free mesh |
Citations (3)
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 |
US20090125282A1 (en) * | 2005-11-07 | 2009-05-14 | Keio University | Numerical structural analysis system based on the load-transfer-path method |
CN101706832A (en) * | 2009-11-25 | 2010-05-12 | 哈尔滨工业大学 | Optimization design method of fibre enhanced composite material marine propeller blade |
-
2010
- 2010-11-29 CN CN 201010565341 patent/CN102479269B/en active Active
Patent Citations (3)
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 |
US20090125282A1 (en) * | 2005-11-07 | 2009-05-14 | Keio University | Numerical structural analysis system based on the load-transfer-path method |
CN101706832A (en) * | 2009-11-25 | 2010-05-12 | 哈尔滨工业大学 | Optimization design method of fibre enhanced composite material marine propeller blade |
Non-Patent Citations (2)
Title |
---|
王骁峰等: "一种新型航天器有效载荷承力结构的优化分析", 《导弹与航天运载技术》 * |
顾王明等: "结构流-固冲击屈曲研究进展", 《力学进展》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106709179A (en) * | 2016-12-21 | 2017-05-24 | 中国燃气涡轮研究院 | Blade aerodynamic distribution data processing method suitable for finite element free mesh |
CN106709179B (en) * | 2016-12-21 | 2020-04-07 | 中国燃气涡轮研究院 | Blade aerodynamic force distribution data processing method |
Also Published As
Publication number | Publication date |
---|---|
CN102479269B (en) | 2013-05-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109902404B (en) | Unified recursion calculation method for structural time-course response integral of different damping forms | |
CN102663229B (en) | Form-finding method of cable-strut tensile structure | |
CN106096119B (en) | Casual acceleration based on large-mass method motivates flowering structure method of topological optimization design | |
CN106021644B (en) | The method for determining mixed dimensional modelling interface constraint equation coefficient | |
CN102446241A (en) | Wing surface structural rigidity simulation method | |
CN102968542A (en) | Method for performing structure analysis on steel pylon by using ANSYS software | |
CN103838913A (en) | Finite element method for curve box girder curved bridge | |
CN104281730A (en) | Great-rotating-deformation plate shell structure dynamic response finite element analysis method | |
CN109459206A (en) | Ground experiment unsteady aerodynamic force loading method | |
CN104807606A (en) | MATLAB-STM32 hybrid power test system and test method thereof | |
CN102479269B (en) | Method for converting fluid load into solid load | |
CN104239588B (en) | A kind of structural-functional analysis method based on ultra-unit construction method | |
CN105160092A (en) | Thermal environment interpolation method suitable for thermal protection system transient temperature field calculation | |
CN102004822A (en) | Modality-based frequency analysis method for space six-freedom-degree parallel motion system | |
CN109783891B (en) | Method for calculating vibration characteristic of spatial pipeline with passive dynamic vibration absorber | |
CN102436514B (en) | Method for updating fluid-solid interaction grid | |
Raju et al. | Nonlinear structural analysis using integrated force method | |
CN104484498B (en) | It is determined that effectively pushing the solar sail transient dynamic analysis method that acceleration is lost | |
CN111008417A (en) | Wind-induced vibration analysis method for continuous structure with large slenderness ratio | |
CN110008543A (en) | A kind of emulation mode for considering neutral axis of the beam and rotating beam dynamic response being influenced | |
Yu et al. | High-fidelity optimization of flapping airfoils for maximum propulsive efficiency | |
Karpel et al. | Rigid-body issues in FFT-based dynamic loads analysis with aeroservoelastic nonlinearities | |
CN104019094B (en) | A kind of Octahedral Truss Units variable geometry truss robot driving lever | |
CN109977498A (en) | A method of FGM ladder beam dynamic response is calculated based on HOC | |
Zhang et al. | Dynamic mesh method based on diffusion equation and nodal rotation for high-aspect-ratio composite wings |
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 |