JP5437879B2 - Stress strain analysis apparatus and stress strain analysis method - Google Patents

Description

The present invention relates to a stress strain analysis apparatus and a stress strain analysis method using a finite element method.

Conventionally, a finite element method is known as an analysis method for an analysis object having a complicated shape or boundary condition. In this finite element method, an analysis target is divided into a plurality of mesh units, and stress is analyzed for each node of the element, thereby analyzing the stress strain of the analysis target. FIG. 1 shows a state before an analysis process to be analyzed by a conventional finite element method algorithm. As shown in the figure, in order to analyze the stress by the finite element method, the analysis object data D is previously divided into a plurality of meshes M such as a quadrangle or a triangle having a plurality of nodes Y (including Gauss points, the same applies hereinafter). It is necessary to divide it.

FIG. 2 shows the result of analyzing the analysis target data D shown in FIG. 1 by an algorithm of the finite element method. In the finite element method algorithm, stresses at the nodes Y1 to Y28 of the elements E-1 to E-18 are analyzed, and the node Y of the element E is moved to the coordinate position after the analysis. As a result, distortion is increased in a portion where stress is high, and the magnitude of stress strain can be visually grasped according to the state of distortion.

The displacement of the node and the displacement of the distance between the nodes can be regarded as the magnitude of the distortion. In the illustrated example, the deformation of the element E is large at the center of the lower surface of the analysis target data D and is small at the center of the upper surface, the bending works as a whole from the deformation shape of each element E, the center of the lower surface of the analysis target data D It can be seen that the stress is concentrated on the part.

The conventional finite element algorithm incorporates a function of creating a contour line of stress based on the calculated stress and displaying it as a contour diagram. When the contour diagram is displayed, the stress distribution can be grasped visually, which increases convenience. In addition, in the conventional numerical analysis method of the finite element method, when a quadrilateral mesh is used, it is possible to analyze stress with high accuracy, and therefore, a method for creating a uniform quadrilateral mesh has been studied (Patent Document 1).

JP-A-5-298409

By the way, in the development of vehicles such as motorcycles and automobiles, the stress of a plurality of nodes in the designated area specified by the user is averaged to obtain stress strain in the designated area, and this is compared with the actual data measured by the strain gauge. This is useful for development.
However, the position and size of the designated area of the user are arbitrary, and the size of the element depends on the size of the mesh M. Therefore, depending on the position and size of the designated area specified by the user and the size of the mesh M, As shown in FIG. 3, a plurality of nodes Y may be relatively biased in the designated area AL designated by the user. In addition, as shown in FIG. 4, the number of nodes Y may be significantly insufficient, and as a result, the reliability of the average value of stress at the nodes Y may be impaired. Further, since the stress value at the node Y also depends on the calculation accuracy of the stress at each node Y, the calculation accuracy may be affected.
Therefore, it is considered that the analysis target data D is subdivided using the mesh M having a small internode distance corresponding to the user's arbitrary region designation, and the calculation of the stress strain in the designated region designated by the user is executed again. It is done.

However, if the subdivision is performed with the mesh M having a small mesh size and the stress strain is recalculated, the number of calculations becomes enormous due to the increase in the number of nodes. There is a problem that takes time.
In addition, the analysis target data D is divided into a plurality of elements so as to be able to cope with the user's arbitrary area designation AL from the beginning without performing such subdivision and recalculation, and the stress strain of the designated area is calculated. However, the problem associated with the increase in the number of nodes cannot be solved.
As another method, it may be possible to perform recalculation by re-dividing only the user's designated area AL. However, this causes a problem that the relationship with other elements cannot be grasped.

The present invention has been devised in view of the above problems, and calculates stress (stress strain) in an arbitrary designated region designated by the user simply and accurately without recalculating from the finite element method algorithm. It is an object of the present invention to provide a stress strain analysis apparatus and a stress strain analysis method capable of performing the above.

In order to solve the above problems, the invention according to claim 1 of the present application regards the analysis target data (D) as the analysis target as an aggregate joined to each other at a finite number of nodes, and the analysis target data (D) In a stress-strain analyzer that divides into a plurality of meshes (M) having a plurality of nodes and analyzes the stress in each element (E) by an algorithm of the finite element method,
A calculation means (17) for calculating stress values of a plurality of nodes (Y) of the analysis object data (D), and a designated area (AL) designated by the user for all areas of the analysis object data (D). Area designation means (18) arbitrarily designated via the input means (3, 4), and a plurality of reference nodes (P1 to P5) for analyzing stress strain in the designated area (AL) AL) on the outline (AL1) and in the vicinity thereof, the reference node generating means (19) for generating in the vicinity of the node (Y) ;
Calculating a stress value of the reference node generating means (19) has generated said plurality of reference nodes (P1 to P5), calculating means for calculating the stress-strain of the specified area on the basis of these stress values (20) With
The stress strain calculation means (20) of the designated area is
Based respective stress values of a plurality of reference nodes (P1 to P4) located on the contour line (AL1), the known stress value of a node of the ambient adjacent to each reference node (P1~P4) (Y) The
According to the distance between each of the reference nodes (P1 to P4) and the adjacent surrounding nodes (Y), it is calculated and calculated by a weighted average in which the weight is increased as the distance is shorter,
And the stress value of the reference node (P5) inside the designated area (AL) is calculated as the stress value at the known node (Y) in the vicinity thereof,
A stress value of the designated region (AL) is calculated by weighting and averaging the obtained stress values of the plurality of reference nodes (P1 to P5) according to positions in the designated region (AL), and based on the stress values. And calculating the stress strain.

According to a second aspect of the present invention, in the first aspect of the present invention, the shape of the designated area (AL) is a rectangle.

According to a third aspect of the present invention, in the second aspect of the present invention, the designated area (AL) has a rectangular shape, and the reference nodes (P1 to P4) include rectangular vertices.

The invention described in claim 4 is the invention described in any one of claims 1 to 3 , characterized in that the reference node is generated at least five locations.

According to the fifth aspect of the present invention, the analysis target data (D) as an analysis target is regarded as an aggregate joined to each other at a finite number of nodes by the stress strain analysis apparatus, and the analysis target is a mesh having a plurality of nodes (M In the stress strain analysis method, the stress in each element (E) is analyzed by a finite element algorithm algorithm.
The stress strain analyzer is
Calculating the stress value calculation means of a plurality of nodes (Y) (17), the stress value of the plurality of nodes (Y) of the analyzed data (D) in the mesh,
The region designation means (18), comprising the steps of: arbitrarily specify specified area (AL) with respect to the total area of the analysis target-flop with data (D),
By reference node generating means (19), the contour line (AL1) on the plurality of reference nodes (P1 to P5) of the designated area for analyzing the stress-strain of the designated area (AL) (AL) and in its interior Generating in the vicinity of the node (Y) ;
The stress distortion calculating means specified area (AL) (20), known around the node (Y) adjacent the respective stress values of the plurality of reference nodes (P1 to P4) to each reference node (P1 to P4) On the basis of the stress value of each of the reference nodes (P1 to P4) and the adjacent surrounding nodes (Y) according to the distance, the shorter the distance, the greater the weighted average,
Calculating a stress value of the reference node (P5) inside the designated area (AL) as a known stress value at the node (Y) in the vicinity thereof ;
The stress values of the plurality of reference nodes (P1 to P5) obtained are weighted and averaged according to the positions in the designated area (AL) to calculate the stress values in the designated area (AL), and the stress values are calculated based on the stress values. Calculating the stress strain of the designated area (AL).

According to the first and fifth aspects of the present invention, it is possible to accurately analyze the stress (stress strain) in an arbitrary designated region designated by the user while using the algorithm of the stress strain analyzer using the conventional finite element method. There is an effect that it can be performed well. Further, it is not necessary to subdivide and calculate the mesh size, and the processing time can be reduced.
Therefore, a stress strain analysis apparatus and a stress strain analysis that can easily and accurately calculate the stress (stress strain) of an arbitrary designated area specified by the user without recalculating from the finite element algorithm. You can get the method.

Further, since the stress in the designated region is calculated by the weighted average in which the weighting is increased as the distance is shorter according to the distance between each reference node (P1 to P4) and the adjacent surrounding nodes (Y), the stress in the designated region is calculated. The value of (stress strain) can be obtained with higher accuracy.

Furthermore, the accuracy of stress analysis in the designated area can be further improved by referring to the value of the reference node inside the designated area.

According to the second and third aspects of the invention, the four vertices of the designated area can be used for positioning of the designated area, so that accurate positioning can be easily performed. Further, since the positioning is easy, the operability is greatly improved.

According to the fourth aspect of the present invention, since the reference node is generated at least five locations, the stress strain can be calculated with high accuracy, and the reliability can be greatly improved.

The figure which shows the state before the analysis processing which analyzes with the algorithm of the conventional finite element method Stress-strain analysis diagram analyzed by conventional finite element algorithm Stress-strain analysis diagram presenting problems of the stress-strain analyzer according to the present invention Stress-strain analysis diagram presenting problems of the stress-strain analyzer according to the present invention Schematic showing a stress strain analyzer according to an embodiment of the present invention The block diagram which shows the structure of the stress analyzer (computer) which concerns on embodiment of this invention Explanatory diagram showing the procedure for obtaining the stress strain of the specified area specified by the user The flowchart which shows the stress analysis program which concerns on embodiment of this invention Figure showing mesh-based result evaluation by this analysis Figure showing the results of geometry-based results analysis

Hereinafter, a stress strain analyzer according to an embodiment of the present invention will be described with reference to the accompanying drawings.
In the present embodiment, components similar to those in the conventional example will be described using basically the same reference numerals (excluding suffixes).
FIG. 5 is a schematic diagram showing a stress-strain analyzer according to an embodiment of the present invention, FIG. 6 is a block diagram showing the configuration of a computer as the stress analyzer, and FIG. FIG. 8 is a flowchart showing a stress analysis program for obtaining the stress at the node of the analysis target data and the stress strain of the designated area designated by the user.

As shown in FIG. 5, the stress strain analysis apparatus 1 according to the present embodiment includes a computer 2, a keyboard 3 and a mouse 4 as input devices for a user to input instructions to the computer 2, and the computer 2. The LCD 5 is a display device that displays data output from the computer 2 on the screen, and the printing device 6 that prints data output from the computer 2. As shown in FIG. 5, the hardware resources of the computer 2 include a hard disk 7 as a storage device for reading and writing data, a memory 8, a backup power source 9 for holding data in the memory 8, and a CPU 10 for processing data. It is configured to include.

The memory 8 holds data referred to by the CPU 10, and the hard disk 7 has a program (hereinafter referred to as a stress analysis program) for causing the hardware resources of the computer 2 to function as the stress strain analysis apparatus 1. Analysis target data D (see FIG. 1) as analysis targets and data necessary for stress strain analysis and display are stored.

The stress analysis program includes a main program and subprograms. The main program displays the hardware resources of the computer 2 and the analysis target data D to be analyzed on the LCD 5, as shown in FIG. Display means 15 to be analyzed, dividing means 16 for dividing the analysis target data D into a plurality of elements E using a rectangular mesh M having a plurality of nodes Y, and stresses of the nodes Y of the elements E are known. It is configured to function as the calculation means 17 for calculating by the finite element method algorithm.

The subprogram uses the area designation means 18 for the user to arbitrarily designate the designated area AL for all areas of the analysis target data D via the mouse 4 or the keyboard 3 as input means. A plurality of reference nodes (see FIG. 7, P1 to P5) for analyzing the stress strain of the designated area AL designated by the area designation means 18 are generated on and inside the contour line AL1 of the designated area AL. The generating means 19 and the stress values of the plurality of generated reference nodes P1 to P4 are respectively determined based on the stress values of the plurality of surrounding nodes Y adjacent to the reference nodes, and the plurality of reference nodes determined. It is configured to function as calculation means 20 that calculates the stress strain of the designated area AL designated by the user based on the value.

Here, the display unit 15, the dividing unit 16, and the calculating unit 17 in the main program are mainly configured by a CPU, a memory, and a hard disk, and display functions, dividing functions, and calculations are performed by exchanging data between them. Realize the function. The display unit 15 displays the analysis target data D on the screen of the LCD 5 according to the procedure of a stress analysis program described later, and the dividing unit 16 selects the mesh M having a shape selected from the plurality of mesh data stored in the hard disk 7. As shown in FIG. 1, the analysis target data D is divided into a plurality of divisions with the number of divisions designated by the user. In the present embodiment, the designated area AL is designated so that the reference node P5 at the approximate center of the designated area AL substantially matches the node Y of the analysis target data D in the designated area AL according to the procedure described later. The number of divisions of the analysis target data D is the same as the conventional one.

The calculation unit 17 of the main program calculates the stress at the plurality of nodes Y of the analysis target data D divided into a plurality by the dividing unit 16 by an algorithm of the finite element method which is one of discrete methods, and calculates these stresses. Create a contour map based on this. Also, the stress and stress strain at each node Y are retrieved as link data and stored in the hard disk 7.

The display means 15 reads the stress of each node Y stored in the hard disk 7 and the data of the contour diagram, arranges each node Y at the coordinate position, and arranges the contour diagram. As a result, the user can visually know the magnitude and direction of the stress from the change in the coordinates of each node Y and the distance between the nodes, and can grasp the stress distribution state from the contour diagram.

In the subprogram, the area designating means 18 creates a designated area AL for the analysis target data D on the LCD 5 screen in cooperation with the CPU 10, mouse button, selection key, confirmation key displayed on the screen, cancel key, etc. Then, this is arranged at a position specified by the user.

The designated area AL may be any of a rectangle, a polygon, and a circle made of an extendable object (frame). In the present embodiment, the designated area AL is created using a rectangular object that expands and contracts in the diagonal direction due to ease of positioning, ease of operation, and high positioning accuracy, and whose area expands and contracts by operation in the diagonal direction. . The CPU 10 detects a user's area designation operation by a mouse button, a selection key, a confirmation key displayed on the screen, or a cancel key, and confirms after the start point and end point of the designated area AL are detected. When the key operation is detected, the user's designated area AL is fixed.

When the movement key is detected, the movement of the designated area AL is detected. When the cancel key operation is detected, the area designation operation by the designated area AL is canceled. Thereby, the designated area AL can be arranged at an arbitrary position in the analysis target data D.

When detecting the operation of the confirmation key, the CPU 10 outputs a reference node generation instruction to the reference node generation means 19. In response to this instruction, the reference node generation means 19 generates reference nodes P1 to P4 on the contour line AL1 of the designated area AL, that is, at each vertex of the rectangular designated area AL, and more preferably within the designated area AL. Generates a reference node P5 at the approximate center point of the designated area AL. When the generation of the reference nodes P1 to P5 is completed, the reference node generation unit 19 outputs a calculation request for the stress (stress strain) of the designated area AL to the calculation unit 20 of the subprogram.

Upon receiving this calculation request, the subprogram calculation means 20 uses the extrapolation method to calculate the values of the reference nodes P1 to P5 from the stress values of the surrounding nodes Y adjacent to the reference nodes P1 to P5. To calculate. Specifically, with respect to the reference node P1 on the left vertex of the upper side shown in FIG. 7, reference is made to the three nodes Y31, Y32, Y33 outside the specified area AL adjacent thereto and the approximate center of the specified area AL. The value of the reference node P1 is used to calculate the stress value of the reference node P1, and the stress value of the node Y34 designated as the stress analysis point by the node P1. Next, the calculation means 20 multiplies the stress values of these nodes Y31, Y32, Y33, and Y34 by the reciprocal of the distance (straight line distance) to the reference node P1, and performs weighted averaging.

For example, as shown in FIG. 7, the stress values of Y31, Y32, Y33, and Y34 are 0.1, 0.1, 0.03, and 0.05, respectively, and the reference nodes are connected from the nodes Y31, Y32, Y33, and Y34, respectively. When the distance to P1 is 2 mm, 4 mm, 3 mm, and 4 mm, the stress value after the weighted average of the reference node P1 is 0.1 × (1/2) (1/2 + 1/4 + 1/3 + 1/4) +0 .1 × (1/4) (1/2 + 1/4 + 1/3 + 1/4) + 0.03 × (1/3) (1/2 + 1/4 + 1/3 + 1/4) + 0.05 × (1/4) (1 /2+1/4+1/3+1/4)=0.073 (MPa).

Similarly, when the stress values after the weighted average of the reference nodes P2, P3, and P4 are obtained, the stress values after the weighted average of the reference nodes P2, P3, and P4 are 0.075 MPa, 0.137 MPa, and 0.135 MPa, respectively. It becomes. On the other hand, if the position of the reference node P5 is arranged at a position that substantially matches the node Y32 when the user designates the region, it can be regarded as a value that approximates the stress value of the node Y32. For this reason, in this example, the stress value of the node Y32 is set to the stress value of the reference node P5 (in this example, 0.1 (MPa)).

Next, the weighted average of each reference node P1, P2, P3, P4, P5 is performed, and this is used as the stress value of the designated area AL designated by the user. The weight is 1/4 at the apex position of the designated area AL and 1 at the center of the designated area AL. For this reason, the weight of the reference node P5 with respect to the weighted average is 1.
Accordingly, the stress in the designated area AL is 0.073 × 1/8 + 0.075 × 1/8 + 0.137 × 1/8 + 0.135 × 1/8 + 0.1 × 1/2 = 0.103 (MPa).
After calculating the stress value of the designated area AL, the stress strain of the designated area AL is calculated based on this stress value.

In this way, the reference node P5 is generated at the center of the designated area AL, and when the area of the designated area AL is designated, the reference node P5 is substantially matched with the node Y32 of the element E, and the outside and inside nodes of the designated area AL are set. Using the values of Y31, Y32, Y33, and Y34, the stress strain of the designated area AL is analyzed using an extrapolation method. These nodes Y31, Y32, Y33, Y34 are accurate values already obtained by the algorithm of the finite element method, and the stress values of the respective reference nodes P1 to P4 are Y31, Y32, Since the weighted average corresponding to the distance is applied to the stress values of Y33 and Y34, the values are highly accurate and reliable. Further, by referring to the value of the reference node P5, which is the value of the approximate center of the designated area AL, the stress analysis accuracy in the designated area AL can be further improved, and the reliability can be improved.
Thereby, the stress value and the stress strain of the designated area AL suitable for various uses such as collation with the measurement value measured with the strain gauge can be obtained.

Next, an analysis procedure of the stress strain analysis apparatus according to the present invention will be described with reference to FIG. When the stress analysis program is activated, the CPU 10 requests the LCD 5 to input analysis conditions in step S1. As the analysis conditions, for example, the material to be analyzed, Young's modulus, load conditions (position, direction, type), constraint conditions, the size of the mesh M used for analysis of stress strain, and the like are input. When the user confirms the input of analysis conditions, the process proceeds to step S2.

In step S2, the CPU 10 functions as the dividing unit 16, and the analysis target data (drawing data) D read in advance from the hard disk 7 or an external storage device (not shown) is converted into a rectangular mesh M (FIG. 1). To be divided into a plurality of parts. When the division ends, the process proceeds to step S3. In step S3, the CPU 10 activates the main program as a calculation means.

When the main program is activated, the CPU 10 proceeds to step S4. In step S4, the CPU 10 functions as the calculation means 17, and analyzes the stress at each node Y of each element E of the divided analysis target data D by a known finite element method algorithm. When the analysis of the stress at each node Y is completed, the process proceeds to step S5.

In step S5, the CPU 10 functions as the calculation means 17, and creates a contour diagram by generating contour lines of stress strain for the entire region of the analysis target data D based on the stress value at each node Y of each element E. To do. When the generation of the contour diagram is completed, the process proceeds to step S6, where the stress value of each node Y and the contour diagram are output to the LCD 5, and the main program is terminated with the output to the LCD 5 held. In step S7, the subprogram is activated.

When the subprogram is activated, the process proceeds to step S8, and a request for inputting a designated area is instructed to the user via the LCD 5. Next, it progresses to step S9 and it is determined whether the designated area | region AL was designated. If the confirm button is pressed in step S9, the process proceeds to step S10.

In step S10, the CPU 10 functions as the calculation means 20, and calculates the stress (stress strain) and the contour diagram of the designated area AL as described with reference to FIGS. When the calculation is completed, the result is output to the LCD 5. After the calculation, when the user operates the mouse release button of the mouse 4, the operation key of the keyboard, or the designation release button displayed on the LCD to input a designation release request to the CPU 10, the CPU 10 designates the designated area AL and the reference node P1. An instruction to erase the display of .about.P5 is output to the LCD 5.

FIG. 9 shows result evaluation on a mesh basis, and FIG. 10 shows result evaluation on a geometry basis. In any case, the arbitrarily designated points all showed the same numerical value (232.021 MPa in this example). Therefore, since any point can be evaluated on a geometry-based display in which a mesh is used as a background and this is not displayed, stress strain analysis can be facilitated.

The present invention is not limited to the above embodiments, and various modifications and applications can be made within the principle of the invention. For example, the region designation is not limited to a rectangle, but is a triangle or five or more corners. Polygons and circles can be set freely.

3: keyboard, 4: mouse, 17: calculation means, 19: reference node generation means, 20: calculation means,
D: Analysis target data, Y: Node, AL: Designated area, AL1: Contour line, P1 to P4: Multiple reference nodes, P5: Reference nodes inside the designated area

Claims (5)

The analysis target data (D) as the analysis target is regarded as an aggregate joined to each other at a finite number of nodes, and the analysis target data (D) is divided into a plurality of elements (M) having a plurality of nodes and each element ( E) In the stress-strain analyzer for analyzing the stress in the finite element algorithm,
A calculation means (17) for calculating stress values of a plurality of nodes (Y) of the analysis object data (D), and a designated area (AL) designated by the user for all areas of the analysis object data (D). Area designation means (18) arbitrarily designated via the input means (3, 4), and a plurality of reference nodes (P1 to P5) for analyzing stress strain in the designated area (AL) AL) on the outline (AL1) and in the vicinity thereof, the reference node generating means (19) for generating in the vicinity of the node (Y) ;
Calculating a stress value of the reference node generating means (19) has generated said plurality of reference nodes (P1 to P5), calculating means for calculating the stress-strain of the specified area on the basis of these stress values (20) With
The stress strain calculation means (20) of the designated area is
Based respective stress values of a plurality of reference nodes (P1 to P4) located on the contour line (AL1), the known stress value of a node of the ambient adjacent to each reference node (P1~P4) (Y) The
According to the distance between each of the reference nodes (P1 to P4) and the adjacent surrounding nodes (Y), it is calculated and calculated by a weighted average in which the weight is increased as the distance is shorter,
And the stress value of the reference node (P5) inside the designated area (AL) is calculated as a known stress value at the nearby node (Y),
A stress value of the designated region (AL) is calculated by weighting and averaging the obtained stress values of the plurality of reference nodes (P1 to P5) according to positions in the designated region (AL), and based on the stress values. stress-strain analysis apparatus and calculates the stress-strain Te The stress strain analysis apparatus according to claim 1, wherein a shape of the designated region is a rectangle. The stress strain analysis apparatus according to claim 2, wherein the reference node includes a rectangular vertex. The stress strain analyzer according to any one of claims 1 to 3, wherein the reference nodes are generated at least five locations. With the stress-strain analyzer , the analysis target data (D) as the analysis target is regarded as an aggregate joined to each other at a finite number of nodes, and the analysis target is divided into a plurality of meshes (M) having a plurality of nodes, In the stress strain analysis method for analyzing the stress in the element (E) by the algorithm of the finite element method,
The stress strain analyzer is
Calculating the stress value calculation means of a plurality of nodes (Y) (17), the stress value of the plurality of nodes (Y) of the analyzed data (D) in the mesh,
The region designation means (18), comprising the steps of: arbitrarily specify specified area (AL) with respect to the total area of the analysis target-flop with data (D),
By reference node generating means (19), the contour line (AL1) on the plurality of reference nodes (P1 to P5) of the designated area for analyzing the stress-strain of the designated area (AL) (AL) and in its interior Generating in the vicinity of the node (Y) ;
The stress distortion calculating means specified area (AL) (20), known around the node (Y) adjacent the respective stress values of the plurality of reference nodes (P1 to P4) to each reference node (P1 to P4) On the basis of the stress value of each of the reference nodes (P1 to P4) and the adjacent surrounding nodes (Y) according to the distance, the shorter the distance, the greater the weighted average,
Calculating a stress value of the reference node (P5) inside the designated area (AL) as a known stress value at the node (Y) in the vicinity thereof ;
The stress values of the plurality of reference nodes (P1 to P5) obtained are weighted and averaged according to the positions in the designated area (AL) to calculate the stress values in the designated area (AL), and the stress values are calculated based on the stress values. Calculating the stress strain of the specified area (AL);
A stress-strain analysis method comprising:

Images