CN104699872A - Numerical analysis method for connection between three-dimensional and two-dimensional finite element units - Google Patents
Numerical analysis method for connection between three-dimensional and two-dimensional finite element units Download PDFInfo
- Publication number
- CN104699872A CN104699872A CN201310660640.0A CN201310660640A CN104699872A CN 104699872 A CN104699872 A CN 104699872A CN 201310660640 A CN201310660640 A CN 201310660640A CN 104699872 A CN104699872 A CN 104699872A
- Authority
- CN
- China
- Prior art keywords
- dimensional
- unit
- connection
- analysis
- rigid
- 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.)
- Pending
Links
Abstract
The invention provides a numerical analysis method for the connection between three-dimensional and two-dimensional finite element units, and aims to save time and reduce analysis cost by replacing part of a test result with numerical simulation. The numerical analysis method is characterized by comprising the following steps: (1) constructing a three-dimensional unit and a two-dimensional unit at actual structure positions respectively, endowing the three-dimensional unit with a simulated rotating freedom degree at the connection between the three-dimensional unit and the two-dimensional unit, and establishing a connection relationship between the two-dimensional unit and a new three-dimensional unit by virtue of a rigid body element; (2) creating a card required by structure analysis according to a connection form and a constraint relationship of the rigid body element; (3) creating a unit rigidity matrix according to the card required by analysis, assembling the unit rigidity matrix with other unit rigidity matrixes of the structure to form an overall rigidity matrix, and creating structure deformation and unit stress distribution according to a structure load and a boundary condition.
Description
Technical field
The invention belongs to technology of numerical simulation, relate to the three-dimensional element model in structural finite element analysis and two-dimensional cell model.
Background technology
In aircaft configuration modeling, in order to power transmission and the details stress state of expression structure better, often can use three-dimensional element (volume elements) and two-dimensional cell (plate unit, Liang Yuan etc.) simultaneously.The node of volume elements only has translational degree of freedom, does not provide the rigidity that rotational freedom is corresponding, and palette and beam element have the rigidity (although Slab element may do not have normal direction rotational stiffness) corresponding with rotational freedom.The inharmonious meeting of this unit firm battle array there will be degree of freedom unmatched problem when causing plate, beam element to connect with body unit.Plate and beam element have the rigidity (although Slab element may do not have normal direction rotational stiffness) corresponding with rotational freedom.So when plate unit or beam is first be connected to volume elements time often produce the unmatched problem of degree of freedom.Existing settling mode all adopts three-dimensional element analysis, and this causes analysis time and cost significantly to increase.
Summary of the invention
Goal of the invention: provide a kind of three-dimensional to be connected numerical analysis method with two dimensional finite element unit, replace part test result with numerical simulation, save time and analysis cost.
Technical scheme: a kind of three-dimensional is connected numerical analysis method with two dimensional finite element unit, comprises
(1) build three-dimensional element and two-dimensional cell respectively in the physical location of structure, give the rotational freedom that three-dimensional element is simulated in the junction of three bit locations and two-dimensional cell, set up the annexation of two-dimensional cell and new three-dimensional element by rigid-body system;
(2) card needed for the type of attachment of rigid-body system and the restriction relation generating structure analysis of rigid-body system;
(3) the firm battle array of card generation unit needed for analysis, and be assembled into overall just battle array with the element stiffness matrix of other unit of structure, the load suffered by integrated structure and boundary condition, the distortion of generating structure and the stress distribution of unit.
Beneficial effect: it is a kind of transition of dissimilar unit that a plate or beam element (two-dimensional cell) same volume elements (three-dimensional element) are linked up.The node that body unit connects only has translation freedoms, and they do not provide the rigidity that rotational freedom is corresponding.The present invention utilizes rigid-body system, by interpolation relation constraint unit, the rotation of connected node can be made organically to associate with the translational degree of freedom of adjacent node, avoid whole two-dimensional cell being transferred to the significantly raising that three-dimensional element calculates analysis time and the cost brought, save time and analysis cost.
Accompanying drawing explanation
Fig. 1 is the connection layout of plate of the present invention unit and volume elements.
Fig. 2 is the connection layout of beam of the present invention unit and volume elements.
Fig. 3 is the Aberration nephogram of structure.
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in further detail.
When setting up finite element model, with rigid-body system, three-dimensional element and two-dimensional cell being coupled together, realizing the object translational degree of freedom of three-dimensional element and the rotational freedom of two-dimensional cell are organically coupled together.
(1) three-dimensional element (volume elements) and two-dimensional cell (plate unit and beam unit) is built respectively in the physical location of structure, give the rotational freedom that three-dimensional element is simulated in the junction of three bit locations and two-dimensional cell, set up the annexation of two-dimensional cell and new three-dimensional element by rigid-body system.
For plate unit and the connection of volume elements, need two rigid-body system, as shown in Figure 1.Rigid-body system (a): passive point is set to node 1, retrains 4,5,6DOF.Initiatively point is set to node 2(retrains 1,2,3DOF), node 3(retrains 1,2,3DOF), node 4(retrains 1,2,3DOF).
Rigid-body system (b): passive point is set to node 2, retrains 4,5,6DOF.Initiatively point is set to node 1(retrains 1,2,3DOF), node 3(retrains 1,2,3DOF), node 4(retrains 1,2,3DOF).
For beam unit and the connection of volume elements, rotational freedom just can associate with translational degree of freedom by a rigid-body system, as shown in Figure 2.
Rigid-body system: passive point is set to node 12, retrains 4,5,6DOF.Initiatively point is set to node 11(retrains 1,2,3DOF), node 13(retrains 1,2,3DOF), node 22(retrains 1,2,3DOF).
(2) card needed for the type of attachment of rigid-body system and the restriction relation generating structure analysis of rigid-body system.
The card that the connection of plate unit and volume elements generates is as shown in table 1:
Table 1
Wherein: RBE3 represents the type of rigid-body system, 901 and 902 represent two the rigid-body system sequence numbers generated; REFGRID represents the node number of initiatively point; REFO represents the degree of freedom of passive constraint; WT1 represents weighting coefficient; C1 represents the degree of freedom of initiatively some constraint; G1,1 and G1,2 node numbers representing initiatively point.
The card that the connection of Liang Yuan and volume elements generates is as shown in table 2:
Table 2
| RBE3 | 903 | 12 | 456 | 1.0 | 123 | 11 | 13 | ||
| 22 |
(3) the firm battle array of card generation unit needed for analysis, and be assembled into overall just battle array with the element stiffness matrix of other unit of structure, the load suffered by integrated structure and boundary condition, the distortion of generating structure, as shown in Figure 3.
When the element stiffness with other unit is assembled into overall just battle array, will in strict accordance with principle corresponding to degree of freedom, namely nodes is by ascending arrangement, each node arranges by degree of freedom 1 ~ 6, to suffer restraints with single-degree-of-freedom at same node and namely in whole matrix, carry out disappearing row relax, simultaneously corresponding in loading matrix row carries out disappearing row relax too.
Claims (1)
1. three-dimensional is connected a numerical analysis method with two dimensional finite element unit, it is characterized in that, comprising:
(1) build three-dimensional element and two-dimensional cell respectively in the physical location of structure, give the rotational freedom that three-dimensional element is simulated in the junction of three bit locations and two-dimensional cell, set up the annexation of two-dimensional cell and new three-dimensional element by rigid-body system;
(2) card needed for the type of attachment of rigid-body system and the restriction relation generating structure analysis of rigid-body system;
(3) the firm battle array of card generation unit needed for analysis, and be assembled into overall just battle array with the element stiffness matrix of other unit of structure, the load suffered by integrated structure and boundary condition, the distortion of generating structure and the stress distribution of unit.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310660640.0A CN104699872A (en) | 2013-12-06 | 2013-12-06 | Numerical analysis method for connection between three-dimensional and two-dimensional finite element units |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310660640.0A CN104699872A (en) | 2013-12-06 | 2013-12-06 | Numerical analysis method for connection between three-dimensional and two-dimensional finite element units |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104699872A true CN104699872A (en) | 2015-06-10 |
Family
ID=53346992
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310660640.0A Pending CN104699872A (en) | 2013-12-06 | 2013-12-06 | Numerical analysis method for connection between three-dimensional and two-dimensional finite element units |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104699872A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105335575A (en) * | 2015-11-25 | 2016-02-17 | 中国航空工业集团公司沈阳飞机设计研究所 | Finite element modeling method of electronic antenna structural integrated laminated board |
| CN106202699A (en) * | 2016-07-07 | 2016-12-07 | 中国飞机强度研究所 | A kind of sensitivity method for solving under many displacement constraints |
| CN109726435A (en) * | 2018-12-04 | 2019-05-07 | 中国航空工业集团公司西安飞机设计研究所 | A kind of method of connector nail group's LOAD FOR |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6968271B1 (en) * | 2000-10-30 | 2005-11-22 | Board Of Trustees Of The University Of Illinois | Method and system for querying in a moving object database |
| CN102622472A (en) * | 2012-02-27 | 2012-08-01 | 西北工业大学 | Method for analyzing load and stress distribution and stress levels of composite mechanical connection pins |
-
2013
- 2013-12-06 CN CN201310660640.0A patent/CN104699872A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6968271B1 (en) * | 2000-10-30 | 2005-11-22 | Board Of Trustees Of The University Of Illinois | Method and system for querying in a moving object database |
| CN102622472A (en) * | 2012-02-27 | 2012-08-01 | 西北工业大学 | Method for analyzing load and stress distribution and stress levels of composite mechanical connection pins |
Non-Patent Citations (3)
| Title |
|---|
| 王家林: "钢筋混凝土结构空间有限元分析的体梁组合单元", 《工程力学》 * |
| 王家林等: "有限元分析中三维体元和梁元的非节点连接", 《重庆交通学院学报》 * |
| 钱仲良: "组合结构有限元分析的连接方法", 《北京水利电力经济管理学院学报》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105335575A (en) * | 2015-11-25 | 2016-02-17 | 中国航空工业集团公司沈阳飞机设计研究所 | Finite element modeling method of electronic antenna structural integrated laminated board |
| CN105335575B (en) * | 2015-11-25 | 2019-06-28 | 中国航空工业集团公司沈阳飞机设计研究所 | A kind of finite element modeling method of the structure-integrated laminate of electronic antenna |
| CN106202699A (en) * | 2016-07-07 | 2016-12-07 | 中国飞机强度研究所 | A kind of sensitivity method for solving under many displacement constraints |
| CN106202699B (en) * | 2016-07-07 | 2019-07-19 | 中国飞机强度研究所 | A kind of sensitivity method for solving under more displacement constraints |
| CN109726435A (en) * | 2018-12-04 | 2019-05-07 | 中国航空工业集团公司西安飞机设计研究所 | A kind of method of connector nail group's LOAD FOR |
| CN109726435B (en) * | 2018-12-04 | 2022-11-22 | 中国航空工业集团公司西安飞机设计研究所 | Method for calculating load of joint nail group |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103612688B (en) | Based on the automobile bottom diskware light weight method of many-body dynamics and topological optimization technology | |
| CN104699872A (en) | Numerical analysis method for connection between three-dimensional and two-dimensional finite element units | |
| CN105787199B (en) | A kind of impeller rotor blade is from hot to the blade profile conversion method of cold conditions | |
| CN102054106A (en) | Structure optimization design method and system | |
| Drager et al. | Smooth distributions are finitely generated | |
| CN106021797B (en) | A kind of design method of crane ring beam | |
| CN105313336A (en) | Thin-shell 3D printing optimizing method | |
| CN106294278B (en) | Adaptive hardware for dynamic reconfigurable array computing system is pre-configured controller | |
| CN103203612B (en) | Auxiliary assembling apparatus | |
| CN103838913A (en) | Finite element method for curve box girder curved bridge | |
| CN103729694A (en) | Method for solving flexible job-shop scheduling problem with improved GA based on polychromatic set hierarchical structure | |
| US9553591B2 (en) | Hybrid architecture for signal processing | |
| CN203133675U (en) | High-density server | |
| CN104376155A (en) | Micro truss structure design platform and method | |
| CN104992006A (en) | Method for establishing CAE linear analysis sub-assembly connection relation | |
| CN103699756A (en) | Loose coupling simulation system and method | |
| CN106096157B (en) | A kind of axial ventilation permanent magnet synchronous motor Calculation Method of Temperature Field | |
| CN109508486B (en) | Bolt-connected finite element high-precision simplified modeling method and system | |
| CN104573178A (en) | Finite element method for calculating integrated impeller strength | |
| CN105426640A (en) | Penalty function based multi-assembly structure system layout optimization design method | |
| CN108009370A (en) | A kind of structural stress sensitivity method for solving | |
| CN111881629A (en) | Pneumatic heat-structure heat conduction coupling nonlinear reduced order model method | |
| EP1943597B8 (en) | Parallel computer architecture for computation of particle interactions | |
| CN103530485B (en) | Full machine beam frame type construction reduces stiffness combine modeling method | |
| CN106991230A (en) | A kind of FEM model method for simplifying of flange arrangement |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150610 |