CN104123400B - Global Local details finite element methods based on force method - Google Patents
Global Local details finite element methods based on force method Download PDFInfo
- Publication number
- CN104123400B CN104123400B CN201310145236.XA CN201310145236A CN104123400B CN 104123400 B CN104123400 B CN 104123400B CN 201310145236 A CN201310145236 A CN 201310145236A CN 104123400 B CN104123400 B CN 104123400B
- Authority
- CN
- China
- Prior art keywords
- load
- analysis
- detail
- model
- point
- 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 belongs to finite element analysis technology, it is related to a kind of efficient, quick, accurate according to the completion partial structurtes details stress analysis of structure aggregate analysis result（Global‑Local）Method.Global Local details finite element methods based on force method are structural internal forces on the detail analysis zone boundary carried out by obtaining needed for overall model, as the outer force boundary condition of the detail analysis model then individually carried out, to carry out detail analysis.Its method relative to tradition using displacement as border, with it is efficient, quick, accurate the characteristics of.
Description
Technical field
The invention belongs to finite element analysis technology, it is related to a kind of efficient, quick, accurate according to structure aggregate analysis result
Complete partial structurtes details stress analysis（Global-Local）Method.
Background technology
Success Determined by Details, embodies especially prominent in Anti-Fatigue design.At this stage, details finite Element Stress point
Analysis turns into the Main Means of Anti fatigue Design.The detail analysis of labyrinth, particularly aircaft configuration are usually using extensive
Volume elements, while considering that geometry, material, the non-linear factor of boundary condition etc. are analyzed to complete.This detail analysis by
In scale of model, finite element solving device solve the limitation such as function generally can not, also It is not necessary to be completed in aggregate analysis model,
Detail analysis model is usually individually set up to complete.
The conventional method of current detail analysis is the finite element linear analysis for the general structure for first completing coarse grid, then thin
Save in model and detail analysis model foundation, and the extraction details mould from aggregate analysis result are carried out to the detailed structure of required analysis
Type boundary condition completes detail analysis.Individually overall and detail analysis can use existing various business CAE softwares more square
Just completion, but must be used from totality to local Global-Local technologies from totality to the transmission of local boundary condition
Come carry out.But there are many restrictions in current tradition Global-Local methods, in accommodation, efficiency, precision at present in work
Lack efficient, quick, accurate method in journey.
Often with Global-Local methods obtained according to aggregate analysis result interpolation using MSC/Patran as the engineering of representative
Detail model boundary displacement field is used as the boundary condition of detail analysis model and carries out detail analysis, this method be primarily present as
Lower problem：
a）Displacement field interpolation is only capable of being directed to translation displacements, and rotation displacement can not interpolation；
b）The displacement field that 1D, 2D unit are produced can not apply on 3D units, and its transition is manually completed, and efficiency is extremely inefficient；
c）Overall model modal displacement interpolation obtain displacement field can not ensure analyzed detailed structure institute it is loaded with it is real
Border institute is loaded consistent.
The content of the invention
The purpose of the present invention：
The present invention provides a kind of efficient, quick, accurate Global-Local methods by boundary condition of boundary force.
Technical scheme：
Global-Local details finite element methods based on force method are to be carried out by obtaining needed for overall model
Detail analysis zone boundary on structural internal force, as the outer force boundary condition of the detail analysis model then individually carried out, come
Carry out detail analysis.Global-Local details finite element methods based on force method, this method comprises the following steps：
A sets up aggregate analysis model：Structural separation gridding is carried out according to the rigidity of structure, correct reflection structure power transmission is set up
The overall model of relation, and the key for determining to need to carry out detail analysis according to overall stress analysis result, design specification requirement
Position；
B determines detail analysis region：According to the position of required detail analysis structure, size, with the annexation on periphery with
And appropriate transition region, the scope and boundary node and unit in detail analysis Local areas are determined in overall model, if overall mould
The mesh generation of over-simplification influences Local selection in type, then overall trellis is adjusted, to ensure suffered by Local areas
Load meets real force-transmitting relation；
C generates boundary element and boundary node：From aggregate analysis model boundary node and border are generated by specified format
Unit information；
D generates the node balance power of boundary node：The load working condition for needing to carry out detail analysis in aggregate analysis is determined, it is complete
Into the finite element analysis of the output of the node balance power containing boundary node of required analysis operating mode；
E obtains Local areas border internal force：The boundary node obtained according to step C is obtained with boundary element information and step D
Aggregate analysis result in node balance power extract Local areas border suffered by structural internal force, the internal force will as then solely
The outer load boundary condition of the vertical detail analysis carried out.
F sets up detail model：According to the structure type of detailed structure, annexation and stand under load mode situation, set up thin
Save analysis model；
G generates detail analysis file：Constrain and contact including outer load, rigid body displacement defined in the model generated in step F,
Annexation, and sequence is solved according to the feature definitions of Solve problems, and generate corresponding detail analysis file；
H generation detail analysis solve file：According to Local areas boundary node coordinate, edge load, detail model, manually
Input load point information, point of load branch link information, are connected using the automatic interconnection technique of edge load point, edge load point branch
Connection technology generates the annexation of border load(ing) point and detail model, and file format is generated as needed for specified finite element solver
Border load(ing) point connection, load card simultaneously merge the final detail calculation Study document of generation with the step G Study documents generated；
I completes detail analysis and solves calculating：Call selected finite element solving device to complete detail analysis to calculate, obtain details point
Analyse result.
The rigidity of structure described in A refers to：The ability that structure resistance is deformed by stress；
Overall model described in A refers to：The coarse grid of use large-size unit comprising full machine or whole part structure
Carry out the FEM model of the discrete foundation of structured grid；
Detail model described in A refers to：For the model of subrange in overall model, it will be considered that the true connection of structure
Relation factor, finite element grid set up using reduced size unit.It is mainly used in obtaining true, the accurate stress of local detail
Strain regime.
Local areas described in B refer to：Overall model Global's corresponding with the follow-up structure for individually carrying out detail analysis
Specific region；
Transition region described in B refers to：The precision of analysis at position is examined for needed for ensureing detail analysis,
The certain limit that Local areas should expand when choosing.Local areas should be that examination region adds transition region, and the scope of transition region is general
For 1.5~2 times of overall model size of mesh opening.
Boundary node described in C refers to：The node of outer load is born with Global boundary interfaces and Local areas in Local areas；
Boundary element described in C refers to：All units being connected in Local areas with boundary node；
Node balance power described in D refers to：All load suffered by individual node are enumerated in finite element analysis, and it includes phase
Unit even is to the active force of the node, external working loads, MPC power, SPC power.Can root accordingly generate result and
The free body FreeBody internal force of the virtual dividing in boundary element information acquisition Local areas of Local areas.
Structural internal force described in E refers to：Virtual cuts complete structure, the interaction between the structure of otch both sides
Power.The internal force will be used as the outer load of the detail analysis model subsequently individually carried out.
Detail model edge load described in E refers to：Effect of the Global areas to Local areas on Local areas boundary node
Structural internal force, and suffered outer in Local areas carry；
Detail model described in F refers to：Consider concrete form, the problem analysis type factor, the finite element of foundation of structure
Grid model；
Detail analysis file described in G refers to complete material property, manually outer load application, rigid body displacement in detail model
Constraint and contact, annexation definition, and complete to solve sequence definition, but the not detail model side including obtaining in step E
The file for submitting the solution analysis of finite element solving device of boundary's load.
The point of load that is manually entered described in H refers to：The load applied manually in H is walked is needed in detail model edge load
Lotus, being manually entered the point of load will not be attached in automatic connection, also not form load card, the point load is in step G
By form needed for analysis, manual definition is applied in detail model；
Point of load branch link information described in H refers to：The branch of load(ing) point when edge load point is connected using branch
Points and branch point coordinate information；
Automatic interconnection technique described in H refers to：Edge load point need to be connected with the node of detail model could complete border
The application of load, this technology is to be automatically performed edge load point and details according to loading point coordinates, detail model node coordinate
The technology of model node connection；
Branch's loading technique described in H refers to：According to the true connection of structure, artificial specified boundary point of load branch
It is connected to multiple spot with detail model, to increase the authenticity of edge load application.Edge load node location is by overall model
Determine, because overall model may carry out simplifying larger to real structure, grid is thicker, causes borderline power transmission and true
Structure has larger difference.
The technique effect of the present invention：
The present invention be directed to the details stress analysis of the particularly aircaft configuration of the structure with complicated annexation, there is provided one
Plant the Global-Local analysis methods using boundary force as boundary condition.Its method relative to tradition using displacement as border, tool
The characteristics of having efficient, quick, accurate.Specific advantage is as follows：
A. ensure that edge load suffered by detail model and institute in aggregate analysis are loaded completely the same, choose enough transition regions
The caused error of boundary condition need not be considered afterwards；
B. there is no any limitation to the cell type of detail model mesh generation, versatility is greatly increased；
C. force method is used, it is directly perceived, succinct, the Global-Local analyses for the CAE platform that collapses can be conveniently realized, this can be abundant
Using the strong point of all kinds of CAE softwares, this has been considerably improved structural analysis and design level.
D. the Global-Local total solutions formed are specification, general.It is easily programmed, so as to greatly improve details
The efficiency of analysis.
E. load(ing) point branch interconnection technique can be substantially reduced simplifies caused edge load point by the transition of overall model
The untrue property of edge load transmission, improves the analysis precision of borderline region.
Brief description of the drawings
Fig. 1 is the techniqueflow chart of the Global-Local details finite element methods based on force method；
Fig. 2 is cantilever beam schematic diagram in example；
Fig. 3 is the aggregate analysis model of 2D units in example；
Fig. 4 is overall model in example（Global）Middle detail analysis region（Local）Position；
Fig. 5 is boundary element in example in example and boundary node definition；
Fig. 6 is the solid element detail model that refines in example；
Fig. 7 is to apply the detail model of all kinds of boundary conditions and annexation in example；
Fig. 8 is overall stress analysis in example（Local areas）Maximum principal stress isogram；
Fig. 9 is details stress analysis maximum principal stress isogram in example.
Embodiment
Global-Local details finite element methods of the invention based on force method are in overall finite element model analysis
As a result basis, detail analysis model is used as using the boundary node internal force of detail analysis zone boundary in overall FEM model
Force boundary condition carries out detail analysis.Detail model only need to be corresponding with the structure of overall model, and in the choosing of detail analysis scope
Ensure enough transitional regions when taking（1.5~2 times of overall model size of mesh opening）It is ensured that interested details position should
Power analysis precision.
It is described in detail with reference to instantiation come the specific steps to the present invention.For convenience of the description selection of method
Simplest example.The example is center drilling cantilever beam structure as shown in Figure 2.Its design parameter is as follows：
It is long：400mm (x to), width 100mm（Y to）, thickness 2mm（Z to）, there is the centre bore of an a diameter of 30mm in plate；Plate is left
Clamped, right-hand member is held to apply the Z load for -50N that makes a concerted effort.
The most dangerous position that centre bore side is the structure can be determined by the overall stress analysis of coarse grid, low cost,
The purpose of detail analysis is exactly accurately to obtain dangerous position stress-strain state by setting up the 3D solid elements of refined net, is
Follow-up all kinds of intensive analysis provide indispensable data.
Global-Local details finite element methods of the invention based on force method solve the specific steps of the problem such as
Under：
Step 1：Set up the plate member aggregate analysis model of coarse grid（Such as Fig. 3）, apply boundary condition, complete to divide needed for details
Analyse the finite element analysis of the generation of power containing node balance of operating mode.Need to be limited in multi-state, large-sized model generation node scope with
Reduce the size of generation destination file；
Step 2:According to the position of required detail analysis structure, size, with the annexation on periphery in overall model really
Determine Local areas（Detail analysis）Scope and boundary node and unit, thick line position is Local areas in the middle of such as Fig. 4；
Step 3:Boundary node and boundary element are generated by Nastran bdf file formats（Thick line as shown in Figure 5 is border
Unit, puts as boundary node）；
Step 4:Node balance power in boundary node and boundary element and aggregate analysis result obtains institute of Local areas
The structural internal force received（Extracted by self-compiling program is automatic from destination file）.
Step 5:According to the concrete condition of detailed structure, detail analysis model corresponding with Fig. 4 thick line regions is set up（Such as
Shown in Fig. 6）；
Step 6:Generation includes outer load, the detail analysis file of rigid body displacement constraint（Do not consider in this other complicated
Nonlinear problem）；
Step 7:The details mould that Local areas boundary node coordinate, load file and the step 6 obtained according to step 3 is completed
The annexation of type file, generation border load(ing) point and detail model, and by analysis model file format generation connection, load card
And merge the final detail analysis file of generation with the Study document that step 6 is generated（Model is as shown in Figure 7）.This step can be by self-programming
Sequence is automatically performed.Because this example is simple in construction, do not applied and branch's loading technique using manual load；
Step 8:Complete detail analysis and solve calculating.
Analysis result is contrasted：
Fig. 8 is the maximum principal stress isogram in Local areas during overall stress is analyzed, and Fig. 9 is Local areas detail analysis
Maximum principal stress isogram, Fig. 8,9 comparative analyses show that edge load applied area in detail analysis model, stress state is multiple
Miscellaneous, this is mainly caused by boundary node application concentrfated load, outside transition region, and influence completely can be ignored.In around-the-hole stress collection
Its result of middle place is more accurate, meets subsequent analysis needs.
It should be pointed out that this example is intended to the explanation to method, it is not embodied to extensive, labyrinth issue handling excellent
Gesture.
Claims (7)
1. the Global-Local details finite element methods based on force method, it is characterized in that, this method comprises the following steps：
A sets up aggregate analysis model：Structural separation gridding is carried out according to the rigidity of structure, correct reflection structure force-transmitting relation is set up
Overall model, and according to overall stress analysis result, design specification requirement determine need carry out detail analysis key position；
B determines detail analysis region：According to the position of required detail analysis structure, size, with the annexation on periphery and suitable
Work as transition region, the scope and boundary node and unit in detail analysis Local areas are determined in overall model, if in overall model
The selection in the mesh generation influence Local areas of over-simplification, then be adjusted to overall trellis, to ensure institute of Local areas stand under load
Lotus meets real force-transmitting relation；
C generates boundary element and boundary node：From aggregate analysis model boundary node and boundary element are generated by specified format
Information；
D generates the node balance power of boundary node：The load working condition for needing to carry out detail analysis in aggregate analysis is determined, institute is completed
The finite element analysis of the output of the node balance power containing boundary node of operating mode need to be analyzed；
E obtains Local areas border internal force：According to the total of step C boundary node and the boundary element information obtained and step D acquisitions
Node balance power in body analysis result extracts the structural internal force suffered by Local areas border, and the internal force will be as then independently entering
The outer load boundary condition of capable detail analysis；
F sets up detail model：According to the structure type of detailed structure, annexation and stand under load mode situation, details point is set up
Analyse model；
G generates detail analysis file：Include outer load, rigid body displacement constraint and contact, connection defined in the model generated in step F
Relation, and sequence is solved according to the feature definitions of Solve problems, and generate corresponding detail analysis file；
H generation detail analysis solve file：According to Local areas boundary node coordinate, edge load, detail model, it is manually entered
Point of load information, point of load branch link information, skill is connected using the automatic interconnection technique of edge load point, edge load point branch
Art generates the annexation of border load(ing) point and detail model, and file format generates border as needed for specified finite element solver
Load(ing) point connection, load card simultaneously merge the final detail calculation Study document of generation with the step G Study documents generated；
I completes detail analysis and solves calculating：Call selected finite element solving device to complete detail analysis to calculate, obtain detail analysis knot
Really.
2. the Global-Local details finite element methods as claimed in claim 1 based on force method, it is characterized in that, A institutes
The overall model stated refers to：The coarse grid of use large-size unit comprising full machine or whole part structure carries out structured grid
The FEM model of discrete foundation.
3. the Global-Local details finite element methods as claimed in claim 1 based on force method, it is characterized in that, B institutes
The Local areas stated refer to：The specific region in overall model Global areas corresponding with the follow-up structure for individually carrying out detail analysis,
Local areas should be that examination region adds transition region；
Transition region described in B refers to：The certain limit that should expand when Local areas choose, the scope of transition region is overall model
1.5~2 times of size of mesh opening.
4. the Global-Local details finite element methods as claimed in claim 1 based on force method, it is characterized in that, C institutes
The boundary node stated refers to：The node of outer load is born with Global areas boundary interface and Local areas in Local areas；
Boundary element described in C refers to：All units being connected in Local areas with boundary node.
5. the Global-Local details finite element methods as claimed in claim 1 based on force method, it is characterized in that, E institutes
The edge load for the detail analysis stated refers to：Global areas are in the structure of the effect in Local areas on Local areas boundary node
Suffered outer load in power, and Local areas.
6. the Global-Local details finite element methods as claimed in claim 1 based on force method, it is characterized in that, F institutes
The detail model stated refers to：Consider concrete form, the problem analysis type factor, the FEM mesh of foundation of structure.
7. the Global-Local details finite element methods as claimed in claim 1 based on force method, it is characterized in that, H institutes
The point of load that is manually entered stated refers to：The load applied manually in H is walked is needed in detail model edge load, is manually entered
The point of load will not be attached in automatic connection, also not form load card, and the point load is by shape needed for analysis in step G
Formula, manual definition is applied in detail model；
Point of load branch link information described in H refers to：Branch's points of load(ing) point when edge load point is connected using branch
And branch point coordinate information；
Automatic interconnection technique described in H refers to：Edge load point is automatically performed according to loading point coordinates, detail model node coordinate
The technology being connected with detail model node；
Branch's interconnection technique described in H refers to：According to the true connection of structure, the artificial specified boundary point of load is branched off into many
Point is connected with detail model, to increase the authenticity of edge load application.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310145236.XA CN104123400B (en) | 2013-04-24 | 2013-04-24 | Global Local details finite element methods based on force method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310145236.XA CN104123400B (en) | 2013-04-24 | 2013-04-24 | Global Local details finite element methods based on force method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104123400A CN104123400A (en) | 2014-10-29 |
| CN104123400B true CN104123400B (en) | 2017-10-20 |
Family
ID=51768811
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310145236.XA Active CN104123400B (en) | 2013-04-24 | 2013-04-24 | Global Local details finite element methods based on force method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104123400B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105335575B (en) * | 2015-11-25 | 2019-06-28 | 中国航空工业集团公司沈阳飞机设计研究所 | A kind of finite element modeling method of the structure-integrated laminate of electronic antenna |
| CN105550464B (en) * | 2016-01-07 | 2018-06-08 | 燕山大学 | A kind of new submodel finite element method based on cut-boundary deformation constraint |
| CN109753687A (en) * | 2018-12-04 | 2019-05-14 | 中国航空工业集团公司西安飞机设计研究所 | A kind of combined type aeroelastic dynamics modeling method |
| CN111462837B (en) * | 2020-03-25 | 2023-06-16 | 华南理工大学 | Image model-oriented abaqus finite element computing Input file automatic generation method |
| CN112149213B (en) * | 2020-09-23 | 2023-04-18 | 中广核工程有限公司 | Method, device and equipment for transmitting finite element model grid data of nuclear island structure |
| CN113821955B (en) * | 2021-09-23 | 2022-08-19 | 中国航空工业集团公司西安飞机设计研究所 | Double-scale finite element iterative analysis method and device for local region of structure |
-
2013
- 2013-04-24 CN CN201310145236.XA patent/CN104123400B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN104123400A (en) | 2014-10-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104123400B (en) | Global Local details finite element methods based on force method | |
| US6799081B1 (en) | Fiber placement and fiber steering systems and corresponding software for composite structures | |
| CN103425842B (en) | A kind of parallel robot rapid development system and method | |
| CN102446241B (en) | Wing surface structural rigidity simulation method | |
| CN107066676A (en) | A kind of finite element automation modeling method based on satellite plate and shell structure | |
| CN102222140B (en) | Design drawing information extraction method of complex bridge structural analysis and modeling | |
| CN112163256A (en) | Method for establishing and analyzing template support model based on finite element analysis and BIM | |
| CN102222121A (en) | Finite element modeling method for modal simulation | |
| CN108090313A (en) | A kind of Complex Rock model of fissuration models recognition methods | |
| CN102663194A (en) | Heading machine collaborative simulation method and model thereof | |
| CN108153931B (en) | Method for calculating and analyzing mechanical properties of pipeline support and hanger of nuclear power plant by VBA (visual basic analysis) in combination with APDL (advanced persistent programming language) | |
| CN103678737A (en) | Lever system virtual assembly method | |
| CN112380650B (en) | Design method of structural part of working device | |
| CN115083541A (en) | Method for calculating stress intensity factor of orthotropic material based on ICVEFG method | |
| CN102375904A (en) | Parameterized computer modeling method for universal main girders of cranes | |
| CN102279126A (en) | Method for determining material performance parameter by combination of testing and CAE simulation | |
| CN103559361B (en) | A kind of optimization method of component strength and stress analysis method thereof | |
| CN105930578B (en) | A kind of dynamic and intelligent structure analysis method of power plant's plant designing | |
| CN105205213B (en) | A kind of lattice material Equivalent Mechanical performance analysis system | |
| CN104318030A (en) | Deformation-equivalence based automatic comprehensive deformation modulus calculating method for arch-dam foundations | |
| CN107239624B (en) | Three-dimensional map method for researching stress characteristics of working device of excavator | |
| Bournival et al. | A mesh-geometry based method for coupling 1D and 3D elements | |
| CN111259575B (en) | Finite element analysis design method for complex steel pipe node integral model | |
| CN103258086A (en) | Analysis method for tool stiffness | |
| Anand et al. | Process automation of simulation using Toolkit/Tool Command Language (TK/TCL) scripting |
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 | ||
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 610091 planning and Development Department of Chengdu aircraft design and Research Institute, 1610 Riyue Avenue, Qingyang District, Chengdu City, Sichuan Province Patentee after: AVIC CHENGDU AIRCRAFT DESIGN & Research Institute Address before: 610091 planning and Development Department of Chengdu aircraft design and Research Institute, 1610 Riyue Avenue, Qingyang District, Chengdu City, Sichuan Province Patentee before: AVIC CHENGDU AIRCRAFT DESIGN & Research Institute |