JP5045347B2 - Data mapping method - Google Patents

Description

  In the present invention, for example, a first analysis model such as a structure analysis model created for a structure having an internal space in which a fluid flows is used, and an analysis method is used for an object different from the structure. The present invention belongs to a technical field related to a data mapping method for assigning data of a second analysis model created differently (for example, a fluid analysis model created for the fluid in the internal space).

  Conventionally, for example, in an engine cylinder head, a structural analysis model composed of a large number of elements (mesh) is created from the three-dimensional design data of the cylinder head, and a large number of fluids such as cooling water flowing in the internal space of the cylinder head are used. A fluid analysis model consisting of the elements (mesh) is created, and each element in a portion corresponding to the surface of the cylinder head that exchanges heat with the fluid (that is, the surface in contact with the internal space (fluid)) in the structural analysis model The data (temperature and heat transfer coefficient) of the element of the fluid analysis model surface (boundary element that is an element corresponding to the surface in contact with the cylinder head of the fluid in the fluid analysis model) is given (data mapping), The structural analysis model performs heat transfer calculations and stress calculations. They are to perform.

The data mapping operation as described above can be performed automatically by a computer as shown in, for example, Patent Document 1, and can be performed in a shorter time compared to the case where the operator performs the operation manually.
JP 2001-82391 A

  By the way, when the data of the fluid analysis model created for the fluid in the internal space is given to the structural analysis model created for the structure having an internal space such as the cylinder head, the internal structure of the structure is included in the structural analysis model. An element in a portion corresponding to a surface in contact with the space is extracted, and data of an element of the fluid analysis model corresponding to the element is assigned to the extracted element. Here, the portion corresponding to the internal space in the structural analysis model and the fluid analysis model have the same shape, and the extracted elements on the surface of the fluid analysis model when the fluid analysis model is superimposed on the portion. The element in contact with is an element of the fluid analysis model corresponding to the extracted element.

  Identification of elements of the fluid analysis model corresponding to each element extracted from the structural analysis model can be automatically performed by a computer. That is, a common three-dimensional coordinate axis is set for both models so that the same part of the structural analysis model corresponding to the internal space and the same part of the fluid analysis model have the same coordinates, and the extracted element Each time, the boundary position at the coordinate position close to the coordinate position of the extracted element is selected by comparing the coordinate position of the element with the coordinate position of all boundary elements on the surface of the fluid analysis model. The boundary element thus identified is identified as an element of the fluid analysis model corresponding to the extracted element.

  However, in the above specific method, for each of the extracted elements, it is necessary to compare the coordinate position of the element with the coordinate positions of all boundary elements on the surface of the fluid analysis model. is there. Particularly, in recent years, the number of elements tends to increase. For this reason, even if a specific operation is automatically performed by a computer, the operation time becomes considerably long.

  The present invention has been made in view of such a point, and an object of the present invention is to create first and second analysis models having different analysis methods such as the structural analysis model and the fluid analysis model. When the data of the second analysis model is given to the element (element to which data is to be given) in the portion corresponding to the predetermined part of the structure in the first analysis model, it corresponds to the predetermined part in the first analysis model An object of the present invention is to be able to specify the element (boundary element) of the second analysis model corresponding to the element in the portion in a short time.

In order to achieve the above object, according to the present invention, the control calculation unit of the computer extracts all elements (elements to which data is to be added) in a portion corresponding to a predetermined part of the structure from the first analysis model, A common three-dimensional coordinate axis is set for the first and second analysis models, and a space in which the first analysis model exists, which is defined by the three-dimensional coordinate axes, is composed of the extracted elements in the first analysis model. A space where the second analysis model exists, which is defined by the three-dimensional coordinate axes, is divided into the plurality of first analysis side regions so that the portion extends over the plurality of first analysis side regions. The boundary element of the second analysis model in the second analysis side region corresponding to the first analysis side region to which the extracted element belongs is divided into a plurality of second analysis side regions so as to correspond to the analysis side regions, respectively. Extracted from the inside And to identify the elements of the second analysis model corresponding to the element.

Specifically, in the first aspect of the present invention, each element in a portion corresponding to a predetermined portion of the structure is in contact with the predetermined portion in the first analysis model composed of a large number of elements created for the structure. The object of the present invention is to provide a data mapping method by a computer that gives data of a second analysis model consisting of a large number of elements, which is created by using an analysis method different from that of the first analysis model for an object different from the structure. To do.

The computer includes an arithmetic processing unit that performs arithmetic processing, and the arithmetic processing unit extracts all elements in a portion corresponding to a predetermined part of the structure from the first analysis model. And the arithmetic processing unit sets a common three-dimensional coordinate axis (X axis, Y axis, and Z axis) for the first and second analysis models, and the first processing unit in a space defined by the three-dimensional coordinate axis. And a coordinate axis setting step for causing the second analysis model to be in the same positional relationship as the structure and the object, and the first analysis model in which the arithmetic processing unit is defined by the three-dimensional coordinate axis. A first division step that divides the space into the plurality of first analysis side regions so that a portion composed of the elements extracted in the extraction step in the first analysis model straddles the plurality of first analysis side regions. When, the processing unit is defined by the three-dimensional coordinate axes, a space where the second analysis model is present, the plurality of second analysis-side region so as to correspond to each of the plurality of first analysis-side region The second dividing step for dividing, and the arithmetic processing unit extracted in the extraction step from boundary elements that are elements in a portion corresponding to the surface of the object that touches the predetermined part in the second analysis model. A specifying step for specifying an element of the second analysis model corresponding to the element, and a data adding step for adding data of the element specified in the specifying step to the element extracted in the extracting step by the arithmetic processing unit. , the specific steps, the arithmetic processing unit, the second analysis in the second analysis-side region corresponding to the first analysis-side region extracted in the extracting step element belongs From among the boundary element Dell, it assumed to be a step of specifying the elements of the second analysis model corresponding to the extracted element.

  According to the above data mapping method, in the extraction step, all elements in the portion corresponding to the predetermined part of the structure, that is, elements to which the data of the second analysis model is to be added are extracted from the first analysis model, and in the coordinate axis setting step A common three-dimensional coordinate axis is set for both models, and in the first division step, the space defined by the three-dimensional coordinate axis in which the first analysis model exists is divided into a plurality of first analysis-side regions, In the two division step, the space defined by the three-dimensional coordinate axis and in which the second analysis model exists is divided into a plurality of second analysis side regions so as to correspond to the plurality of first analysis side regions, respectively. Then, in the specifying step, the element of the second analysis model corresponding to the extracted element is specified from the boundary elements of the second analysis model, and in the data providing step, the specified element is included in the extracted element. Is given.

  In the specifying step, the second analysis model corresponding to the extracted element is selected from the boundary elements of the second analysis model in the second analysis side area corresponding to the first analysis side area to which the extracted element belongs. Elements are identified. As a result, it is not necessary to compare the coordinate positions of the extracted elements with the coordinate positions of all the boundary elements of the second analysis model, and the number of comparison target elements of the extracted one element can be reduced. Therefore, it becomes possible to specify the elements of the second analysis model extracted from the first analysis model and corresponding to the elements to which the data is to be added in a short time, thereby reducing the work time of the entire data mapping. can do.

According to a second aspect of the present invention, in the first aspect of the invention, the specifying step includes, for each element extracted by the arithmetic processing unit, the coordinate position of the element and a first analysis side region to which the element belongs. The boundary of the second analysis model at the coordinate position closest to the coordinate position of the extracted element is compared with the coordinate positions of all the boundary elements of the second analysis model in the second analysis side region corresponding to It is assumed that this is a step of selecting an element and specifying the selected boundary element as an element of the second analysis model corresponding to the extracted element.

  That is, when the size of the element of the second analysis model is smaller than the element of the first analysis model and there are a plurality of elements of the second analysis model corresponding to each element extracted from the first analysis model, the plurality of these It is conceivable to add the extracted element data to the extracted element by, for example, averaging the data. However, in this case, an extra work time such as averaging is required. However, in the present invention, only one boundary element at the coordinate position closest to the coordinate position of the extracted element is specified as the element of the second analysis model corresponding to the extracted element. It is not necessary to take extra work time such as averaging and the elements of the second analysis model corresponding to the extracted elements can be easily specified. Therefore, the work time of the entire data mapping can be further reduced. On the other hand, even if only one element data is given to each extracted element, for example, the result of the heat conduction analysis or the thermal stress analysis of the structure is not greatly affected, and no major problem occurs. .

According to a third aspect of the present invention, in the first or second aspect of the invention, the first division step includes a center of a portion configured by the element extracted by the extraction step in the first analysis model by the arithmetic processing unit. The space where the first analysis model is found is determined with respect to the center, the first first analysis side region on the X axis positive side, the Y axis positive side, and the Z axis positive side, and the X axis positive side. , A second first analysis side region on the Y axis positive side and the Z axis negative side, a third first analysis side region on the X axis positive side, the Y axis negative side and the Z axis positive side, and the X axis A fourth first analysis side region on the positive side, the Y-axis negative side, and the Z-axis negative side; a fifth first analysis-side region on the X-axis negative side, the Y-axis positive side, and the Z-axis positive side; Sixth first analysis side region on the X-axis negative side, Y-axis positive side, and Z-axis negative side, and a seventh first analysis-side region on the X-axis negative side, Y-axis negative side, and Z-axis positive side X axis negative side, Y axis A step of dividing into eight first analysis-side region of the eighth of the first analysis-side region on the side and the Z-axis negative side, the second division step, the arithmetic processing unit, the second analysis model The center of the portion composed of the surface elements is obtained, and the space where the second analysis model exists is located on the X axis positive side, the Y axis positive side, and the Z axis positive side with respect to the center. The first second analysis side region corresponding to the first first analysis side region, and the second first analysis side region on the X axis positive side, the Y axis positive side, and the Z axis negative side correspond to the second first analysis side region. A second second analysis side region, a third second analysis side region corresponding to the third first analysis side region on the X axis positive side, the Y axis negative side, and the Z axis positive side, and X A fourth second analysis side region corresponding to the fourth first analysis side region, an X axis negative side, a Y axis positive side, and a Z axis on the positive axis side, the Y axis negative side, and the Z axis negative side Positive The fifth second analysis side region corresponding to the fifth first analysis side region, and the sixth first analysis side on the X axis negative side, the Y axis positive side, and the Z axis negative side A sixth second analysis side region corresponding to the region, and a seventh second analysis side corresponding to the seventh first analysis side region on the X axis negative side, the Y axis negative side, and the Z axis positive side Eight second analysis side regions including the region and the eighth second analysis side region corresponding to the eighth first analysis side region on the X axis negative side, the Y axis negative side, and the Z axis negative side It is assumed that the step is a division step.

  As a result, the number of boundary elements of the second analysis model in each second analysis side region can be made substantially uniform, and the number of comparison target elements for any of the extracted elements can be set to all boundary elements of the second analysis model. Can be reduced to about 1/8 of In addition, the space where the first and second analysis models exist can be easily divided. Therefore, it is possible to significantly reduce the time required for specifying the element of the second analysis model corresponding to the extracted element while suppressing an increase in the division time.

As described above, according to the data mapping method of the present invention , the first and second analysis models defined by the common three-dimensional coordinate axes set by the control calculation unit of the computer for the first and second analysis models are obtained. The existing space is divided into a plurality of first and second analysis side regions corresponding to each other, and is in the second analysis side region corresponding to the first analysis side region to which the element extracted from the first analysis model belongs. By specifying the element of the second analysis model corresponding to the extracted element from the boundary elements of the second analysis model, it is necessary to specify the element of the second analysis model corresponding to the extracted element. The time can be shortened, and the work time of the entire data mapping can be shortened.

  Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings. In the present embodiment, the first analysis model is a structure analysis model created for a structure having an internal space through which fluid flows, and the second analysis model is a fluid in the internal space of the structure (as a predetermined part of the structure). And a fluid analysis model created for an object that is in contact with a surface in contact with the internal space of the structure and is separate from the structure), and conducts heat conduction analysis and thermal stress analysis of the structure.

  FIG. 1 shows an operation procedure when conducting heat conduction analysis or thermal stress analysis of the structure. In the present embodiment, the internal space has an open part that opens to the outside of the structure, and the fluid flows into and out of the internal space through the open part. It may be what was done. The steps of the following flowcharts are basically processes performed by a computer, but also include processes performed by an operator. When an operator performs the process, the fact is clearly indicated.

  First, three-dimensional (3D) design data (for example, CAD data) of a structure having an internal space through which fluid flows is input from a storage device that stores the three-dimensional design data (step S1).

  Subsequently, a fluid analysis model including a large number of elements (tetramesh, hexamesh, etc.) is created for the fluid in the internal space from the 3D design data (step S2). The fluid analysis model is created using existing general-purpose finite element software. Then, the operator sets the boundary condition of the fluid (step S3), executes the fluid calculation (thermal fluid analysis) based on the fluid analysis model and the boundary condition (step S4), and shows the calculation result. Post processing for outputting data (temperature and heat transfer coefficient) for each element is performed (step S5).

  On the other hand, in parallel with the thermal fluid analysis for the fluid, a structural analysis model composed of a large number of elements (tetramesh, hexamesh, etc.) is created for the structure from the 3D design data (step S6). Similar to the fluid analysis model, this structural analysis model is also performed using existing general-purpose finite element method software. Note that the size of the element of the structural analysis model is larger than the size of the element of the fluid analysis model.

  Subsequently, all elements in the portion corresponding to the surface in contact with the internal space of the structure, that is, elements to which the data of the fluid analysis model should be added (hereinafter referred to as boundary condition setting target elements) are extracted from the structural analysis model. Then, the extracted all boundary condition setting target elements (including reference elements described later) are output (step S8).

  Next, data mapping is performed to assign the fluid analysis model data to the extracted boundary condition setting target element (the fluid analysis model data is set as the boundary condition) (step S9), and the above-described assignment is given to the structural analysis model. The obtained data is taken into an input file for heat transfer calculation (or thermal stress calculation) (step S10), and heat transfer calculation (or thermal stress calculation) is executed in the structural analysis model (step S11). Post processing for outputting data (temperature or thermal stress) for each element is performed (step S12).

  In the present embodiment, the extraction of the boundary condition setting target element in step S7 is specifically performed as shown in the flowchart of FIG.

  That is, first, an operator designates an element in a portion corresponding to the outer surface of the structure in the structural analysis model (that is, an element on the outer surface of the structural analysis model) (step S21). The operator can select any element in the portion corresponding to the vicinity of the open portion from among all the elements in the portion corresponding to the surface in contact with the internal space of the structure in the structural analysis model (visible to the operator). Element) is set as a reference element (step S22). This reference element is the boundary condition setting target element.

  Then, in the structural analysis model, the element in the portion corresponding to the surface in contact with the internal space adjacent to the reference element is extracted as a boundary condition setting target element. Specifically, the surface of the structure (the surface in contact with the outer surface and the inner space) adjacent to the reference element (the boundary condition setting target element extracted in step S25 described later when returning from step S28 described later). ) Is selected (step S23), and it is determined whether or not the selected element is an element in the portion corresponding to the outer surface of the structure designated in step S21 (step S23). Step S24).

  When the selected element is not an element in a portion corresponding to the outer surface, that is, an element in a portion corresponding to a surface in contact with the internal space (when the determination in step S24 is NO), the selected element is selected. Elements are extracted as boundary condition setting target elements (step S25), and then the number of boundary condition setting target elements is counted up (step S26), and then the process proceeds to step S28.

  On the other hand, when the selected element is an element in a portion corresponding to the outer surface (when the determination in step S24 is YES), the selected element is not extracted as a boundary condition setting target element (step S27), the process proceeds to step S28.

  In step S28, it is determined whether or not the count number of the boundary condition setting target element is increased (step S28). When the count number is increased (when the determination in step S28 is YES), the above step is performed. Returning to S23, this time, an element adjacent to the extracted boundary condition setting target element is selected, and the processes of steps S24 to S28 are repeated. On the other hand, when the count number does not increase (when the determination in step S28 is NO), that is, extraction of all the boundary condition setting target elements is completed, and only the elements in the portion corresponding to the outer surface are selected. Sometimes, the boundary condition setting target element extraction process is stopped.

  The element extracted as the boundary condition setting target element is stored in the storage device of the computer. However, in this embodiment, in order to simplify the processing, the selected element is extracted as the boundary condition setting target element. The element which does not do is not memorized. For this reason, although the element extracted as the boundary condition setting target element is not selected again, the element that has been selected but not extracted as the boundary condition setting target element may be selected again. Even if it does in this way, since it is not extracted as a boundary condition setting object element, there is no problem.

  In the above extraction method, elements in the portion corresponding to the surface in contact with the internal space of the structure in the structural analysis model can be extracted easily and in a short time. That is, the element in the part corresponding to the surface in contact with the internal space of the structure is different from the element in the part corresponding to the outer surface and is inside the structural analysis model. Therefore, it is necessary to extract by sequentially removing elements from the outer surface in the structural analysis model. On the other hand, the reference element is an element in a portion corresponding to the surface in contact with the internal space, but is an element in a portion corresponding to the vicinity of the open portion. The setting can be easily made by the operator himself / herself while specifying the structural analysis model. An element adjacent to the reference element and corresponding to the surface of the structure (an element in which at least one surface of the element is not in contact with another element) can be automatically extracted by a computer. Here, since the reference element is an element in a portion corresponding to the surface in contact with the internal space in the structural analysis model, the element in the portion corresponding to the surface of the structure adjacent to the reference element one after another is Usually, the element is in a portion corresponding to the surface in contact with the internal space, but may be an element in a portion corresponding to the outer surface of the structure. Since the element in the part corresponding to the outer surface can be easily identified, it is determined whether the element is in the part corresponding to the outer surface, and the part corresponding to the outer surface of the structure. If it is not an element, it can be determined that the element is in a part corresponding to the surface in contact with the internal space, and as described above, the element is extracted and the element in the part corresponding to the outer surface of the structure Do not extract. In this way, it is possible to easily and quickly extract elements in a portion corresponding to the surface in contact with the internal space of the structure in the structural analysis model.

  Here, the reference element may be one or plural, and it is preferable to set a plurality of reference elements particularly when a plurality of open portions are provided. In this case, it is preferable that the plurality of reference elements are elements in portions corresponding to the vicinity of different open portions in the structural analysis model. In this way, the boundary condition setting target element extraction process can be completed in a shorter time.

  If the internal space does not have an open portion that opens to the outside of the structure (when the internal space is closed), the extraction method cannot be used. A method of extracting by sequentially removing from the surface elements may be used. Further, even when the internal space has an open portion as in the present embodiment, a method of extracting by sequentially removing elements from the outer surface may be used.

  Specifically, the data mapping in step S9 is performed as shown in the flowchart of FIG.

First, as the same portion of the part and the fluid dynamics model, which corresponds to the internal space of the upper Symbol structural analysis model are the same coordinates, common three-dimensional coordinate axes for the both models (X axis, Y axis and Z Axis) is set (step S51). That is, in the space defined by the three-dimensional coordinate axis, the positional relationship between the two models is made the same as the positional relationship between the structure and the fluid. Here, the portion corresponding to the internal space in the structural analysis model has the same shape as the fluid analysis model, and in the space defined by the three-dimensional coordinate axis, the fluid analysis is performed on the portion corresponding to the internal space in the structural analysis model. The models just overlap .

  Subsequently, a space defined by the three-dimensional coordinate axis in which the structural analysis model exists is divided into a plurality of structural analysis side regions (first analysis side). The region is divided into a plurality of structural analysis side regions so as to straddle the region. In the present embodiment, first, the coordinates of the center of the portion constituted by the boundary condition setting target element in the structural analysis model (hereinafter referred to as the structural analysis model center) by the three-dimensional coordinate axis are obtained (step S52). Specifically, the maximum value and the minimum value of the X coordinate are extracted from the center coordinate data of each of the boundary condition setting target elements by extracting the maximum value and the minimum value of the X coordinate, the Y coordinate, and the Z coordinate. Is the X coordinate of the center of the structural analysis model, the average value of the maximum and minimum values of the Y coordinate is the Y coordinate of the center of the structural analysis model, and the average value of the maximum and minimum values of the Z coordinate. Is the Z coordinate of the center of the structural analysis model.

  Then, the space where the structural analysis model exists is defined with respect to the first structural analysis side region F1 on the X-axis positive side, the Y-axis positive side, and the Z-axis positive side with respect to the center of the structural analysis model, Side, Y-axis positive side and Z-axis negative side second structural analysis side region F2, X-axis positive side, Y-axis negative side and Z-axis positive side third structural analysis-side region F3, X A fourth structural analysis side region F4 on the positive axis side, the negative Y axis side, and the negative side of the Z axis, and a fifth structural analysis side region F5 on the negative side of the X axis, the positive side of the Y axis, and the positive side of the Z axis; The sixth structural analysis side region F6 on the X axis negative side, the Y axis positive side, and the Z axis negative side, and the seventh structural analysis side region on the X axis negative side, the Y axis negative side, and the Z axis positive side The region is divided into eight structural analysis side regions including F7 and an eighth structural analysis side region F8 on the X axis negative side, the Y axis negative side, and the Z axis negative side (step S53). Thereby, the part comprised by the said boundary condition setting object element in a structural analysis model straddles a plurality (in this embodiment, eight) structure analysis side area | regions.

  After the division, the boundary condition setting target elements are classified into F1 to F8 structure analysis side regions (step S54).

  Next, a plurality of fluid analysis side regions (second analysis) are defined so that the space where the fluid analysis model exists, defined by the three-dimensional coordinate axes, corresponds to the plurality of structural analysis side regions (F1 to F8), respectively. Side area). In the present embodiment, first, the coordinates of the center of the portion composed of the surface elements in the fluid analysis model (hereinafter referred to as the fluid analysis model center) by the three-dimensional coordinate axis are obtained (step S55). Specifically, the maximum value and the minimum value of the X coordinate are extracted from the center coordinate data of each element of the fluid analysis model by extracting the maximum value and the minimum value of the X coordinate, the Y coordinate, and the Z coordinate, respectively. Is the X coordinate of the center of the fluid analysis model, the average value of the maximum and minimum values of the Y coordinate is the Y coordinate of the center of the fluid analysis model, and the average value of the maximum and minimum values of the Z coordinate. Is the Z coordinate of the center of the fluid analysis model.

  Then, the space where the fluid analysis model exists is located on the X-axis positive side, the Y-axis positive side, and the Z-axis positive side with respect to the fluid analysis model center. One fluid analysis side region C1, a second fluid analysis side region C2 corresponding to the second structural analysis side region F2 on the X axis positive side, the Y axis positive side, and the Z axis negative side, and the X axis The third fluid analysis side region C3 corresponding to the third structural analysis side region F3 on the positive side, the Y axis negative side, and the Z axis positive side, the X axis positive side, the Y axis negative side, and the Z axis negative side A fourth fluid analysis side region C4 corresponding to the fourth structural analysis side region F4 on the side, and the fifth structural analysis side on the X axis negative side, the Y axis positive side, and the Z axis positive side. Corresponds to the fifth fluid analysis side region C5 corresponding to the region F5 and the sixth structural analysis side region F6 on the X axis negative side, Y axis positive side and Z axis negative side. A sixth fluid analysis side region C6, a seventh fluid analysis side region C7 corresponding to the seventh structural analysis side region F7 on the X axis negative side, the Y axis negative side, and the Z axis positive side; The fluid is divided into eight fluid analysis side regions, which are the eighth fluid analysis side region C8 corresponding to the eighth structural analysis side region F8 on the negative axis side, the negative Y axis side, and the negative Z axis side (step S56). ).

  The coordinates of the center of the structural analysis model and the coordinates of the center of the fluid analysis model are basically the same (in this embodiment, the data of the center coordinates of each boundary condition setting target element and each of the fluid analysis model as described above. Since the average value of the maximum value and the minimum value of the center coordinate data of the element is used, there is a possibility of slight deviation), and the structural analysis side region and the fluid analysis side region corresponding to each other are in the above three-dimensional coordinate axes. It will be located in the same part in the defined space.

  After the division, the fluid analysis model elements are classified into C1 to C8 fluid analysis side regions (step S57).

  Next, the element of the fluid analysis model corresponding to the boundary condition setting target element is identified from the boundary elements that are elements on the surface of the fluid analysis model (elements corresponding to the surface in contact with the predetermined portion of the fluid) To do.

  Specifically, it corresponds to the boundary condition setting target element from the elements (boundary elements) on the surface of the fluid analysis model in the fluid analysis side area corresponding to the structural analysis side area to which the boundary condition setting target element belongs. Identify the elements of the fluid analysis model.

  In the present embodiment, first, the element of the fluid analysis model corresponding to the boundary condition setting target element in the first structure analysis side region F1 is specified (steps S58 and S59). That is, the element of the fluid analysis model corresponding to the boundary condition setting target element is specified from the elements on the surface of the fluid analysis model in the first fluid analysis side area C1 corresponding to the first structure analysis side area F1. . Specifically, for each boundary condition setting target element in the first structural analysis side region F1, the coordinate position of the boundary condition setting target element and all the elements on the surface of the fluid analysis model in the first fluid analysis side region C1 Is compared with the coordinate position of the element in the fluid analysis model at the coordinate position closest to the coordinate position of the boundary condition setting target element (that is, in the structure analysis model in the space defined by the three-dimensional coordinate axis). When the fluid analysis model is overlaid on the part corresponding to the internal space, the fluid analysis model surface element in contact with the boundary condition setting target element is selected, and the selected element corresponds to the boundary condition setting target element. Identified as an element of the fluid analysis model.

  Subsequently, the element data (temperature and heat transfer coefficient) of the identified fluid analysis model surface are assigned to each boundary condition setting target element in the first structural analysis side region F1 (step S60). And similarly about the boundary condition setting object element in the next 2nd structural analysis side area | region F2, the element of the fluid analysis model corresponding to the said boundary condition setting object element is specified similarly, 2nd structural analysis side area | region The element data of the specified fluid analysis model surface is assigned to the boundary condition setting target element in F2. This is repeated up to the eighth structural analysis side region F8 (steps S59 to S62), the elements on the surface of the fluid analysis model corresponding to the boundary condition setting target elements in all the structural analysis side regions are specified, and the boundary condition setting is performed. Data of the element on the surface of the identified fluid analysis model is assigned to the target element.

  The above-mentioned data may be added every time an element on the surface of the fluid analysis model is specified, and the elements on the surface of the fluid analysis model corresponding to the boundary condition setting target elements in all the structural analysis side regions are specified. You may do it all together later.

  Next, the data given to the boundary condition setting target element is output in the format for structural analysis in order to take in the input file for heat transfer calculation (or thermal stress calculation) in step S10 (step S63). .

  The structural analysis side region and the fluid analysis side region are not limited to F1 to F8 and C1 to C8, and a portion constituted by boundary condition setting target elements in the structural analysis model is divided into a plurality of structural analysis side regions. As long as it straddles, the number and the shape of the structure analysis side region may be any, and the fluid analysis side region corresponds to the structure analysis side region in the same part in the space defined by the three-dimensional coordinate axes. It only has to be located.

  Here, the case where a structure is a cylinder head of an engine is shown as an example.

  FIG. 4 shows a structural analysis model M1 made of a large number of elements (tetramesh in this example) created for the cylinder head. Reference numeral 1 denotes a portion corresponding to the cooling water inlet to the water jacket. The cooling water inlet is provided on a contact surface of the cylinder head with the cylinder block so that the cooling water flows from the cylinder block to the water jacket. Although not visible in FIG. 4, there is a portion corresponding to the cooling water outlet on the left side surface of FIG. 4 in the structural analysis model M1 (in FIG. 6, the portion corresponding to the cooling water outlet is denoted by reference numeral 2). ing). Thus, the cylinder head has a water jacket as an internal space through which cooling water (fluid) flows, and this water jacket is provided with a cooling water inlet and outlet as an open portion that opens to the outside of the cylinder head. The cooling water flows into and out of the water jacket through the cooling water inlet and outlet.

  When extracting the boundary condition setting target element (the element in the portion corresponding to the surface in contact with the water jacket) from the structural analysis model M1, the operator looks at the structural analysis model M1 as shown in FIG. Among all the elements in the part corresponding to the surface in contact with the element, the element in the part corresponding to the vicinity of the cooling water inlet, which is a visible element (for example, the cooling water inlet in the structural analysis model of FIG. 4). An element E1) shown in FIG. 5 with an enlarged portion is set as a reference element. The orientation of the structural analysis model M1 can be changed so that an operator can see the side of the cooling water outlet, and in this way, an element in a portion corresponding to the vicinity of the cooling water outlet is set as a reference element. can do.

  The reference element set in this way is used as a core, and the elements in the portion corresponding to the surface in contact with the water jacket that are adjacent to each other are automatically extracted as the boundary condition setting target elements. Here, an element in a portion corresponding to the cylinder block contact surface provided with the cooling water inlet (for example, an element E2 shown in FIG. 5) or an element in a portion corresponding to the side surface provided with the cooling water outlet is provided. However, these elements are selected (see step S23 above), but these elements are elements that are pre-designated by the operator and are in a portion corresponding to the outer surface of the cylinder head. It is never extracted. Then, all the elements in the portion corresponding to the surface in contact with the water jacket are extracted, and there is no extracted element (the element (element E2 etc.) in the portion corresponding to the outer surface of the cylinder head is selected). The extraction process ends.

  FIG. 6 shows a portion composed of the extracted boundary condition setting target elements in the structural analysis model M1. Note that FIG. 6 is different from FIG. 4 in the viewing direction, and the portion 2 corresponding to the cooling water outlet is visible.

  FIG. 7 is composed of elements (boundary elements) on the surface of a fluid analysis model M2 made of a large number of elements (in this example, hexamesh) created for cooling water that is a fluid flowing through the water jacket of the cylinder head. This part is shown. 6 and 7, common three-dimensional coordinate axes (X axis, Y axis, and Z axis) are set for the structural analysis model M1 and the fluid analysis model M2.

  In the space defined by the three-dimensional coordinate axis, when the fluid analysis model M2 is superimposed on a portion corresponding to the internal space in the structural analysis model M1, the boundary analysis setting target element is extracted from the structural analysis model M1. The portion and the surface of the fluid analysis model M2 just touch each other. While the number of elements on the surface of the fluid analysis model M2 is about 110,000, the number of boundary condition setting target elements extracted from the structural analysis model M1 is about 80,000. This is because the size of the element of the structural analysis model M1 is larger than the size of the element of the fluid analysis model M2. For this reason, although there are a plurality of elements (boundary elements) on the surface of the fluid analysis model M2 in contact with one boundary condition setting target element, in the present embodiment, as elements of the fluid analysis model M2 corresponding to one boundary condition setting target element Specifies only one element (the element at the coordinate position closest to the coordinate position of the boundary condition setting target element). Even if it does in this way, it will not have a big influence on the result of the heat conduction analysis or thermal stress analysis of a cylinder head, and a big problem will not arise.

  A plurality of elements may be specified as elements of the fluid analysis model M2 corresponding to one boundary condition setting target element. In this case, for example, the element (s) on the surface of the fluid analysis model M2 at a coordinate position within a predetermined distance from the coordinate position of the boundary condition setting target element is used as an element of the fluid analysis model M2 corresponding to the boundary condition setting target element. What is necessary is just to specify and give to the said boundary condition setting object element by averaging the data of these some elements.

  In this embodiment, when specifying the element of the fluid analysis model M2 corresponding to the boundary condition setting target element, the coordinate position of the boundary condition setting target element is not compared with the coordinate positions of all elements on the surface of the fluid analysis model M2. The space where the structural analysis model M1 and the fluid analysis model M2 exist, defined by the three-dimensional coordinate axes, is divided into eight structural analysis side regions F1 to F8 and fluid analysis side regions C1 to C8 corresponding to each other. Then, the coordinate position of the boundary condition setting target element is compared with the coordinate position of the surface of the fluid analysis model M2 in the fluid analysis side area corresponding to the structural analysis side area to which the boundary condition setting target element belongs. The number of comparison target elements of one boundary condition setting target element can be reduced to about 1/8 (110,000 / 8) compared with all elements on the surface of the fluid analysis model M2. Kill.

  Therefore, in this embodiment, all elements in the part corresponding to the surface in contact with the internal space of the structure are extracted from the structural analysis model as boundary condition setting target elements, and a common three-dimensional coordinate axis is set for both models. The space defined by the three-dimensional coordinate axes in which the structural analysis model exists, the plurality of the plurality of structural analysis side regions in the structural analysis model so that the portion constituted by the boundary condition setting target elements extends over the plurality of structural analysis side regions. The space where the fluid analysis model exists, which is divided into structural analysis side regions and defined by the three-dimensional coordinate axes, is divided into a plurality of fluid analysis side regions so as to correspond to the plurality of structural analysis side regions, respectively. The boundary condition setting target from the fluid analysis model surface elements (boundary elements) in the fluid analysis side area corresponding to the structural analysis side area to which the boundary condition setting target element belongs Since the element of the fluid analysis model corresponding to the element is specified, the element of the fluid analysis model corresponding to the boundary condition setting target element can be specified in a short time. The work time can be shortened.

  In addition, since the element at the coordinate position closest to the coordinate position of the boundary condition setting target element is specified as the element of the fluid analysis model corresponding to the boundary condition setting target element, the data of a plurality of elements is averaged In addition, it is possible to easily identify the elements of the fluid analysis model corresponding to the boundary condition setting target elements, and to further reduce the work time of the entire data mapping. it can.

  Further, the space defined by the three-dimensional coordinate axis in which the structural analysis model exists is divided into the eight structural analysis-side regions F1 to F8, and the space defined by the three-dimensional coordinate axis and in which the fluid analysis model exists exists. Since the fluid analysis model is divided into the eight fluid analysis side regions C1 to C8, the number of elements on the surface of the fluid analysis model in each fluid analysis side region can be made substantially uniform. The number can be reduced to about 1/8 compared with all elements on the surface of the fluid analysis model, and the space where the structural analysis model and the fluid analysis model exist can be easily divided. Therefore, it is possible to greatly reduce the time required for specifying the element of the fluid analysis model corresponding to the boundary condition setting target element while suppressing an increase in the division time.

  In the above embodiment, the case where the fluid flows in the internal space of the structure has been described. However, the present invention can be applied even when the fluid flows in contact with a predetermined portion such as the outer surface of the structure. it can. The second analysis model is not limited to the fluid analysis model, and may be an analysis model created for an object different from the structure. For example, when powder flows in the internal space of a structure, create a powder analysis model for the powder and assign the data of the powder analysis model to the structure analysis model created for the structure. You just have to do it. Further, the first analysis model is not limited to the structural analysis model, and the first and second analysis models may be analysis models having different analysis methods.

  In the first analysis model created for a structure, the present invention relates to an object that is in contact with the predetermined part of each element in a part corresponding to the predetermined part of the structure and that is different from the structure. The analysis model is useful for a data mapping method for assigning data of a second analysis model created by using a different analysis method. In particular, the structure analysis model created for a structure having an internal space through which a fluid flows This is useful for a data mapping method for assigning data of a fluid analysis model created for a fluid in an internal space.

It is a flowchart which shows the operation | movement procedure in the case of performing the heat conduction analysis and the thermal stress analysis of the structure which has the internal space where a fluid distribute | circulates. It is a flowchart which shows the detailed procedure of extraction of a boundary condition setting object element. It is a flowchart which shows the detailed procedure of a data mapping. It is a perspective view which shows the structural analysis model created about the cylinder head. It is a figure which expands and shows the part corresponded to the cooling water inlet in the structural analysis model of FIG. FIG. 5 is a perspective view showing a part composed of boundary condition setting target elements in the structural analysis model of FIG. 4. It is a perspective view which shows the part comprised by the element in the surface in a fluid analysis model.

1 part corresponding to the cooling water inlet in the structural analysis model 2 part E1 corresponding to the cooling water outlet in the structural analysis model element E2 in the part corresponding to the vicinity of the cooling water inlet in the surface in contact with the water jacket in the structural analysis model In the part corresponding to the outer surface of the cylinder head in

Claims (3)

In the first analysis model composed of a number of elements created for a structure, each element in a portion corresponding to the predetermined portion of the structure is in contact with the predetermined portion, and the object different from the structure is described above. The first analysis model is a data mapping method by a computer that assigns data of a second analysis model composed of a large number of elements, created using a different analysis method,
The computer has an arithmetic processing unit that performs arithmetic processing.
An extraction step in which the arithmetic processing unit extracts all elements in a portion corresponding to a predetermined part of the structure from the first analysis model;
The arithmetic processing unit sets a common three-dimensional coordinate axis (X axis, Y axis, and Z axis) for the first and second analysis models, and the first and second in a space defined by the three-dimensional coordinate axis. (2) a coordinate axis setting step for causing the analysis model to have the same positional relationship as the structure and the object;
The arithmetic processing unit includes a plurality of first analyzes in which a portion defined by the three-dimensional coordinate axes and including the elements extracted in the extraction step in the first analysis model is a plurality of first analysis models. A first dividing step of dividing the plurality of first analysis side regions so as to straddle the side regions;
The arithmetic processing unit divides a space defined by the three-dimensional coordinate axis in which the second analysis model exists into a plurality of second analysis side regions so as to correspond to the plurality of first analysis side regions, respectively. A second division step;
The arithmetic processing unit, the out of the boundary element is an element in the portion corresponding to the surface in contact with the predetermined portion of the object in the second analysis model, the second analysis model that correspond to elements extracted in the extraction step Specific steps to identify the elements of
The arithmetic processing unit includes a data providing step for giving data of the element specified in the specifying step to the element extracted in the extracting step,
In the specifying step, the arithmetic processing unit includes the boundary element of the second analysis model in the second analysis side region corresponding to the first analysis side region to which the element extracted in the extraction step belongs. A data mapping method characterized by being a step of specifying an element of a second analysis model corresponding to an extracted element. The data mapping method according to claim 1, wherein
In the specifying step, for each element extracted by the extraction step , the calculation processing unit includes the second analysis side region corresponding to the coordinate position of the element and the first analysis side region to which the element belongs. By comparing the coordinate positions of all the boundary elements of the analysis model, the boundary element of the second analysis model at the coordinate position closest to the coordinate position of the extracted element is selected, and the selected boundary element is A data mapping method characterized by being a step of specifying as an element of a second analysis model corresponding to an extracted element. The data mapping method according to claim 1 or 2,
In the first division step, the arithmetic processing unit obtains a center of a portion composed of the elements extracted in the extraction step in the first analysis model, and the space where the first analysis model exists is determined as the center. In contrast, the first first analysis side region on the X axis positive side, the Y axis positive side, and the Z axis positive side, and the second one on the X axis positive side, the Y axis positive side, and the Z axis negative side. 1 analysis side region, a third first analysis side region on the X axis positive side, Y axis negative side and Z axis positive side, and a fourth on the X axis positive side, Y axis negative side and Z axis negative side The first analysis side region, the fifth first analysis side region on the X-axis negative side, the Y-axis positive side, and the Z-axis positive side, and the X-axis negative side, the Y-axis positive side, and the Z-axis negative side Sixth first analysis side region, seventh first analysis side region on the X-axis negative side, Y-axis negative side, and Z-axis positive side, X-axis negative side, Y-axis negative side, and Z-axis negative side 8 with the first analysis side area in A step of dividing the first analysis-side region,
In the second division step, the arithmetic processing unit obtains a center of a portion constituted by surface elements in the second analysis model, and a space where the second analysis model exists is defined with respect to the center. A first second analysis side region corresponding to the first first analysis side region on the X-axis positive side, the Y-axis positive side, and the Z-axis positive side, the X-axis positive side, the Y-axis positive side, and the Z-axis A second second analysis side region corresponding to the second first analysis side region on the negative axis side, and a third third side on the X axis positive side, the Y axis negative side, and the Z axis positive side. A third second analysis side region corresponding to one analysis side region and a fourth second corresponding to the fourth first analysis side region on the X axis positive side, the Y axis negative side and the Z axis negative side. 2 analysis side regions, a fifth second analysis side region corresponding to the fifth first analysis side region on the X axis negative side, the Y axis positive side, and the Z axis positive side, the X axis negative side, Y axis positive side and Z A sixth second analysis side region corresponding to the sixth first analysis side region on the negative side, and the seventh first side on the X axis negative side, the Y axis negative side, and the Z axis positive side. A seventh second analysis side region corresponding to the analysis side region, and an eighth second corresponding to the eighth first analysis side region on the X axis negative side, the Y axis negative side, and the Z axis negative side. A data mapping method characterized by being a step of dividing into eight second analysis side regions with the analysis side region.

Images