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 PDF

Info

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
Application number
CN201510191027.8A
Other languages
Chinese (zh)
Other versions
CN104778325A (en
Inventor
陈志福
陈璞
任燕翔
李敏
段斌
屈凯峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Yun Dao Zhi Zao Science And Technology Ltd
Original Assignee
Beijing Yun Dao Zhi Zao Science And Technology Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Yun Dao Zhi Zao Science And Technology Ltd filed Critical Beijing Yun Dao Zhi Zao Science And Technology Ltd
Priority to CN201510191027.8A priority Critical patent/CN104778325B/en
Publication of CN104778325A publication Critical patent/CN104778325A/en
Application granted granted Critical
Publication of CN104778325B publication Critical patent/CN104778325B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

Face load processing method and processing device based on surface cell
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.
CN201510191027.8A 2015-04-21 2015-04-21 Face load processing method and processing device based on surface cell Active CN104778325B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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