CN104778325A - Surface element-based surface load processing method and surface element-based surface load processing device - Google Patents
Surface element-based surface load processing method and surface element-based surface load processing device Download PDFInfo
- Publication number
- CN104778325A CN104778325A CN201510191027.8A CN201510191027A CN104778325A CN 104778325 A CN104778325 A CN 104778325A CN 201510191027 A CN201510191027 A CN 201510191027A CN 104778325 A CN104778325 A CN 104778325A
- Authority
- CN
- China
- Prior art keywords
- surface cell
- load
- solid element
- preset
- 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.)
- Granted
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a surface element-based surface load processing method and a surface element-based surface load processing device. The surface element-based surface load processing method includes the following steps: surface elements are created on a solid element on which surface load acts; according to a preset rule, the incidence relation between the surface elements and the surface load is generated; according to the incidence relation, preset models of the surface elements are determined; according to the preset models, the number of nodes is determined; according to the preset models and the number of the nodes, preprocessed data are generated; the preprocessed data are solved, so that a general algebraic equation is obtained; the general algebraic equation is solved to obtain node displacements, so that a corresponding stress and strain field can be calculated according to the node displacements. By utilizing the surface elements to process the load, the invention can bring the following benefits: (1) the processing of the surface load can be completed as long as all the surface elements are circulated, so the process of processing the surface load is simplified; (2) all different types of solid elements can share one or a small number of surface elements to carry out the processing of boundary conditions, the multiplexing of surface element function codes is realized, and the maintenance and development costs of a calculation program are reduced.
Description
Technical field
The present invention relates to face load process field, in particular to based on the face load disposal route of surface cell and device.
Background technology
Finite element simulation is widely used in industrial design field, as machinery, the energy, Aero-Space etc.Dissimilar unit (Elements) is the core component that finite element program calculates, it is used for physical characteristics in description short space region and geometrical property, and is the coefficient calculating the algebraic equation needed by dissimilar load abbreviation.Unit is the basis realizing the physical field calculating such as stress-strain field, electromagnetic field, temperature field.
Calculating for three-D space structure realizes mainly through 3D solid unit (as tetrahedron element, hexahedral element, prism elements etc.) and Shell Finite Element.The face load of solid element can be pressure, temperature, electromotive force and current density etc.Need the node load these face load being converted into equivalence in computation process, finally all the node load of equivalence is converted into the right-hand vector of overall algebraic equation.
In conventional finite element simulation calculation flow process, be placed on unit inside to the process of elemental area load to process, finite-element preprocessing program needs given side load to act on which face of certain unit, in computation process, unit needs the shape according to different face load type and different faces, realizes different computing functions.Because face load may have Various Complex situation, the process of face load is part comparatively loaded down with trivial details during Elementary Function realizes.Wherein, the complex situations of face load comprise: according to the difference of face loading area, and face load can be divided into loads in plane and face external applied load.Different according to the space distribution of face magnitude of load, uniformly distributed load and non-uniformly distributed load can be divided into again.And if it is not the shape difference of cell surface (such as triangle, quadrilateral), identical to the computing method of solid element yet.
Inventor finds under study for action, carries out differentiation process bring many troubles to unit calculating in prior art to dissimilar face load in a unit.One is store the face sequence number that unit number that load information should store load effect stores load effect again, needs the geological information required for these information searchings in computation process, calculate and search procedure cumbersome, easily make mistakes.On the other hand, the face load handling procedure of the unit (as tetrahedron element, hexahedral element, prism elements etc.) of different geometries has the part of many repetitions but can not be multiplexing, increases the difficulty that program code is safeguarded.
Summary of the invention
The object of the present invention is to provide the face load disposal route based on surface cell and device, the process that opposite load carries out processing can be simplified, and code is multiplexing, the very big Maintenance Development cost reducing calculation procedure.
First aspect, embodiments provides a kind of face load disposal route based on surface cell, comprising:
The effect obtained at sliding-model control has on the solid element of face load sets up surface cell;
According to preset rules, generate the incidence relation of surface cell and face load;
According to the incidence relation of surface cell and face load, determine the preset model of surface cell;
According to preset model, determine the nodal point number of node surface cell;
According to preset model and the nodal point number of surface cell, generate pre-processing data;
According to the constraint condition preset, pre-processing data is solved, obtains overall algebraic equation;
Overall algebraic equation is solved, obtains displacement of joint, so that the follow-up stress-strain field corresponding according to displacement of joint computational entity unit.
In conjunction with first aspect, embodiments provide the first possible embodiment of first aspect, wherein, the effect obtained at sliding-model control has on the solid element of face load sets up surface cell, comprising:
Determine that effect that sliding-model control obtains has the preset kind of the solid element of face load;
According to preset kind, solid element is decomposed, obtain preset model solid element;
According to preset order, set up the surface cell that each preset model solid element is corresponding.
In conjunction with the first possible embodiment of first aspect, embodiments provide the embodiment that the second of first aspect is possible, wherein, according to preset rules, generate the incidence relation of surface cell and face load, comprising:
According to preset rules, generate the incidence relation of surface cell and preset model solid element;
According to the incidence relation of surface cell and preset model solid element, generate the incidence relation of surface cell and face load.
In conjunction with the embodiment that the second of first aspect is possible, embodiments provide the third possible embodiment of first aspect, wherein, according to the constraint condition preset, pre-processing data solved, obtains overall algebraic equation and comprise:
Circulate all solid elements, calculates the stiffness matrix that each solid element is corresponding;
Circulation all surface unit, is all converted to equivalent node load by face load corresponding for each surface cell;
Stiffness matrix and equivalent node load are carried out general assembly, obtains overall algebraic equation.
In conjunction with the third possible embodiment of first aspect, embodiments provide the 4th kind of possible embodiment of first aspect, wherein, circulation all surface unit comprises:
According to surface cell and the incidence relation of solid element and the predetermined number of solid element, effects on surface element number;
According to numbering the surface cell number obtained, circulation all surface unit.
In conjunction with the 4th kind of possible embodiment of first aspect, embodiments provide the 5th kind of possible embodiment of first aspect, wherein, face load corresponding for each surface cell be all converted to equivalent node load and comprise:
According to formula
force per unit area on surface cell will be acted on
be converted into equivalent node force f
e;
Wherein, f
erepresent the equivalent node force that equivalent node load is corresponding;
represent the distributed force acting on load above surface cell, and
three components of presentation surface distributed force under rectangular coordinate system;
the area of space of presentation surface unit, area of space comprise in curved surface and plane any one; N presentation surface unit form function matrix; N
1, N
2, N
3..., N
mfor the shape function of each node corresponding, 4≤m≤8, and
In conjunction with first aspect, embodiments provide the 6th kind of possible embodiment of first aspect, wherein, before the effect obtained at sliding-model control has the solid element of face load sets up surface cell, also comprise:
According to physical property and the geometric areas of tested object, set up the geometric model of tested object;
Sliding-model control is carried out to geometric model, obtains the solid element that multiple effect has face load.
Second aspect, the embodiment of the present invention additionally provides a kind of face load-handling device based on surface cell, comprising:
First sets up unit, and the effect obtained at sliding-model control has on the solid element of face load sets up surface cell;
First generation unit, for according to preset rules, generates first and sets up the surface cell of unit foundation and the incidence relation of face load;
First determining unit, for according to the surface cell of the first generation unit generation and the incidence relation of face load, determines the preset model of surface cell;
Second determining unit, for the preset model determined according to the first determining unit, determines the nodal point number of node surface cell;
Second generation unit, for the nodal point number that the preset model of surface cell determined according to the first determining unit and the second determining unit are determined, generates pre-processing data;
First solves unit, for according to the constraint condition preset, solves, obtain overall algebraic equation to the pre-processing data that the second generation unit generates;
Second solves unit, solving, obtaining displacement of joint, so that the follow-up stress-strain field corresponding according to displacement of joint computational entity unit for solving to first the overall algebraic equation that unit obtains.
In conjunction with second aspect, embodiments provide the first possible embodiment of second aspect, wherein, first sets up unit comprises:
Determine subelement, determine that effect that sliding-model control obtains has the preset kind of the solid element of face load;
Decomposing subelement, according to determining that solid element decomposes by the preset kind that subelement is determined, obtaining preset model solid element;
Set up subelement, according to preset order, set up each and decompose surface cell corresponding to the preset model solid element that obtains of subelement.
In conjunction with the first possible embodiment of second aspect, embodiments provide the embodiment that the second of second aspect is possible, wherein, the first generation unit comprises:
First generates subelement, for according to preset rules, generates the incidence relation of surface cell and preset model solid element;
Second generates subelement, for generating the surface cell of subelement generation and the incidence relation of preset model solid element according to first, generates the incidence relation of surface cell and face load.
The face load disposal route based on surface cell that the embodiment of the present invention provides and device, adopt the effect obtained at sliding-model control to have on the solid element of face load and set up surface cell, according to preset rules, generate the incidence relation of surface cell and face load, according to the incidence relation of surface cell and face load, determine the preset model of surface cell, according to preset model, determine the nodal point number of node surface cell, according to preset model and the nodal point number of surface cell, generate pre-processing data, according to the constraint condition preset, pre-processing data is solved, obtains overall algebraic equation, overall algebraic equation is solved, obtain displacement of joint, so that the follow-up stress-strain field corresponding according to displacement of joint computational entity unit, with of the prior art, dissimilar face load to be calculated to unit and brings compared with many troubles carrying out differentiation process in a unit, it utilizes surface cell to carry out load process, utilize surface cell to carry out load process and can bring following benefit: first, as long as store the surface cell number of load effect when storing load information, the all surface unit that circulates in computation process can complete the process of opposite load, simplify the process that opposite load carries out processing, the second, all dissimilar solid elements can share the process that a kind of or a few surface cell can complete most boundary condition, realize the multiplexing of surface cell function code most possibly, reduce the Maintenance Development cost of calculation procedure.
For making above-mentioned purpose of the present invention, feature and advantage become apparent, preferred embodiment cited below particularly, and coordinate appended accompanying drawing, be described in detail below.
Accompanying drawing explanation
In order to be illustrated more clearly in the technical scheme of the embodiment of the present invention, be briefly described to the accompanying drawing used required in embodiment below, be to be understood that, the following drawings illustrate only some embodiment of the present invention, therefore the restriction to scope should be counted as, for those of ordinary skill in the art, under the prerequisite not paying creative work, other relevant accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 shows the process flow diagram of a kind of face load disposal route based on surface cell that the embodiment of the present invention provides;
Fig. 2 shows another kind that the embodiment of the present invention the provides process flow diagram based on the face load disposal route of surface cell;
Fig. 3 shows another kind that the embodiment of the present invention the provides process flow diagram based on the face load disposal route of surface cell;
Fig. 4 shows another kind that the embodiment of the present invention the provides process flow diagram based on the face load disposal route of surface cell;
Fig. 5 shows another kind that the embodiment of the present invention the provides process flow diagram based on the face load disposal route of surface cell;
Fig. 6 shows the structural representation of a kind of face load-handling device based on surface cell that the embodiment of the present invention provides;
It is a kind of based on the first structural representation setting up unit in the face load-handling device of surface cell that Fig. 7 shows that the embodiment of the present invention provides;
Fig. 8 shows a kind of structural representation based on the first generation unit in the face load-handling device of surface cell that the embodiment of the present invention provides;
It is a kind of based on the first structural representation solving unit in the face load-handling device of surface cell that Fig. 9 shows that the embodiment of the present invention provides;
Figure 10 shows another kind that the embodiment of the present invention provides based on the first structural representation solving unit in the face load-handling device of surface cell;
Figure 11 shows another kind that the embodiment of the present invention provides based on the structural representation in the face load-handling device of surface cell;
It is a kind of based on the distribution schematic diagram of pressure at node place in the face load disposal route of surface cell that Figure 12 shows that the embodiment of the present invention provides.
Embodiment
Below in conjunction with accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.The assembly of the embodiment of the present invention describing and illustrate in usual accompanying drawing herein can be arranged with various different configuration and design.Therefore, below to the detailed description of the embodiments of the invention provided in the accompanying drawings and the claimed scope of the present invention of not intended to be limiting, but selected embodiment of the present invention is only represented.Based on embodiments of the invention, the every other embodiment that those skilled in the art obtain under the prerequisite not making creative work, all belongs to the scope of protection of the invention.
Finite element analysis (FEA, Finite Element Analysis) utilizes the method for mathematical approach to simulate actual physical system (geometry and load working condition).Also utilize simple and interactional element, i.e. unit, just can remove with the unknown quantity of limited quantity the real system approaching unlimited unknown quantity.At present, finite element simulation is widely used in industrial design field, as machinery, the energy, Aero-Space etc.Dissimilar unit (Elements) is the core component (solving the discretize in territory) that finite element program calculates, it is used for physical characteristics in description short space region and geometrical property, and is the coefficient calculating the algebraic equation needed by dissimilar load abbreviation.Unit is the basis realizing the physical field calculating such as stress-strain field, electromagnetic field, temperature field.
In conventional finite element simulation calculation flow process, unit inside is placed on to the process of elemental area load and processes.Finite-element preprocessing program needs given side load to act on which face of certain unit (namely will know the unit inner face sequence number of the solid element numbering that load acts on and load institute acting surface) simultaneously; In computation process, unit needs according to different face load type (i.e. mesh shape) and different face shape (in this way plane or curved surface), realizes different computing functions.Because face load may have Various Complex situation, the process of face load is part comparatively loaded down with trivial details during Elementary Function realizes.Wherein, the complex situations of face load comprise: according to the difference of face loading area, and face load can be divided into loads in plane and face external applied load.Different according to the space distribution of face magnitude of load, uniformly distributed load and non-uniformly distributed load can be divided into again.And if the shape difference of cell surface (such as triangle, quadrilateral) computing method are not identical yet.
Concrete, differentiation process is carried out in a unit to dissimilar face load and brings many troubles to unit calculating.One is store the face sequence number that unit number that load information should store load effect stores load effect again, needs the geological information required for these information searchings, and carry out computing according to these geological information opposite load in computation process; Wherein, the calculating in said process and search procedure cumbersome, easily make mistakes.On the other hand, the unit of different geometries is (as tetrahedron element, hexahedral element and prism elements etc.) face load handling procedure have the part of many repetitions but can not be multiplexing, as tetrahedron and hexahedron are all made up of multiple quadrilateral, identical quadrangular portion is included in its face load handling procedure, and calling of same section program cannot be realized in unit, again can only write identical code, increase the difficulty that staff's program code is safeguarded.
The present invention proposes the method unifying treated side load with the layer of surface unit being attached to 3D solid surface, and the process of face load is no longer placed on unit inside and carries out.Wherein, so-called surface cell refers to that thickness is zero, and be used for face load to be converted into a class unit of equivalent node load, its shape can be triangle, quadrilateral etc., further, the geometric properties that surface cell is concrete is determined by the surface characteristics of accompanying solid element.In order to process complicated boundary condition in computation process, coordinate solid element, simplify load treatment scheme.
The disposal route of face load and device in Finite Element provided by the invention, benefit that load process brings has to utilize surface cell to carry out: first, as long as store the surface cell number of load effect when storing load information, the all surface unit that circulates in computation process can complete the process of opposite load, the geological information searching needs need not be entered in solid element again, simplify the process that opposite load carries out processing; The second, all dissimilar solid elements can share the process that a kind of or a few surface cell can complete most boundary condition, realize the multiplexing of surface cell function code most possibly, reduce the Maintenance Development cost of calculation procedure.
Specific embodiments of the invention are as follows:
With reference to figure 1, a kind of face load disposal route based on surface cell, comprising:
101, the effect obtained at sliding-model control has on the solid element of face load sets up surface cell.
Concrete, in finite element analysis, first pre-treatment modeling is carried out to tested object, and sliding-model control is carried out to the geometric model of tested object, obtain the solid element acting on face load; Wherein, this solid element is multiple, obtains as which solid element to each solid element numbering; In the present embodiment, solid element is preferably polygon, when solid element is polygon (as triangle and quadrilateral), each face of this solid element is all numbered, obtain solid element face sequence number, then determine the numbering of solid element of face load effect and the face sequence number of the solid element of this numbering, obtain the incidence relation of solid element and face load.
Then the solid element obtained at above-mentioned sliding-model control sets up surface cell, so that follow-up calculating and the process carrying out face load according to this surface cell.
Wherein, so-called surface cell refers to that thickness is zero, and be used for face load to be converted into a class unit of equivalent node load, its shape can be triangle, quadrilateral etc., further, the geometric properties that surface cell is concrete is determined by the surface characteristics of accompanying solid element.
102, according to preset rules, the incidence relation of surface cell and face load is generated.
Concrete, preset rules can preset, and such as makes some surface cells replace certain face of certain solid element, and makes another solid element replace another face etc. of this solid element.This preset rules specifically can be arranged as required, and the present invention does not do concrete restriction to this.According to preset rules, the incidence relation of the surface cell that record is set up and solid element, and the incidence relation of surface cell and face load is generated according to the incidence relation of this surface cell and solid element, said process acts on which number face of certain solid element by face load, be converted to and face load is acted on surface cell, and replace which number face of certain solid element with this surface cell.
103, according to the incidence relation of surface cell and face load, the preset model of surface cell is determined.
Concrete, the preset model of surface cell is which number face for replacing which solid element, therefore this preset model is under the effect of himself, and its shape can be arranged arbitrarily, as quadrilateral, triangle etc.
104, according to preset model, the nodal point number of node surface cell is determined.
Concrete, determining under above-mentioned preset model, be quadrangular surface unit with surface cell in the present embodiment for example is described, as known that surface cell is quadrilateral, and it is No. 1 face for representing No. 1 solid element, then can according to above-mentioned preset model, whether No. 1 face determining No. 1 solid element there is face load, if have, determining that this face load acts on the nodal point number on quadrangular surface unit.
Wherein, this nodal point number is preferably greater than or equal to polygonal basic angle point number, and as when surface cell is quadrangular surface unit, nodal point number is more than or equal to 4; When surface cell is triangular surface unit, the node that nodal point number is more than or equal to 3 surface cells depends on the nodal point number on solid element surface.For single order solid element, Surface Nodes number equals surperficial polygonal limit number, and for high-order solid element, Surface Nodes number is greater than surperficial polygonal limit number.The nodal point number of surface cell is corresponding.
105, according to preset model and the nodal point number of surface cell, pre-processing data is generated.
Concrete, the site position on the number of pre-processing data and surface cell, each surface cell on node number, each surface cell and face load quantity etc. corresponding to each surface cell.
In the present embodiment, be which number face that instead of certain solid element with surface cell, therefore according to the incidence relation of this surface cell and face load, namely can generate above-mentioned etc. pre-processing data, carry out solving calculating according to this pre-processing data so that follow-up.
Wherein, the preset model of surface cell determines according to the preset kind of solid element, and namely the surface cell in the present embodiment represents which number face of certain solid element, then the type in which number face of this solid element is exactly the preset kind of surface cell.
106, according to the constraint condition preset, pre-processing data is solved, obtains overall algebraic equation.
Concrete, namely the constraint condition preset process degree of freedom restrained in geometric model; In fact, solving in computation process in reality, is need some variablees to retrain, and using these variablees as constant coefficient, namely replaces complicated problem with simple question, is solving simple question, finally could solve complicated problem; Namely this process is be defined as constant coefficient to the design variables in geometric model, and according to fixed constant coefficient, pre-processing data is solved, and the result solved is carried out general assembly, then overall algebraic equation is obtained, follow-up this overall algebraic equation to be solved, required result can be obtained---i.e. displacement of joint.
107, overall algebraic equation is solved, obtain displacement of joint, so that the follow-up stress-strain field corresponding according to displacement of joint computational entity unit.
Concrete, namely this process is that the overall algebraic equation of assembling in step 104 solves, and namely the overall algebraic equation after assembling is the equation only containing displacement of joint variable, can calculate displacement of joint according to mathematical operation.
Then, according to the displacement of joint of trying to achieve and in conjunction with the general-purpose algorithm in industry, can calculate stress-strain field corresponding for the solid element of correspondence, so that staff can according to the tensile strength of this stress-strain field determination tested object.Wherein, the distribution situation of this stress-strain field and stress value and all stress values.
The disposal route of face load in Finite Element provided by the invention, benefit that load process brings has to utilize surface cell to carry out: first, as long as store the surface cell number of load effect when storing load information, the all surface unit that circulates in computation process can complete the process of opposite load, the geological information searching needs need not be entered in solid element again, simplify the process that opposite load carries out processing; The second, all dissimilar solid elements can share the process that a kind of or a few surface cell can complete most boundary condition, realize the multiplexing of surface cell function code most possibly, reduce the Maintenance Development cost of calculation procedure.
With reference to figure 2, in the present invention, consider follow-up convenience of calculation and the problem that can be used to repetition program code, in this Finite Element face load disposal route in, the effect obtained at sliding-model control in step 101 has on the solid element of face load sets up surface cell, specifically comprises the steps:
201, determine that effect that sliding-model control obtains has the preset kind of the solid element of face load.
Concrete, the preset kind of solid element comprises: 3D solid unit is (as tetrahedron element, hexahedral element, prism elements etc.) and Shell Finite Element, first determine the particular type of this solid element, on this solid element, set up surface cell according to this type determined so that follow-up.
202, according to preset kind, solid element is decomposed, obtain preset model solid element.
Concrete, after determining the type of solid element, the different decomposition method corresponding according to all kinds decomposes whole solid element; As to quadrilateral, namely according to six of quadrilateral, the face load that whether mesh shape, every one side of the face load of every one side are planes, this face is corresponding to be face load corresponding to loads in plane or face external applied load and this face be even load also non-uniform load is decomposed solid element; As diabolo, namely according to leg-of-mutton four sides, the face load that whether mesh shape, every one side of the face load of every one side are planes, this face is corresponding to be face load corresponding to loads in plane or face external applied load and this face be even load also non-uniform load is decomposed solid element.
Wherein, all corresponding preset model solid element of each model that above-mentioned decomposition obtains, writes down the quantity of this preset model solid element, so that the surface cell that follow-up foundation is identical with the preset model solid element of this quantity.
Certainly, concrete isolation is not limited to said method, other be in solid element carry out in decomposition range all in protection scope of the present invention.
203, according to preset order, the surface cell that each preset model solid element is corresponding is set up.
Concrete, after decomposition terminates, obtain multiple preset model, set up each preset model solid element to the surface cell of correspondence.Wherein, the quantity of surface cell is identical with the quantity of preset model solid element, and the quantity of surface cell can be not identical with preset model solid element yet.When not identical, acting on the preset model solid element of face load as determined in preset model solid element, then having on the preset model solid element of face load in effect and setting up surface cell.
With reference to figure 3, in the present invention, the specific implementation of step 103 is as follows:
301, according to preset rules, the incidence relation of surface cell and preset model solid element is generated.
302, according to the incidence relation of surface cell and preset model solid element, the incidence relation of surface cell and face load is generated.
Concrete, incidence relation first according to preset rules determination surface cell and solid element in the embodiment of the present invention, because solid element is relevant with face load known in advance, therefore the incidence relation of surface cell and face load can be obtained by both incidence relations above-mentioned.
The processing procedure of opposite of the present invention load is by surface cell agent entity unit, by surface cell process complicated boundary condition in computation process, and makes it coordinate solid element, simplify load treatment scheme, with reference to figure 4, wherein, use the concrete processing mode of face load as follows:
401, circulate all solid elements, calculates the stiffness matrix that each solid element is corresponding.
In the present embodiment, the mode of calculated rigidity matrix is identical with the mode of existing computational entity unit, and this step be not described.
402, circulation all surface unit, is all converted to equivalent node load by face load corresponding for each surface cell.
Concrete, for the ease of circulation all surface unit, first according to surface cell and the incidence relation of solid element and the predetermined number of solid element, effects on surface unit is numbered, and obtains surface cell number; Then according to numbering the surface cell number obtained, circulation all surface unit, and be equivalent node load by face load transfer corresponding for the surface cell determined.
Circulate in the present embodiment surface cell, is equivalent node load according to node coordinate on surface cell by certain algorithm by face load transfer.Compare the solid element numbering and the face sequence number of acting surface of searching and having face load effect, judge the shape in face according to face sequence number and search the geometric coordinate of node on face.Being the traditional process of equivalent node load according to face shape and node coordinate by calculating face load transfer, eliminating the trouble of the face of searching and node information.
403, stiffness matrix and equivalent node load are carried out general assembly, obtain overall algebraic equation.
Concrete, by the left end item of the stiffness matrix calculated each solid element as algebraic equation, using the replacer of above-mentioned solid element---multiple equivalent node information that the face LOAD FOR that each surface cell is corresponding obtains are as the right-hand vector of algebraic equation, obtain multiple left end item and right-hand vector, make that multiple left end item is corresponding with multiple right-hand vector always dresses up overall algebraic equation.
Be quadrangular surface unit with surface cell in the present embodiment be example pair, the concrete grammar that face load corresponding for each surface cell is all converted to equivalent node load be described, comprise:
According to formula
force per unit area on surface cell will be acted on
be converted into equivalent node force f
e;
Wherein, f
erepresent the equivalent node force that equivalent node load is corresponding;
represent the distributed force acting on load above surface cell, and
three components of presentation surface distributed force under rectangular coordinate system;
the area of space of presentation surface unit, area of space comprise in curved surface and plane any one; N presentation surface unit form function matrix; N
1, N
2, N
3..., N
mfor the shape function of each node corresponding, 4≤m≤8, and
In the present embodiment, surface cell is described for quadrangular surface unit, and the present invention proposes with the quadrangular surface unit being attached to 3D solid surface, is used for processing the face load diatibution power in Structure Calculation, as pressure.The process of face load is no longer placed on 3D solid unit inside and carries out.The thickness of surface cell is zero, and it is only used for face load to be converted into equivalent node load.
In addition, with reference to figure 5, before the effect obtained at sliding-model control has the solid element of face load sets up surface cell, also comprise:
501, according to physical property and the geometric areas of tested object, the geometric model of tested object is set up.
502, sliding-model control is carried out to geometric model, obtain the solid element that multiple effect has face load.
Below for quadrangular surface unit, brief description is carried out to the process that this quadrangular surface unit matching solid element opposite load processes, with reference to Figure 12:
The shape of quadrangular surface unit can be quadrilateral or space curved surface quadrilateral in plane; The nodal point number of quadrangular surface unit can be 4,5,6,7,8, and wherein 1-4 node is limit end points, and No. 5-8 is node in limit, as shown in figure 12 (wherein, Q7 and Q is not shown in Figure 12).
Figure 12. quadrangular surface unit, node Q represents, the size and Orientation of distribution pressure at node place represents with arrow shown in P.
Quadrangular surface unit is according to formula
force per unit area on surface cell will be acted on
be converted into equivalent node force f
e.
Wherein,
Each variable implication is as follows:
---the area of space (curved surface or plane) of surface cell;
N---surface cell form function matrix; N
1, N
2, N
3..., N
mfor the shape function of each node corresponding, 4≤m≤8 and
And
---three components of force per unit area under rectangular coordinate system.
Below for linear static calculation process, illustrate how quadrangular surface unit coordinates solid element to complete its computing function.As Figure 1-Figure 5.Based in the calculation process of surface cell, on the corresponding quadrangular surface unit of pre-treatment program definition face load effect.First circulate when solving solid element, complete Rigidity Calculation and the stiffness matrix assembling of solid element, then circulation quadrangular surface unit, the distributed load acted on quadrilateral curved surface is treated to corresponding equivalent node load, equivalent node load joins in integrated solution equation by solver, solving equation.Obtain the displacement of unit node according to integrated solution non trivial solution, finally calculate the stress-strain field be concerned about according to unit displacement of joint.
The face load disposal route based on surface cell that the embodiment of the present invention provides, with of the prior art, dissimilar face load to be calculated to unit and brings compared with many troubles carrying out differentiation process in a unit, it utilizes surface cell to carry out load process, utilize surface cell to carry out load process and can bring following benefit: first, as long as store the surface cell number of load effect when storing load information, the all surface unit that circulates in computation process can complete the process of opposite load, simplifies the process that opposite load carries out processing; The second, all dissimilar solid elements can share the process that a kind of or a few surface cell can complete most boundary condition, realize the multiplexing of surface cell function code most possibly, reduce the Maintenance Development cost of calculation procedure.
With reference to figure 6, the invention provides a kind of face load-handling device based on surface cell, comprising:
First sets up unit 11, and the effect obtained at sliding-model control has on the solid element of face load sets up surface cell;
First generation unit 12, for according to preset rules, generates first and sets up the surface cell of unit 11 foundation and the incidence relation of face load;
First determining unit 13, for according to the surface cell of the first generation unit 12 generation and the incidence relation of face load, determines the preset model of surface cell;
Second determining unit 14, for the preset model determined according to the first determining unit 13, determines the nodal point number of node surface cell;
Second generation unit 15, for the nodal point number that the preset model of surface cell determined according to the first determining unit 13 and the second determining unit 14 are determined, generates pre-processing data;
First solves unit 16, for according to the constraint condition preset, solves, obtain overall algebraic equation to the pre-processing data that the second generation unit 15 generates;
Second solves unit 17, solving, obtaining displacement of joint, so that the follow-up stress-strain field corresponding according to displacement of joint computational entity unit for solving to first the overall algebraic equation that unit 16 obtains.
Further, with reference to figure 7, should based in the face load-handling device of surface cell, first sets up unit 11 comprises:
Determine subelement 111, for determining that effect that sliding-model control obtains has the preset kind of the solid element of face load.
Decomposing subelement 112, for according to determining that solid element decomposes by the preset kind that subelement 111 is determined, obtaining preset model solid element;
Set up subelement 113, for according to preset order, set up each and decompose surface cell corresponding to preset model solid element that subelement 112 obtains.
Further, with reference to figure 8, in this Finite Element face load treating apparatus in, the first generation unit 12 comprises:
First generates subelement 121, for according to preset rules, generates the incidence relation of surface cell and preset model solid element;
Second generates subelement 122, for generating the surface cell of subelement 121 generation and the incidence relation of preset model solid element according to first, generates the incidence relation of surface cell and face load.
Further, with reference to figure 9, the treating apparatus of face load in this Finite Element, first solves unit 16 comprises:
Computation subunit 161, for all solid elements that circulates, calculates the stiffness matrix that each solid element is corresponding.
Conversion subelement 162, for all surface unit that circulates, is all converted to equivalent node load by face load corresponding for each surface cell;
General assembly subelement 163, the equivalent node load that stiffness matrix and conversion subelement 162 for computation subunit 161 being calculated are changed carries out general assembly, obtains overall algebraic equation.
Further, with reference to Figure 10, in this Finite Element face load treating apparatus in, first solves unit 16 also comprises:
Numbering subelement 164, for according to the incidence relation of surface cell and solid element and the predetermined number of solid element, effects on surface element number;
Circulation subelement 165, for numbering the surface cell number obtained according to numbering subelement 164, circulation all surface unit.
Further, in this Finite Element face load treating apparatus in, conversion subelement 162 specifically for, according to formula
force per unit area on surface cell will be acted on
be converted into equivalent node force f
e;
Wherein, f
erepresent the equivalent node force that equivalent node load is corresponding;
represent the distributed force acting on load above surface cell, and
three components of presentation surface distributed force under rectangular coordinate system;
the area of space of presentation surface unit, area of space comprise in curved surface and plane any one; N presentation surface unit form function matrix; N
1, N
2, N
3..., N
mfor the shape function of each node corresponding, 4≤m≤8, and
Further, with reference to Figure 11, in this Finite Element face load treating apparatus in also comprise:
Second sets up unit 18, for according to the physical property of tested object and geometric areas, sets up the geometric model of tested object;
Sliding-model control unit 19, carries out sliding-model control for the geometric model setting up unit 18 foundation to second, obtains the solid element that multiple effect has face load.
The face load-handling device based on surface cell that the embodiment of the present invention provides, adopts the effect obtained at sliding-model control to have on the solid element of face load and sets up surface cell, according to preset rules, generate the incidence relation of surface cell and face load, according to the incidence relation of surface cell and face load, determine the preset model of surface cell, according to preset model, determine the nodal point number of node surface cell, according to preset model and the nodal point number of surface cell, generate pre-processing data, according to the constraint condition preset, pre-processing data is solved, obtains overall algebraic equation, overall algebraic equation is solved, obtain displacement of joint, so that the follow-up stress-strain field corresponding according to displacement of joint computational entity unit, with of the prior art, dissimilar face load to be calculated to unit and brings compared with many troubles carrying out differentiation process in a unit, it utilizes surface cell to carry out load process, utilize surface cell to carry out load process and can bring following benefit: first, as long as store the surface cell number of load effect when storing load information, the all surface unit that circulates in computation process can complete the process of opposite load, simplify the process that opposite load carries out processing, the second, all dissimilar solid elements can share the process that a kind of or a few surface cell can complete most boundary condition, realize the multiplexing of surface cell function code most possibly, reduce 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, also can be that the independent physics of unit exists, also can two or more unit in a unit integrated.
If described function using the form of SFU software functional unit realize and as independently production marketing or use time, can be stored in a computer read/write memory medium.Based on such understanding, the part of the part that technical scheme of the present invention contributes to prior art in essence in other words or this technical scheme can embody with the form of software product, this computer software product is stored in a storage medium, comprising some instructions in order to make a computer equipment (can be personal computer, server, or the network equipment etc.) perform all or part of step of method described in each embodiment of the present invention.And aforesaid storage medium comprises: USB flash disk, portable hard drive, ROM (read-only memory) (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disc or CD etc. various can be program code stored medium.
The above; be only the specific embodiment of the present invention, but protection scope of the present invention is not limited thereto, is anyly familiar with those skilled in the art in the technical scope that the present invention discloses; change can be expected easily or replace, all should be encompassed within protection scope of the present invention.Therefore, protection scope of the present invention should described be as the criterion with the protection domain of claim.
Claims (10)
1., based on a face load disposal route for surface cell, it is characterized in that, comprising:
The effect obtained at sliding-model control has on the solid element of face load sets up surface cell;
According to preset rules, generate the incidence relation of described surface cell and described load;
According to the incidence relation of described surface cell and described load, determine the preset model of described surface cell;
According to described preset model, determine the nodal point number of described node surface cell;
According to preset model and the nodal point number of surface cell, generate pre-processing data;
According to the constraint condition preset, described pre-processing data is solved, obtains overall algebraic equation;
Described overall algebraic equation is solved, obtains displacement of joint, calculate stress-strain field corresponding to described solid element so that follow-up according to described displacement of joint.
2. the face load disposal route based on surface cell according to claim 1, is characterized in that, the described effect obtained at sliding-model control has on the solid element of face load sets up surface cell, comprising:
Determine that effect that sliding-model control obtains has the preset kind of the solid element of face load;
According to described preset kind, described solid element is decomposed, obtain preset model solid element;
According to preset order, set up the surface cell that described in each, preset model solid element is corresponding.
3. the face load disposal route based on surface cell according to claim 2, is characterized in that, described according to preset rules, generates the incidence relation of described surface cell and described load, comprising:
According to preset rules, generate the incidence relation of described surface cell and described preset model solid element;
According to the incidence relation of described surface cell and described preset model solid element, generate the incidence relation of described surface cell and face load.
4. the face load disposal route based on surface cell according to claim 3, is characterized in that, the constraint condition that described basis is preset, and solves, obtain overall algebraic equation and comprise described pre-processing data:
Circulate all described solid elements, calculates the stiffness matrix that solid element described in each is corresponding;
Circulate all described surface cells, and face load corresponding for surface cell described in each is all converted to equivalent node load;
Described stiffness matrix and described equivalent node load are carried out general assembly, obtains overall algebraic equation.
5. the face load disposal route based on surface cell according to claim 4, is characterized in that, all described surface cells of described circulation comprise:
According to described surface cell and the incidence relation of described solid element and the predetermined number of described solid element, described surface cell is numbered;
According to numbering the surface cell number obtained, circulate all described surface cells.
6. the face load disposal route based on surface cell according to claim 5, is characterized in that, describedly face load corresponding for surface cell described in each is all converted to equivalent node load comprises:
According to formula
force per unit area on surface cell will be acted on
be converted into equivalent node force f
e;
Wherein, f
erepresent the equivalent node force that equivalent node load is corresponding;
represent the distributed force acting on load above surface cell, and
three components of presentation surface distributed force under rectangular coordinate system;
the area of space of presentation surface unit, described area of space comprise in curved surface and plane any one; N presentation surface unit form function matrix; N
1, N
2, N
3..., N
mfor the shape function of each node corresponding, 4≤m≤8, and
7. the face load disposal route based on surface cell according to claim 1, is characterized in that, before the described effect obtained at sliding-model control has the solid element of face load sets up surface cell, also comprises:
According to physical property and the geometric areas of tested object, set up the geometric model of described tested object;
Sliding-model control is carried out to described geometric model, obtains the solid element that multiple effect has face load.
8., based on a face load-handling device for surface cell, it is characterized in that, comprising:
First sets up unit, and the effect obtained at sliding-model control has on the solid element of face load sets up surface cell;
First generation unit, for according to preset rules, generates described first and sets up the surface cell of unit foundation and the incidence relation of described load;
First determining unit, for according to the described described surface cell of the first generation unit generation and the incidence relation of described load, determines the preset model of described surface cell;
Second determining unit, for the preset model determined according to described first determining unit, determines the nodal point number of described node surface cell;
Second generation unit, for the described nodal point number that the preset model of surface cell determined according to the first determining unit and described second determining unit are determined, generates pre-processing data;
First solves unit, for according to the constraint condition preset, solves, obtain overall algebraic equation to the described pre-processing data that described second generation unit generates;
Second solves unit, solving, obtaining displacement of joint, calculating stress-strain field corresponding to described solid element so that follow-up according to described displacement of joint for solving to described first the described overall algebraic equation that unit obtains.
9. the face load-handling device based on surface cell according to claim 8, is characterized in that, the described effect obtained at sliding-model control has on the solid element of face load sets up surface cell, and described first sets up unit comprises:
Determine subelement, determine that effect that sliding-model control obtains has the preset kind of the solid element of face load;
Decompose subelement, according to the described described preset kind determining that subelement is determined, described solid element is decomposed, obtain preset model solid element;
Set up subelement, according to preset order, set up and described in each, to decompose surface cell corresponding to described preset model solid element that subelement obtains.
10. the face load-handling device based on surface cell according to claim 9, is characterized in that, described first generation unit comprises:
First generates subelement, for according to preset rules, generates the incidence relation of described surface cell and described preset model solid element;
Second generates subelement, for generating the described surface cell of subelement generation and the incidence relation of described preset model solid element according to described first, generates the incidence relation of described 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 true CN104778325A (en) | 2015-07-15 |
| CN104778325B 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) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105825035A (en) * | 2016-05-16 | 2016-08-03 | 中国航空工业集团公司西安飞机设计研究所 | Equivalent treating method for surface distribution force generated when wing supporting poles are axially pressed |
| CN105868511A (en) * | 2016-05-18 | 2016-08-17 | 苏州智设信息技术有限公司 | Method and device for determining airplane stiffened plate unit and loads |
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 |
|---|
| 张志华: "动车组铝合金车体结构强度分析", 《中国优秀硕士学位论文全文数据库-工程科技Ⅱ辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105825035A (en) * | 2016-05-16 | 2016-08-03 | 中国航空工业集团公司西安飞机设计研究所 | Equivalent treating method for surface distribution force generated when wing supporting poles are axially pressed |
| CN105825035B (en) * | 2016-05-16 | 2019-03-22 | 中国航空工业集团公司西安飞机设计研究所 | The equivalent way of surface distributed force when a kind of wing strut axial compression |
| CN105868511A (en) * | 2016-05-18 | 2016-08-17 | 苏州智设信息技术有限公司 | Method and device for determining airplane stiffened plate unit and loads |
| CN105868511B (en) * | 2016-05-18 | 2019-03-08 | 苏州数设科技有限公司 | A kind of method and apparatus of determining aircraft reinforcement plate unit and load |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104778325B (en) | 2018-05-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lipnikov et al. | A framework for developing a mimetic tensor artificial viscosity for Lagrangian hydrocodes on arbitrary polygonal meshes | |
| Murthy et al. | Nonparametric stochastic modeling of structural uncertainty in rotordynamics: unbalance and balancing aspects | |
| Larsson et al. | Fast and robust approximation of smallest enclosing balls in arbitrary dimensions | |
| CN104778325A (en) | Surface element-based surface load processing method and surface element-based surface load processing device | |
| Behera et al. | Fuzzy finite element based solution of uncertain static problems of structural mechanics | |
| Videla et al. | Explicit integration of the stiffness matrix of a four‐noded‐plane‐elasticity finite element | |
| Chen | Discrete element analysis methods of generic differential quadratures | |
| Rao | Simulation based engineering in solid mechanics | |
| Zhang et al. | Study on the stress-stiffening effect and modal synthesis methods for the dynamics of a spatial curved beam | |
| Do | Boundary control of slender beams under deterministic and stochastic loads | |
| Ryan et al. | Automating the derivation of the equations of motion of a multibody dynamic system with uncertainty using polynomial chaos theory and variational work | |
| Mirsajedi et al. | A multizone moving mesh algorithm for simulation of flow around a rigid body with arbitrary motion | |
| CN104750946A (en) | Method and device for treating surface load in finite element method | |
| CN104156606B (en) | The method for solving for nuclein migration governing equation based on mixing computing architecture | |
| Sapre et al. | Finite element mesh smoothing using cohort intelligence | |
| Kulikov et al. | Hydrodynamic modeling of self-gravitating astrophysical objects on tetrahedral meshes | |
| Oliveto et al. | 3D finite-deformation beam model with viscous damping: Computational aspects and applications | |
| Azevedo | Fracture mechanics using the natural neighbour radial point interpolation method | |
| Kolberg et al. | Parallel verified linear system solver for uncertain input data | |
| Veron et al. | Large displacements of light thin flexible structures coupled with heavy fluids using co-simulation between finite element and fast boundary element solvers | |
| US11775710B1 (en) | Methods and systems for simulating multistage cyclic symmetry assemblies | |
| Gaidzhurov et al. | Application of the Double Approximation Method for Constructing Stiffness Matrices of Volumetric Finite Elements | |
| Farina et al. | An inverse preconditioner for a free surface ocean circulation model | |
| Ganguli | Isospectral Vibrating Systems | |
| Berezyuk et al. | Test example of calculation of thin-walled structures with truncated conical shell |
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 |