CN112685936B - Modeling method for shell mother-of-pearl microstructure finite element analysis - Google Patents
Modeling method for shell mother-of-pearl microstructure finite element analysis Download PDFInfo
- Publication number
- CN112685936B CN112685936B CN202011565362.7A CN202011565362A CN112685936B CN 112685936 B CN112685936 B CN 112685936B CN 202011565362 A CN202011565362 A CN 202011565362A CN 112685936 B CN112685936 B CN 112685936B
- Authority
- CN
- China
- Prior art keywords
- grid
- node
- splitting
- serial number
- finite element
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to a modeling method for shell mother-of-pearl microstructure finite element analysis, which realizes the high-efficiency and high-precision modeling of a finite element, can be widely used for various types of finite element analysis work, and solves the defects of low efficiency and low precision of the existing modeling method. The method has good applicability when Cohesive element finite element modeling, parametric analysis finite element modeling with scale effect structure and large-batch grid finite element modeling are carried out.
Description
Technical Field
The invention relates to a modeling method for finite element analysis of a shell mother-of-pearl microstructure, which is used for establishing a finite element simulation analysis model. Applicable objects comprise Cohesive unit finite element modeling, parametric analysis finite element modeling with scale effect structure and large-batch grid finite element modeling.
Background
The shell nacre microstructure finite element analysis method is commonly used for structural analysis, complex unit modeling, scale effect structure parametric analysis modeling and large-batch grid modeling are required in the shell nacre microstructure finite element analysis process, and the traditional modeling method is high in modeling difficulty, low in efficiency and incapable of meeting the calculation requirements in precision. The existence of these problems makes finite element analysis more difficult, and the traditional modeling is mainly the following method:
the modeling is carried out through CAE software, a geometric model is firstly established, then grid division is carried out, and the part is complex in work, takes more time and is low in efficiency;
when a large batch of grids are established through CAE software, the grid precision is low, the grid rule degree cannot meet the calculation requirement, a large amount of repairing work needs to be carried out, and the modeling efficiency is low;
when a specific unit is modeled, such as a Cohesive unit, the difficulty in modeling through CAE software is high; at present, for the Cohesive unit modeling, the built units are generally processed by a program, and Cohesive units are inserted among the units to realize the Cohesive unit modeling; this method can be used for small-batch unit modeling, which is inefficient when modeling large-batch units.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the defects of the prior art are overcome, and the method for modeling the complex structure and parameterizing the modeling is provided. The method has the advantages that high-efficiency and high-precision modeling of the complex structure is realized, high-efficiency and high-precision parametric analysis of the finite element model is realized, the efficiency and precision of the finite element analysis are greatly improved, and the requirement of finite element analysis of the shell mother-of-pearl microstructure at the present stage is met.
The invention is realized by the following technical scheme: a modeling method for finite element analysis of shell mother-of-pearl microstructures comprises the following steps:
(1) Carrying out finite element meshing according to the geometric shape of the shell mother-of-pearl microstructure, specifically simulating a soft structure in the shell mother-of-pearl microstructure according to a zero-thickness Cohesive unit, and simulating a hard structure in the shell mother-of-pearl microstructure by using a four-node unit or a hexahedral solid unit;
(2) Determining grid node serial numbers in grids and grid serial number generation modes, wherein the grid node serial numbers and the grid serial number generation modes are the same; determining grid node serial numbers and grid serial number generation functions according to the generation mode;
(3) Determining the node serial numbers and the grid serial numbers of the shell mother-of-pearl microstructure grid according to the grid node serial numbers and the grid serial number generation function in combination with the finite element grid divided in the step (1);
(4) Determining grid node data and grid data of the shell mother-of-pearl microstructure according to the grid node serial number and the grid serial number generation mode;
the grid data are grid serial numbers and grid node serial numbers under the same grid; the grid node data is a grid node serial number and a coordinate value corresponding to the grid node serial number;
(5) And sorting the generated grid node data, the grid data, the corresponding material attributes, loads, boundary conditions and simulation parameters according to a calculation file format required by finite element analysis to complete modeling.
Preferably, when the finite element mesh divided in the step (1) is a two-dimensional mesh, the mesh node serial numbers and the mesh serial number generation mode are both expressed by two digits, wherein the ten digit is a vertical direction serial number, and the one digit is a horizontal direction serial number.
Preferably, the grid node sequence number and the grid sequence number generating function expression are as follows: f (x) = x, F (y) = y, x represents a horizontal direction number, and y represents a vertical direction number.
Preferably, the processing mode of the coincident grid nodes for the zero-thickness four-node coherent unit is to split the coincident grid nodes into three directions, namely an X axis, a Y axis and a 45-degree direction, from the original grid nodes; the Cohesive unit splitting direction in the X direction is the X-axis direction, the Cohesive unit splitting mode in the Y direction is the Y-axis direction, and the Cohesive unit splitting mode in the Y direction is a two-node splitting mode; the node splitting directions at the intersection of the X and the Y are the X axis, the Y axis and the 45-degree direction, and the four-node splitting mode is adopted.
Preferably, the mesh node data and the mesh data are determined in step (4) by:
assuming that the serial number of any grid node is YX, the coordinates are (a X (X-1) and a X (Y-1)), wherein a is the side length of the grid where the grid node is located;
assume that any grid node has a grid serial number of YX, and grid node serial numbers of YX, Y (X + 1), (Y + 1) (X + 1), and (Y + 1) X.
Preferably, when the finite element mesh divided in the step (1) is a three-dimensional mesh, the mesh node and the mesh sequence number generation mode are both expressed by three digits, wherein ten digits are sequence numbers in the vertical direction, one digit is a sequence number in the horizontal direction, and hundred digits are sequence numbers in the direction perpendicular to the horizontal direction and the vertical direction.
Preferably, the grid node and grid sequence number generating function expression is as follows: f (x) = x, F (y) = y, F (z) = z, x represents a horizontal direction number, y represents a vertical direction number, and z represents a vertical direction number perpendicular to the horizontal and vertical directions.
Preferably, the processing mode of the coincident grid nodes for the zero-thickness eight-node coherent unit is that the coincident nodes are formed by splitting original grid nodes, and the splitting direction has four directions, namely an X-axis direction, a Y-axis direction, a 45-degree direction and a Z-axis direction;
the edge region of the coherent unit adopts a four-node splitting mode, and the splitting direction is other two directions except the direction of the edge; the intersection point of the edge and the edge adopts an eight-node splitting mode, the splitting directions are the X-axis direction, the Y-axis direction, the Z-axis direction and the 45-degree direction, the other areas adopt a two-node splitting mode, and the splitting direction is the normal direction of the surface of the area.
Preferably, the mesh node data and the mesh data are determined in step (4) by:
assuming that the serial number of any grid node is YXZ, the coordinate is (a X (X-1), a X (Y-1), a X (Z-1)), wherein a is the side length of the grid where the grid node is located;
assuming that any mesh has a mesh serial number of YXZ, Y (X + 1) Z, (Y + 1) XZ, YX (Z + 1), Y (X + 1) (Z + 1), (Y + 1) (X + 1) (Z + 1), and Y + 1) X (Z + 1).
Preferably, the calculation file format is Ls-dyna calculation file format or Abaqus calculation file format.
Compared with the prior art, the invention has the beneficial effects that: the invention realizes the high-efficiency and high-precision modeling of the finite element parametric analysis, solves the problems of difficult modeling and low modeling efficiency and precision of the existing parametric analysis, and can be widely used for the finite element parametric analysis. For some complex shapes, if the shape function is known, node and grid generation can be carried out according to the method. The method can greatly improve the modeling efficiency, and takes Cohesive three-dimensional unit modeling as an example, when large-batch grid modeling (about a million grids) is carried out, the modeling time of the traditional insertion method is 2-3 days, the modeling time of the method is 2-3 minutes, and the modeling efficiency is greatly improved.
Drawings
FIG. 1 is an example of two-dimensional grid data;
FIG. 2 is a two-dimensional node data example;
FIG. 3 illustrates a two-dimensional planar rectangular node and a grid ordering method;
FIG. 4 shows a shell element on the left side and a 0-thickness 4-node coherent element on the right side;
FIG. 5 is a two-dimensional node splitting diagram;
FIG. 6 shows solid elements on the left side and 8-node coherent elements with 0 thickness on the right side;
FIG. 7 is a schematic diagram of a three-dimensional Cohesive unit grid node splitting area;
FIG. 8 is a three-dimensional node splitting diagram;
FIG. 9 is a flow chart of a method;
FIG. 10 models the program interface.
Detailed Description
The invention is further described in detail below with reference to the drawings and examples. By analyzing the finite element calculation file, it can be known that the finite element model mainly consists of grid node data and grid data, and generally, the finite element model structure has a specific rule, such as a cube, a cuboid, or a regular model aggregate. And (3) directly generating finite element node data and grid data by analyzing the structural rule of the finite element model, and finishing the establishment of the finite element model. The invention discloses a modeling method for finite element analysis of a shell mother-of-pearl microstructure, which comprises the following steps as shown in figure 9:
(1) Carrying out finite element mesh division according to the geometric shape of the shell mother-of-pearl microstructure, specifically simulating a soft structure in the shell mother-of-pearl microstructure according to a Cohesive unit, and simulating a hard structure in the shell mother-of-pearl microstructure by using a four-node unit (two-dimensional) or a hexahedral solid unit (three-dimensional);
(2) Determining grid node and grid serial number generation modes in the grid, wherein the grid node and grid serial number generation modes are the same; determining a grid node and a grid sequence number generating function according to the generating mode; and determining the number of generating functions according to the freedom degree of the coordinate system.
When the finite element grid divided in the step (1) is a two-dimensional grid, the grid node and the grid serial number generation mode are both expressed by two digits, wherein the ten digit is a serial number in the vertical direction, and the one digit is a serial number in the horizontal direction; the corresponding grid node and grid sequence number generation function expression is as follows: f (x) = x, F (y) = y, x represents a horizontal direction number, and y represents a vertical direction number.
As in the two-dimensional grid of fig. 3, the decimal number is a node serial number, the tens number of the decimal number is a node serial number in the vertical direction, and the ones number is a node serial number in the horizontal direction; the large number is a grid serial number, the tens digit of the large number is a grid serial number in the vertical direction, and the single digit of the large number is a grid serial number in the horizontal direction; the side length of the grid is a =1, and the node 11 is a zero point;
when the finite element mesh divided in the step (1) is a three-dimensional mesh, the generation modes of mesh nodes and mesh serial numbers are expressed by three digits, wherein the ten digit is a serial number in the vertical direction, the one digit is a serial number in the horizontal direction, and the hundred digit is a serial number in the direction perpendicular to the horizontal direction and the vertical direction; the finite element node and grid sequence number generating function expression is as follows: f (x) = x, F (y) = y, F (z) = z, x represents a horizontal direction number, y represents a vertical direction number, and z represents a vertical direction number perpendicular to the horizontal and vertical directions.
(3) Determining the node serial number and the grid serial number of the shell mother-of-pearl microstructure grid according to the grid node and grid serial number generation function in combination with the finite element grid divided in the step (1);
(4) Determining grid node data and grid data of the shell mother-of-pearl microstructure according to a grid node sequence number generation mode;
the grid data are grid serial numbers and grid node serial numbers under the same grid are shown in figure 1; the grid node data is a grid node serial number and a coordinate value corresponding to the grid node serial number, as shown in fig. 2;
when the divided finite element mesh is a two-dimensional mesh, the node and mesh generation function is F (x) = x, F (y) = y;
according to the grid node serial number, the coordinates of any grid node YX can be determined, wherein the coordinates of any grid node YX are a (X-1, Y-1), and the coordinates of 23 nodes are 1 (3-1, 2-1) = (2, 1);
from the grid serial numbers, the grid node serial numbers can be determined, and for any grid YX, the grid node serial numbers are YX, Y (X + 1), (Y + 1) X, such as 25 grids, the grid node serial numbers are 25, 26, 36 and 35.
The processing mode of the coincident grid nodes aiming at the Cohesive unit is to split the coincident grid nodes into the original grid nodes, wherein the splitting direction comprises three directions, namely an X axis, a Y axis and a 45-degree direction; as shown in fig. 4 and 5, horizontal and vertical lines in the drawings are a coherence unit region with zero thickness and four nodes, a coherence unit splitting direction in the X direction is the X-axis direction, a coherence unit splitting mode in the Y direction is the Y-axis direction, and the coherence unit splitting mode is a two-node splitting mode; the node splitting directions at the intersection of the X and the Y are the X axis, the Y axis and the 45-degree direction, and the four-node splitting mode is adopted.
And (2) when the finite element grid divided in the step (1) is a three-dimensional grid, expressing the grid node and the grid serial number generation mode by three digits, wherein the ten digits are serial numbers in the vertical direction, the ones digits are serial numbers in the horizontal direction, and the hundred digits are serial numbers in the direction perpendicular to the horizontal direction and the vertical direction. The grid node and grid sequence number generating function expression is as follows: f (x) = x, F (y) = y, F (z) = z, x represents a horizontal direction number, y represents a vertical direction number, and z represents a vertical direction number perpendicular to the horizontal and vertical directions.
Aiming at the processing mode of coincident grid nodes of the coherent unit, the coincident nodes are assumed to be formed by splitting original grid nodes, and the splitting direction comprises four directions, namely an X-axis direction, a Y-axis direction, a 45-degree direction and a Z-axis direction; as shown in fig. 6, the right graph in the figure shows a Cohesive cell area with zero thickness and eight nodes; as shown in fig. 7 and 8, the circle portion in fig. 7 is a two-node splitting method, and the splitting direction is the normal direction (X-axis, Y-axis, Z-axis direction) of the surface; the area shown by the edge is in a four-node splitting mode, and the splitting direction of the area is other two directions (an X axis and a Y axis, an X axis and a Z axis, and a Y axis and a Z axis) except the direction of the edge; the region shown by the dots is in an eight-node splitting mode, and the splitting directions are the X-axis, the Y-axis, the Z-axis and the 45-degree direction.
Assuming that the serial number of any grid node is YXZ, the coordinate is (a X (X-1), a X (Y-1), a X (Z-1)), wherein a is the side length of the grid where the grid node is located;
assuming that any mesh has a mesh serial number of YXZ, Y (X + 1) Z, (Y + 1) XZ, YX (Z + 1), Y (X + 1) (Z + 1), (Y + 1) (X + 1) (Z + 1), and Y + 1) X (Z + 1).
(5) And sorting the generated grid node data, the generated grid data, the corresponding material attribute, the load, the boundary condition and the simulation parameter according to a calculation file format required by finite element analysis to complete modeling.
The invention can directly generate grid node data and grid data through a program, introduces a calculation file obtained after modeling into LS-dyna and other calculation software for calculation, combines with the input of corresponding loads to obtain a series of mechanical properties related to the shell mother-of-pearl microstructure, and provides an effective modeling method for the mechanical property research, such as the research on tensile property, compression property, impact resistance and the like.
Example 1
The invention is further illustrated below in conjunction with the Cohesive two-dimensional element modeling example.
As shown in fig. 4, a square structure with a side length of 80mm is taken as an example, and is composed of square shell units with a side length of 4 × 4 mm, and a four-node and zero-thickness Cohesive unit is arranged between squares. The concrete modeling steps are as follows:
1) Determining a grid node and grid serial number calibration mode, namely sequentially generating grid nodes and grid data in X and Y directions;
2) Determining a grid node sequence number generating function; as can be seen from fig. 4 and 5, the Cohesive cell has two coincident mesh nodes, the two-dimensional four-node cell junction has two coincident mesh nodes, and the two-dimensional four-node cell junction has four coincident mesh nodes; the coincident grid nodes may be assumed to be split from the original nodes in three directions, the X-axis, Y-axis, and 45-degree directions, as shown in fig. 5. When the grid node serial number calibration mode is determined, the part of grid node labels need to be specially considered. The mesh node sequence number generation function is:
f (x) = x, F (y) = y, normal mesh node
F (X) = X + M, F (y) = y, X-axis direction split grid node
F (x) = x, F (Y) = Y + M, Y-axis direction split grid node
F (x) = x + M, F (y) = y + M, 45-degree direction split grid node
M is a constant (set according to modeling conditions on the premise of not influencing other grid node labels, and is prevented from being overlapped with other grid node labels);
3) Determining a grid sequence number generating function; because the Cohesive unit exists in the grid, when the grid sequence number generating function is determined, the sequence number of the Cohesive unit needs to be specially considered. The grid order generation function is:
f (x) = x, F (y) = y, case unit
F (X) = X + M, F (y) = y, X-direction coherent unit
F (x) = x, F (Y) = Y + M, Y-direction complex unit
M is a constant (set according to modeling conditions and prevented from being overlapped with other grid labels on the premise of not influencing other unit labels);
for example, for the grid node data of the splitting part, the coordinates can be determined according to the grid node serial number, and any grid node Y (X + M) has the coordinates of a · (X-1, Y-1), such as 2 (3 + 6) node, and has the coordinates of 1 · (3-1, 2-1) = (2, 1);
for the grid data at the splitting position, the grid node sequence number can be determined according to the grid sequence number, for any grid Y (X + M), the grid node sequence numbers are YX, Y (X +1+ M), (Y + 1) (X +1+ M), and (Y + 1) X, such as 2 (3 + 6) grid, the grid node sequence numbers are 23, 2 (3 +1+ 6), (2 +1+ 6) 3, and (2 + 1) 3, namely 23, 210, 93, 33.
4) Calculating a file format through Ls-dyna, and directly outputting grid node data and grid data;
5) Arranging and outputting corresponding material attributes, loads, boundary conditions, simulation parameters and the like according to a finite element analysis file format according to the Ls-dyna calculation file format;
6) Finally outputting a calculation file to complete finite element modeling;
7) Modeling program (written in Python language), fig. 10 is a modeling program interface, and the brick structure in the interface corresponds to the hard structure of the shell mother-of-pearl microstructure.
Example 2
The invention is further explained by combining the Cohesive three-dimensional unit modeling embodiment.
As shown in fig. 8, a 80mmx80mmx20mm rectangular parallelepiped structure is taken as an example, and is composed of a 20mm × 20mmx2mm rectangular parallelepiped with a side length, a 5mmx5mmx5mm solid unit with a grid size, and eight-node and zero-thickness Cohesive units are arranged between squares.
The three-dimensional model in the modeling step needs to add one more degree of freedom, and compared with the two-dimensional model generation method, the method is mainly different in the steps 2) and 3), and the three-dimensional model is more complex.
Determining a node sequence number generating function; as can be seen from fig. 6, the three-dimensional model is complex, and the generated area is divided into four parts, namely, the uppermost layer, the lowermost layer, the middle odd-numbered layer and the middle even-numbered layer; each layer can be divided into a two-node splitting area, a four-node splitting area and an eight-node splitting area, wherein the two-node splitting area is provided with two coincident nodes, the four-node splitting area is provided with four coincident nodes, and the eight-node splitting area is provided with eight coincident nodes; the coincident nodes can be assumed to be formed by splitting the original nodes, and the splitting directions include four directions, namely X-axis, Y-axis, 45-degree and Z-axis directions. A schematic diagram of a split node is shown in fig. 8. The node sequence number generation function is:
f (x) = x, F (y) = y, F (z) = z, normal node
F (X) = X + M, F (y) = y, F (z) = z, X-axis direction split node
F (x) = x, F (Y) = Y + M, F (z) = z, Y-axis direction split node
F (x) = x + M, F (y) = y + M, F (z) = z,45 degree direction split node
F (x) = x, F (y) = y, F (Z) = Z + M, Z-axis direction split node
M is a constant (set according to modeling conditions and prevented from being overlapped with other grid labels on the premise of not influencing other node labels);
3) Determining a grid sequence number generating function; because the grid has the coherent unit, when the grid sequence number generating function is determined, the sequence number of the coherent unit needs to be specially considered. The grid sequence number generation function is:
f (x) = x, F (y) = y, F (z) = z, physical unit
F (X) = X + M, F (y) = y, F (z) = z, X-direction coherent unit
F (x) = x, F (Y) = Y + M, F (z) = z, Y-direction coherent cell
F (x) = x, F (y) = y, F (Z) = Z + M, cohesive cell in Z direction
M is a constant (set according to the modeling condition without affecting other unit labels)
For example, for the grid node data of the splitting part, the coordinates can be determined according to the grid node number, and any grid node ZY (X + M) has the coordinates of a · (X-1, Y-1, Z-1), such as 12 (3 + 6) node, and has the coordinates of 1 · (3-1, 2-1, 1-1) = (2, 1, 0);
for the grid data of the splitting position, the grid node number can be determined according to the grid number, for any grid ZY (X + M), the grid node number is ZYX, ZY (X +1+ M), Z (Y + 1) X, (Z + 1) YX, (Z + 1) Y (X +1+ M), (Z + 1) (Y + 1) X is 12 (3 + 6) grid, the grid nodes are numbered 123, 12 (3 +1+ 6), 1 (2 +1+ 6) 3, 1 (2 + 1) 3, (1 + 1) 23, (1 + 1) 2 (3 +1+ 6), (1 + 1) (2 +1+ 6) 3, (1 + 1) (2 + 1) 3, namely 123, 1210, 193, 133, 223, 2210, 293, 233.
The remaining steps were processed as in example 1.
The invention has not been described in detail in part in the common general knowledge of a person skilled in the art.
Claims (10)
1. A modeling method for finite element analysis of a shell mother-of-pearl microstructure is characterized by comprising the following steps:
(1) Carrying out finite element meshing according to the geometric shape of the shell mother-of-pearl microstructure, specifically simulating a soft structure in the shell mother-of-pearl microstructure according to a zero-thickness Cohesive unit, and simulating a hard structure in the shell mother-of-pearl microstructure by using a four-node unit or a hexahedral solid unit;
(2) Determining grid node serial numbers in grids and grid serial number generation modes, wherein the grid node serial numbers and the grid serial number generation modes are the same; determining grid node serial numbers and grid serial number generation functions according to the generation mode;
(3) Determining the node serial numbers and the grid serial numbers of the shell pearl shell microstructure grid according to the grid node serial numbers and the grid serial number generation function in combination with the finite element grid divided in the step (1);
(4) Determining grid node data and grid data of the shell mother-of-pearl microstructure according to the grid node serial number and the grid serial number generation mode;
the grid data are grid serial numbers and grid node serial numbers under the same grid; the grid node data is a grid node serial number and a coordinate value corresponding to the grid node serial number;
(5) And sorting the generated grid node data, the generated grid data, the corresponding material attribute, the load, the boundary condition and the simulation parameter according to a calculation file format required by finite element analysis to complete modeling.
2. A modeling method as claimed in claim 1, characterized in that: and (2) when the finite element grid divided in the step (1) is a two-dimensional grid, expressing the grid node serial number and the grid serial number generation mode by two digits, wherein the ten digit is a serial number in the vertical direction, and the one digit is a serial number in the horizontal direction.
3. A modeling method as claimed in claim 2, characterized in that: the grid node sequence number and the grid sequence number generation function expression are as follows: f (x) = x, F (y) = y, x represents a horizontal direction number, and y represents a vertical direction number.
4. A modeling method in accordance with claim 3, characterized by: the processing mode of the coincident grid nodes aiming at the zero-thickness four-node coherent unit is to split the coincident grid nodes into three directions, namely an X axis direction, a Y axis direction and a 45-degree direction, from the original grid nodes; the splitting direction of the coherent unit in the X direction is the X-axis direction, the splitting mode of the coherent unit in the Y direction is the Y-axis direction, and the splitting mode is a two-node splitting mode; the node splitting directions at the intersection of the X and the Y are the X axis, the Y axis and the 45-degree direction, and the four-node splitting mode is adopted.
5. A modeling method as claimed in claim 3, characterized in that: in the step (4), the grid node data and the grid data are determined in the following way:
assuming that the serial number of any grid node is YX, the coordinates are (a X (X-1) and a X (Y-1)), wherein a is the side length of the grid where the grid node is located;
assume that any grid node has a grid serial number of YX, and grid node serial numbers of YX, Y (X + 1), (Y + 1) (X + 1), and (Y + 1) X.
6. The modeling method of claim 1, wherein: and (2) when the finite element grid divided in the step (1) is a three-dimensional grid, expressing the grid node and the grid serial number generation mode by three digits, wherein the ten digits are serial numbers in the vertical direction, the ones digits are serial numbers in the horizontal direction, and the hundred digits are serial numbers in the direction perpendicular to the horizontal direction and the vertical direction.
7. The modeling method of claim 6, wherein: the grid node and grid sequence number generating function expression is as follows: f (x) = x, F (y) = y, F (z) = z, x represents a horizontal direction number, y represents a vertical direction number, and z represents a vertical direction number perpendicular to the horizontal and vertical directions.
8. The modeling method of claim 7, wherein: the processing mode of the coincident grid nodes aiming at the zero-thickness eight-node coherent unit is that the coincident nodes are formed by splitting original grid nodes, and the splitting direction has four directions, namely X-axis, Y-axis, 45-degree and Z-axis directions;
the edge region of the coherent unit adopts a four-node splitting mode, and the splitting direction is other two directions except the direction of the edge; the intersection point of the edge and the edge adopts an eight-node splitting mode, the splitting directions are the X-axis direction, the Y-axis direction, the Z-axis direction and the 45-degree direction, the other areas adopt a two-node splitting mode, and the splitting direction is the normal direction of the surface of the area.
9. The modeling method of claim 7, wherein: in the step (4), the grid node data and the grid data are determined in the following way:
assuming that the serial number of any grid node is YXZ, the coordinate is (a X (X-1), a X (Y-1), a X (Z-1)), wherein a is the side length of the grid where the grid node is located;
assuming that any mesh number is YXZ, the mesh node numbers thereof are YXZ, Y (X + 1) Z, (Y + 1) XZ, YX (Z + 1), Y (X + 1) (Z + 1), (Y + 1) (X + 1) (Z + 1), and Y + 1) X (Z + 1).
10. The modeling method of claim 1, wherein: the calculation file format is Ls-dyna calculation file format or Abaqus calculation file format.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011565362.7A CN112685936B (en) | 2020-12-25 | 2020-12-25 | Modeling method for shell mother-of-pearl microstructure finite element analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011565362.7A CN112685936B (en) | 2020-12-25 | 2020-12-25 | Modeling method for shell mother-of-pearl microstructure finite element analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112685936A CN112685936A (en) | 2021-04-20 |
| CN112685936B true CN112685936B (en) | 2022-10-28 |
Family
ID=75453356
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011565362.7A Active CN112685936B (en) | 2020-12-25 | 2020-12-25 | Modeling method for shell mother-of-pearl microstructure finite element analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112685936B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113283146B (en) * | 2021-06-09 | 2022-12-30 | 深圳大学 | Two-dimensional coherent unit global embedding method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104281726A (en) * | 2014-05-27 | 2015-01-14 | 北京宇航系统工程研究所 | Unfolding finite element modeling method for complex cylinder structural plane |
| CN104679955A (en) * | 2015-02-15 | 2015-06-03 | 北京宇航系统工程研究所 | Method for carrying out finite element parametric modeling on triangular grid stiffener cylinder structure |
| CN107491616A (en) * | 2017-08-24 | 2017-12-19 | 北京航空航天大学 | A kind of structure finite element parametric modeling method suitable for grid configuration rudder face |
| CN108959786A (en) * | 2018-07-12 | 2018-12-07 | 华东交通大学 | A kind of parametric modeling method of imitative mother-of-pearl laminated composite plate |
-
2020
- 2020-12-25 CN CN202011565362.7A patent/CN112685936B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104281726A (en) * | 2014-05-27 | 2015-01-14 | 北京宇航系统工程研究所 | Unfolding finite element modeling method for complex cylinder structural plane |
| CN104679955A (en) * | 2015-02-15 | 2015-06-03 | 北京宇航系统工程研究所 | Method for carrying out finite element parametric modeling on triangular grid stiffener cylinder structure |
| CN107491616A (en) * | 2017-08-24 | 2017-12-19 | 北京航空航天大学 | A kind of structure finite element parametric modeling method suitable for grid configuration rudder face |
| CN108959786A (en) * | 2018-07-12 | 2018-12-07 | 华东交通大学 | A kind of parametric modeling method of imitative mother-of-pearl laminated composite plate |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112685936A (en) | 2021-04-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Pandolfi et al. | An efficient adaptive procedure for three-dimensional fragmentation simulations | |
| CN107562976B (en) | Finite element modeling and analysis of crack propagation in multiple planes of a structure | |
| CN110399684B (en) | Data extraction method, system and computer readable storage medium for ABAQUS crack propagation | |
| CN113297650B (en) | BIM technology-based unit type glass curtain wall construction method and system | |
| CN112685936B (en) | Modeling method for shell mother-of-pearl microstructure finite element analysis | |
| Jaworski et al. | Implementation features of composite materials effective mechanical characteristics finding method based on microlevel cellular structural models | |
| CN103106313B (en) | Roll consequent order reconstructing method | |
| CN112231863B (en) | Solar wing cell array substrate modeling method, device, equipment and storage medium | |
| CN107145630B (en) | Plate shell structure design and analysis integrated method based on CAD (computer-aided design) cutting curved surface | |
| CN110321571A (en) | A kind of mechanics parameter numerical value extracting method of honeycomb plate and shell structure | |
| CN115630542B (en) | Reinforced layout optimization method of thin-wall reinforced structure | |
| RU2689803C1 (en) | Method for automatized construction of three-dimensional model of heterogeneous structure of composite material with fibers | |
| EP3432171B1 (en) | Analysis model creation assistance device and analysis model creation assistance method | |
| CN113722965B (en) | Fracture simulation method based on integral-generalized finite difference numerical discrete operator | |
| CN105243238A (en) | Integrated quick product iteration forming device and method | |
| CN107832537B (en) | Numerical simulation method for residual stress in complex form | |
| CN111859268B (en) | Magnetic tensor abnormal spatial domain fast forward modeling method based on grid point lattice | |
| KR101327964B1 (en) | Fem modeling method for vibration analysis of commercial vessel in shipbuilding industry | |
| CN117010094B (en) | Single-machine batch modeling method and system for satellite structure finite element model | |
| CN106251401A (en) | A kind of generation system and method for hexahedron 20 node unit grid | |
| CN118070622B (en) | Rapid modeling method and system for realizing unstructured grid stretching based on FLAC software | |
| CN110868325B (en) | Uniform grid division method capable of reducing rigidity matrix construction difficulty | |
| CN116090113B (en) | Variable node suspension unit grid finite element model of block structure | |
| JP2005259043A (en) | Three-dimensional mesh generation method | |
| CN117240717B (en) | Method for processing grid of slicing structure, optimizing method, storage medium and equipment |
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 |