CN104778325B - Face load processing method and processing device based on surface cell - Google Patents
Face load processing method and processing device based on surface cell Download PDFInfo
- Publication number
- CN104778325B CN104778325B CN201510191027.8A CN201510191027A CN104778325B CN 104778325 B CN104778325 B CN 104778325B CN 201510191027 A CN201510191027 A CN 201510191027A CN 104778325 B CN104778325 B CN 104778325B
- Authority
- CN
- China
- Prior art keywords
- surface cell
- load
- solid element
- face
- unit
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provides a kind of face load processing method and processing devices based on surface cell, and the effect of being included in has establishes surface cell on the solid element of face load；According to preset rules, generation surface cell and the incidence relation of face load；The preset model of surface cell is determined according to above-mentioned incidence relation；Nodal point number is determined according to preset model；According to preset model and nodal point number generation pre-processing data；Pre-processing data is solved to obtain overall algebraic equation；Overall algebraic equation is solved, obtains displacement of joint to calculate corresponding stress-strain field according to displacement of joint, load processing is carried out using surface cell, following benefit can be brought：1st, the processing of opposite load can be completed in Xun Huan all surface unit, simplifies the processing procedure of opposite load；2nd, all different types of solid elements can share the processing that a kind of or a few surface cell completes boundary condition, realize the multiplexing of surface cell function code, reduce the Maintenance Development cost of calculation procedure.
Description
Technical field
The present invention relates to face load process field, in particular to the face load processing method based on surface cell and
Device.
Background technology
Finite element simulation is widely used in industrial design field, such as machinery, the energy, aerospace.Different types of list
First (Elements) is the core component that finite element program calculates, it is used for describing the physical characteristic in a short space region
And geometrical property, and the coefficient for the algebraic equation that different types of load abbreviation is needed for calculating.Unit is to realize stress
The basis of the physics field computation such as strain field, electromagnetic field, temperature field.
Mainly pass through 3D solid unit (such as tetrahedron element, hexahedral element, rib for the calculating of three-D space structure
Pole unit etc.) and Shell Finite Element realize.The face load of solid element can be pressure, temperature, potential and current density etc..
It needs these face load being converted into equivalent node load in calculating process, is finally all converted into equivalent node load always
The right-hand vector of body algebraic equation.
In conventional finite element simulation calculation flow process, what is handled inside unit is placed on to the processing of elemental area load,
Finite-element preprocessing program needs given side load to act on which face of some unit, and unit needs in calculating process
According to different face load types and the shape in different faces, different computing functions is realized.Due to there are many face load possibility
Complex situations, the processing of face load is part relatively complicated during Elementary Function is realized.Wherein, the complex situations bag of face load
It includes：According to the difference of face loading area, face load can be divided into loads in plane knead dough external applied load.According to face magnitude of load
Spatial distribution is different, and can be divided into uniform load and non-uniform load.And the if shape of cell surface different (such as three
Angular, quadrangle), the computational methods of solid element are also differed.
Inventor has found under study for action, carries out differentiation in a unit to different types of face load in the prior art
It handles and brings many troubles to unit calculating.One be store load information should store load effect unit number store again
The face sequence number of load effect, needs, according to the required geological information of these information searchings, to calculate and searched in calculating process
Journey is cumbersome, easily malfunctions.On the other hand, unit (such as tetrahedron element, hexahedral element, the prism of different geometries
Unit etc.) face load processing routine there are many repeat part cannot but be multiplexed, increase program code maintenance difficulty.
The content of the invention
It is an object of the invention to provide the face load processing method and processing devices based on surface cell, can simplify opposite load
The process that lotus is handled, and realize the multiplexing of code, the very big Maintenance Development cost for reducing calculation procedure.
In a first aspect, an embodiment of the present invention provides a kind of face load processing method based on surface cell, including：
Surface cell is established on the solid element that the effect obtained in sliding-model control has face load；
According to preset rules, generation surface cell and the incidence relation of face load；
According to the incidence relation of surface cell and face load, the preset model of surface cell is determined；
According to preset model, the nodal point number of node surface cell is determined；
According to the preset model and nodal point number of surface cell, pre-processing data is generated；
According to default constraints, pre-processing data is solved, obtains overall algebraic equation；
Overall algebraic equation is solved, obtains displacement of joint, so as to subsequently according to displacement of joint computational entity unit
Corresponding stress-strain field.
With reference to first aspect, an embodiment of the present invention provides the first possible embodiment of first aspect, wherein,
The effect that sliding-model control obtains has establishes surface cell on the solid element of face load, including：
Determine that effect that sliding-model control obtains has the preset kind of solid element of face load；
Solid element is decomposed according to preset kind, obtains preset model solid element；
According to preset order, each corresponding surface cell of preset model solid element is established.
The possible embodiment of with reference to first aspect the first, an embodiment of the present invention provides second of first aspect
Possible embodiment, wherein, according to preset rules, generation surface cell and the incidence relation of face load, including：
According to preset rules, the incidence relation of generation surface cell and preset model solid element；
According to surface cell and the incidence relation of preset model solid element, generation surface cell is associated with face load
System.
Second of possible embodiment with reference to first aspect, an embodiment of the present invention provides the third of first aspect
Possible embodiment, wherein, according to default constraints, pre-processing data is solved, obtains overall algebraic equation
Including：
All solid elements are cycled, calculate the corresponding stiffness matrix of each solid element；
All surface unit is cycled, the corresponding face load of each surface cell is converted into equivalent node load；
Stiffness matrix and equivalent node load are subjected to general assembly, obtain overall algebraic equation.
The third possible embodiment with reference to first aspect, an embodiment of the present invention provides the 4th kind of first aspect
Possible embodiment, wherein, Xun Huan all surface unit includes：
According to the predetermined number of the incidence relation and solid element of surface cell and solid element, surface cell is compiled
Number；
According to the surface cell number that number obtains, all surface unit is cycled.
The 4th kind of possible embodiment with reference to first aspect, an embodiment of the present invention provides the 5th kind of first aspect
Possible embodiment, wherein, the corresponding face load of each surface cell is converted into equivalent node load includes：
According to formulaIt will act at the force per unit area on surface cellIt is converted into equivalent knot
Point power fe；
Wherein, feRepresent the corresponding equivalent node force of equivalent node load；Expression acts on load above surface cell
Distributed force, and Represent three components of the force per unit area under rectangular coordinate system；
Represent the area of space of surface cell, area of space includes any one in curved surface and plane；N represents surface cell shape function
Matrix；N1,N2,N3,...,NmFor the shape function of each node of correspondence, 4≤m≤8, and
With reference to first aspect, an embodiment of the present invention provides the 6th kind of possible embodiment of first aspect, wherein,
The effect that sliding-model control obtains, which has, establishes on the solid element of face load before surface cell, further includes：
According to the physical property and geometric areas of test object, the geometrical model of test object is established；
Sliding-model control is carried out to geometrical model, obtains multiple solid elements for acting on and having face load.
Second aspect, the embodiment of the present invention additionally provide a kind of face load-handling device based on surface cell, including：
First establishes unit, and surface cell is established on the solid element that the effect obtained in sliding-model control has face load；
First generation unit, for according to preset rules, the surface cell of unit foundation and face load to be established in generation first
Incidence relation；
First determination unit, for the surface cell that is generated according to the first generation unit and the incidence relation of face load, really
Determine the preset model of surface cell；
Second determination unit for the preset model determined according to the first determination unit, determines the knot of node surface cell
Points；
Second generation unit, the preset model and second of the surface cell for being determined according to the first determination unit determine list
The definite nodal point number of member, generates pre-processing data；
First solves unit, for according to default constraints, to the pre-processing data of the second generation unit generation into
Row solves, and obtains overall algebraic equation；
Second solves unit, and the overall algebraic equation for being obtained to the first solution unit solves, and obtains node position
It moves, so as to subsequently according to the corresponding stress-strain field of displacement of joint computational entity unit.
With reference to second aspect, an embodiment of the present invention provides the first possible embodiment of second aspect, wherein,
One, which establishes unit, includes：
Determination subelement determines that effect that sliding-model control obtains has the preset kind of solid element of face load；
Subelement is decomposed, the preset kind determined according to determination subelement decomposes solid element, obtains default mould
Type solid element；
Subelement is established, according to preset order, each is established and decomposes the preset model solid element pair that subelement obtains
The surface cell answered.
With reference to the first possible embodiment of second aspect, an embodiment of the present invention provides second of second aspect
Possible embodiment, wherein, the first generation unit includes：
First generation subelement, for according to preset rules, generation surface cell to be associated with preset model solid element
Relation；
Second generation subelement, for the surface cell and preset model solid element generated according to the first generation subelement
Incidence relation, the incidence relation of generation surface cell and face load.
Face load processing method and processing device provided in an embodiment of the present invention based on surface cell, using in sliding-model control
Obtained effect has establishes surface cell on the solid element of face load；According to preset rules, generation surface cell and face load
Incidence relation；According to the incidence relation of surface cell and face load, the preset model of surface cell is determined；According to default mould
Type determines the nodal point number of node surface cell；According to the preset model and nodal point number of surface cell, pre-processing data is generated；Root
According to default constraints, pre-processing data is solved, obtains overall algebraic equation；Overall algebraic equation is asked
Solution, obtains displacement of joint, so as to subsequently according to the corresponding stress-strain field of displacement of joint computational entity unit, and in the prior art
The differentiation that carried out to different types of face load in a unit handle and to unit calculating many troubles brought to compare, profit
Load processing is carried out with surface cell, following benefit can be brought by carrying out load processing using surface cell：First, store load
As long as storing the surface cell number of load effect during information, opposite can be completed in Xun Huan all surface unit in calculating process
The processing of load simplifies the process that opposite load is handled；Second, all different types of solid elements can share one
The processing of most boundary conditions can be completed in kind or a few surface cell, realizes surface cell function most possibly
The multiplexing of code reduces the Maintenance Development cost of calculation procedure.
For the above objects, features and advantages of the present invention is enable to be clearer and more comprehensible, preferred embodiment cited below particularly, and coordinate
Appended attached drawing, is described in detail below.
Description of the drawings
It in order to illustrate the technical solution of the embodiments of the present invention more clearly, below will be to needed in the embodiment attached
Figure is briefly described, it should be understood that the following drawings illustrates only certain embodiments of the present invention, therefore is not construed as pair
The restriction of scope, for those of ordinary skill in the art, without creative efforts, can also be according to this
A little attached drawings obtain other relevant attached drawings.
Fig. 1 shows a kind of flow for face load processing method based on surface cell that the embodiment of the present invention is provided
Figure；
Fig. 2 shows the flow for another face load processing method based on surface cell that the embodiment of the present invention is provided
Figure；
Fig. 3 shows the flow for another face load processing method based on surface cell that the embodiment of the present invention is provided
Figure；
Fig. 4 shows the flow for another face load processing method based on surface cell that the embodiment of the present invention is provided
Figure；
Fig. 5 shows the flow for another face load processing method based on surface cell that the embodiment of the present invention is provided
Figure；
Fig. 6 shows that a kind of structure for face load-handling device based on surface cell that the embodiment of the present invention is provided is shown
It is intended to；
Fig. 7 shows in a kind of face load-handling device based on surface cell that the embodiment of the present invention is provided that first builds
The structure diagram of vertical unit；
Fig. 8 shows the first life in a kind of face load-handling device based on surface cell that the embodiment of the present invention is provided
Into the structure diagram of unit；
Fig. 9 shows in a kind of face load-handling device based on surface cell that the embodiment of the present invention is provided that first asks
Solve the structure diagram of unit；
Figure 10 shows in another face load-handling device based on surface cell that the embodiment of the present invention provided the
One solves the structure diagram of unit；
Figure 11 is shown in another face load-handling device based on surface cell that the embodiment of the present invention is provided
Structure diagram；
Figure 12 shows pressure in a kind of face load processing method based on surface cell that the embodiment of the present invention is provided
Distribution schematic diagram at node.
Specific embodiment
Below in conjunction with attached drawing in the embodiment of the present invention, the technical solution in the embodiment of the present invention is carried out clear, complete
Ground describes, it is clear that described embodiment is only part of the embodiment of the present invention, instead of all the embodiments.Usually exist
The component of the embodiment of the present invention described and illustrated in attached drawing can configure to arrange and design with a variety of herein.Cause
This, the detailed description of the embodiment of the present invention to providing in the accompanying drawings is not intended to limit claimed invention below
Scope, but it is merely representative of the selected embodiment of the present invention.Based on the embodiment of the present invention, those skilled in the art are not doing
Go out all other embodiments obtained on the premise of creative work, belong to the scope of protection of the invention.
Finite element analysis (FEA, Finite Element Analysis) is the method using mathematical approach to actual physical
System (geometry and load working condition) is simulated.Also utilize simple and interaction element, i.e. unit, it is possible to limited
The unknown quantity of quantity goes to approach the real system of unlimited unknown quantity.At present, finite element simulation is widely used in industrial design field,
Such as machinery, the energy, aerospace.Different types of unit (Elements) is that the core component that finite element program calculates (is asked
Solve the discretization in domain), it is used for describing physical characteristic and geometrical property in a short space region, and by different type
Load abbreviation for calculate need algebraic equation coefficient.Unit is to realize the physics such as stress-strain field, electromagnetic field, temperature field
The basis of field computation.
In conventional finite element simulation calculation flow process, what is handled inside unit is placed on to the processing of elemental area load.
Finite-element preprocessing program, which needs given side load to act on which face of some unit, (will know that load is made simultaneously
Solid element is numbered and the unit inner face sequence number of load institute acting surface)；In calculating process, unit is needed according to different
Face load type (i.e. mesh shape) and different face shape (plane or curved surface in this way), realize different computing functions.By
May be there are many complex situations in face load, the processing of face load is part relatively complicated during Elementary Function is realized.Wherein, face
The complex situations of load include：According to the difference of face loading area, face load can be divided into loads in plane knead dough external applied load.
It is different according to the spatial distribution of face magnitude of load, and uniform load and non-uniform load can be divided into.And if cell surface
Different (such as the triangle, quadrangle) computational methods of shape also differ.
Specifically, carrying out in a unit differentiation to different types of face load and handling many is brought to unit calculating
Trouble.One is to store load information to store the face sequence number that the unit number that load acts on stores load effect again, is calculated
It is needed in the process according to the required geological information of these information searchings, and is calculated according to these geological information opposite load
Processing；Wherein, the calculating in the above process and search procedure are cumbersome, easily malfunction.On the other hand, different geometries
The face load processing routine of unit (such as tetrahedron element, hexahedral element and prism elements etc.) there are many part repeated but not
It can be multiplexed, if tetrahedron and hexahedron are all made of multiple quadrangles, identical four are included in the load processing routine of face
Side shape part, and the calling of same section program can not be realized in unit, identical code can only be write again, increase work
The difficulty of automation program's code maintenance.
The present invention proposes that the method for face load, face is uniformly processed with the layer of surface unit for being attached to 3D solid surface
The processing of load is no longer placed on inside unit and carries out.Wherein, so-called surface cell refers to that thickness is zero, for face load is converted
For a kind of unit of equivalent node load, shape can be triangle, quadrangle etc., also, the specific geometry of surface cell
Feature is determined by the surface characteristics of accompanying solid element.To handle complicated boundary condition in calculating process, cooperation is real
Body unit simplifies load process flow.
The processing method and processing device of face load in Finite Element provided by the invention is carried out using surface cell at load
Managing the benefit brought has：First, as long as storing the surface cell number of load effect during storage load information, in calculating process
The processing of opposite load can be completed in cycle all surface unit, believes without the geometry of searching needs is entered back into solid element
Breath, simplifies the process that opposite load is handled；Second, all different types of solid elements can share a kind of or a small number of
The processing of most boundary conditions can be completed in several surface cells, realizes answering for surface cell function code most possibly
With reducing the Maintenance Development cost of calculation procedure.
Specific embodiments of the present invention are as follows：
With reference to figure 1, a kind of face load processing method based on surface cell, including：
101st, surface cell is established on the solid element that the effect obtained in sliding-model control has face load.
Specifically, in finite element analysis, pre-treatment modeling is carried out to test object first, and to the geometry of test object
Model carries out sliding-model control, and the effect of obtaining has the solid element of face load；Wherein, the solid element is multiple, to each reality
Body unit numbers that such as which solid element obtained；In the present embodiment, solid element is preferably polygon, when solid element is polygon
During shape (such as triangle and quadrangle), each face of the solid element is numbered, obtains solid element face sequence number, so
The number of the solid element of face load effect and the face sequence number of the solid element of the number are determined afterwards, are obtained solid element and are carried with face
The incidence relation of lotus.
Then surface cell is established on the solid element obtained in above-mentioned sliding-model control, so as to single subsequently according to the surface
The calculating and processing of first carry out face load.
Wherein, so-called surface cell refers to that thickness is zero, single for face load to be converted into one kind of equivalent node load
Member, shape can be triangle, quadrangle etc., also, the specific geometric properties of surface cell are by accompanying solid element
Surface characteristics determine.
102nd, according to preset rules, generation surface cell and the incidence relation of face load.
Specifically, preset rules can be preset, for example some surface cell is made to replace some solid element
Some face and another solid element is made to replace another face etc. of the solid element.The preset rules specifically can root
According to needing to be configured, the present invention is not particularly limited this.According to preset rules, the surface cell of foundation and entity list are recorded
The incidence relation of member, and according to the surface cell and the incidence relation generation surface cell of solid element and associating for face load
Face load is acted on which number face of some solid element, is converted to face load acting on surface by system, the above process
On unit, and which number face of some solid element is replaced with the surface cell.
103rd, according to the incidence relation of surface cell and face load, the preset model of surface cell is determined.
Specifically, the preset model of surface cell is for which number face of which solid element to be replaced, therefore this is default
Under the action of its own, shape can arbitrarily be set model, such as quadrangle, triangle.
104th, according to preset model, the nodal point number of node surface cell is determined.
Specifically, in the case where above-mentioned preset model is determined, it is by quadrangular surface unit of surface cell in the present embodiment
Example illustrates, and is such as knowing that surface cell is quadrangle, and it is intended to indicate that No. 1 face of No. 1 solid element, then can be with
According to above-mentioned preset model, determine whether there is face load on No. 1 face of No. 1 solid element, if so, determining face load work
With the nodal point number on quadrangular surface unit.
Wherein, which is preferably greater than or equal to the basic angle point number of polygon, such as when surface cell is quadrangular surface
During unit, nodal point number is more than or equal to 4；When surface cell is triangular surface unit, nodal point number is more than or equal to 3 surface cells
The nodal point number on solid element surface is depended on node.For single order solid element, Surface Nodes number is equal to surface polygon
Side number, for high-order solid element, Surface Nodes number is more than the side number of surface polygon.The nodal point number of surface cell is in contrast
It should.
105th, according to the preset model and nodal point number of surface cell, pre-processing data is generated.
Specifically, node number, each surface are single on the number of pre-processing data, that is, surface cell, each surface cell
Corresponding face load quantity in site position and each surface cell in member etc..
It is which number face with surface cell instead of some solid element in the present embodiment, therefore according to the surface cell
With the incidence relation of face load, you can to generate above-mentioned etc. pre-processing data, subsequently to be carried out according to the pre-processing data
It solves and calculates.
Wherein, the preset model of surface cell is determined according to the preset kind of solid element, in the present embodiment
Surface cell is which number face for representing some solid element, then the type in which number face of the solid element is exactly surface cell
Preset kind.
106th, according to default constraints, pre-processing data is solved, obtains overall algebraic equation.
Specifically, default constraints is to handle degree of freedom restrained in geometrical model；In fact, it is asked in actual
It solves in calculating process, is to need to constrain some variables, using these variables as constant coefficient, i.e., is replaced with simple question
The problem of the problem of complicated, solves to simple question, ability final solving complexity；The process is to geometrical model
In design variables be defined as constant coefficient, and pre-processing data is solved according to fixed constant coefficient, and by the result of solution
General assembly is carried out, overall algebraic equation is then obtained, subsequently the totality algebraic equation is solved, you can obtain required knot
Fruit --- i.e. displacement of joint.
107th, overall algebraic equation is solved, obtains displacement of joint, so as to subsequently according to displacement of joint computational entity
The corresponding stress-strain field of unit.
It is solved specifically, the process is the overall algebraic equation assembled in step 104, the overall algebraically after assembling
Equation is the equation for containing only displacement of joint variable, and displacement of joint can be calculated according to mathematical operation.
Then, according to the general-purpose algorithm in the displacement of joint and combination industry acquired, you can be calculated corresponding reality
The corresponding stress-strain field of body unit, so that staff can determine that the tension of test object is strong according to the stress-strain field
Degree.Wherein, the distribution situation of the stress-strain field, that is, stress value, that is, all stress values.
The processing method of face load in Finite Element provided by the invention carries out load processing using surface cell and brings
Benefit have：First, as long as storing the surface cell number of load effect during storage load information, institute is cycled in calculating process
There is the processing that opposite load can be completed in surface cell, without the geological information of searching needs is entered back into solid element, simplify
The process that opposite load is handled；Second, all different types of solid elements can share a kind of or a few table
The processing of most boundary conditions can be completed in face unit, realizes the multiplexing of surface cell function code, drop most possibly
The Maintenance Development cost of low calculation procedure.
With reference to figure 2, in the present invention, it is contemplated that subsequent convenience of calculation and ask repeat that program code can be used
It inscribes, in the Finite Element in the processing method of face load, the effect obtained in step 101 in sliding-model control has face load
Surface cell is established on solid element, is specifically comprised the following steps：
201st, determine that effect that sliding-model control obtains has the preset kind of solid element of face load.
Specifically, the preset kind of solid element includes：3D solid unit (such as tetrahedron element, hexahedral element, rib
Pole unit etc.) and Shell Finite Element, it is first determined the concrete type of the solid element, so as to subsequently according to the definite type at this
Surface cell is established on solid element.
202nd, solid element is decomposed according to preset kind, obtains preset model solid element.
Specifically, after the type of solid element is determined, according to the corresponding different decomposition method of various types to whole
A solid element is decomposed；It is such as to quadrangle, i.e., the mesh shape of the face load per one side, every according to six faces of quadrangle
Whether it is on one side that plane, the corresponding face load in the face are that the corresponding face load of loads in plane or face external applied load and the face is equal
Even load also non-uniform load decomposes solid element；Such as to triangle, i.e., according to the four sides of triangle, the face per one side
Whether the mesh shape of load per one side is that plane, the corresponding face load in the face are loads in plane or face external applied load and should
The corresponding face load in face is that uniform load also non-uniform load decomposes solid element.
Wherein, each model that above-mentioned decomposition obtains all corresponds to a preset model solid element, writes down the default mould
The quantity of type solid element, subsequently to establish the surface cell identical with the preset model solid element of the quantity.
Certainly, specific isolation is not limited to the above method, other, which are in, carries out in decomposition range solid element
All within the scope of the present invention.
203rd, according to preset order, each corresponding surface cell of preset model solid element is established.
Specifically, after decomposition, multiple preset models are obtained, establish each preset model solid element to correspondence
Surface cell.Wherein, the quantity of surface cell is identical with the quantity of preset model solid element, and the quantity of surface cell also may be used
To be differed with preset model solid element.In the case that different, determine that effect has such as in preset model solid element
Then the preset model solid element of face load establishes surface cell i.e. on the preset model solid element for having face load is acted on
It can.
With reference to figure 3, in of the invention, the specific implementation of step 103 is as follows：
301st, according to preset rules, the incidence relation of generation surface cell and preset model solid element.
302nd, according to surface cell and the incidence relation of preset model solid element, surface cell and the pass of face load are generated
Connection relation.
Specifically, the incidence relation of surface cell and solid element is first determined in the embodiment of the present invention according to preset rules,
Since solid element is relevant with face load previously known, thus by both above-mentioned incidence relation can obtain surface cell with
The incidence relation of face load.
The processing procedure of opposite load of the present invention is by surface cell agent entity unit, passes through table in calculating process
Face cell processing complicated boundary condition, and it is made to coordinate solid element, simplify load process flow, with reference to figure 4, wherein, it uses
The specific processing mode of face load is as follows：
401st, all solid elements are cycled, calculate the corresponding stiffness matrix of each solid element.
In the present embodiment, the mode of calculated rigidity matrix is identical with the mode of existing computational entity unit, this step does not exist
It illustrates.
402nd, all surface unit is cycled, the corresponding face load of each surface cell is converted into equivalent node carries
Lotus.
Specifically, for the ease of cycling all surface unit, first according to surface cell and the incidence relation of solid element
And the predetermined number of solid element, surface cell is numbered, obtains surface cell number；Then the table obtained according to number
Face unit number cycles all surface unit, and is equivalent node load by the definite corresponding face load transfer of surface cell.
Surface cell is cycled in the present embodiment, is turned face load by certain algorithm according to node coordinate on surface cell
It is changed to equivalent node load.Compared to the face sequence number for searching the solid element number for having face load effect and acting surface, according to face
Sequence number judges the geometric coordinate of node in the shape and lookup face in face.According to face shape and node coordinate by calculating face load
The traditional process of equivalent node load is converted to, eliminates the trouble of lookup face and node information.
403rd, stiffness matrix and equivalent node load are subjected to general assembly, obtain overall algebraic equation.
Specifically, will be to the left end item of stiffness matrix that each solid element is calculated as algebraic equation, it will
The replacer of above-mentioned solid element --- multiple equivalent node information that the corresponding face LOAD FOR of each surface cell obtains
As the right-hand vector of algebraic equation, multiple left end items and right-hand vector are obtained, makes the general assembly corresponding with multiple right-hand vectors of multiple left end items
Into overall algebraic equation.
By surface cell to be right exemplified by quadrangular surface unit in the present embodiment, the corresponding face of each surface cell is carried
The specific method that lotus is converted into equivalent node load illustrates, including：
According to formulaIt will act at the force per unit area on surface cellIt is converted into equivalent knot
Point power fe；
Wherein, feRepresent the corresponding equivalent node force of equivalent node load；Expression acts on load above surface cell
Distributed force, and Represent three components of the force per unit area under rectangular coordinate system；
Represent the area of space of surface cell, area of space includes any one in curved surface and plane；N represents surface cell shape function
Matrix；N1,N2,N3,...,NmFor the shape function of each node of correspondence, 4≤m≤8, and
In the present embodiment, surface cell is illustrated by taking quadrangular surface unit as an example, and the present invention is proposed with being attached to
The quadrangular surface unit on 3D solid surface, the face load distributed force in being calculated for processing structure, such as pressure.Face load
Processing is no longer placed on inside 3D solid unit and carries out.The thickness of surface cell is zero, it is only used for face load being converted into
Imitate node load.
In addition, with reference to established on the solid element that figure 5, the effect obtained in sliding-model control have face load surface cell it
Before, it further includes：
501st, according to the physical property and geometric areas of test object, the geometrical model of test object is established.
502nd, sliding-model control is carried out to geometrical model, obtains multiple solid elements for acting on and having face load.
Below by taking quadrangular surface unit as an example, at the quadrangular surface unit matching solid element opposite load
The process of reason is briefly described, with reference to figure 12：
The shape of quadrangular surface unit can be quadrangle or space curved surface quadrangle in plane；Quadrangular surface unit
Nodal point number can be 4,5,6,7,8, wherein 1-4 nodes are side end point, and No. 5-8 is node in side, (its as shown in figure 12
In, Q7 and Q are not shown in FIG. 12).
Figure 12 quadrangular surface units, node are represented with Q, are distributed size and Orientation of the pressure at node arrow shown in P
Head represents.
Quadrangular surface unit is according to formulaIt will act at the force per unit area on surface cell
It is converted into equivalent node force fe。
Wherein,
Each variable meaning is as follows：
--- the area of space (curved surface or plane) of surface cell；
N --- surface cell form function matrix；N1,N2,N3,...,NmFor the shape function of each node of correspondence, 4≤m≤8
And And--- force per unit area
Three components under rectangular coordinate system.
Below by taking linear static calculation process as an example, illustrate how quadrangular surface unit coordinates solid element to complete
Its computing function.As Figure 1-Figure 5.In the calculation process based on surface cell, pre-treatment program defines face load effect
On corresponding quadrangular surface unit.Solid element is cycled during solution first, completes the Rigidity Calculation of solid element and rigidity square
Battle array assembling, then cycles quadrangular surface unit, and the distributed load processing that will act on quadrangle curved surface is corresponding equivalent
Equivalent node load is added in integrated solution equation by node load, solver, solves equation.According to integrated solution equation
Solution obtains the displacement of unit node, finally calculates stress-strain field of concern according to unit displacement of joint.
Face load processing method provided in an embodiment of the present invention based on surface cell, with of the prior art to inhomogeneity
The face load of type carried out in a unit differentiation handle calculating to unit many troubles is brought to compare, utilize surface cell
Load processing is carried out, following benefit can be brought by carrying out load processing using surface cell：First, as long as during storage load information
The surface cell number of load effect is stored, the place that opposite load can be completed in all surface unit is cycled in calculating process
Reason, simplifies the process that opposite load is handled；Second, all different types of solid elements can share a kind of or a small number of
The processing of most boundary conditions can be completed in several surface cells, realizes answering for surface cell function code most possibly
With reducing the Maintenance Development cost of calculation procedure.
With reference to figure 6, the present invention provides a kind of face load-handling device based on surface cell, including：
First establishes unit 11, and surface list is established on the solid element that the effect obtained in sliding-model control has face load
Member；
First generation unit 12, for according to preset rules, surface cell and the face of the foundation of unit 11 to be established in generation first
The incidence relation of load；
First determination unit 13, for associating according to the surface cell that the first generation unit 12 generates and face load
System determines the preset model of surface cell；
Second determination unit 14 for the preset model determined according to the first determination unit 13, determines node surface cell
Nodal point number；
Second generation unit 15, for the preset model according to the definite surface cell of the first determination unit 13 and second really
14 definite nodal point number of order member, generates pre-processing data；
First solves unit 16, for according to default constraints, to the preprocessing number of the second generation unit 15 generation
According to being solved, overall algebraic equation is obtained；
Second solves unit 17, and the overall algebraic equation for being obtained to the first solution unit 16 solves, and is tied
Point displacement, so as to subsequently according to the corresponding stress-strain field of displacement of joint computational entity unit.
Further, should be based in the face load-handling device of surface cell with reference to figure 7, first, which establishes unit 11, includes：
Determination subelement 111, for determining that effect that sliding-model control obtains has the default class of the solid element of face load
Type.
Subelement 112 is decomposed, the preset kind for being determined according to determination subelement 111 decomposes solid element,
Obtain preset model solid element；
Subelement 113 is established, for according to preset order, establishing each and decomposing the preset model that subelement 112 obtains
The corresponding surface cell of solid element.
Further, with reference to figure 8, in the Finite Element in the processing unit of face load, the first generation unit 12 includes：
First generation subelement 121, for according to preset rules, the pass of generation surface cell and preset model solid element
Connection relation；
Second generation subelement 122, for according to the surface cell that the first generation subelement 121 generates and preset model reality
The incidence relation of body unit, generation surface cell and the incidence relation of face load.
Further, with reference to figure 9, the processing unit of face load in the Finite Element, the first solution unit 16 includes：
Computation subunit 161 for cycling all solid elements, calculates the corresponding stiffness matrix of each solid element.
For cycling all surface unit, the corresponding face load of each surface cell is converted for conversion subunit 162
For equivalent node load；
General assembly subelement 163, what stiffness matrix and conversion subunit 162 for computation subunit 161 to be calculated were converted
Equivalent node load carries out general assembly, obtains overall algebraic equation.
Further, with reference to figure 10, in the Finite Element in the processing unit of face load, the first solution unit 16 also wraps
It includes：
Number subelement 164, for the default of the incidence relation according to surface cell and solid element and solid element
Number numbers surface cell；
Subelement 165 is cycled, for numbering obtained surface cell number according to number subelement 164, cycles all surface
Unit.
Further, in the Finite Element in the processing unit of face load, conversion subunit 162 is specifically used for, according to
FormulaIt will act at the force per unit area on surface cellIt is converted into equivalent node force fe；
Wherein, feRepresent the corresponding equivalent node force of equivalent node load；Expression acts on load above surface cell
Distributed force, and Represent three components of the force per unit area under rectangular coordinate system；
Represent the area of space of surface cell, area of space includes any one in curved surface and plane；N represents surface cell shape function
Matrix；N1,N2,N3,...,NmFor the shape function of each node of correspondence, 4≤m≤8, and
Further, with reference to figure 11, further included in the Finite Element in the processing unit of face load：
Second establishes unit 18, for the physical property and geometric areas according to test object, establishes the several of test object
What model；
Sliding-model control unit 19, the geometrical model for establishing the foundation of unit 18 to second carry out sliding-model control, obtain
To multiple solid elements for acting on and having face load.
Face load-handling device provided in an embodiment of the present invention based on surface cell, using what is obtained in sliding-model control
Effect has establishes surface cell on the solid element of face load；According to preset rules, generation surface cell is associated with face load
Relation；According to the incidence relation of surface cell and face load, the preset model of surface cell is determined；According to preset model, determine
The nodal point number of node surface cell；According to the preset model and nodal point number of surface cell, pre-processing data is generated；According to default
Constraints solves pre-processing data, obtains overall algebraic equation；Overall algebraic equation is solved, is tied
Point displacement, so as to subsequently according to the corresponding stress-strain field of displacement of joint computational entity unit, with of the prior art to difference
The face load of type carries out differentiation and handles many troubles to unit calculating to be brought to compare in a unit, single using surface
Member carries out load processing, and following benefit can be brought by carrying out load processing using surface cell：First, store load information when only
The surface cell number of load effect is stored, the place that opposite load can be completed in all surface unit is cycled in calculating process
Reason, simplifies the process that opposite load is handled；Second, all different types of solid elements can share a kind of or a small number of
The processing of most boundary conditions can be completed in several surface cells, realizes answering for surface cell function code most possibly
With reducing the Maintenance Development cost of calculation procedure.
In addition, each functional unit in each embodiment of the present invention can be integrated in a processing unit, it can also
That unit is individually physically present, can also two or more units integrate in a unit.
If the function is realized in the form of SFU software functional unit and is independent production marketing or in use, can be with
It is stored in a computer read/write memory medium.Based on such understanding, technical scheme is substantially in other words
The part contribute to the prior art or the part of the technical solution can be embodied in the form of software product, the meter
Calculation machine software product is stored in a storage medium, is used including some instructions so that a computer equipment (can be
People's computer, server or network equipment etc.) perform all or part of the steps of the method according to each embodiment of the present invention.
And foregoing storage medium includes：USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic disc or CD.
The above description is merely a specific embodiment, but protection scope of the present invention is not limited thereto, any
Those familiar with the art in the technical scope disclosed by the present invention, can readily occur in change or replacement, should all contain
Lid is within protection scope of the present invention.Therefore, protection scope of the present invention described should be subject to the protection scope in claims.
Claims (9)
1. a kind of face load processing method based on surface cell, which is characterized in that including：
Surface cell is established on the solid element that the effect obtained in sliding-model control has face load；
According to preset rules, the incidence relation of the surface cell and the face load is generated；
According to the incidence relation of the surface cell and the face load, the preset model of the surface cell is determined；
According to the preset model, the nodal point number of the surface cell is determined；
According to the preset model and nodal point number of surface cell, pre-processing data is generated；
According to default constraints, the pre-processing data is solved, obtains overall algebraic equation；
The overall algebraic equation is solved, obtains displacement of joint, so as to subsequently according to calculating the displacement of joint
The corresponding stress-strain field of solid element；
It is described according to default constraints, the pre-processing data is solved, obtaining overall algebraic equation includes：Xun Huan
All solid elements calculate the corresponding stiffness matrix of each described solid element；All surface cells are cycled, it will
The corresponding face load of each described surface cell is converted into equivalent node load；By the stiffness matrix and the equivalent knot
Point load carries out general assembly, obtains overall algebraic equation.
2. the face load processing method according to claim 1 based on surface cell, which is characterized in that described in discretization
The effect that processing obtains has establishes surface cell on the solid element of face load, including：
Determine that effect that sliding-model control obtains has the preset kind of solid element of face load；
The solid element is decomposed according to the preset kind, obtains preset model solid element；
According to preset order, the corresponding surface cell of each described preset model solid element is established.
3. the face load processing method according to claim 2 based on surface cell, which is characterized in that the basis is preset
Rule generates the incidence relation of the surface cell and the face load, including：
According to preset rules, the incidence relation of the surface cell and the preset model solid element is generated；
According to the surface cell and the incidence relation of the preset model solid element, the surface cell and face load are generated
Incidence relation.
4. the face load processing method according to claim 3 based on surface cell, which is characterized in that the Xun Huan is all
The surface cell includes：
According to the surface cell and the incidence relation of the solid element and the predetermined number of the solid element, to described
Surface cell is numbered；
According to the surface cell number that number obtains, all surface cells are cycled.
5. the face load processing method according to claim 4 based on surface cell, which is characterized in that described by each
The corresponding face load of the surface cell, which is converted into equivalent node load, to be included：
According to formulaIt will act at the force per unit area on surface cellIt is converted into equivalent node
Power fe；
Wherein, feRepresent the corresponding equivalent node force of equivalent node load；Expression acts on the distribution of load above surface cell
Power, and Represent three components of the force per unit area under rectangular coordinate system；Represent surface
The area of space of unit, the area of space include any one in curved surface and plane；N represents surface cell shape function square
Battle array；N1, N2, N3..., NmFor the shape function of each node of correspondence, 4≤m≤8, and
6. the face load processing method according to claim 1 based on surface cell, which is characterized in that described in discretization
The effect that processing obtains, which has, establishes on the solid element of face load before surface cell, further includes：
According to the physical property and geometric areas of test object, the geometrical model of the test object is established；
Sliding-model control is carried out to the geometrical model, obtains multiple solid elements for acting on and having face load.
7. a kind of face load-handling device based on surface cell, which is characterized in that including：
First establishes unit, and surface cell is established on the solid element that the effect obtained in sliding-model control has face load；
First generation unit, for according to preset rules, the surface cell of unit foundation and the face to be established in generation described first
The incidence relation of load；
First determination unit, for associating according to the surface cell that first generation unit generates and the face load
Relation determines the preset model of the surface cell；
Second determination unit for the preset model determined according to first determination unit, determines the knot of the surface cell
Points；
Second generation unit, the preset model and described second of the surface cell for being determined according to the first determination unit determine list
The definite nodal point number of member, generates pre-processing data；
First solves unit, for according to default constraints, to the preprocessing number of second generation unit generation
According to being solved, overall algebraic equation is obtained；The first solution unit includes：Computation subunit, for cycling all entities
Unit calculates the corresponding stiffness matrix of each solid element；Conversion subunit, will be each for cycling all surface unit
The corresponding face load of a surface cell is converted into equivalent node load；General assembly subelement, for calculate computation subunit
Stiffness matrix and the equivalent node load of conversion subunit conversion carry out general assembly, obtain overall algebraic equation；
Second solves unit, for being solved to the overall algebraic equation that the described first solution unit obtains, is tied
Point displacement, subsequently to calculate the corresponding stress-strain field of the solid element according to the displacement of joint.
8. the face load-handling device according to claim 7 based on surface cell, which is characterized in that described in discretization
The effect that processing obtains, which has, establishes surface cell on the solid element of face load, described first, which establishes unit, includes：
Determination subelement determines that effect that sliding-model control obtains has the preset kind of solid element of face load；
Subelement is decomposed, the solid element is decomposed, obtained by the preset kind determined according to the determination subelement
To preset model solid element；
Subelement is established, according to preset order, establishes each described described preset model entity list for decomposing subelement and obtaining
The corresponding surface cell of member.
9. the face load-handling device according to claim 8 based on surface cell, which is characterized in that first generation
Unit includes：
First generation subelement, for according to preset rules, generating the surface cell and the preset model solid element
Incidence relation；
Second generation subelement, for the surface cell generated according to the described first generation subelement and the preset model
The incidence relation of solid element generates the incidence relation of the surface cell and face load.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510191027.8A CN104778325B (en) | 2015-04-21 | 2015-04-21 | Face load processing method and processing device based on surface cell |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510191027.8A CN104778325B (en) | 2015-04-21 | 2015-04-21 | Face load processing method and processing device based on surface cell |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104778325A CN104778325A (en) | 2015-07-15 |
| CN104778325B true CN104778325B (en) | 2018-05-25 |
Family
ID=53619786
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510191027.8A Active CN104778325B (en) | 2015-04-21 | 2015-04-21 | Face load processing method and processing device based on surface cell |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104778325B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105825035B (en) * | 2016-05-16 | 2019-03-22 | 中国航空工业集团公司西安飞机设计研究所 | The equivalent way of surface distributed force when a kind of wing strut axial compression |
| CN105868511B (en) * | 2016-05-18 | 2019-03-08 | 苏州数设科技有限公司 | A kind of method and apparatus of determining aircraft reinforcement plate unit and load |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101561840A (en) * | 2009-06-02 | 2009-10-21 | 肖锋 | Numerical value simulation method for spot welding connection and failure thereof |
| US7844421B2 (en) * | 2007-01-15 | 2010-11-30 | Seiko Epson Corporation | Hybrid finite element method for traveling surface acoustic waves with thickness effect |
-
2015
- 2015-04-21 CN CN201510191027.8A patent/CN104778325B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7844421B2 (en) * | 2007-01-15 | 2010-11-30 | Seiko Epson Corporation | Hybrid finite element method for traveling surface acoustic waves with thickness effect |
| CN101561840A (en) * | 2009-06-02 | 2009-10-21 | 肖锋 | Numerical value simulation method for spot welding connection and failure thereof |
Non-Patent Citations (1)
| Title |
|---|
| 动车组铝合金车体结构强度分析;张志华;《中国优秀硕士学位论文全文数据库-工程科技Ⅱ辑》;20080915(第09期);第5-12、60-62页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104778325A (en) | 2015-07-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Roget et al. | Robust and efficient overset grid assembly for partitioned unstructured meshes | |
| CN108090963A (en) | A kind of numerical computation method of fibre reinforced composites thermal residual strain at low temperature | |
| Chen et al. | Performance of the MLPG method for static shakedown analysis for bounded kinematic hardening structures | |
| US11119466B2 (en) | Implicit method and an algorithm for flexible functionally tailorable slicing for additive manufacturing | |
| Liu et al. | A programming framework for large scale numerical simulations on unstructured mesh | |
| Knap et al. | A computational framework for scale‐bridging in multi‐scale simulations | |
| Kumar et al. | Synthesis of c0 path-generating contact-aided compliant mechanisms using the material mask overlay method | |
| Rama et al. | Real-time modelling of diastolic filling of the heart using the proper orthogonal decomposition with interpolation | |
| Avery et al. | A computationally tractable framework for nonlinear dynamic multiscale modeling of membrane woven fabrics | |
| JP7221062B2 (en) | Fluid analysis system, fluid analysis method, and fluid analysis program | |
| Fadel et al. | Packing optimization of free-form objects in engineering design | |
| Mesri et al. | Dynamic parallel adaption for three dimensional unstructured meshes: Application to interface tracking | |
| CN104778325B (en) | Face load processing method and processing device based on surface cell | |
| Duh et al. | Discretization of non-uniform rational B-spline (NURBS) models for meshless isogeometric analysis | |
| Fernandes et al. | ALE incompressible fluid–shell coupling based on a higher-order auxiliary mesh and positional shell finite element | |
| de la Fuente | An efficient procedure to obtain exact solutions in random vibration analysis of linear structures | |
| Khan et al. | Case based reasoning support for adaptive finite element analysis: mesh selection for an integrated system | |
| CN104750946B (en) | The processing method and processing device of face load in Finite Element | |
| Sapre et al. | Finite element mesh smoothing using cohort intelligence | |
| Sanfui et al. | Symbolic and numeric kernel division for graphics processing unit-based finite element analysis assembly of regular meshes with modified sparse storage formats | |
| Carrington et al. | Three-dimensional ALE-FEM method for fluid flow in domains with moving boundaries part 1: algorithm description | |
| CN114297877A (en) | Multi-working-condition simulation automation system and method for rod structure metamaterial structure | |
| Kumar | Object-Oriented Finite Element Programming for Engineering Analysis in C++. | |
| Jamroz et al. | Asynchronous communication in spectral-element and discontinuous Galerkin methods for atmospheric dynamics–a case study using the High-Order Methods Modeling Environment (HOMME-homme_dg_branch) | |
| Do et al. | Transient analysis of functionally graded plates using extreme gradient boosting |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |