CN107145651B - ABAQUS 3D Infinite Element Boundary Rapid Modeling Method Based on INP File - Google Patents

ABAQUS 3D Infinite Element Boundary Rapid Modeling Method Based on INP File Download PDF

Info

Publication number
CN107145651B
CN107145651B CN201710270976.4A CN201710270976A CN107145651B CN 107145651 B CN107145651 B CN 107145651B CN 201710270976 A CN201710270976 A CN 201710270976A CN 107145651 B CN107145651 B CN 107145651B
Authority
CN
China
Prior art keywords
infinite
inp
boundary
nodes
node
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.)
Expired - Fee Related
Application number
CN201710270976.4A
Other languages
Chinese (zh)
Other versions
CN107145651A (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN201710270976.4A priority Critical patent/CN107145651B/en
Publication of CN107145651A publication Critical patent/CN107145651A/en
Application granted granted Critical
Publication of CN107145651B publication Critical patent/CN107145651B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)

Abstract

The invention discloses an ABAQUS three-dimensional infinite element boundary rapid modeling method based on an INP file, which belongs to the field of infinite unit boundary modeling. The method is simple to operate, the ABAQUS infinite element boundary model can be generated in a short time, and the modeling rate of the ABAQUS is greatly improved.

Description

ABAQUS three-dimensional infinite element boundary rapid modeling method based on INP file
Technical Field
The invention relates to an infinite unit boundary modeling method, in particular to an ABAQUS three-dimensional infinite element boundary rapid modeling method based on an INP file.
Background
The advent and rapid development of computers have provided powerful tools for engineering analysis, making it possible to solve many large-scale engineering practical problems using numerical simulations. The finite element method is used as the most common numerical simulation means and is widely applied to the engineering fields of rock and soil, structures, earthquakes, water conservancy and the like. The problem of using finite element method is how to use finite model to simulate real infinite area, the common solution is to neglect the influence of infinite area boundary according to Saint-Venn principle, use truncation boundary to take large enough model size to divide finite element mesh of geometric size, and apply corresponding approximate constraint boundary condition on artificial boundary. In actual calculation, the size identification of the model which is large enough is difficult, the calculation cost is low when the area is small, but the result precision is poor; when the area is large, certain precision requirements are met, but the calculation cost is high. Especially for the wave computation problem, the finite element computation results are often distorted by reflections and scattering of the wave at the mesh boundaries.
The reasonable artificial boundary is established by introducing the boundary unit, so that the reflection of the finite element boundary on the boundary of various fluctuations in the foundation caused by artificial truncation can be reduced as much as possible, the number of units is greatly reduced, and the calculation efficiency is improved. Currently, artificial boundaries such as viscous boundaries, paraxial approximation boundaries, transmission boundaries, viscoelastic boundaries, infinite element boundaries, and the like are common.
As a mode for simulating an infinite area, the appearance of an infinite element provides an important way for overcoming the defects of a finite element calculation method and solving the problem of infinite area boundary simulation. Infinite elements are an extension of finite elements in concept, and the main idea is to simulate an infinite physical field by geometrically bounding a "finite" element of infinite size. Also, since the infinite element must reflect the boundary characteristics of the near field or combine with finite elements that model the near field, it tends to be infinite in one direction only, and is therefore also referred to as a semi-infinite element. Infinite elements still fall into the category of finite elements in the broad finite element concept. In summary, infinite element is proposed to overcome the problem of infinite domain in finite element, and is often used to solve more complicated infinite problem with the conventional finite element, which is a supplement to finite element method, so it is more advantageous than other numerical methods such as boundary element to solve infinite domain problem.
A large number of experience in solving the infinite domain problem has shown that: the finite element and infinite element coupling model has wide practicability in solving the practical engineering problem. In particular, in the research of the problem of endogenous fluctuation, such as train vibration and explosion, and the seismic dynamics problem of exogenous vibration, infinite elements show obvious superiority in simulating and approximately simulating infinite domain problems. In addition, the finite element and infinite element coupling model is widely applied to the fields of electromagnetism, thermodynamics, acoustics and the like, and obtains good simulation effect. A plurality of infinite unit types are improved in large-scale finite element calculation software ABAQUS and ANSYS, and the infinite domain problem of finite element calculation is simulated. The ABAQUS is widely applied to power calculation in the field of geotechnical engineering by the powerful power nonlinear calculation capability.
ABAQUS provides first-order and second-order infinite elements including plane strain, plane stress, axisymmetric and three-dimensional infinite element elements, which are based on static calculation components such as ZienkiewiczThe simulation method is developed by dynamic response analysis such as Lysmer and the like, can be used for solving the problem of static infinite domain and solving the problem of local source vibration in the domain as an infinite element dynamic artificial boundary, namely, the simulation of external traveling waves which penetrate through the artificial boundary from the finite domain to the infinite domain is effective, and the problem of external source incidence is successfully solved by carrying out secondary development on the basis of ABAQUS infinite elements. Three-dimensional infinite units common to ABAQUS infinite element power artificial boundaries include CIN3D8, CIN3D12R(S)And CIN3D18R (S), etc., which may be combined with standard finite elements to simulate the near field region with finite elements and the far field region with infinite elements.
Although the ABAQUS boundary infinite unit provides convenience for static infinite element and dynamic artificial boundary simulation. However, since the node number in the infinite element ensures that the first surface of the element is the interface of the finite element and the infinite element to ensure the directivity of the infinite element, the extending direction of the element is from the near field to the far field. An infinite unit cannot be directly defined in the ABAQUS, only an infinite element boundary part can be preset during modeling, other unit types are adopted for distinguishing, the infinite element boundary part is exported to an ABAQUS input file (. inp), and then the comparison model manually modifies the unit node sequence of the infinite element boundary of the ABAQUS input file (. inp). Particularly for a three-dimensional boundary infinite unit, the modification process is complicated, more manpower and time are consumed especially when the analysis task amount is large, and a situation that manual operation errors cause grid operation errors often occurs.
The invention content is as follows:
the invention aims to solve the technical problem of providing an ABAQUS three-dimensional infinite element boundary rapid modeling method based on an INP file, overcoming the defects of complicated modeling process and low modeling efficiency of the conventional boundary infinite element model of the ABAQUS, considering the characteristics of ABAQUS infinite elements and the characteristics of ABAQUS input INP file model definition, having simple operation and high automation degree, generating infinite element grids in a short time and greatly improving the modeling efficiency of the ABAQUS.
In order to solve the technical problems, the invention adopts the technical scheme that:
according to a unit node definition rule in an ABAQUS three-dimensional infinite unit (3D infinite element), a unit and node definition method of an ABAQUS input file (. INP) is utilized, a data line defined by the ABAQUS infinite unit is output by means of Python, an input file defined by a model is completed according to a syntax rule of an INP file, and the ABAQUS/CAE is imported to generate a boundary infinite unit.
The technical scheme adopted by the invention for solving the technical problem comprises the following steps:
establishing a finite element model, establishing a specified Part boundary node set as a near field node set of an infinite element, and exporting an ABAQUS input file (Job-1. inp);
modifying the INP file to define a far-field node of an infinite unit by using a keyword (NCOPY) and storing the remote node as a new INP file (Job-2. INP);
Abaqus/CAE imports the modified INP file (Job-2.INP), looks up the generated infinite element far-field node and exports a new INP file (Job-3. INP);
an InP file (Job-3.INP) is specified, a keyword line ([ ELEMENT ]) is inserted into a position to define an infinite unit, and the infinite unit is stored as Job-4. INP;
python programming constructs infinite element boundary unit definition data lines and inputs the infinite element boundary unit definition data lines into an input file (Job-4.INP) of ABAQUS or a specified INP file;
the INP file (Job-4.INP) is opened in a mode of importing the ABAQUS model, and the ABAQUS designates that Part three-dimensional infinite unit boundaries are automatically generated.
Further, the specific steps of establishing an infinite element near-field node set in the finite element model and exporting an ABAQUS input file (Job-1.inp) comprise:
establishing an ABAQUS/CAE three-dimensional finite element model, and establishing a node set on a boundary surface of a designated Part as a near-field node set of a corresponding boundary infinite element by using a set tool under a Part module according to a research problem, namely four nodes on a first surface of the infinite element;
the analysis task Job-1 was created under the Job Module for the model of the problem and the ABAQUS Input file (Job-1.inp) was exported on Job Manager via Write Input.
Modifying the INP file to generate an infinite element far-field node by using the keyword ([ NCOPY ]), and storing the infinite element far-field node as the main content of a new INP file (Job-2.INP) comprises the following steps:
modifying the INP file Job-1.INP, and inserting a plurality of lines of keywords ([ NCOPY ]) defined by nodes in a specified Part definition line according to different boundary regions;
copying and shifting the near field nodes by using a keyword (NCOPY) to generate the far field nodes of different area boundary infinite elements;
storing the modified INP file as a new INP file (Job-2. INP);
the steps of importing the modified INP file Job-2.INP by Abaqus/CAE, checking the generated infinite element far-field node and exporting a new INP file Job-3.INP comprise:
importing the modified INP File Job-2.INP in a File → Import → Model manner through a graphical user interface GUI of Abaqus/CAE, wherein infinite elements of the far-field node are generated but cannot be displayed temporarily;
under the condition of a Part module, checking generated infinite element far-field nodes in different areas by clicking a Show node labels of Show node labels under View → Part Display Options → Mesh;
establishing a Job-3 task by the imported model in a Job module, and exporting a new INP file Job-3.INP through WriteInput on Job Manager;
the INP file Job-3.INP specifies the position to insert the keyword line ELEMENTS to define an infinite unit, and further includes for Job-4. INP:
editing an INP file Job-3.INP, and inserting a keyword line ELEMENT in a designated line to define an infinite unit;
and the INP file inserted with the ELEMENT keyword line and the data line is saved as a new INP file Job-4. INP.
Further, the main contents of the far-field node for generating infinite elements of different zone boundaries by copying and shifting the near-field node by using the keyword (. times.NCOPY) comprise:
the key word lines and data lines defined by the nodes must be edited according to the grammar rules of the INP file of ABAQUS;
new nodes are defined by NCOPY key rows and data rows.
Further, the key line and data line must be edited according to the syntax rules of the INP file of ABAQUS by the following steps:
1) the far-field node defines a Keyword Line (Keyword Line) as:
*NCOPY,CHANGE NUMBER＝N,OLD SET＝near_field_NodeSetName,NEW SET＝far_field_NodeSetName,SHIFT
2) the far-field node definition Data Line (Data Line) is:
Shifting-x,Shifting-y,Shifting-z
the key line and data line insertion position defined by the above nodes must be located before the definition End line End Part key of Part where the infinite element boundary is located, and behind the definition data line of the near field node set.
Defining a new node by NCOPY key line and data line mainly refers to:
copying and offsetting a certain distance through a near-field node, defining an infinite element boundary far-field node, and adding an integer value N to the node number of a new node on the basis of the original node number;
the near-field nodes and far-field nodes in the infinite element direction of the same infinite unit are in one-to-one correspondence, and the unit numbers have a difference of N;
the data line defines coordinate offset values of the newly created far-field node on the basis of the original near-field node.
Further, editing the INP file (Job-3.INP) and inserting a keyword (, ELEMENT) at a specific position to define details of an infinite ELEMENT means:
defining a new infinite cell by an ELEMENT key row and a data row;
the key word lines and data lines defined by the infinite unit must be edited according to the grammar rules of the INP file of ABAQUS;
the key line and data line insertion positions of the element definition must be located before the End-of-definition line End Part key of Part where the infinite element boundary is located, and the far-field node set definition data line and the finite element definition data line are located after.
Further, defining the specific content of the boundary infinite unit by the ELEMENT key word line and the data line includes:
the definition of the ABAQUS unit and the node is mainly to create a geometric shape and divide a unit grid through a preprocessor ABAQUS/CAE to form a unit node, wherein the ABAQUS infinite unit is defined in an INP file by manually editing the INP file, inserting an ELEMENT key word, assigning a unit number to the infinite unit, defining the unit through the specified unit node number, and simultaneously grouping the boundary infinite unit.
Considering that the unit number assigned on the ASSEMBLY (ASSEMBLY) is not unique, the combination Part (Part) is needed to distinguish, and the unit number assigned on the Part is unique, so the definition of the boundary infinite element is only defined on the designated Part;
an infinite unit is defined mainly through an ELEMENT keyword, the unit type is specified to be CIN3D8, and an infinite boundary ELEMENT set on a corresponding boundary area is established;
the data line defined by the infinite unit is a unit number and a node number formed by the units;
the data lines defined by the ABAQUS infinite cells can also be read from other INP files containing cell definition data;
the above-described node number arrangement order defined by the data line cell satisfies the node order defined by the three-dimensional infinite cell (CIN3D8) in the ABAQUS help document.
The key word line and data line defined by the infinite unit must be edited according to the grammar rule of INP file of ABAQUS, which specifically means that:
1) there are two ways for a Keyword Line (Keyword Line) defined by an infinite unit at an infinite unit boundary, which are respectively:
①*Element,type＝CIN3D8,Elset＝Infinite_Element_SetName
②*Element,type＝CIN3D8,Elset＝Infinite_Element_SetName,INPUT＝filename.inp
2) the corresponding infinite unit definition Data Line (Data Line) can correspond to the above key Line in two ways:
① the data line defined by the infinite element is directly inserted into the next line of the position of the infinite element definition key word of the corresponding area, which mainly comprises an infinite element number and a composition node number, namely:
Infinite_element_Num,NodeP1,NodeP2,NodeP3,NodeP4,NodeP5,NodeP6,NodeP7,NodeP8
② Inp, the file content is simply the line of data containing the current bounding area Infinity definition of the previous entry-including Infinity Unit number and component node number, i.e.:
Infinite_element_Num,NodeP1,NodeP2,NodeP3,NodeP4,NodeP5,NodeP6,NodeP7,NodeP8。
further, the data line Infinite _ element _ Num, node p1, node p2, node p3, node p4, node p5, node p6, node p7, node p8 defined by Infinite element include the following related conventions:
the node numbers of the four far field nodes of claims 5 and 6 are greater than the node number value of the corresponding near field node (node on finite element boundary) on the infinite direction element edge by N, that is:
Figure BDA0001277372020000061
according to the infinite unit definition rule, four nodes NodeP1, NodeP2, NodeP3 and NodeP4 in the near field of the first surface of the infinite unit are required to be arranged in a sequence of the four nodes in the near field in a counterclockwise direction when the four nodes are observed from the far field to the near field of the infinite unit;
the unit number designated by the infinite unit definition is determined by the number of finite units on the current Part, the number of infinite units in the defined boundary area and the order of the infinite unit definition in the area.
Further, the construction of the infinite element boundary unit definition data line by using Python programming and inputting the infinite element boundary unit definition data line into an input file (Job-4.INP) of the ABAQUS or a specified INP file specifically comprises the following steps:
1) reading a node number and a corresponding coordinate, a unit number and a component node number in an input file (Job-4.inp) of the ABAQUS by utilizing Python programming, and storing the node number and the corresponding coordinate, the unit number and the component node number in a specified list;
2) reading a key line number defined by an infinite unit in an INPUT file (Job-4.INP) of ABAQUS by utilizing Python programming and storing the key line number in a specified variable, or reading an INPUT parameter written by an infinite unit definition data line-the INP file.
3) Traversing the node list of the finite element model according to the characteristics of the node coordinates, searching nodes on different boundaries according to the characteristics of the coordinate values of the nodes of the boundaries, and adding the nodes to the specified boundary node list;
4) searching boundary units by traversing the unit list of the finite element model according to the relationship between the unit composition node coordinates and the characteristics of the node coordinates on the boundary, and adding the boundary units to the specified boundary unit list;
5) traversing the appointed boundary node list, taking out any four nodes A, B, C and D, and judging whether the four nodes are four nodes on the first surface of the boundary infinite unit near field;
6) if the four nodes A, B, C and D are four nodes on the first surface of the boundary infinite unit near field, solving a vector included angle relation according to the quantity product of the vectors to determine the relative position relation of the four nodes;
7) determining the arrangement sequence of the four nodes in the counterclockwise direction as Q1, Q2, Q3 and Q4 by using the vector product among the vectors;
8) the infinite unit serial number is determined by a Part unit number a specified in the finite model, a defined boundary infinite unit number b and an infinite unit definition sequence c in the region, and then the current infinite unit number is Q ═ a + b + c;
9) given a data line output specifying a boundary cell infinite element definition:
Q,Q1,Q2,Q3,Q4,Q1+N,Q2+N,Q3+N,Q4+N
wherein: q is an infinite element number, Q1, Q2, Q3, Q4 are four near-field nodes, Q1+ N, Q2+ N, Q3+ N, and Q4+ N are four far-field nodes;
10) two infinite unit definition data line output modes, wherein the former mainly adopts a key line ① and a data line ① aiming at infinite unit definition, namely the unit definition data line is directly inserted into the next line of the key definition position in Job-4.inp and is stored;
the latter mainly adopts a key line ② and a data line ② for the infinite unit definition, writes the data line of the unit definition into an INPUT parameter file.
Further, the specific steps of traversing the specified boundary node list, taking out any four nodes, and judging whether the node is four nodes on the first surface of the infinite meta-model near field are as follows:
by using the member test function of Python, if four nodes taken out from the boundary node list are four of the boundary unit composition nodes, the four nodes are considered as four near-field nodes on the first surface of the boundary infinite unit.
If the four nodes A, B, C and D are four nodes on the first surface of the boundary infinite unit near field, solving two vector included angles according to the vector quantity product to determine the relative position relation of the four nodes, and the specific steps comprise:
with point a as the first point of the quadrilateral, the sizes of ∠ BAC, &lttttransition = & &ltt/t &gttbad and ∠ CAD are respectively calculated using the vector quantity product and summed:
Figure BDA0001277372020000081
let SUM ∠ BAC + ∠ BAD + ∠ CAD:
1) if SUM is 2 ∠ BAC, point a and point D, point B and point C are diagonal nodes, and the counterclockwise direction of the nodes is a, B, D, C or a, C, D, B;
2) if SUM is 2 ∠ BAD, point a and point C, point B and point D are diagonal nodes to each other, and the counterclockwise direction of the nodes is sequentially a, B, C, D or a, D, C, B;
3) if SUM is 2 ∠ CAD, point a and point B, point C and point D are diagonal nodes to each other, and the nodes are in order a, D, B, C or a, C, B, D counterclockwise.
The specific steps of determining the arrangement sequence of the four nodes in the counterclockwise direction by using the vector product among the vectors comprise:
after the relative position relations of four near-field nodes P1, P2, P3 and P4 in the first plane of the boundary infinite unit are determined, two vectors are formed by connecting any three adjacent nodes in the four nodes end to end, and assuming that P1 and P3 and P2 and P4 are corner points of each otherTaking the point P1 as a starting point, the adjacent nodes are P2 and P4, the anticlockwise sequence of the four nodes is P1, P2, P3, P4 or P1, P4, P3 and P2, and a vector is calculated
Figure BDA0001277372020000082
And
Figure BDA0001277372020000083
cross product of
Figure BDA0001277372020000084
Figure BDA0001277372020000085
If it is
Figure BDA0001277372020000091
Normal vector to point P1 on the current boundary surface
Figure BDA0001277372020000092
When the included angle is an acute angle, namely:
Figure BDA0001277372020000093
the points P1, P2 and P3 are arranged in a counterclockwise direction, and the counterclockwise arrangement sequence of the four nodes is P1, P2, P3 and P4;
if it is
Figure BDA0001277372020000094
Normal vector to point P1 on the current boundary surface
Figure BDA0001277372020000095
When the included angle is an obtuse angle, namely:
Figure BDA0001277372020000096
the points P1, P2 and P3 are arranged clockwise, and the four nodes are arranged counterclockwise in the order of P1, P4, P3 and P2.
Further, the INP file (Job-4.INP) is opened in an ABAQUS model importing mode, and the specific steps of the ABAQUS that the Part three-dimensional infinite unit boundary is automatically generated comprise:
opening ABAQUS/CAE, opening an INP File (Job-4.INP) under a target folder through a Model Import mode (Import → Model) under a File menu, completing the modeling process of a boundary infinite unit, and checking the infinite element modeling effect under a specified Part.
Compared with the prior art, the invention has the beneficial effects that: the method mainly comprises the steps of utilizing a unit and a node definition method of an ABAQUS input file (.inp) according to a unit node definition rule in an ABAQUS three-dimensional infinite unit, completing automatic output of an infinite unit definition data line by certain Python programming, completing an input file of model definition according to a grammar rule of an INP file, and importing the ABAQUS/CAE to generate a boundary infinite unit. The method and the device have the advantages of avoiding the complicated process of manually modifying the INP file in the traditional method, having high automation degree, generating the infinite unit grid in a short time and greatly improving the modeling efficiency of the ABAQUS, along with simple operation.
Drawings
FIG. 1 is a flow chart of the technical solution of the present invention.
Figure 2 is an ABAQUS three-dimensional finite element model.
Fig. 3 is an infinite element far-field node generated using NCOPY.
FIG. 4 is a diagram of generating infinite element effects.
FIG. 5 is a model bottom surface bounding region infinite cell.
FIG. 6 is a model left boundary region infinite cell.
FIG. 7 is a model right border area infinite cell.
FIG. 8 is a model rear boundary region infinite element.
FIG. 9 is a model front edge boundary region infinite unit.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and specific embodiments.
As shown in fig. 1, the ABAQUS three-dimensional infinite element boundary fast modeling method based on the INP file in the embodiment of the present invention includes:
the method comprises the steps that a three-dimensional model is divided based on a three-dimensional unit grid, the rear central point of the three-dimensional model is located on the origin of a coordinate system, and the X axis, the Y axis and the Z axis of the coordinate system are the symmetry axes of the three-dimensional model;
101. establishing a finite element model, establishing a specified Part boundary node set as a near field node set of an infinite element, and exporting an ABAQUS input file (Job-1. inp);
establishing an ABAQUS/CAE three-dimensional finite element model as shown in FIG. 2, establishing a node set on a boundary surface of a designated Part as a near-field node (four nodes on a first surface of an infinite unit) set of a corresponding boundary infinite element by using a set tool under a Part module according to a research problem, wherein the node set conditions of different boundary areas are shown in Table 1; the analysis task Job-1 was created under the Job Module for the model of the problem and the ABAQUS Input file (Job-1.inp) was exported on Job Manager via Write Input.
102. Modifying the INP file to define a far-field node of an infinite unit by using a keyword (NCOPY) and storing the remote node as a new INP file (Job-2. INP);
modifying the INP file Job-1.INP, inserting a plurality of rows of node-defined keywords ([ NCOPY ]) in the specified Part definition line according to different boundary regions, and commanding as shown below
*Ncopy,Change Number＝2000,Old Set＝Node-Y1,New Set＝Node-Y2,Shift
0.0,-40.0,0.0
*Ncopy,Change Number＝4000,Old Set＝Node-X1,New Set＝Node-X3,Shift
-40.0,0.0,0.0
*Ncopy,Change Number＝6000,Old Set＝Node-X2,New Set＝Node-X4,Shift
40.0,0.0,0.0
*Ncopy,Change Number＝8000,Old Set＝Node-Z1,New Set＝Node-Z3,Shift
0.0,0.0,40.0
*Ncopy,Change Number＝10000,Old Set＝Node-Z2,New Set＝Node-Z4,Shift
0.0,0.0,-40.0
*Ncopy,Change Number＝12000,Old Set＝Node-X1Z1,New Set＝Node-X3Z3,Shift
-40.0,0.0,40.0
*Ncopy,Change Number＝14000,Old Set＝Node-X1Z2,New Set＝Node-X3Z4,Shift
-40.0,0.0,-40.0
*Ncopy,Change Number＝16000,Old Set＝Node-X2Z1,New Set＝Node-X4Z3,Shift
40.0,0.0,40.0
*Ncopy,Change Number＝18000,Old Set＝Node-X2Z2,New Set＝Node-X4Z4,Shift
40.0,0.0,-40.0
Copying the near field nodes by the specified offset distance by using the key word (NCOPY) to generate a far field node set with infinite elements at different zone boundaries as shown in Table 2, and storing the modified INP file as a new INP file (Job-2. INP).
The infinite element boundary far-field nodes are defined by copying and offsetting the near-field nodes for a certain distance, the node number of a new node is added with an integer value N on the basis of the original node number, one-to-one correspondence between the near-field nodes and the far-field nodes in the infinite element direction of the same infinite element is realized, the unit number difference N is obtained, the N value is selected in combination with the total number of nodes in the model, and the node number of the original finite element is ensured not to be repeated with the newly generated node number. The data line defines coordinate offset values of the newly created far-field node on the basis of the original near-field node.
TABLE 1 different boundary regions Infinite element near field and far field node sets
Figure BDA0001277372020000111
103, importing the modified INP file (Job-2.INP) by Abaqus/CAE, checking the generated infinite element far-field node and exporting a new INP file (Job-3. INP);
importing the modified INP File (Job-2.INP) in a File → Import → Model mode through a GUI (graphical user interface) of Abaqus/CAE, generating an infinite element far-field Node, under the condition of a Part module, clicking a Show Node display unit Node number under View → Part display Options → Mesh to View the generated infinite element far-field Node in different areas, as shown in FIG. 3, placing the imported Model in a Job module, creating a Job-3 task, and exporting a new INP File (Job-3.INP) through Write Input on JobManager.
Inserting a keyword line ([ ELEMENT ]) into a designated position of an INP file (Job-3.INP) to define an infinite unit, and storing the infinite unit as Job-4. INP;
the INP file (Job-3.INP) is edited according to the syntax rules of the INP file of ABAQUS, a new ELEMENT is defined by ELEMENT key lines and data lines, and the key line and data line insertion position defined by the ELEMENT must be located before the End-of-definition line End Part key of Part where the infinite ELEMENT boundary is located, and the far-field node set definition data line and the finite ELEMENT definition data line follow.
An INP file of ABAQUS mainly defines an infinite unit through an ELEMENT keyword, specifies the unit type as CIN3D8, and establishes an infinite boundary ELEMENT set on a corresponding boundary area; the data line is the unit number and the node number formed by the unit, the definition rule of the data line is given to the three-dimensional infinite unit ABAQUS help document, and the data line can also be read from other INP files containing unit definition data.
The above-described node number arrangement order defined by the data line cell satisfies the node order defined by the three-dimensional infinite cell (CIN3D8) in the ABAQUS help document.
The INP file of ABAQUS uses ELEMENT to perform infinite ELEMENT definition in two ways:
1) the data line defined by the infinite unit is directly inserted into the next line of the position of the infinite element definition key word of the corresponding area, and the command is as follows;
*Element,type＝CIN3D8,Elset＝bottomy
*Element,type＝CIN3D8,Elset＝leftx
*Element,type＝CIN3D8,Elset＝rightx
*Element,type＝CIN3D8,Elset＝frontz
*Element,type＝CIN3D8,Elset＝backz
*Element,type＝CIN3D8,Elset＝leftfrontcorner
*Element,type＝CIN3D8,Elset＝leftbackcorner
*Element,type＝CIN3D8,Elset＝rightfrontcorner
*Element,type＝CIN3D8,Elset＝rightbackcorner
*Element,type＝CIN3D8,Elset＝frontleftcorner
*Element,type＝CIN3D8,Elset＝frontrightcorner
*Element,type＝CIN3D8,Elset＝backleftcorner
*Element,type＝CIN3D8,Elset＝backrightcorner
2) inp, the content of which is simply the data line containing the current bounding area infinite unit definition that has previously completed the input-including the infinite unit number and the component node number, the command is as follows.
*Element,type＝CIN3D8,Elset＝bottomy,INPUT＝bottomy_element1.inp
*Element,type＝CIN3D8,Elset＝leftx,INPUT＝leftx_element1.inp
*Element,type＝CIN3D8,Elset＝rightx,INPUT＝rightx_element1.inp
*Element,type＝CIN3D8,Elset＝frontz,INPUT＝frontz_element1.inp
*Element,type＝CIN3D8,Elset＝backz,INPUT＝backz_element1.inp
*Element,type＝CIN3D8,Elset＝leftfrontcorner,INPUT＝leftfrontcorner_element.inp
*Element,type＝CIN3D8,Elset＝leftbackcorner,INPUT＝leftbackcorner_element.inp
*Element,type＝CIN3D8,Elset＝rightfrontcorner,INPUT＝rightfrontcorner_element.inp
*Element,type＝CIN3D8,Elset＝rightbackcorner,INPUT＝rightbackcorner_element.inp
*Element,type＝CIN3D8,Elset＝frontleftcorner,INPUT＝frontleftcorner_element.inp
*Element,type＝CIN3D8,Elset＝frontrightcorner,INPUT＝frontrightcorner_element.inp
*Element,type＝CIN3D8,Elset＝backleftcorner,INPUT＝backleftcorner_element.inp
*Element,type＝CIN3D8,Elset＝backrightcorner,INPUT＝backrightcorner_element.inp
Python programming constructs an infinite element boundary unit definition data line and inputs the infinite element boundary unit definition data line into an input file (Job-4.INP) of ABAQUS or a specified INP file, and the specific process is as follows:
1) reading a node number and a corresponding coordinate, a unit number and a component node number in an input file (Job-4.inp) of the ABAQUS by utilizing Python programming, and storing the node number and the corresponding coordinate, the unit number and the component node number in a specified list;
2) reading a key line number defined by an infinite unit in an INPUT file (Job-4.inp) of ABAQUS by utilizing Python programming and storing the key line number in a specified variable, or reading an INPUT parameter filename. inp written by an infinite unit definition data line;
3) traversing the node list of the finite element model according to the characteristics of the node coordinates, searching nodes on different boundaries according to the characteristics of the coordinate values of the nodes of the boundaries, and adding the nodes to the specified boundary node list;
4) searching boundary units by traversing the unit list of the finite element model according to the relationship between the unit composition node coordinates and the characteristics of the node coordinates on the boundary, and adding the boundary units to the specified boundary unit list;
5) traversing the appointed boundary node list, taking out any four nodes A, B, C and D, and judging whether the four nodes are four nodes on the first surface of the boundary infinite unit near field;
6) if the four nodes A, B, C and D are four nodes on the first surface of the boundary infinite unit near field, solving a vector included angle relation according to the quantity product of the vectors to determine the relative position relation of the four nodes;
7) determining the arrangement sequence of the four nodes in the counterclockwise direction as Q1, Q2, Q3 and Q4 by using the vector product among the vectors;
8) the infinite unit serial number is determined by a Part unit number a specified in the finite model, a defined boundary infinite unit number b and an infinite unit definition sequence c in the region, and then the current infinite unit number is Q ═ a + b + c;
9) given a data line output specifying a boundary cell infinite element definition:
Q,Q1,Q2,Q3,Q4,Q1+N,Q2+N,Q3+N,Q4+N
wherein: q is an infinite element number, Q1, Q2, Q3, Q4 are four near-field nodes, Q1+ N, Q2+ N, Q3+ N, and Q4+ N are four far-field nodes;
10) two infinite unit definition data line output modes, wherein the former mainly adopts a key line ① and a data line ① aiming at infinite unit definition, namely the unit definition data line is directly inserted into the next line of the key definition position in Job-4.inp and is stored;
the latter mainly adopts a key line ② and a data line ② for the infinite unit definition, writes the data line of the unit definition into an INPUT parameter file.
106. The INP file (Job-4.INP) is opened in a mode of importing the ABAQUS model, and the ABAQUS designates that Part three-dimensional infinite unit boundaries are automatically generated.
Opening ABAQUS/CAE, opening an INP File (Job-4.INP) under a target folder through a Model Import mode (Import → Model) under a File menu, completing a modeling process of a boundary infinite unit, and viewing an infinite element modeling effect under a designated Part, as shown in FIG. 4. Fig. 5-9 are diagrams of infinite elements of the generation of different bounding regions.

Claims (10)

1.一种基于INP文件的ABAQUS三维无限元边界快速建模方法，其特征在于，根据ABAQUS三维无限单元中单元节点定义规则，利用ABAQUS输入文件.inp的单元和节点定义方法，借助Python编程实现ABAQUS无限单元定义数据行的自动输出，按照INP文件的语法规则完成模型定义的输入文件，导入ABAQUS/CAE生成边界无限单元；其具体步骤包括：1. a kind of ABAQUS three-dimensional infinite element boundary rapid modeling method based on INP file, it is characterized in that, according to the element node definition rule in ABAQUS three-dimensional infinite element, utilize the element and the node definition method of ABAQUS input file .inp, realize with the aid of Python programming The ABAQUS infinite element defines the automatic output of the data line, completes the input file of the model definition according to the grammar rules of the INP file, and imports the ABAQUS/CAE to generate the boundary infinite element; the specific steps include: 建立有限元模型并创建指定Part边界节点集作为无限元的近场节点集，导出该ABAQUS输入文件Job-1.inp；Establish a finite element model and create a specified Part boundary node set as a near-field node set of an infinite element, and export the ABAQUS input file Job-1.inp; 修改INP文件，利用关键字*NCOPY定义无限单元的远场节点，并另存为新的INP文件Job-2.inp；Modify the INP file, use the keyword *NCOPY to define the far-field node of the infinite element, and save it as a new INP file Job-2.inp; 上述远场节点定义的关键字行和数据行插入位置必须位于无限元边界所在Part的定义结束行*End Part关键字之前，近场节点集合定义数据行之后；The insertion position of the keyword row and data row defined by the above far-field node must be located before the definition end row *End Part keyword of the Part where the infinite element boundary is located, and after the data row defined by the near-field node set; ABAQUS/CAE导入修改后INP文件Job-2.inp，查看生成的无限元远场节点并导出新INP文件Job-3.inp；ABAQUS/CAE imports the modified INP file Job-2.inp, checks the generated infinite element far-field node and exports the new INP file Job-3.inp; INP文件Job-3.inp指定位置插入关键字行*ELEMENT定义无限单元，并另存为Job-4.inp；The INP file Job-3.inp inserts the keyword line *ELEMENT to define an infinite unit at the specified position and saves it as Job-4.inp; 所述指定位置位于无限元边界所在Part的定义结束行*End Part关键字之前，远场节点集合定义数据行和有限元单元定义数据行之后；The specified position is located before the definition end row *End Part keyword of Part where the infinite element boundary is located, and after the far-field node set definition data row and the finite element element definition data row; Python编程构建无限元边界单元定义数据行并输入到ABAQUS的输入文件Job-4.inp或者指定的INP文件；Python programming constructs the infinite element boundary element definition data line and inputs it to the input file Job-4.inp of ABAQUS or the specified INP file; 通过ABAQUS模型导入的方式打开所述指定的INP文件或Job-4.inp，ABAQUS指定Part三维无限单元边界自动生成完毕。Open the specified INP file or Job-4.inp by importing the ABAQUS model, and the ABAQUS specified Part three-dimensional infinite element boundary is automatically generated. 2.根据权利要求1所述的基于INP文件的ABAQUS三维无限元边界快速建模方法，其特征在于，2. the ABAQUS three-dimensional infinite element boundary rapid modeling method based on INP file according to claim 1, is characterized in that, 建立有限元模型中无限元近场节点集的并导出ABAQUS输入文件Job-1.inp包括：Create an infinite element near-field node set in the finite element model and export the ABAQUS input file Job-1.inp including: 建立ABAQUS/CAE三维有限元模型，根据研究问题在Part模块下利用set工具建立指定Part的边界面上的节点集作为对应边界无限元的近场节点集，即无限单元第一个面上的四个节点；Establish an ABAQUS/CAE three-dimensional finite element model, and use the set tool under the Part module to establish the node set on the boundary surface of the specified Part as the near-field node set corresponding to the boundary infinite element, that is, the four nodes on the first surface of the infinite element. nodes; 在Job模块下为所述研究问题的模型建立分析任务Job-1，并在Job Manager上通过Write Input导出该ABAQUS输入文件Job-1.inp；Establish an analysis task Job-1 for the model of the research problem under the Job module, and export the ABAQUS input file Job-1.inp through Write Input on the Job Manager; 修改INP文件，利用关键字*NCOPY生成无限单元的远场节点，并另存为新的INP文件Job-2.inp包括：Modify the INP file, use the keyword *NCOPY to generate the far-field node of the infinite element, and save it as a new INP file Job-2.inp including: 修改INP文件Job-1.inp，根据不同的边界区域在指定Part定义行内插入多行节点定义的关键字*NCOPY；Modify the INP file Job-1.inp, and insert the keyword *NCOPY defined by multi-line nodes in the specified Part definition line according to different boundary areas; 利用关键字*NCOPY成不同区域边界无限元的远场节点；Use the keyword *NCOPY to form the far-field nodes of infinite elements at the boundaries of different regions; 将修改后的INP文件另存为新的INP文件Job-2.inp；Save the modified INP file as a new INP file Job-2.inp; ABAQUS/CAE导入修改后INP文件Job-2.inp，查看生成的无限元远场节点并导出新INP文件Job-3.inp包括：ABAQUS/CAE imports the modified INP file Job-2.inp, view the generated infinite element far-field node and export the new INP file Job-3.inp including: 通过ABAQUS/CAE的图形用户界面GUI以File→Import→Model的方式导入修改后的INP文件Job-2.inp，无限元的远场节点生成，但暂时无法显示；Import the modified INP file Job-2.inp through the graphical user interface GUI of ABAQUS/CAE in the way of File→Import→Model, and the infinite element far-field node is generated, but it cannot be displayed temporarily; 在Part模块条件下，通过View→Part Display Options→Mesh下点选Show NodeLabels显示单元节点号，查看生成的不同区域无限元远场节点；Under the condition of Part module, click Show NodeLabels under View→Part Display Options→Mesh to display the element node number, and view the generated infinite element far-field nodes in different regions; 导入模型在Job模块内，新建Job-3任务，通过Job Manager上的Write Input导出新的INP文件Job-3.inp；Import the model in the Job module, create a new Job-3 task, and export the new INP file Job-3.inp through Write Input on the Job Manager; INP文件Job-3.inp指定位置插入关键字行*ELEMENT定义无限单元，并另存为Job-4.inp包括：The INP file Job-3.inp inserts the keyword line *ELEMENT to define an infinite unit at the specified location, and saves it as Job-4.inp including: 编辑INP文件Job-3.inp，在指定行插入关键字行*ELEMENT进行无限单元定义；Edit the INP file Job-3.inp, insert the keyword line *ELEMENT in the specified line to define infinite units; 插入*ELEMENT关键字行和数据行后的INP文件保存为新的INP文件Job-4.inp。The INP file after inserting the *ELEMENT keyword line and data line is saved as a new INP file Job-4.inp. 3.根据权利要求2所述的基于INP文件的ABAQUS三维无限元边界快速建模方法，其特征在于，利用关键字*NCOPY生成不同区域边界无限元的远场节点包括：3. the ABAQUS three-dimensional infinite element boundary fast modeling method based on INP file according to claim 2, is characterized in that, utilizes the keyword *NCOPY to generate the far-field node of different area boundary infinite elements and comprises: 节点定义的关键字行和数据行必须按照ABAQUS的INP文件的语法规则编辑；The keyword line and data line defined by the node must be edited according to the syntax rules of the INP file of ABAQUS; 通过*NCOPY关键字行和数据行定义新节点。New nodes are defined by *NCOPY keyword lines and data lines. 4.根据权利要求3所述的基于INP文件的ABAQUS三维无限元边界快速建模方法，其特征在于，关键字行和数据行必须按照ABAQUS的INP文件的语法规则编辑包括：4. the ABAQUS three-dimensional infinite element boundary rapid modeling method based on INP file according to claim 3, is characterized in that, keyword row and data row must edit according to the grammar rule of the INP file of ABAQUS: 1)远场节点定义关键字行Keyword Line为：1) The far-field node defines the Keyword Line as: *NCOPY,CHANGE NUMBER＝N,OLD SET＝near_field_NodeSetName,NEW SET＝far_field_NodeSetName,SHIFT*NCOPY,CHANGE NUMBER=N,OLD SET=near_field_NodeSetName,NEW SET=far_field_NodeSetName,SHIFT 2)远场节点定义数据行Data Line为：2) The far-field node defines the data line Data Line as: Shifting-x,Shifting-y,Shifting-zShifting-x, Shifting-y, Shifting-z 所述节点定义的关键字行和数据行插入位置必须位于无限元边界所在Part的定义结束行*End Part关键字之前，近场节点集合定义数据行之后；The keyword row of the node definition and the data row insertion position must be located before the definition end row *End Part keyword of the Part where the infinite element boundary is located, and after the near field node set defines the data row; 通过*NCOPY关键字行和数据行定义新节点包括：Defining new nodes via *NCOPY keyword lines and data lines includes: 通过近场节点复制偏移指定距离，定义无限元边界远场节点，新节点的节点号在原有节点号基础上增加整数值N；The specified distance is copied and offset by the near-field node, and the far-field node of the infinite element boundary is defined, and the node number of the new node is increased by an integer value N based on the original node number; 实现了同一无限单元延无限元方向上近场节点和远场节点一一对应，且单元号相差N；A one-to-one correspondence between near-field nodes and far-field nodes in the direction of the same infinite element extending from the infinite element is realized, and the element numbers differ by N; 数据行定义了新创建的远场节点在原有近场节点基础上的坐标偏移值。The data line defines the coordinate offset value of the newly created far-field node based on the original near-field node. 5.根据权利要求2所述的基于INP文件的ABAQUS三维无限元边界快速建模方法，其特征在于，编辑INP文件Job-3.inp，在指定位置插入关键字*ELEMENT行进行无限单元定义包括：5. the ABAQUS three-dimensional infinite element boundary rapid modeling method based on INP file according to claim 2, is characterized in that, edit INP file Job-3.inp, insert keyword *ELEMENT row at specified position and carry out infinite element definition including : 通过*ELEMENT关键字行和数据行定义无限单元；Define infinite elements through *ELEMENT keyword lines and data lines; 无限单元定义的关键字行和数据行必须按照ABAQUS的INP文件的语法规则编辑；The keyword line and data line of the infinite element definition must be edited according to the syntax rules of the INP file of ABAQUS; 所述节点定义的关键字行和数据行插入位置必须位于无限元边界所在Part的定义结束行*End Part关键字之前，远场节点集合定义数据行和有限元单元定义数据行之后。The keyword line and data line insertion position of the node definition must be located before the definition end line *End Part keyword of Part where the boundary of the infinite element is located, and after the far-field node set definition data line and the finite element element definition data line. 6.根据权利要求3所述的基于INP文件的ABAQUS三维无限元边界快速建模方法，其特征在于，通过*ELEMENT关键字行和数据行定义边界无限单元包括：6. the ABAQUS three-dimensional infinite element boundary rapid modeling method based on INP file according to claim 3, is characterized in that, by *ELEMENT keyword row and data row definition boundary infinite element comprises: ABAQUS的单元和节点的定义是通过其前处理器ABAQUS/CAE，创建几何形状并划分单元网格，形成单元节点，这里的ABAQUS无限单元在INP文件中定义是通过手动编辑INP文件，插入*ELEMENT关键字，给无限单元分配单元号，并通过指定单元节点号定义单元，同时给边界无限元分组；The definition of elements and nodes of ABAQUS is through its preprocessor ABAQUS/CAE, which creates geometric shapes and divides the element mesh to form element nodes. The ABAQUS infinite element here is defined in the INP file by manually editing the INP file and inserting *ELEMENT keyword, assign a unit number to the infinite element, define the element by specifying the element node number, and group the boundary infinite elements at the same time; 考虑到在装配件ASSEMBLY上分配的单元号不唯一，需要结合部件Part区分，而在部件上分配的单元号唯一，所以边界无限元的定义只在指定的Part定义；Considering that the unit number assigned on the assembly ASSEMBLY is not unique, it needs to be distinguished by the part Part, and the unit number assigned on the part is unique, so the definition of the boundary infinite element is only defined in the specified Part; 通过*ELEMENT关键字定义无限单元，指定单元类型为CIN3D8，并建立相应边界区域上无限边界元集合；Define an infinite element through the *ELEMENT keyword, specify the element type as CIN3D8, and establish an infinite boundary element set on the corresponding boundary area; 无限单元定义的数据行是指定单元号和单元组成的节点号；The data line defined by the infinite element is the specified element number and the node number composed of the element; ABAQUS无限单元定义的数据行或从其他包含单元定义数据的INP文件中读取；ABAQUS infinite element definition data lines or read from other INP files containing element definition data; 数据行单元定义的节点号排列顺序满足ABAQUS帮助文档中三维无限单元CIN3D8定义的节点次序；The order of the node numbers defined by the data line element satisfies the node order defined by the three-dimensional infinite element CIN3D8 in the ABAQUS help document; 无限单元定义的关键字行和数据行必须按照ABAQUS的INP文件的语法规则编辑包括：The keyword line and data line of the infinite element definition must be edited according to the syntax rules of the INP file of ABAQUS including: 1)无限元边界无限单元定义的关键字行KeywordLine有两种方式，分别为：1) There are two ways to define the keyword line KeywordLine of the infinite element boundary, namely: ①*Element,type＝CIN3D8,Elset＝Infinite_Element_SetName①*Element, type=CIN3D8, Elset=Infinite_Element_SetName ②*Element,type＝CIN3D8,Elset＝Infinite_Element_SetName,INPUT＝filename.inp②*Element, type=CIN3D8, Elset=Infinite_Element_SetName, INPUT=filename.inp 2)相应的无限单元定义数据行DataLine由两种方式，分别与关键字行对应：2) The corresponding infinite unit defines the data line DataLine in two ways, corresponding to the keyword line respectively: ①无限单元定义的数据行直接插入对应区域无限元定义关键字所在位置的下一行，包括无限单元号和组成节点号，即：①The data line defined by the infinite element is directly inserted into the next line where the keyword of the infinite element definition in the corresponding region is located, including the infinite element number and the constituent node number, namely: Infinite_element_Num,NodeP1,NodeP2,NodeP3,NodeP4,NodeP5,NodeP6,NodeP7,NodeP8Infinite_element_Num,NodeP1,NodeP2,NodeP3,NodeP4,NodeP5,NodeP6,NodeP7,NodeP8 ②无限单元定义的数据行由filename.inp文件给定，该filename.inp文件内容只是包含预先完成输入的当前边界区域无限单元定义的数据行—包括无限单元号和组成节点号，即：② The data line defined by the infinite element is given by the filename.inp file, and the content of the filename.inp file is only the data line defined by the infinite element in the current boundary area of the pre-completed input - including the infinite element number and the constituent node number, namely: Infinite_element_Num,NodeP1,NodeP2,NodeP3,NodeP4,NodeP5,NodeP6,NodeP7,NodeP8。Infinite_element_Num,NodeP1,NodeP2,NodeP3,NodeP4,NodeP5,NodeP6,NodeP7,NodeP8. 7.根据权利要求6所述的基于INP文件的ABAQUS三维无限元边界快速建模方法，其特征在于，无限单元定义的数据行Infinite_element_Num,NodeP1,NodeP2,NodeP3,NodeP4,NodeP5,NodeP6,NodeP7,NodeP8包括：7. ABAQUS three-dimensional infinite element boundary fast modeling method based on INP file according to claim 6, is characterized in that, the data row Infinite_element_Num of infinite element definition, NodeP1, NodeP2, NodeP3, NodeP4, NodeP5, NodeP6, NodeP7, NodeP8 include: 远场四个节点的节点号比对应沿无限方向单元边上近场节点，即有限元边界上节点，的节点号值大N，即：The node numbers of the four nodes in the far field are larger by N than the node numbers corresponding to the near field nodes along the edge of the element along the infinite direction, that is, the nodes on the finite element boundary, namely:
Figure FDA0002370126550000051
Figure FDA0002370126550000051
 根据无限单元定义规则，无限单元第一面内近场四个节点NodeP1、NodeP2、NodeP3和NodeP4，要满足从无限单元远场向近场观测时近场四个节点为逆时针方向的顺序排列；According to the definition rules of infinite elements, the four nodes NodeP1, NodeP2, NodeP3 and NodeP4 in the near field of the first plane of the infinite element must be arranged in a counterclockwise order when observing from the far field of the infinite element to the near field; 无限单元定义指定的单元号是由当前Part上有限单元个数，已定义边界区域无限单元个数以及区域内无限单元定义的次序决定。The element number specified by the infinite element definition is determined by the number of finite elements on the current Part, the number of infinite elements in the defined boundary region, and the order of the definition of infinite elements in the region. 8.根据权利要求7所述的基于INP文件的ABAQUS三维无限元边界快速建模方法，其特征在于，利用Python编程构建无限元边界单元定义数据行并输入到ABAQUS的输入文件Job-4.inp或者指定的INP文件包括以下步骤：8. the ABAQUS three-dimensional infinite element boundary fast modeling method based on INP file according to claim 7, is characterized in that, utilizes Python programming to construct infinite element boundary element definition data row and input to the input file Job-4.inp of ABAQUS Or the specified INP file includes the following steps: 1)利用Python编程读取ABAQUS的输入文件Job-4.inp中的节点号及相应坐标，单元号及组成节点号，并存储在指定的列表中；1) Use Python programming to read the node number and corresponding coordinates, unit number and component node number in the input file Job-4.inp of ABAQUS, and store them in the specified list; 2)利用Python编程读取ABAQUS的输入文件Job-4.inp中无限单元定义的关键字行号并存储在指定的变量中，或者读取无限单元定义数据行写入的INPUT参数—ABAQUS读取的INP文件filename.inp；2) Use Python programming to read the keyword line number defined by the infinite element in the input file Job-4.inp of ABAQUS and store it in the specified variable, or read the INPUT parameter written by the infinite element definition data line - ABAQUS read INP file filename.inp; 3)根据节点坐标特点遍历所述有限元模型的节点列表，根据各边界节点坐标值特点查找不同边界上节点，并添加到指定边界节点列表；3) traverse the node list of the finite element model according to the node coordinate characteristics, search for nodes on different boundaries according to the coordinate value characteristics of each boundary node, and add to the specified boundary node list; 4)通过遍历有限元模型的单元列表，依据单元组成节点坐标与边界上节点坐标特点的关系，查找边界单元，并添加到指定的边界单元列表；4) By traversing the element list of the finite element model, according to the relationship between the node coordinates of the element and the coordinates of the nodes on the boundary, find the boundary element and add it to the specified boundary element list; 5)遍历指定边界节点列表，并取出其中任意四个节点A、B、C和D，判断这四个节点是否为边界无限单元近场第一个面上的四个节点；5) Traverse the specified boundary node list, and take out any four nodes A, B, C and D among them, and judge whether these four nodes are the four nodes on the first surface of the near field of the boundary infinite element; 6)若A、B、C和D四个节点是边界无限单元近场第一个面上的四个节点，则根据向量的数量积求向量夹角关系确定四个节点的相对位置关系；6) If the four nodes A, B, C and D are the four nodes on the first surface of the near field of the boundary infinite element, the relative position relationship of the four nodes is determined by calculating the angle relationship between the vectors according to the quantity product of the vectors; 7)利用向量间的向量积确定四个节点的逆时针方向的排列顺序为Q1、Q2、Q3、Q4；7) Use the vector product between the vectors to determine that the counterclockwise order of the four nodes is Q1, Q2, Q3, Q4; 8)无限单元序号由有限模型中指定Part单元个数a，已定义边界无限单元个数为b以及该区域内无限单元定义次序c决定，则当前无限单元号为Q＝a+b+c；8) The number of the infinite element is determined by the number of Part elements a specified in the finite model, the number of infinite elements in the defined boundary is b and the definition order of infinite elements in this area is determined by c, then the current infinite element number is Q=a+b+c; 9)给出指定边界单元无限元定义的数据行输出：9) Give the data line output for the infinite element definition of the specified boundary element: Q,Q1,Q2,Q3,Q4,Q1+N,Q2+N,Q3+N,Q4+NQ,Q1,Q2,Q3,Q4,Q1+N,Q2+N,Q3+N,Q4+N 其中：Q为无限元单元号，Q1,Q2,Q3,Q4为近场四个节点，Q1+N,Q2+N,Q3+N,Q4+N为远场四个节点；Among them: Q is the infinite element unit number, Q1, Q2, Q3, Q4 are four nodes in the near field, Q1+N, Q2+N, Q3+N, Q4+N are four nodes in the far field; 10)两种无限单元定义数据行输出方式，前者是针对无限单元定义采用关键字行①和数据行①，即单元定义数据行直接插入在Job-4.inp中关键字定义位置的下一行并保存；10) Two infinite unit definition data line output methods, the former uses keyword line ① and data line ① for infinite unit definition, that is, the unit definition data line is directly inserted into the next line of the keyword definition position in Job-4.inp and save; 后者是针对无限单元定义采用关键字行②和数据行②，将单元定义的数据行写入无限单元定义行中的INPUT参数filename.inp文件中，ABAQUS在导入模型时会自动读取filename.inp文件中的数据行，生成无限单元。The latter uses keyword line ② and data line ② for the infinite element definition, and writes the data line of the element definition into the INPUT parameter filename.inp file in the infinite element definition line. ABAQUS will automatically read the filename when importing the model. Data lines in the inp file, generating infinite cells. 9.根据权利要求8所述的基于INP文件的ABAQUS三维无限元边界快速建模方法，其特征在于，遍历指定边界节点列表，并取出其中任意四个节点，判断该节点是否为无限元模型近场第一个面上的四个节点包括：9. the ABAQUS three-dimensional infinite element boundary fast modeling method based on INP file according to claim 8, is characterized in that, traverse the specified boundary node list, and take out any four nodes wherein, judge whether this node is the infinite element model near. The four nodes on the first face of the field include: 利用Python的成员测试功能，若边界节点列表上取出的四个节点为边界单元组成节点的其中四个节点，则这四个节点为边界无限单元第一个面上的四个近场节点；Using Python's member test function, if the four nodes taken out from the boundary node list are four of the nodes that constitute the boundary element, these four nodes are the four near-field nodes on the first surface of the boundary infinite element; 若A、B、C和D四个节点是边界无限单元近场第一个面上的四个节点，根据向量数量积求两个向量夹角确定四个节点的相对位置关系包括：If the four nodes A, B, C and D are the four nodes on the first surface of the near field of the boundary infinite element, the relative position relationship of the four nodes is determined by calculating the angle between the two vectors according to the product of the vector quantities: 以点A作为四边形的第一个点，利用向量数量积分别计算∠BAC,∠BAD和∠CAD的大小并求和：Taking point A as the first point of the quadrilateral, calculate the magnitudes of ∠BAC, ∠BAD and ∠CAD separately by using the vector quantity product and sum them up:
Figure FDA0002370126550000061
Figure FDA0002370126550000061
 令SUM＝∠BAC+∠BAD+∠CAD：Let SUM=∠BAC+∠BAD+∠CAD: 1)若SUM＝2∠BAC，则点A和点D，点B和点C互为对角节点，则节点逆时针方向顺序为A、B、D、C或者A、C、D、B；1) If SUM=2∠BAC, then point A and point D, point B and point C are diagonal nodes to each other, then the counterclockwise order of nodes is A, B, D, C or A, C, D, B; 2)若SUM＝2∠BAD，则点A和点C，点B和点D互为对角节点，则节点逆时针方向顺序为A、B、C、D或者A、D、C、B；2) If SUM=2∠BAD, then point A and point C, point B and point D are diagonal nodes to each other, then the counterclockwise order of nodes is A, B, C, D or A, D, C, B; 3)若SUM＝2∠CAD，则点A和点B，点C和点D互为对角节点，则节点逆时针方向顺序为A、D、B、C或者A、C、B、D；3) If SUM=2∠CAD, then point A and point B, point C and point D are diagonal nodes to each other, then the counterclockwise order of nodes is A, D, B, C or A, C, B, D; 利用向量间的向量积确定四个节点的逆时针方向的排列顺序包括：Using the vector product between vectors to determine the counterclockwise order of the four nodes includes: 确定边界无限单元第一面内的四个近场节点P1、P2、P3、P4相对位置关系后，通过四个节点中的任意三个相邻节点首尾相接形成两个向量，假设P1与P3、P2与P4互为角点，以P1点为起始点，相邻节点必为P2和P4，则四个节点的逆时针顺序为P1、P2、P3、P4或者P1、P4、P3、P2，计算向量
Figure FDA0002370126550000071
与
Figure FDA0002370126550000072
的向量积
Figure FDA0002370126550000073
After determining the relative positional relationship of the four near-field nodes P1, P2, P3, and P4 in the first surface of the boundary infinite element, two vectors are formed by connecting any three adjacent nodes of the four nodes end to end, assuming that P1 and P3 , P2 and P4 are the corners of each other, with P1 as the starting point, the adjacent nodes must be P2 and P4, then the counterclockwise order of the four nodes is P1, P2, P3, P4 or P1, P4, P3, P2, compute vector
Figure FDA0002370126550000071
and
Figure FDA0002370126550000072
the vector product of
Figure FDA0002370126550000073
Figure FDA0002370126550000074
Figure FDA0002370126550000074
 若
Figure FDA0002370126550000075
与当前边界面上点P1的法向量
Figure FDA0002370126550000076
夹角为锐角时，即：like
Figure FDA0002370126550000075
with the normal vector of point P1 on the current boundary surface
Figure FDA0002370126550000076
When the included angle is an acute angle, that is:
Figure FDA0002370126550000077
Figure FDA0002370126550000077
 则点P1、P2、P3为逆时针方向排列，四个节点的逆时针排列顺序为P1、P2、P3、P4；Then the points P1, P2, P3 are arranged counterclockwise, and the counterclockwise order of the four nodes is P1, P2, P3, P4; 若
Figure FDA0002370126550000078
与当前边界面上点P1的法向量
Figure FDA0002370126550000079
夹角为钝角时，即：like
Figure FDA0002370126550000078
with the normal vector of point P1 on the current boundary surface
Figure FDA0002370126550000079
When the included angle is an obtuse angle, that is:
Figure FDA00023701265500000710
Figure FDA00023701265500000710
 则点P1、P2、P3为顺时针方向排列，四个节点的逆时针排列顺序为P1、P4、P3、P2。Then the points P1, P2, and P3 are arranged clockwise, and the counterclockwise order of the four nodes is P1, P4, P3, and P2. 10.根据权利要求7所述的基于INP文件的ABAQUS三维无限元边界快速建模方法，其特征在于，通过ABAQUS模型导入的方式打开所述INP文件或Job-4.inp，ABAQUS指定Part三维无限单元边界自动生成完毕包括：10. the ABAQUS three-dimensional infinite element boundary rapid modeling method based on INP file according to claim 7, is characterized in that, open described INP file or Job-4.inp by the mode of ABAQUS model import, ABAQUS specifies Part three-dimensional infinite. The cell boundaries are automatically generated, including: 打开ABAQUS/CAE，通过File菜单下的模型导入方式Import→Model，打开目标文件夹下的INP文件Job-4.inp，完成边界无限单元的建模过程，在指定Part下查看无限元建模效果。Open ABAQUS/CAE, use the model import method Import→Model under the File menu, open the INP file Job-4.inp in the target folder, complete the modeling process of the boundary infinite element, and view the infinite element modeling effect under the specified Part. .
CN201710270976.4A 2017-04-24 2017-04-24 ABAQUS 3D Infinite Element Boundary Rapid Modeling Method Based on INP File Expired - Fee Related CN107145651B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710270976.4A CN107145651B (en) 2017-04-24 2017-04-24 ABAQUS 3D Infinite Element Boundary Rapid Modeling Method Based on INP File

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710270976.4A CN107145651B (en) 2017-04-24 2017-04-24 ABAQUS 3D Infinite Element Boundary Rapid Modeling Method Based on INP File

Publications (2)

Publication Number Publication Date
CN107145651A CN107145651A (en) 2017-09-08
CN107145651B true CN107145651B (en) 2020-04-28

Family

ID=59774761

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710270976.4A Expired - Fee Related CN107145651B (en) 2017-04-24 2017-04-24 ABAQUS 3D Infinite Element Boundary Rapid Modeling Method Based on INP File

Country Status (1)

Country Link
CN (1) CN107145651B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107679341B (en) * 2017-10-30 2021-03-02 南京理工大学 A finite element parametric modeling method of barrel structure
CN109446734B (en) * 2018-12-17 2022-09-27 西安建筑科技大学 Method for adding spring unit in steel pipe concrete finite element model
CN111046608B (en) * 2019-12-19 2023-07-18 东南大学 A finite element model modeling method of basalt fiber asphalt mortar
CN111209694B (en) * 2019-12-30 2023-12-15 北京工业大学 A structural identification method for truss structure stiffness and axial force
CN111159951B (en) * 2019-12-31 2023-08-18 哈尔滨工业大学（深圳） A Coupling Method of Finite Element and Boundary Element Based on ABAQUS

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102004039982A1 (en) * 2004-08-12 2006-02-23 Basf Ag Microstructure-based modeling of mechanical characteristics of block copolymer comprises generating hypothetical morphologies and cross-linking finite elements, producing several domains and simulating deformation of cross-link of elements
CN104361631A (en) * 2014-10-24 2015-02-18 哈尔滨工业大学 Method for building vibration response three-dimensional numerical models of railway roadbeds in seasonally frozen ground regions in parameterized manner

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102004039982A1 (en) * 2004-08-12 2006-02-23 Basf Ag Microstructure-based modeling of mechanical characteristics of block copolymer comprises generating hypothetical morphologies and cross-linking finite elements, producing several domains and simulating deformation of cross-link of elements
CN104361631A (en) * 2014-10-24 2015-02-18 哈尔滨工业大学 Method for building vibration response three-dimensional numerical models of railway roadbeds in seasonally frozen ground regions in parameterized manner

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
山区特殊地质条件下地基-箱基-上部结构地震响应分析;尚军;《天津城建大学学报》;20140831;第20卷(第3期);全文 *
背后空洞引起高速铁路隧道衬砌裂缝形态的FEM对比分析;雷波等;《铁道标准设计》;20150930;第59卷(第9期);全文 *

Also Published As

Publication number Publication date
CN107145651A (en) 2017-09-08

Similar Documents

Publication Publication Date Title
CN107145651B (en) ABAQUS 3D Infinite Element Boundary Rapid Modeling Method Based on INP File
Nielson On marching cubes
CN109359351B (en) Three-dimensional modeling method for pipeline and pipeline design method
Blacker Meeting the challenge for automated conformal hexahedral meshing
JP5045853B2 (en) Calculation data generation apparatus, calculation data generation method, and calculation data generation program
CN110765523A (en) BIM technology-based rapid construction method for deep foundation pit support structure
CN104091027A (en) Visual workshop layout optimization method and system
CN110569580B (en) Urban street space vitality simulation method, system, computer equipment and medium
Chen et al. Triangulated manifold meshing method preserving molecular surface topology
CN109858161A (en) A kind of Abaqus Meshing Method based on Midas modeling and Matlab conversion
CN108763625A (en) A kind of method and system of FLAC3D softwares pre-treatment
CN103559339A (en) Offshore oil underwater production virtual simulation system
CN106294572A (en) A kind of elevation map rapid generation based on AutoCAD data source and ArcGIS instrument
Jabi Parametric spatial models for energy analysis in the early design stages
US10846449B1 (en) Conversion of block model-based circuit designs into circuit implementations
CN115358123A (en) EDA model Mesh classified display method
CN114494584A (en) Method, device, device and medium for generating three-dimensional white mold based on shapefile
CN113486429A (en) Automatic construction method of space intersection structure based on interpolation algorithm
CN113220199A (en) Method for assisting in establishing finite element model
CN107194997B (en) A 3D Visualization Method of Block Theoretical Equilibrium Region Map
US7086030B1 (en) Incremental netlisting
US9141734B2 (en) System and method of refining a topological indexed mesh
Gaither et al. The National Grid Project: a system overview
Zhang Research on 3D Variable Design of Corrugated Box Based on SolidWorks
Goanță Comparative study of" revolve" command in the average-class design software taught in the higher education

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200428

Termination date: 20210424