Analysis apparatus and analysis method

Abstract

To reduce calculation cost in finite element analysis.SOLUTION: An analyzer comprises: a specification unit that, based on an analysis result of finite element analysis using a first model (model M1) including a first base part 21 and a specific shape part (cusp part 31), specifies a first stiffness matrix that is the stiffness matrix of the first model; a derivation unit that subtracts a second stiffness matrix that is the stiffness matrix of the first base part 21 from the first stiffness matrix to derive a third stiffness matrix that is the stiffness matrix of the specific shape part; and a creation unit that creates a second model in which a shell element having the third stiffness matrix set as its stiffness matrix is attached to a second base part.SELECTED DRAWING: Figure 3

Description

[Technical Field]

[0001] This disclosure relates to an analytical device and an analytical method. [Background technology]

[0002] In product development, there is a need to predict the stress and deformation behavior that occurs in components. As a method for predicting the stress and deformation behavior that occurs in components, finite element analysis, which is an analysis using the finite element method, is used, for example, as disclosed in Patent Document 1. [Prior art documents] [Patent Documents]

[0003] [Patent Document 1] Japanese Patent Publication No. 2019-160198 [Overview of the project] [Problems that the invention aims to solve]

[0004] Finite element analysis requires dividing a model into multiple meshes (i.e., multiple elements). As the model's geometry becomes more complex, the number of meshes increases exponentially, leading to higher computational costs. Therefore, reducing the computational cost in finite element analysis is highly desirable.

[0005] The purpose of this disclosure is to provide an analysis device and analysis method that can reduce the computational cost in finite element analysis. [Means for solving the problem]

[0006] In order to solve the above problems, the analysis device of the present disclosure includes a specifying unit that specifies a first stiffness matrix, which is the stiffness matrix of the first model, based on the analysis result of finite element analysis using a first model including a first base part and a specific shape part, a deriving unit that subtracts a second stiffness matrix, which is the stiffness matrix of the first base part, from the first stiffness matrix to derive a third stiffness matrix, which is the stiffness matrix of the specific shape part, and a generating unit that generates a second model in which shell elements with the third stiffness matrix set as the stiffness matrix are attached to the second base part.

[0007] The specifying unit may specify the second stiffness matrix based on the analysis result of finite element analysis using a third model obtained by removing the specific shape part from the first model.

[0008] The first base part may have the basic shape of the workpiece, and the specific shape part may have a shape more complex than the first base part.

[0009] The specific shape part may have the shape of a machining mark formed on the workpiece.

[0010] The stiffness matrix may include a matrix that defines the relationship between in-plane stress and in-plane strain.

[0011] The stiffness matrix may include a matrix that defines the relationship between moment and curvature.

[0012] The stiffness matrix may include a matrix that defines the relationship between out-of-plane shear stress and out-of-plane shear strain.

[0013] In order to solve the above problems, the analysis method of the present disclosure includes a step of specifying a first stiffness matrix, which is the stiffness matrix of the first model, based on the analysis result of finite element analysis using a first model including a first base part and a specific shape part, a step of subtracting a second stiffness matrix, which is the stiffness matrix of the first base part, from the first stiffness matrix to derive a third stiffness matrix, which is the stiffness matrix of the specific shape part, and a step of generating a second model in which shell elements with the third stiffness matrix set as the stiffness matrix are attached to the second base part.

Advantages of the Invention

[0014] According to this disclosure, the computational cost in finite element analysis can be reduced. [Brief explanation of the drawing]

[0015] [Figure 1] Figure 1 is a block diagram showing an example of the functional configuration of an analysis device according to an embodiment of this disclosure. [Figure 2] Figure 2 is a flowchart showing an example of the processing flow performed by the analysis device according to the embodiment of this disclosure. [Figure 3] Figure 3 is a schematic diagram showing a model used in the processing performed by the analysis apparatus according to the embodiment of this disclosure, the model including a first base portion and a cusp portion. [Figure 4] Figure 4 is a schematic diagram showing a model used in the processing performed by the analysis apparatus according to the embodiment of this disclosure, which is the model obtained by removing the cusp portion from the model in Figure 3. [Figure 5] Figure 5 is a schematic diagram showing a model used in the processing performed by the analysis apparatus according to the embodiment of this disclosure, in which a shell element, whose third stiffness matrix is ​​set as a stiffness matrix, is attached to the second base portion. [Modes for carrying out the invention]

[0016] Embodiments of this disclosure will be described below with reference to the attached drawings. The dimensions, materials, and other specific numerical values ​​shown in the embodiments are merely examples for the purpose of facilitating understanding and do not limit this disclosure unless otherwise specified. In this specification and drawings, elements having substantially the same function or configuration are denoted by the same reference numerals to avoid redundant explanations, and elements not directly related to this disclosure are omitted from the illustrations.

[0017] Figure 1 is a block diagram showing an example of the functional configuration of the analysis device 1 according to this embodiment. The analysis device 1 performs finite element analysis to predict the stress and deformation behavior occurring in a component. The analysis device 1 includes, for example, a central processing unit (CPU), a ROM storing programs, and RAM as a work area. As shown in Figure 1, the analysis device 1 includes, for example, a generation unit 11, an analysis unit 12, a specification unit 13, and a derivation unit 14. The functions of the analysis device 1 described below may be realized by a single device or may be divided among multiple devices.

[0018] The generation unit 11 generates a model that reproduces the shape of the part to be analyzed using finite element analysis. For example, the generation unit 11 generates a model in response to user operations performed on the analysis device 1. Models used in finite element analysis are divided into multiple meshes (i.e., multiple elements). The generation of meshes in the model is also performed by the generation unit 11.

[0019] The analysis unit 12 performs finite element analysis using the model generated by the generation unit 11. In the finite element analysis performed by the analysis unit 12, governing equations relating to mechanics are set for each element of the model, and a solution that satisfies all of the set governing equations is sought. Specifically, in the finite element analysis performed by the analysis unit 12, the governing equations set for each element are equations that define the relationship between deformation and stress.

[0020] For example, equation (1) below is used as the governing equation set for each element. Equation (1) is an equation that defines the relationship between in-plane stress and in-plane strain, as well as the relationship between moment, curvature and torsion. In the following, we will explain assuming that the first direction and the second direction are orthogonal to each other in the plane of the calculation target. Furthermore, the direction orthogonal to the first and second directions will be explained as the third direction. Note that in Figures 3 to 5 described later, the X direction and Y direction correspond to the first direction and the second direction, respectively, and the Z direction corresponds to the third direction.

[0021]

number

[0022] N in Equation (1) 11 , N 22 , N 12 is the in-plane stress. N 11 is the normal stress in the first direction. N 22 is the normal stress in the second direction. N 12 is the shear stress on the plane parallel to the first and second directions. M in Equation (1) 11 , M 22 , M 12 is the moment. M 11 is the bending moment about the axis in the first direction. M 22 is the bending moment about the axis in the second direction. M 12 is the torsional moment.

[0023] ε in Equation (1) 11 , ε 22 , ε 12 is the in-plane strain. ε 11 is the normal strain in the first direction. ε 22 is the normal strain in the second direction. ε 12 is the shear strain on the plane parallel to the first and second directions. κ in Equation (1) 11 , κ 22 , κ 12 are the curvature and the twist rate. κ 11 is the curvature in the first direction. κ 22 is the curvature in the second direction. κ 12 is the twist rate.

[0024] K in Equation (l) is the stiffness matrix that defines the relationship between the in-plane stress and the in-plane strain, as well as the relationship between the moment and the curvature and the twist rate. The stiffness matrix K is represented by the following Equation (2).

[0025]

Equation

[0026] In equation (2), A is the in-plane stiffness matrix. In equation (2), B is the coupling stiffness matrix. In equation (2), D is the bending stiffness matrix. The in-plane stiffness matrix A gives N in equation (1). 11 , N 22 , N 12 and ε 11 , ε 22 , ε 12 The relationship between the in-plane stress and the in-plane strain is defined by the bending stiffness matrix D, which determines the relationship between M in equation (1). 11 M 22 M 12 and κ 11 κ 22 κ 12 The relationship between moment, curvature, and torsion is defined.

[0027] For example, in addition to equation (1) above, equation (3) below is also used as a governing equation set for each element. Equation (3) is an equation that defines the relationship between out-of-plane shear stress and out-of-plane shear strain.

[0028]

number

[0029] S in equation (3) 13 S 23 This is the out-of-plane shear stress. 13 S is the shear stress in the plane parallel to the first and third directions. 23 γ in equation (3) is the shear stress in the plane parallel to the second and third directions. 13 gamma 23 This is the out-of-plane shear strain. γ 13 This is the shear strain in the planes parallel to the first and third directions. γ 23 This is the shear strain in the planes parallel to the second and third directions. In equation (3), E is the stiffness matrix that defines the relationship between out-of-plane shear stress and out-of-plane shear strain. The stiffness matrix E gives S in equation (3). 13 S 23 and gamma 13 gamma23 The relationship between out-of-plane shear stress and out-of-plane shear strain is defined.

[0030] The following describes an example in which equations (1) and (3) are used as governing equations in the finite element analysis performed by the analysis unit 12. However, as will be discussed later, the equations used as governing equations are not limited to this example.

[0031] The identification unit 13 identifies the stiffness matrix of the model used in the finite element analysis based on the analysis results of the finite element analysis. The derivation unit 14 uses the information obtained by the identification unit 13 to derive the stiffness matrix of a specific shape part (for example, the cusp part 31 in Figure 3, which will be described later). In the analysis device 1, the computational cost in finite element analysis is reduced mainly by the processing performed by the identification unit 13 and the derivation unit 14.

[0032] Figure 2 is a flowchart showing an example of the processing flow performed by the analysis device 1 according to this embodiment. The processing flow shown in Figure 2 is executed, for example, when a predetermined operation is performed on the analysis device 1 by the user.

[0033] The following describes examples in which Model M1 shown in Figure 3, Model M2 shown in Figure 4, and Model M3 shown in Figure 5 are used in the processing performed by the analysis device 1. However, the models used in the processing performed by the analysis device 1 are set appropriately according to the shape of the part to be analyzed, so this example is not exhaustive. As will be described later, Model M1 shown in Figure 3 corresponds to an example of the first model. Model M2 shown in Figure 4 corresponds to an example of the third model. Model M3 shown in Figure 5 corresponds to an example of the second model.

[0034] When the processing flow shown in Figure 2 begins, in step S101, the generation unit 11 creates the model M1 shown in Figure 3. Model M1 is a model that reproduces a part of the component to be analyzed. As shown in Figure 3, model M1 includes a first base part 21 and a cusp part 31.

[0035] The cusp portion 31 corresponds to an example of a specific shape portion having a specific shape formed on the surface of the part to be analyzed. The specific shape is a complex shape compared to a flat, simple shape. For example, the specific shape has at least a recess or a convex portion. The cusp portion 31 has the shape of a machining mark formed on the workpiece. In the example in Figure 3, the cusp portion 31 is a portion in which multiple recesses extending in the X direction are arranged in the Y direction. The first base portion 21 is the portion of the part to be analyzed from which the specific shape portion is formed, excluding the specific shape portion. In other words, the first base portion 21 has a flat, simple shape. Specifically, the first base portion 21 has the basic shape of the workpiece. The specific shape portion has a more complex shape than the first base portion 21. In the example in Figure 3, the first base portion 21 has a flat plate shape. In Figure 3, the boundary between the first base portion 21 and the cusp portion 31 in model M1 is shown by a dashed line.

[0036] As shown in Figure 3, the model M1, when viewed as a whole, is formed in a substantially flat plate shape that extends in the X and Y directions and has thickness in the Z direction. The first base portion 21 and the cusp portion 31 also extend in the X and Y directions, respectively, and have thickness in the Z direction. The first base portion 21 and the cusp portion 31 are adjacent to each other in the Z direction. Note that the X, Y, and Z directions are orthogonal to each other. The model M1 is divided into multiple elements. Specifically, as shown in Figure 3, the model M1 is divided in the X, Y, and Z directions.

[0037] Following step S101 in Figure 2, in step S102, the analysis unit 12 performs a finite element analysis using model M1. Specifically, the analysis unit 12 sets equations (1) and (3) for each element of model M1 and finds a solution that satisfies all the set governing equations. The stiffness matrices K and E set for each element are predetermined. Note that the stiffness matrices K and E may differ between elements in model M1.

[0038] Following step S102, in step S103, the identification unit 13 identifies the first stiffness matrix, which is the stiffness matrix of model M1. Specifically, the identification unit 13 identifies the first stiffness matrix based on the analysis results of a finite element analysis using model M1. The first stiffness matrix is ​​the stiffness matrix K and stiffness matrix E of the entire model M1, assuming that the entire model M1 is a single element.

[0039] Following step S103, in step S104, the generation unit 11 creates the model M2 shown in Figure 4. As shown in Figure 4, model M2 is a model obtained by removing the cusp portion 31 from model M1 in Figure 3. In other words, model M2 is a model that reproduces only the first base portion 21 which has a flat and simple shape. As shown in Figure 4, model M2 is divided in the X, Y, and Z directions. In the example in Figure 4, the meshing method of the first base portion 21 in model M2 is different from the meshing method of the first base portion 21 in model M1 in Figure 3. However, the meshing method of the first base portion 21 in model M2 may be the same as the meshing method of the first base portion 21 in model M1.

[0040] Following step S104 in Figure 2, in step S105, the analysis unit 12 performs a finite element analysis using model M2. Specifically, the analysis unit 12 sets equations (1) and (3) for each element of model M2 and finds a solution that satisfies all the set governing equations. The stiffness matrices K and E set for each element are predetermined. Note that the stiffness matrices K and E may differ among elements in model M2.

[0041] Following step S105, in step S106, the identification unit 13 identifies the second stiffness matrix, which is the stiffness matrix of model M2. Since model M2 is a model that reproduces only the first base part 21, the second stiffness matrix is ​​the stiffness matrix of the first base part 21. Specifically, the identification unit 13 identifies the second stiffness matrix based on the analysis results of a finite element analysis using model M2. The second stiffness matrix is ​​the stiffness matrix K and stiffness matrix E of the entire first base part 21 when the entire first base part 21 is assumed to be a single element.

[0042] Following step S106, in step S107, the derivation unit 14 derives a third stiffness matrix, which is the stiffness matrix of the cusp portion 31. Specifically, the derivation unit 14 derives the third stiffness matrix by subtracting the second stiffness matrix specified in step S106 from the first stiffness matrix specified in step S103. The third stiffness matrix is ​​the stiffness matrix K and stiffness matrix E of the entire cusp portion 31, assuming the entire cusp portion 31 as a single element.

[0043] For example, the derivation unit 14 derives the stiffness matrix obtained by subtracting the stiffness matrix K of the first base portion 21, which was identified as the second stiffness matrix in step S106, from the stiffness matrix K of model M1, which was identified as the first stiffness matrix in step S103, as the stiffness matrix K of the cusp portion 31 among the third stiffness matrices. For example, the derivation unit 14 derives the stiffness matrix obtained by subtracting the stiffness matrix E of the first base portion 21, which was identified as the second stiffness matrix in step S106, from the stiffness matrix E of model M1, which was identified as the first stiffness matrix in step S103, as the stiffness matrix E of the cusp portion 31 among the third stiffness matrices.

[0044] Following step S107, in step S108, the generation unit 11 creates the model M3 shown in Figure 5. Model M3 is a model that reproduces a part of the component under analysis that is different from the part reproduced by model M1. On the surface of the part of the component under analysis that is reproduced by model M3, a cusp portion 31 is actually formed, similar to the part reproduced by model M1.

[0045] As shown in Figure 5, in model M3, shell elements 41 are attached to the surface of the second base portion 22. Shell elements 41 are a type of element used in finite element analysis and are planar elements. Shell elements 41 are created on a surface. Although shell elements 41 theoretically have rigidity equal to the plate thickness, they are not divided in the plate thickness direction. Specifically, multiple shell elements 41 are attached to the surface of the second base portion 22 so as to cover the entire surface of the second base portion 22. Note that in Figure 5, for ease of understanding, the second base portion 22 and the shell elements 41 are shown separated. The second base portion 22 has a flat, simple shape, similar to the first base portion 21 of model M1. In the example in Figure 5, the second base portion 22 has a flat plate shape. Note that the shape and dimensions of the second base portion 22 are fundamentally different from those of the first base portion 21, as the parts reproduced by model M1 and the parts reproduced by model M3 are different.

[0046] The shell element 41 is a shell element in which the third stiffness matrix derived in step S107 is set as the stiffness matrix. In other words, the shell element 41 is an element used in place of a solid element that reproduces the shape of the cusp portion 31. Multiple such shell elements 41 are attached to the surface of the second base portion 22 so as to cover the entire surface of the second base portion 22. In other words, model M3 is a model in which shell elements 41 in which the third stiffness matrix is ​​set as the stiffness matrix are attached to the second base portion 22.

[0047] Following step S108 in Figure 2, in step S109, the analysis unit 12 performs a finite element analysis using model M3, and the processing flow shown in Figure 2 is completed. Specifically, the analysis unit 12 sets equations (1) and (3) for each element of model M3 and finds a solution that satisfies all the set governing equations. The stiffness matrices K and E set for each element of the second base unit 22 in model M3 are predetermined. The stiffness matrices K and E derived as the third stiffness matrices in step S107 are set for each shell element 41 in model M3.

[0048] As explained above, in the analysis device 1, the identification unit 13 identifies the first stiffness matrix, which is the stiffness matrix of the first model, based on the analysis results of a finite element analysis using a first model (model M1 in the above example) that includes a first base unit 21 and a specific shape part (cusp unit 31 in the above example). The derivation unit 14 derives the third stiffness matrix, which is the stiffness matrix of the specific shape part, by subtracting the second stiffness matrix, which is the stiffness matrix of the first base unit 21, from the first stiffness matrix. The generation unit 11 generates a second model (model M3 in the above example) by attaching shell elements 41, for which the third stiffness matrix is ​​set as the stiffness matrix, to the second base unit 22. As a result, when performing finite element analysis on a part of the component to be analyzed in which a specific shape part is formed on the surface, finite element analysis can be performed without using a model that reproduces the shape of the specific shape part. In a model that reproduces the shape of the specific shape part (for example, model M1 above), the number of meshes becomes enormous as the complexity of the model shape increases. On the other hand, in the analysis device 1, by using shell elements 41, finite element analysis can be performed without using a model that reproduces the shape of a specific part, thus suppressing the increase in the number of meshes that occurs as the shape of the model becomes more complex. Therefore, the computational cost in finite element analysis can be reduced.

[0049] The shell elements 41 generated by the processing flow shown in Figure 2 are used each time finite element analysis is performed on other parts of the component being analyzed where the cusp portion 31 is formed on the surface, after the processing flow shown in Figure 2 is completed. Specifically, in the finite element analysis of other parts of the component being analyzed where the cusp portion 31 is formed on the surface, the generation unit 11 generates a model by attaching the shell elements 41 generated by the above process to the base portion corresponding to the other part. Then, the analysis unit 12 performs finite element analysis using the obtained model. This makes it possible to sequentially perform finite element analysis on multiple parts of the component being analyzed where the cusp portion 31 is formed on the surface, while reducing computational costs.

[0050] In particular, the specific unit 13 identifies the second stiffness matrix based on the analysis results of a finite element analysis using a third model (model M2 in the above example) obtained by removing a specific shape part (cusp part 31 in the above example) from the first model (model M1 in the above example). This allows the stiffness matrix of the entire first base part 21, assuming the entire first base part 21 as a single element, to be appropriately identified as the second stiffness matrix. However, the specific unit 13 may identify the second stiffness matrix without relying on the analysis results of a finite element analysis. For example, if the second stiffness matrix is ​​known, the specific unit 13 can identify the second stiffness matrix without performing a finite element analysis. Since the target of the calculation of the second stiffness matrix is ​​a simple plate shape, the second stiffness matrix may be obtained by a manual calculation based on material mechanics, which is the mechanics of elastic bodies. For example, the specific unit 13 may identify the second stiffness matrix by a calculation other than finite element analysis. In these cases, steps S104 and S105 may be omitted.

[0051] In particular, the first base portion 21 has the basic shape of the workpiece, while the specific shape portion (cusp portion 31 in the above example) has a more complex shape than the first base portion 21. This reduces the computational cost in finite element analysis when the workpiece is the object of analysis. However, the object of analysis is not limited to the workpiece.

[0052] In particular, the specific shape is a cusp portion 31 having the shape of a machining mark formed on the workpiece. This makes it possible to suppress the increase in the number of meshes due to the complexity of the model shape when performing finite element analysis on the portion of the part to be analyzed in which the cusp portion 31 is formed on the surface. Therefore, the computational cost in finite element analysis of a part in which the cusp portion 31 is formed on the surface can be reduced. However, the specific shape is not limited to a cusp portion 31 and may have various shapes.

[0053] In particular, the stiffness matrix includes a matrix that defines the relationship between in-plane stress and in-plane strain (stiffness matrix K in equation (1) in the example above). This reduces the computational cost in finite element analysis for predicting in-plane stress and strain occurring in a part.

[0054] In particular, the stiffness matrix includes a matrix that defines the relationship between moment, curvature, and torsion (stiffness matrix K in equation (1) in the example above). This reduces the computational cost in finite element analysis for predicting the moment, curvature, and torsion occurring in a component.

[0055] In particular, the stiffness matrix includes a matrix that defines the relationship between out-of-plane shear stress and out-of-plane shear strain (stiffness matrix E in equation (3) in the example above). This reduces the computational cost in finite element analysis for predicting out-of-plane shear stress and out-of-plane shear strain occurring in a component.

[0056] While embodiments of this disclosure have been described above with reference to the attached drawings, it goes without saying that this disclosure is not limited to such embodiments. It will be obvious to those skilled in the art that various modifications or alterations can be conceived within the scope of the claims, and these will naturally also fall within the technical scope of this disclosure.

[0057] The above describes an example in which equations (1) and (3) are used as governing equations in the finite element analysis performed by the analysis unit 12. However, the equations used as governing equations are not limited to the above example. For example, equation (3) may not be used, and only equation (1) may be used as the governing equation. For example, instead of equation (1), the in-plane stress N may be used as the governing equation. 11 , N 22 , N 12 and in-plane distortion ε 11 , ε 22 , ε 12 An equation that defines only the relationship with [the given element] may also be used.

[0058] This disclosure can contribute, for example, to Sustainable Development Goal (SDG) 12, "Ensure sustainable consumption and production patterns." [Explanation of symbols]

[0059] 1 Analysis device 11 Generation part 13 Specific section 14 Derivation part 21. First base section 22. Second base section 31 Cusp part (particularly shaped part) 41 Shell elements M1 Model (First Model) M2 Model (3rd Model) M3 Model (2nd Model)

Claims

1. A specific unit identifies a first stiffness matrix, which is the stiffness matrix of the first model, based on the analysis results of a finite element analysis using a first model including a first base unit and a specific shape unit. A derivation unit that derives a third stiffness matrix, which is the stiffness matrix of the specific shape portion, by subtracting the second stiffness matrix, which is the stiffness matrix of the first base portion, from the first stiffness matrix, A generation unit that generates a second model in which the shell elements, set as the third stiffness matrix, are attached to the second base portion, Equipped with, Analysis device.

2. The specified part determines the second stiffness matrix based on the analysis results of a finite element analysis using a third model obtained by removing the specified shape part from the first model. The analysis apparatus according to claim 1.

3. The first base portion has the basic shape of the workpiece, The aforementioned specific shaped portion has a more complex shape than the first base portion. The analysis apparatus according to claim 1 or 2.

4. The aforementioned specific shaped portion has the shape of a processing mark formed on the workpiece. The analysis apparatus according to claim 3.

5. The stiffness matrix includes a matrix that defines the relationship between in-plane stress and in-plane strain. The analysis apparatus according to claim 1 or 2.

6. The stiffness matrix includes a matrix that defines the relationship between moment, curvature, and torsion. The analysis apparatus according to claim 1 or 2.

7. The stiffness matrix includes a matrix that defines the relationship between out-of-plane shear stress and out-of-plane shear strain. The analysis apparatus according to claim 1 or 2.

8. A step of identifying a first stiffness matrix, which is the stiffness matrix of the first model, based on the analysis results of a finite element analysis using a first model including a first base portion and a specific shape portion. The steps include: subtracting the second stiffness matrix, which is the stiffness matrix of the first base portion, from the first stiffness matrix to derive the third stiffness matrix, which is the stiffness matrix of the specific shaped portion; The steps include generating a second model by attaching shell elements, for which the third stiffness matrix is ​​set as a stiffness matrix, to the second base portion, including, Analysis method.

Images