Engine cylinder temperature field finite element result curve automatic extraction method and medium

Abstract

The present application relates to a kind of engine cylinder temperature field finite element result curve automatic extraction method and medium;Including: making engine cylinder axial and circumferential temperature field finite element calculation result curve extraction standardization rule;Based on TCL script language, finite element calculation result curve extraction is programmed and programmed;Adjust geometric coordinate system;Engine cylinder geometry is pretreated;Engine cylinder is modeled by finite element;Carrying out engine cylinder temperature field finite element simulation;Using hyperview software reads finite element simulation calculation result file;Calling the engine cylinder axial and circumferential temperature field finite element calculation result curve automatic extraction program prepared;The present application uses TCL script language, custom finite element calculation result curve extraction standardization rule, realizes the engine cylinder axial and circumferential temperature field finite element calculation result curve extraction automation, compared with traditional manual work cycle, work efficiency is shortened by more than 95%, improves work efficiency.

Description

Technical Field

[0001] This invention relates to a method and medium for automatically extracting finite element result curves of engine cylinder temperature field, specifically to a method and computer-readable storage medium for automatically extracting finite element calculation result curves of engine cylinder axial and circumferential temperature fields. Background Technology

[0002] In the development and design of automotive engines, finite element simulation of the engine block's temperature field is a mandatory analysis item. This process requires extracting and evaluating the finite element calculation results curves of the engine cylinder's axial and circumferential temperature fields. However, extracting these curves involves numerous factors, including extraction location definition, geometric processing requirements, mesh generation requirements, and curve extraction itself. The data processing is complex, time-consuming, and prone to errors. Furthermore, inconsistent processing methods among different engineers lead to inconsistent analysis results. Existing technologies focus on better engine block structure design, enabling the evaluation of engine block temperature performance, but they cannot effectively avoid the aforementioned problems. To address these issues, this paper utilizes the TCL scripting language, combined with standardized rules and design flowcharts for extracting the finite element calculation results curves of the engine cylinder's axial and circumferential temperature fields, to automate the extraction process, significantly improving work efficiency and the consistency of calculation results. Summary of the Invention

[0003] The technical problem to be solved by the present invention is to overcome the above-mentioned problems existing in the prior art and to provide an automatic extraction method for the finite element calculation results curves of the axial and circumferential temperature fields of engine cylinders.

[0004] It should be noted that, in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article, or apparatus.

[0005] To solve the above-mentioned technical problems, the present invention is implemented using the following technical solution:

[0006] An automatic extraction method for finite element method results curves of engine cylinder temperature field includes the following steps:

[0007] Step 1: Extract standardized rules and define a coordinate system based on the finite element calculation results of the axial and circumferential temperature fields of the engine cylinder.

[0008] Step 2: Preprocess the engine cylinder geometry;

[0009] Step 3: Perform finite element modeling on the engine cylinder;

[0010] Step 4: Perform finite element simulation of the engine block temperature field;

[0011] Step 5: Use HyperView software to read the finite element simulation calculation result file;

[0012] Step 6: Calculate the location of key nodes and obtain their numbers;

[0013] Step 7: Extract the curves from the finite element calculation of the axial temperature field of the engine cylinder.

[0014] Step 8: Extract the curve of the finite element calculation result of the circumferential temperature field of the engine cylinder.

[0015] Furthermore, the definition of the coordinate system described in step one specifically includes:

[0016] The point where the center line of cylinder bore 1 intersects the horizontal plane containing the center line of main bearing journal is defined as the origin of the rectangular coordinate system. The direction from cylinder 1 to cylinder 2 along the crankshaft is defined as the X-axis of the rectangular coordinate system, and the direction from the main bearing cap to the cylinder block is defined as the Z-axis of the rectangular coordinate system. The Y-axis of the rectangular coordinate system is determined according to the right-hand rule.

[0017] Furthermore, the preprocessing of the engine cylinder geometry described in step two specifically includes the following steps:

[0018] S1. The distance from the top dead center of the engine piston to the top of the cylinder is L, and the piston stroke is S. Then the distance from the bottom dead center of the engine piston to the top of the cylinder is L+S. Using the XY plane as the cutting plane, the engine cylinder geometry is cut at the top dead center and bottom dead center positions of the engine piston.

[0019] S2. Using the XZ plane and YZ plane as cutting planes respectively, cut the engine cylinder geometry through the engine cylinder geometry center point.

[0020] Furthermore, the finite element modeling of the engine cylinder in step three specifically includes: placing the finite element mesh of the cylinder surface into the same group and naming it F-LINER: the distribution of the finite element mesh nodes on the engine cylinder surface satisfies that the distribution of nodes in each XY plane is completely consistent, excluding intermediate nodes.

[0021] Furthermore, step four involves performing finite element simulation of the engine block temperature field; specifically, this includes defining the material properties of the finite element model, applying loads and boundary conditions to the finite element model, and performing finite element analysis.

[0022] Furthermore, step six, which involves calculating the locations of key nodes and obtaining their numbers, specifically includes operating on the finite element mesh group F-LINER on the cylinder surface:

[0023] S1. Obtain the coordinates of the node with the maximum temperature of each cylinder barrel (XTM) j YTM j ZTM j ), j is an integer of 1, 2, 3...n, where n is the number of cylinders;

[0024] S2. Calculate the location of the critical node at the top of the cylinder and obtain the critical node number; select the node with the largest Z-value and a Y-value equal to 0 (i.e., Z = Z). max There are 2n nodes (x, y = 0), where n is the number of cylinders. These 2n nodes are arranged in ascending order of their x-coordinates, and the coordinate of the i-th node is (x, y = 0). i ,0,Z max Let i be an integer from 1 to 2n; then the top radius of the j-th cylinder is RT. j =(XT 2j -XT 2j-1 ) / 2, where j is an integer from 1 to n; then the coordinates of the critical node at the top of the j-th cylinder are KP. topj-1 (XT 2j ,0,Z max ), KP topj-2 (XT 2j -RT j RT j Z max ), KP topj-3 (XT 2j-1 ,0,Z max ), KP topj-4 (XT 2j -RT j -RT j Z max Based on the coordinates of the critical node at the top of the j-th cylinder, find and obtain the critical node KP at the top of the j-th cylinder. topj-1 KP topj-2 KP topj-3 KP topj-4 The node numbers are KPN in sequence. topj-1 KPN topj-2 KPN topj-3 KPN topj-4 ;

[0025] S3. Calculate the location of key nodes at the bottom of the cylinder and obtain their numbers; select the node with the smallest Z-value and a Y-value equal to 0 (i.e., Z = Z). minThere are 2n nodes (x, y = 0), where n is the number of cylinders. These 2n nodes are arranged in ascending order of their x-coordinates, and the coordinate of the i-th node is (x, y = 0). i ,0,Z min Let i be an integer from 1 to 2n; then the bottom radius of the j-th cylinder is RB. j =(XB 2j -XB 2j-1 ) / 2, where j is an integer from 1 to n; then the coordinates of the critical node at the bottom of the j-th cylinder are KP. botj-1 (XB 2j ,0,Z min ), KP botj-2 (XB 2j -RB j RB j Z min ), KP botj-3 (XB 2j-1 ,0,Z min ), KP botj-4 (XB 2j -RB j -RB j Z min Based on the coordinates of the critical node at the bottom of the j-th cylinder, find and obtain the critical node KP at the bottom of the j-th cylinder. botj-1 KP botj-2 KP botj-3 KP botj-4 The node numbers are KPN in sequence. botj-1 KPN botj-2 KPN botj-3 KPN botj-4 ;

[0026] S4. Calculate the location of the critical node at the top dead center of the cylinder piston and obtain the critical node number; select the node coordinate value Z as Z. max -L, Y value equals 0 (i.e. Z = Z) max There are 2n nodes in total, where n is the number of cylinders; these 2n nodes are arranged in ascending order of their X-coordinates, and the coordinate of the i-th node is (XTP). i ,0,Z max -L), where i is an integer from 1 to 2n; then the radius of the piston top dead center position of the j-th cylinder is RTP. j =(XTP) 2j -XTP 2j-1 ) / 2, where j is an integer from 1 to n; then the coordinates of the key node at the top dead center of the piston in the j-th cylinder are:

[0027] KP ptopj-1 (XTP 2j ,0,Z max -L), KP ptopj-2 (XTP 2j -RTP j RTP j Z max -L), KP ptopj-3 (XTP 2j-1 ,0,Z max -L), KP ptopj-4 (

[0028] XTP 2j -RTP j -RTP j Z max -L); Based on the coordinates of the key node at the top dead center position of the piston in the j-th cylinder, find and obtain the key node KP at the top dead center position of the piston in the j-th cylinder. ptopj-1 KP ptopj-2 KP ptopj-3 KP ptopj-4 The node numbers are KPN in sequence. ptopj-1 KPN ptopj-2 KPN ptopj-3 KPN ptopj-4 ;

[0029] S5. Calculate the location of the critical node at the bottom dead center of the cylinder piston and obtain the critical node number; select the node coordinate value Z as Z. max -LS, Y value equals 0 (i.e. Z = Z) max There are 2n nodes in total, where n is the number of cylinders; these 2n nodes are arranged in ascending order of their X-coordinates, and the coordinate of the i-th node is (XBP). i ,0,Z max -LS), where i is an integer from 1 to 2n; then the radius of the piston bottom dead center position of the j-th cylinder is RBP. j =(XBP) 2j -XBP 2j-1 ) / 2, where j is an integer from 1 to n; then the coordinates of the key node at the bottom dead center of the piston in the j-th cylinder are:

[0030] KP pbomj-1 (XBP 2j ,0,Z max -LS), KP pbomj-2 (XBP 2j -RBP j RBP j Z max -LS), KPpbomj-3 (XBP 2j-1 ,0,Z max -LS),

[0031] KP pbomj-4 (XBP 2j -RBP j -RBP j Z max -LS); Based on the coordinates of the key node at the bottom dead center of the piston in the j-th cylinder, find and obtain the key node KP at the bottom dead center of the piston in the j-th cylinder. pbomj-1 KP pbomj-2 KP pbomj-3 KP pbomj-4 The node numbers are KPN in sequence. pbomj-1 KPN pbomj-2 KPN pbomj-3 KPN pbomj-4 ;

[0032] S6. Calculate the location of the critical node at the highest cylinder temperature and obtain the critical node number; for the j-th cylinder, select the node coordinate Z value as ZTM. j The Y value is equal to 0 (i.e., Z = ZTM). j There are 2n nodes (x, y = 0), where n is the number of cylinders and j is an integer from 1 to n. These 2n nodes are arranged in ascending order of their x-coordinates, and the coordinate of the i-th node is (x, y = 0). i ,0,ZTM j Let i be an integer from 1 to 2n; then the radius of the highest temperature point in the j-th cylinder is RTMP. j =(XTMP) 2j -XTMP 2j-1 ) / 2, where j is an integer from 1 to n; then the coordinates of the key node at the highest temperature point of the j-th cylinder are KP. tmpj-1 (XTMP 2j ,0,ZTM j ), KP tmpj-2 (XTMP 2j -RTMP j RTMP j ZTM j ), KP tmpj-3 (XTMP 2j-1 ,0,ZTM j ), KP tmpj-4 (XTMP 2j -RTMP j -RTMP j ZTM jBased on the coordinates of the key node at the highest temperature point of the j-th cylinder, find and obtain the key node KP at the highest temperature point of the j-th cylinder. tmpj-1 KP tmpj-2 KP tmpj-3 KP tmpj-4 The node numbers are KPN in sequence. tmpj-1 KPN tmpj-2 KPN tmpj-3 KPN tmpj-4 .

[0033] Furthermore, the extraction of the finite element calculation result curve of the axial temperature field of the engine cylinder in step seven, for the j-th cylinder (where j is an integer from 1 to n), includes the following steps:

[0034] S1. In Hyperview software, access the annotations-measures module, add measure groups, select node path, and enter KPN sequentially. topj-1 and KPN botj-1 Select "create curves", select "value" for Y axis and "distance" for X axis. Select "new plot" for "place on", select "ok" to generate curves, change "header" to "LINERj", change "label" in "curvelabel" to "0deg", change "Y axis label" to "T℃", and "X axis label" to "distance" (unit: mm). Use "file-export-curves", select "ok" to save the curve data, and use "capture frame area" to save the curve plot.

[0035] S2. In Hyperview software, access the annotations-measures module, add measure groups, select node path, and enter KPN sequentially. topj-2 and KPN botj-2 Select "create curves", select "value" for Y axis and "distance" for X axis. Select "new plot" for "place on", select "ok" to generate curves, change "header" to "LINERj", change "label" in "curvelabel" to "90deg", change "Y axis label" to "T℃", and change "X axis label" to "distance" in mm. Use "file-export-curves", select "ok" to save the curve data, and use "capture frame area" to save the curve plot.

[0036] S3. In Hyperview software, access the annotations-measures module, add measure groups, select node path, and enter KPN sequentially. topj-3 and KPN botj-3 Select "create curves", select "value" for Y axis and "distance" for X axis. Select "new plot" for "place on", select "ok" to generate curves, change "header" to "LINERj", change "label" in "curvelabel" to "180deg", change "Y axis label" to "T℃", and change "X axis label" to "distance" in mm. Use "file-export-curves", select "ok" to save the curve data, and use "capture framearea" to save the curve plot.

[0037] S4. In Hyperview software, access the annotations-measures module, add measure groups, select node path, and enter KPN sequentially. topj-4 and KPN botj-4 Select "Create Curves", choose "Value" for Y axis and "Distance" for X axis. Select "New Plot" for "Place On", and click "OK" to generate the curve. Change the header to "LINERj", change the label item in "Curve Label" to "270deg", change the Y axis label to "T℃", and change the X axis label to "Distance" in mm. Use "File-Export-Curves", click "OK" to save the curve data, and use "Capture Framearea" to save the curve plot.

[0038] Further, the finite element calculation result curve of the circumferential temperature field of the engine cylinder described in step eight is extracted; for the j-th cylinder, the coordinates of the center of the piston top dead center position are (XTO). j =(XT 2j +XT 2j-1 ) / 2, 0, Z max -L), the coordinates of the center of the cylinder piston bottom dead center position are (XBO) j =(XB 2j +XB 2j-1 ) / 2, 0, Z max -LS), the coordinates of the center of the cylinder's highest temperature point are (XTO). j =(XTMP) 2j +XTMP2j-1 ) / 2, 0, ZTM j (), where j is an integer of 1, 2, 3...n, specifically including the following steps:

[0039] S1. In Hyperview software, access the annotations-measures module, add measure groups, select node path, and enter KPN sequentially. ptopj-1 KPN ptopj-2 KPN ptopj-3 KPN ptopj-4 and KPN ptopj-1 Arrange the nodes in the order of the path, and retrieve the node numbers and coordinates of all nodes on the path. There are m nodes in total, and the coordinates of the k-th node are KTP. k (XKTP k YKTP k Z max -L), where k is an integer from 1 to m, then the angle between the straight line connecting the kth node and the center of the circle at the top dead center of the cylinder piston and the positive x-axis is α. top-k ;

[0040] Select "Create Curves," choose "Value" for the Y-axis and "Distance" for the X-axis. Choose "New Plot" for "Place On," and select "OK" to generate the curve. Change the header to "LINERj," and in the curve label section, change the label to "TDC," the Y-axis label to "T℃," and the X-axis label to "Deg." In the data values ​​section, replace the X-values ​​with the angle α between the node and the X-axis, according to the node order. top-k Use file-export-curves, select OK to save the curve data, and use captureframe area to save the curve graph;

[0041] S2. In Hyperview software, access the annotations-measures module, add measure groups, select node path, and enter KPN sequentially. pbomj-1 KPN pbomj-2 KPN pbomj-3 KPN pbomj-4 and KPN pbomj-1 Arrange the nodes in the order of the path, retrieve the node numbers and coordinates of all nodes on the path. There are m nodes in total, and the coordinates of the k-th node are KBP. k (XKBP k YKBP k Z maxLet k be an integer from 1 to m, and k be the number of nodes. Then the angle between the straight line connecting the kth node to the center of the circle at the bottom dead center of the cylinder piston and the positive x-axis is α. bom-k ;

[0042] Select "Create Curves," choose "Value" for the Y-axis and "Distance" for the X-axis. Choose "New Plot" for "Place On," and select "OK" to generate the curve. Change the header to "LINERj," and in the curve label section, change the label to "BDC," the Y-axis label to "T℃," and the X-axis label to "Deg." In the data values ​​section, replace the X-values ​​with the angle α between the node and the X-axis, according to the node order. bom-k Use file-export-curves, select OK to save the curve data, and use capture frame area to save the curve graph;

[0043] S3. In Hyperview software, access the annotations-measures module, add measure groups, select node path, and enter KPN sequentially. tmpj-1 KPN tmpj-2 KPN tmpj-3 KPN tmpj-4 and KPN tmpj-1 Arrange the nodes in the order of the path, retrieve the node numbers and coordinates of all nodes on the path. There are m nodes in total, and the coordinates of the k-th node are KMP. k (XKMP k YKMP k ZTM j Let k be an integer from 1 to m. Then, the angle between the straight line connecting the k-th node and the center of the circle where the cylinder reaches its highest temperature point and the positive x-axis is α. max-k ;

[0044] Select "Create Curves," choose "Value" for the Y-axis and "Distance" for the X-axis. Choose "New Plot" for "Place On," and select "OK" to generate the curve. Change the header to "LINERj," and in the curve label section, change the label to "MAXTEMP." Change the Y-axis label to "T℃" and the X-axis label to "Deg." In the data values ​​section, replace the X-values ​​with the angle α between the node and the X-axis, according to the node order. max-kUse file-export-curves, select OK to save the curve data, and use capture frame area to save the curve graph.

[0045] Furthermore, the angle with the positive x-axis is α. top-k ;

[0046] ①When XKTP k -XTO j ≥0 and YKTP k When ≥0, α top-k =arctan[YKTP k / (XKTP k -XTO j )] / π*180;

[0047] ②When XKTP k -XTO j ≥0 and YKTP k When <0, α top-k ={arctan[YKTP] k / (XKTP k -XTO j )]+2π} / π*180;

[0048] ③ When XKTP k -XTO j <0 and YKTP k ≥0 or XKTP k -XTO j <0 and YKTP k When <0, α top-k ={arctan[YKTP] k / (XKTP k -XTO j )]+π} / π*180;

[0049] The angle between the x-axis and the positive x-axis is α. bom-k ;

[0050] ① When XKBP k XBO j ≥0 and YKBP k When ≥0, α bom-k =arctan[YKBP k / (XKBP k XBO j )] / π*180;

[0051] ②When XKBP k XBO j ≥0 and YKBP k When <0, αbom-k ={arctan[YKBP k / (XKBP k XBO j )]+2π} / π*180;

[0052] ③ When XKBP k XBO j <0 and YKBP k ≥0 or XKBP k XBO j <0 and YKBP k When <0, α bom-k ={arctan[YKBP k / (XKBP k XBO j )]+π} / π*180;

[0053] The angle between the x-axis and the positive x-axis is α. max-k :

[0054] ① When XKMP k -XTO j ≥0 and YKMP k When ≥0, α max-k =arctan[YKMP k / (XKMP k -XTO j )] / π*180;

[0055] ② When XKMP k -XTO j ≥0 and YKMP k When <0, α max-k ={arctan[YKMP k / (XKMP k -XTO j )]+2π} / π*180;

[0056] ③ When XKMP k -XTO j <0 and YKMP k ≥0 or XKMP k -XTO j <0 and YKMP k When <0, α max-k ={arctan[YKMP k / (XKMP k -XTO j )]+π} / π*180.

[0057] A computer-readable storage medium having a computer program stored thereon that implements the method described above.

[0058] Compared with the prior art, the beneficial effects of the present invention are:

[0059] This invention uses the TCL scripting language to customize standardized rules for extracting the finite element calculation results curves of the axial and circumferential temperature fields of the engine cylinder. Through programming, the extraction of the finite element calculation results curves of the axial and circumferential temperature fields of the engine cylinder is automated. Compared with the traditional manual work cycle, the work efficiency is shortened by more than 95%, which greatly improves the work efficiency and the consistency of the calculation results. Attached Figure Description

[0060] The invention will now be further described with reference to the accompanying drawings:

[0061] Figure 1 Standardization rules and flowcharts for extracting finite element calculation results curves of axial and circumferential temperature fields of engine cylinders;

[0062] Figure 2a Define the origin of the rectangular coordinate system Figure 1 ;

[0063] Figure 2b Figure 2 defines the origin of the rectangular coordinate system;

[0064] Figure 3 To cut the cylinder geometry at the top dead center and bottom dead center positions of the engine piston;

[0065] Figure 4 The cylinder geometry is cut using the XZ and YZ planes as cutting planes and the cylinder geometric center point as the reference point.

[0066] Figure 5a Finite element modeling of engine cylinders Figure 1 ;

[0067] Figure 5b Figure 2 shows the finite element model of the engine cylinder.

[0068] Figure 6 This is a schematic diagram of the node with the maximum temperature.

[0069] Figure 7a Diagram showing the location and numbering of key nodes at the top of the cylinder. Figure 1 ;

[0070] Figure 7b Schematic diagram 2 showing the location and numbering of key nodes at the top of the cylinder;

[0071] Figure 8 This is a schematic diagram showing the node numbers and coordinates of all nodes at the top dead center position of the cylinder. Detailed Implementation

[0072] To make the objectives, technical solutions, and advantages of this invention clearer, the technical solutions of the embodiments of this invention will be described in more detail below with reference to the accompanying drawings. In the drawings, the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The described embodiments are some, but not all, embodiments of this invention. The embodiments described below with reference to the accompanying drawings are exemplary and intended to explain this invention, and should not be construed as limiting the invention. All other embodiments obtained by those skilled in the art based on the embodiments of this invention without creative effort are within the scope of protection of this invention. The embodiments of this invention will be described in detail below with reference to the accompanying drawings.

[0073] S1. Establish standardized rules for extracting finite element calculation result curves of the engine cylinder bore axial and circumferential temperature fields, and design the program flowchart as follows: Figure 1 As shown.

[0074] S11. Define the coordinate system: Define the point where the horizontal plane containing the center line of cylinder bore 1 and the center line of the main journal is located as the origin of the rectangular coordinate system, such as... Figure 2a , Figure 2b As shown, the direction from cylinder 1 to cylinder 2 along the crankshaft is defined as the X-axis of the rectangular coordinate system, the direction from the main bearing cap to the cylinder block is defined as the Z-axis of the rectangular coordinate system, and the Y-axis of the rectangular coordinate system is determined according to the right-hand rule.

[0075] S12. Preprocess the engine cylinder geometry:

[0076] S121. The distance from the top dead center of the engine piston to the top of the cylinder is L, and the piston stroke is S. Therefore, the distance from the bottom dead center of the engine piston to the top of the cylinder is L+S. Using the XY plane as the cutting plane, the engine cylinder geometry is cut at the top dead center and bottom dead center positions of the engine piston. Figure 3 As shown;

[0077] S122. Using the XZ and YZ planes as cutting planes respectively, cut the engine cylinder geometry through the engine cylinder geometry center point, as follows: Figure 4 As shown.

[0078] S13. Perform finite element modeling on the engine cylinder bore, placing the finite element meshes on the cylinder bore surface into the same group, named F-LINER: the node distribution of the finite element meshes on the engine cylinder bore surface must satisfy the condition that the node distribution in each XY plane (i.e., the node coordinates Z values ​​are the same) is completely consistent (excluding intermediate nodes), such as... Figure 5a , Figure 5b As shown.

[0079] S14. Perform finite element simulation of engine cylinder block temperature field: This includes defining the material properties of the finite element model, applying loads and boundary conditions to the finite element model, and performing finite element analysis.

[0080] S15. Use HyperView software to read the finite element simulation calculation result file.

[0081] S16. Calculate the location of key nodes and obtain their numbers: Perform the following operations on the finite element mesh group F-LINER on the cylinder surface.

[0082] S161. Obtain the coordinates of the node with the maximum temperature of each cylinder barrel (XTM). j YTM j ZTM j ), j is an integer of 1, 2, 3...n, and n is the number of cylinders.

[0083] S162. Calculate the location of the critical node at the top of the cylinder and obtain the critical node number. Select the node with the largest Z-value and a Y-value equal to 0 (i.e., Z = Z). max There are 2n nodes (x, y = 0), where n is the number of cylinders. These 2n nodes are arranged in ascending order of their x-coordinates, and the coordinate of the i-th node is (x, y = 0). i ,0,Z max Let i be an integer from 1 to 2n; then the top radius of the j-th cylinder is RT. j =(XT 2j -XT 2j-1 ) / 2, where j is an integer from 1 to n; then the coordinates of the critical node at the top of the j-th cylinder are KP. topj-1 (XT 2j ,0,Z max ), KP topj-2 (XT 2j -RT j RT j Z max ), KP topj-3 (XT 2j-1 ,0,Z max ), KP topj-4 (XT 2j -RT j -RT j Z max Based on the coordinates of the critical node at the top of the j-th cylinder, find and obtain the critical node KP at the top of the j-th cylinder. topj-1 KP topj-2 KP topj-3 KP topj-4 The node numbers are KPN in sequence. topj-1 KPNtopj-2 KPN topj-3 KPN topj-4 .

[0084] S163. Calculate the location of key nodes at the bottom of the cylinder and obtain their numbers. Select nodes with the smallest Z-value and a Y-value equal to 0 (i.e., Z = Z). min There are 2n nodes (x, y = 0), where n is the number of cylinders. These 2n nodes are arranged in ascending order of their x-coordinates, and the coordinate of the i-th node is (x, y = 0). i ,0,Z min Let i be an integer from 1 to 2n; then the bottom radius of the j-th cylinder is RB. j =(XB 2j -XB 2j-1 ) / 2, where j is an integer from 1 to n; then the coordinates of the critical node at the bottom of the j-th cylinder are KP. botj-1 (XB 2j ,0,Z min ), KP botj-2 (XB 2j -RB j RB j Z min ), KP botj-3 (XB 2j-1 ,0,Z min ), KP botj-4 (XB 2j -RB j -RB j Z min Based on the coordinates of the critical node at the bottom of the j-th cylinder, find and obtain the critical node KP at the bottom of the j-th cylinder. botj-1 KP botj-2 KP botj-3 KP botj-4 The node numbers are KPN in sequence. botj-1 KPN botj-2 KPN botj-3 KPN botj-4 .

[0085] S164. Calculate the location of the critical node at the top dead center of the cylinder piston and obtain the critical node number. Select the node coordinate value Z as Z. max -L, Y value equals 0 (i.e. Z = Z) max There are 2n nodes in total, where n is the number of cylinders; these 2n nodes are arranged in ascending order of their X-coordinates, and the coordinate of the i-th node is (XTP). i ,0,Z max-L), where i is an integer from 1 to 2n; then the radius of the piston top dead center position of the j-th cylinder is RTP. j =(XTP) 2j -XTP 2j-1 ) / 2, where j is an integer from 1 to n; then the coordinates of the key node at the top dead center of the piston in the j-th cylinder are:

[0086] KP ptopj-1 (XTP 2j ,0,Z max -L), KP ptopj-2 (XTP 2j -RTP j RTP j Z max -L), KP ptopj-3 (XTP 2j-1 ,0,Z max -L), KP ptopj-4 (

[0087] XTP 2j -RTP j -RTP j Z max -L); Based on the coordinates of the key node at the top dead center position of the piston in the j-th cylinder, find and obtain the key node KP at the top dead center position of the piston in the j-th cylinder. ptopj-1 KP ptopj-2 KP ptopj-3 KP ptopj-4 The node numbers are KPN in sequence. ptopj-1 KPN ptopj-2 KPN ptopj-3 KPN ptopj-4 .

[0088] S165. Calculate the location of the critical node at the bottom dead center of the cylinder piston and obtain the critical node number. Select the node coordinate value Z as Z. max -LS, Y value equals 0 (i.e. Z = Z) max There are 2n nodes in total, where n is the number of cylinders; these 2n nodes are arranged in ascending order of their X-coordinates, and the coordinate of the i-th node is (XBP). i ,0,Z max -LS), where i is an integer from 1 to 2n; then the radius of the piston bottom dead center position of the j-th cylinder is RBP. j =(XBP) 2j -XBP 2j-1 ) / 2, where j is an integer from 1 to n; then the coordinates of the key node at the bottom dead center of the piston in the j-th cylinder are:

[0089] KP pbomj-1 (XBP 2j ,0,Z max -LS), KP pbomj-2 (XBP 2j -RBP j RBP j Z max -LS), KP pbomj-3 (XBP 2j-1 ,0,Z max -LS),

[0090] KP pbomj-4 (XBP 2j -RBP j -RBP j Z max -LS); Based on the coordinates of the key node at the bottom dead center of the piston in the j-th cylinder, find and obtain the key node KP at the bottom dead center of the piston in the j-th cylinder. pbomj-1 KP pbomj-2 KP pbomj-3 KP pbomj-4 The node numbers are KPN in sequence. pbomj-1 KPN pbomj-2 KPN pbomj-3 KPN pbomj-4 .

[0091] S166. Calculate the location of the critical node at the highest cylinder temperature and obtain the critical node number. For the j-th cylinder, select the node coordinate Z value as ZTM. j The Y value is equal to 0 (i.e., Z = ZTM). j There are 2n nodes (x, y = 0), where n is the number of cylinders and j is an integer from 1 to n. These 2n nodes are arranged in ascending order of their x-coordinates, and the coordinate of the i-th node is (x, y = 0). i ,0,ZTM j Let i be an integer from 1 to 2n; then the radius of the highest temperature point in the j-th cylinder is RTMP. j =(XTMP) 2j -XTMP 2j-1 ) / 2, where j is an integer from 1 to n; then the coordinates of the key node at the highest temperature point of the j-th cylinder are KP. tmpj-1 (XTMP 2j ,0,ZTM j ), KP tmpj-2 (XTMP 2j -RTMP j RTMP j ZTMj ), KP tmpj-3 (XTMP 2j-1 ,0,ZTM j ), KP tmpj-4 (XTMP 2j -RTMP j -RTMP j ZTM j Based on the coordinates of the key node at the highest temperature point of the j-th cylinder, find and obtain the key node KP at the highest temperature point of the j-th cylinder. tmpj-1 KP tmpj-2 KP tmpj-3 KP tmpj-4 The node numbers are KPN in sequence. tmpj-1 KPN tmpj-2 KPN tmpj-3 KPN tmpj-4 .

[0092] S17. Extraction of finite element calculation results curves of axial temperature field of engine cylinder: For the j-th cylinder, j is an integer from 1 to n.

[0093] S171. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. topj-1 and KPN botj-1 Select "create curves", select "value" for the Y axis and "distance" for the X axis. Select "new plot" for "place on", select "ok" to generate the curve, change the header to "LINERj", change the label item in the curve label to "0deg", change the Y axis label to "T (°C)" and the X axis label to "distance (mm)", use "file-export-curves", select "ok" to save the curve data, and use "captureframe area" to save the curve plot.

[0094] S172. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. topj-2 and KPN botj-2Select "create curves", select "value" for Y axis and "distance" for X axis. Select "new plot" for "place on", select "ok" to generate curves, change "header" to "LINERj", change "label" in "curve label" to "90deg", change "Y axis label" to "T (°C)" and "X axis label" to "distance (mm)", use "file-export-curves", select "ok" to save the curve data, and use "captureframe area" to save the curve plot.

[0095] S173. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. topj-3 and KPN botj-3 Select "create curves", select "value" for Y axis and "distance" for X axis. Select "new plot" for "place on", select "ok" to generate curves, change "header" to "LINERj", change "label" in "curve label" to "180deg", change "Y axis label" to "T (°C)" and "X axis label" to "distance (mm)", use "file-export-curves", select "ok" to save the curve data, and use "captureframe area" to save the curve plot.

[0096] S174. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. topj-4 and KPN botj-4 Select "create curves", choose "value" for Y axis and "distance" for X axis. Select "new plot" for "place on", select "ok" to generate curves, change "header" to "LINERj", change the "label" item in "curve label" to "270deg", change "Y axis label" to "T (°C)" and "X axis label" to "distance (mm)", use "file-export-curves", select "ok" to save the curve data, and use "captureframe area" to save the curve plot.

[0097] S18. Extraction of finite element calculation results of the circumferential temperature field of the engine cylinder: For the j-th cylinder, the coordinates of the center of the piston top dead center are (XTO). j =(XT 2j +XT 2j-1 ) / 2, 0, Z max -L), the coordinates of the center of the cylinder piston bottom dead center position are (XBO) j =(XB 2j +XB 2j-1 ) / 2, 0, Z max -LS), the coordinates of the center of the cylinder's highest temperature point are (XTO). j =(XTMP) 2j +XTMP 2j-1 ) / 2, 0, ZTM j ), where j is an integer of 1, 2, 3...n;

[0098] S181. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. ptopj-1 KPN ptopj-2 KPN ptopj-3 KPN ptopj-4 and KPN ptopj-1 Arrange the nodes in the order of the path, and retrieve the node numbers and coordinates of all nodes on the path. There are m nodes in total, and the coordinates of the k-th node are KTP. k (XKTP k YKTP k Z max -L), k is an integer from 1, 2, 3...m, then the angle α between the straight line connecting the kth node and the center of the circle at the top dead center of the cylinder piston and the positive x-axis. top-k For ① when XKTP k -XTO j ≥0 and YKTP k When ≥0, α top-k =arctan[YKTP k / (XKTP k -XTO j )] / π*180;②When XKTP k -XTO j ≥0 and YKTP k When <0, α top-k ={arctan[YKTP] k / (XKTP k -XTO j )]+2π} / π*180;③When XKTP k -XTOj <0 and YKTP k ≥0 or XKTP k -XTO j <0 and YKTP k When <0, α top-k ={arctan[YKTP] k / (XKTP k -XTO j )]+π} / π*180. Select "Create Curves", choose "Value" for Y-axis and "Distance" for X-axis. Select "New Plot" for "Place On", select "OK" to generate the curve. Change the header to "LINERj", change the label in "Curve Label" to "TDC", change the Y-axis label to "T (°)", and change the X-axis label to "Deg (°)". In "Data Values", replace the X-values ​​with the angle α between the node and the X-axis, according to the node order. top-k Use file-export-curves, select OK to save the curve data, and use capture frame area to save the curve graph;

[0099] S182. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. pbomj-1 KPN pbomj-2 KPN pbomj-3 KPN pbomj-4 and KPN pbomj-1 Arrange the nodes in the order of the path, and retrieve the node numbers and coordinates of all nodes on the path. There are m nodes in total, and the coordinates of the kth node are KBP. k (XKBP k YKBP k Z max -LS), k is an integer from 1 to m, then the angle α between the straight line connecting the kth node and the center of the circle at the bottom dead center of the cylinder piston and the positive x-axis. bom-k For ① when XKBP k XBO j ≥0 and YKBP k When ≥0, α bom-k =arctan[YKBP k / (XKBP k XBO j )] / π*180;②When XKBP k XBO j ≥0 and YKBPk When <0, α bom-k ={arctan[YKBP k / (XKBP k XBO j )]+2π} / π*180;③When XKBP k XBO j <0 and YKBP k ≥0 or XKBP k XBO j <0 and YKBP k When <0, α bom-k ={arctan[YKBP k / (XKBP k XBO j )]+π} / π*180. Select "Create Curves", choose "Value" for Y-axis and "Distance" for X-axis. Select "New Plot" for "Place On", select "OK" to generate the curve. Change the header to "LINERj", change the label in "Curve Label" to "BDC", change the Y-axis label to "T (°)", and change the X-axis label to "Deg (°)". In "Data Values", replace the X-values ​​with the angle α between the node and the X-axis, according to the node order. bom-k Use file-export-curves, select OK to save the curve data, and use capture frame area to save the curve graph;

[0100] S183. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. tmpj-1 KPN tmpj-2 KPN tmpj-3 KPN tmpj-4 and KPN tmpj-1 Arrange the nodes in the order of the path, retrieve the node numbers and coordinates of all nodes on the path. There are m nodes in total, and the coordinates of the k-th node are KMP. k (XKMP k YKMP k ZTM j Let k be an integer from 1 to m. Then, the angle α between the straight line connecting the k-th node and the center of the circle where the cylinder reaches its highest temperature and the positive x-axis is calculated. max-k For ① when XKMP k -XTO j ≥0 and YKMP k When ≥0, αmax-k =arctan[YKMP k / (XKMP k -XTO j )] / π*180;②When XKMP k -XTO j ≥0 and YKMP k When <0, α max-k ={arctan[YKMP k / (XKMP k -XTO j )]+2π} / π*180;③When XKMP k -XTO j <0 and YKMP k ≥0 or XKMP k -XTO j <0 and YKMP k When <0, α max-k ={arctan[YKMP k / (XKMP k -XTO j )]+π} / π*180. Select "Create Curves", choose "Value" for Y-axis and "Distance" for X-axis. Select "New Plot" for "Place On", select "OK" to generate the curve. Change the header to "LINERj", change the label in "Curve Label" to "MAX TEMP", change the Y-axis label to "T" (°C), and change the X-axis label to "Deg" (°). In "Data Values", replace the X-values ​​with the angle α between the node and the X-axis, according to the node order. max-k Use file-export-curves, select OK to save the curve data, and use capture frame area to save the curve graph;

[0101] S2. Programming and programming of processes S16-S18 based on TCL scripting language.

[0102] S3. Adjust the geometric coordinate system to satisfy the coordinate system defined in S11: The point where the horizontal plane containing the center line of cylinder bore 1 and the center line of the main journal intersects is defined as the origin of the rectangular coordinate system, such as... Figure 2a , Figure 2b As shown, the X-axis of the rectangular coordinate system is located along the crankshaft from cylinder 1 to cylinder 2, and the Z-axis is located along the main bearing cap to the cylinder block. The Y-axis of the rectangular coordinate system is determined according to the right-hand rule.

[0103] S4. Preprocess the engine cylinder geometry:

[0104] S41. The distance from the top dead center of the engine piston to the top of the cylinder is L, and the piston stroke is S. Therefore, the distance from the bottom dead center of the engine piston to the top of the cylinder is L+S. Using the XY plane as the cutting plane, cut the engine cylinder geometry at the top dead center and bottom dead center positions of the engine piston. Figure 3 As shown;

[0105] S42. Using the XZ and YZ planes as cutting planes respectively, cut the engine cylinder geometry through the engine cylinder geometry center point, as follows: Figure 4 As shown.

[0106] S5. Perform finite element modeling on the engine cylinder bore, placing the finite element meshes on the cylinder bore surface into the same group, named F-LINER: the node distribution of the finite element meshes on the engine cylinder bore surface must satisfy the condition that the node distribution in each XY plane (i.e., the node coordinates Z values ​​are the same) is completely consistent (excluding intermediate nodes), such as... Figure 5a , Figure 5b As shown.

[0107] S6. Perform finite element simulation of engine cylinder block temperature field: This includes defining the material properties of the finite element model, applying loads and boundary conditions to the finite element model, and performing finite element analysis.

[0108] S7. Use HyperView software to read the finite element simulation calculation result file.

[0109] S8. Call the automatic extraction program for the finite element calculation results of the engine cylinder axial and circumferential temperature fields, which was compiled in S2.

[0110] S81. Use the program to automatically calculate the location of key nodes and obtain the key node number.

[0111] S82. The finite element calculation results curve of the axial temperature field of the engine cylinder is automatically extracted using the program.

[0112] S83. The program automatically extracts the finite element calculation result curve of the circumferential temperature field of the engine cylinder.

[0113] A readable storage medium includes the method for automatically extracting the finite element result curve of the engine cylinder temperature field. Specific implementation examples:

[0115] For engine structures with two cylinders

[0116] S1. Establish standardized rules for extracting finite element calculation result curves of the engine cylinder bore axial and circumferential temperature fields, and design the program flowchart as follows: Figure 1 As shown.

[0117] S11. Define the coordinate system: Define the point where the horizontal plane containing the center line of cylinder bore 1 and the center line of the main journal is located as the origin of the rectangular coordinate system, such as... Figure 2a , Figure 2b As shown, the direction from cylinder 1 to cylinder 2 along the crankshaft is defined as the X-axis of the rectangular coordinate system, the direction from the main bearing cap to the cylinder block is defined as the Z-axis of the rectangular coordinate system, and the Y-axis of the rectangular coordinate system is determined according to the right-hand rule.

[0118] S12. Preprocess the engine cylinder geometry:

[0119] S121. The distance from the top dead center of the engine piston to the top of the cylinder is L = 10mm, and the piston stroke is S = 40mm. Therefore, the distance from the bottom dead center of the engine piston to the top of the cylinder is L + S = 50mm. Using the XY plane as the cutting plane, cut the engine cylinder geometry at the top dead center and bottom dead center positions of the engine piston. Figure 3 As shown;

[0120] S122. Using the XZ and YZ planes as cutting planes respectively, cut the engine cylinder geometry through the engine cylinder geometry center point, as follows: Figure 4 As shown.

[0121] S13. Perform finite element modeling on the engine cylinder bore, placing the finite element meshes on the cylinder bore surface into the same group, named F-LINER: the node distribution of the finite element meshes on the engine cylinder bore surface must satisfy the condition that the node distribution in each XY plane (i.e., the node coordinates Z values ​​are the same) is completely consistent (excluding intermediate nodes), such as... Figure 5a , Figure 5b As shown.

[0122] S14. Perform finite element simulation of engine cylinder block temperature field: This includes defining the material properties of the finite element model, applying loads and boundary conditions to the finite element model, and performing finite element analysis.

[0123] S15. Taking HyperView software as an example, read the finite element simulation calculation result file.

[0124] S16. Calculate the locations of key nodes and obtain their numbers: See [link / reference] Figure 6 The following operations were performed on the finite element mesh group F-LINER on the cylinder surface:

[0125] S161. Obtain the coordinates of the node with the maximum temperature of each cylinder. The coordinates of the node with the maximum temperature of the first cylinder are (XTM1 = -15, YTM1 = -25.98, ZTM1 = 73.33) and the coordinates of the node with the maximum temperature of the second cylinder are (XTM2 = 85, YTM2 = -25.98, ZTM2 = 86.67).

[0126] S162. Calculate the location of the critical node at the top of the cylinder and obtain the critical node number. Select the node with the largest Z-value and a Y-value equal to 0 (i.e., Z = Z). max There are 4 nodes with X=110, Y=0; these 4 nodes are arranged in ascending order of their X-coordinates. The first node has coordinates (X=110, Y=0). max =110), the coordinates of the second node are (XT2=30, 0, Z). max =110), the coordinates of the third node are (XT3 = 40, 0, Z = 110), max =110), the coordinates of the 4th node are (XT4 = 100, 0, Z = 110), max =110); then the top radius of the first cylinder is RT1 = (XT2 - XT1) / 2 = 30, and the top radius of the second cylinder is RT2 = (XT4 - XT3) / 2 = 30; then the coordinates of the key node at the top of the first cylinder are KP. top1-1 (XT2 = 30, 0, Z) max =110), KP top1-2 (XT2-RT1=0,RT1=30,Z max =110), KP top1-3 (XT1 = -30, 0, Z) max =110), KP top1-4 (XT2-RT1=0,-RT1=-30,Z max =110), then the coordinates of the critical node at the top of the second cylinder are KP. top2-1 (XT4 = 100, 0, Z) max =110), KP top2-2 (XT4-RT2=70,RT2=30,Z max =110), KP top2-3 (XT3 = 40, 0, Z) max =110), KP top2-4 (XT4-RT2=70,-RT2=-30,Z max =110); Based on the coordinates of the key node at the top of the first cylinder, find and obtain the key node at the top of the first cylinder.

[0127] KP top1-1 KP top1-2 KP top1-3 KP top1-4 The node numbers are KPN in sequence. top1-1 =1, KPN top1-2 =2, KPN top1-3 =3, KPN top1-4=4, based on the coordinates of the critical node at the top of the second cylinder, find and obtain the critical node KP at the top of the second cylinder. top2-1 KP top2-2 KP top2-3 KP top2-4 The node numbers are KPN in sequence. top2-1 =5, KPN top2-2 =6, KPN top3-3 =7, KPN top2-4 =8, such as Figure 7a , Figure 7b As shown.

[0128] S163. Calculate the location of key nodes at the bottom of the cylinder and obtain their numbers. Select nodes with the smallest Z-value and a Y-value equal to 0 (i.e., Z = Z). min There are 4 nodes with X=50, Y=0. These 4 nodes are arranged in ascending order of their X-coordinates. The first node has coordinates (XB1 = ​​-30, 0, ZB1 = 0). min =50), the coordinates of the second node are (XB2=30, 0, Z =50). min =50), the coordinates of the third node are (XB3 = 40, 0, Z = 50), min =50), the coordinates of the 4th node are (XB4 = 100, 0, Z = 50). min =50); then the bottom radius of the first cylinder is RB1 = (XB2 - XB1) / 2 = 30, then the bottom radius of the second cylinder is RB2 = (XB4 - XB3) / 2 = 30; then the coordinates of the key node at the bottom of the first cylinder are KP. bot1-1 (XB2 = 30, 0, Z) min =50), KP bot1-2 (XB2-RB1=0,RB1=30,Z min =50), KP bot1-3 (XB1 = ​​-30, 0, Z) min =50), KP bot1-4 (XB2-RB1=0,-RB1=-30,Z min =50), then the coordinates of the critical node at the bottom of the second cylinder are KP. bot2-1 (XB4 = 100, 0, Z) min =50), KP bot2-2 (XB4-RB2=70,RB2=30,Z min =50), KP bot2-3 (XB3 = 40, 0, Z) min =50), KP bot2-4 (XB4-RB2=70,-RB2=-30,Z min=50); Based on the coordinates of the key node at the bottom of the first cylinder, find and obtain the key node KP at the bottom of the first cylinder. bot1-1 KP bot1-2 KP bot1-3 KP bot1-4 The node numbers are KPN in sequence. bot1-1 =9, KPN bot1-2 =10, KPN bot1-3 =11, KPN bot1-4 =12, based on the coordinates of the critical node at the bottom of the second cylinder, find and obtain the critical node KP at the bottom of the second cylinder. bot2-1 KP bot2-2 KP bot2-3 KP bot2-4 The node numbers are KPN in sequence. bot2-1 =13, KPN bot2-2 =14, KPN bot2-3 =15, KPN bot2-4 =16.

[0129] S164. Calculate the location of the critical node at the top dead center of the cylinder piston and obtain the critical node number. Select the node coordinate value Z as Z. max -L = 100, Y value equals 0 (i.e. Z = Z) max For nodes with L = 100 and Y = 0, the processing method is the same as in S162 and S163, yielding the critical node KP at the top dead center position of the piston in the first cylinder. ptop1-1 KP ptop1-2 KP ptop1-3 KP ptop1-4 The node coordinates are (30, 0, 100), (0, 30, 100), (-30, 0, 100), (0, -30, 100), and the node number is KPN. ptop1-1 =17, KPN ptop1-2 =18, KPN ptop1-3 =19, KPN ptop1-4 =20, the critical node KP at the top dead center of the piston in the second cylinder. ptop2-1 KP ptop2-2 KP ptop2-3 KP ptop2-4 The node coordinates are (100, 0, 100), (70, 30, 100), (40, 0, 100), (70, -30, 100), and the node numbers are KPN. ptop2-1 =21, KPN ptop2-2 =22, KPN ptop2-3 =23, KPN ptop2-4 =24.

[0130] S165. Calculate the location of the critical node at the bottom dead center of the cylinder piston and obtain the critical node number. Select the node coordinate value Z as Z. max -LS = 60, Y value equals 0 (i.e. Z = Z) max For nodes with -LS=60, Y=0, the processing method is the same as for S162 and S163, yielding the critical node KP at the bottom dead center position of the piston in the first cylinder. pbom1-1 KP pbom1-2 KP pbom1-3 KP pbom1-4 The node coordinates are (30, 0, 60), (0, 30, 60), (-30, 0, 60), (0, -30, 60), and the node number is KPN. pbom1-1 =25, KPN pbom1-2 =26, KPN pbom1-3 =27, KPN pbom1-4 =28, the critical node KP at the top dead center position of the piston in the second cylinder. pbom2-1 KP pbom2-2 KP pbom2-3 KP pbom2-4 The node coordinates are (100, 0, 60), (70, 30, 60), (40, 0, 60), (70, -30, 60), and the node numbers are KPN. pbom2-1 =29, KPN pbom2-2 =30, KPN pbom2-3 =31, KPN pbom2-4 =32.

[0131] S166. Calculate the location of the critical node at the highest cylinder temperature and obtain the critical node number. For the first cylinder, select the node with Z-coordinate value ZTM1 = 73.33 and Y-coordinate value equal to 0 (i.e., Z = ZTM1 = 73.33, Y = 0). The processing method is the same as in S162 and S163, thus obtaining the critical node KP at the highest temperature point of the first cylinder. tmp1-1 KP tmp1-2 KP tmp1-3 KP tmp1-4 The node coordinates are (30, 0, 73.33), (0, 30, 73.33), (-30, 0, 73.33), (0, -30, 73.33), and the node number is KPN. tmp1-1 =33, KPN tmp1-2 =34, KPN tmp1-3 =35, KPN tmp1-4=36; For the second cylinder, select the node with Z coordinate value ZTM1 = 86.67 and Y value equal to 0 (i.e., Z = ZTM1 = 86.67, Y = 0), and process it in the same way as S162 and S163 to obtain the key node KP at the highest temperature point of the second cylinder. tmp2-1 KP tmp2-2 KP tmp2-3 KP tmp2-4 The node coordinates are (100, 0, 86.67), (70, 30, 86.67), (40, 0, 86.67), (70, -30, 86.67), and the node number is KPN. tmp2-1 =37, KPN tmp2-2 =38, KPN tmp2-3 =39, KPN tmp2-4 =40.

[0132] S17. Extraction of curves from finite element calculation results of the axial temperature field of the engine cylinder:

[0133] For the first cylinder

[0134] S171. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. top1-1 =1 and KPN bot1-1 =9, select Create Curves, select Value for Y axis, select Distance for X axis, select New Plot for Place On, select OK to generate curves, change the header to LINER1, change the label item in Curve Label to 0deg, change the Y axis label to T (°C), change the X axis label to Distance (mm), use File-Export-Curves, select OK to save the curve data, and use CaptureFrame Area to save the curve plot.

[0135] S172. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. top1-2 =2 and KPN bot1-2=10, select Create Curves, select Value for Y axis, select Distance for X axis, select New Plot for Place On, select OK to generate curves, change the header to LINER1, change the label item in Curve Label to 90deg, change the Y axis label to T (°C), change the X axis label to Distance (mm), use File-Export-Curves, select OK to save the curve data, and use CaptureFrame Area to save the curve plot.

[0136] S173. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. top1-3 =3 and KPN bot1-3 =11, select Create Curves, select Value for Y axis, select Distance for X axis, select New Plot for Place On, select OK to generate curves, change the header to LINER1, change the label item in Curve Label to 180deg, change the Y axis label to T (°C), change the X axis label to Distance (mm), use File-Export-Curves, select OK to save the curve data, and use CaptureFrame Area to save the curve plot.

[0137] S174. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. top1-4 =4 and KPN bot1-4 =12, select Create Curves, select Value for Y axis, select Distance for X axis, select New Plot for Place On, select OK to generate curves, change the header to LINER1, change the label item in Curve Label to 270deg, change the Y axis label to T (°C), and change the X axis label to Distance (mm), use File-Export-Curves, select OK to save the curve data, and use CaptureFrame Area to save the curve plot.

[0138] For the second cylinder

[0139] S171. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. top2-1 =5 and KPN bot2-1 =13, select Create Curves, select Value for Y axis, select Distance for X axis, select New Plot for Place On, select OK to generate curves, change the header to LINER2, change the label item in Curve Label to 0deg, change the Y axis label to T (°C), change the X axis label to Distance (mm), use File-Export-Curves, select OK to save the curve data, and use CaptureFrame Area to save the curve plot.

[0140] S172. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. top2-2 =6 and KPN bot2-2 =14, select Create Curves, select Value for Y axis, select Distance for X axis, select New Plot for Place On, select OK to generate curves, change the header to LINER2, change the label item in Curve Label to 90deg, change the Y axis label to T (°C), change the X axis label to Distance (mm), use File-Export-Curves, select OK to save the curve data, and use CaptureFrame Area to save the curve plot.

[0141] S173. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. top3-3 =7 and KPN bot2-3=15, select Create Curves, select Value for Y axis, select Distance for X axis, select New Plot for Place On, select OK to generate curves, change the header to LINER2, change the label item in Curve Label to 180deg, change the Y axis label to T (°C), change the X axis label to Distance (mm), use File-Export-Curves, select OK to save the curve data, and use CaptureFrame Area to save the curve plot.

[0142] S174. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. top2-4 =8 and KPN bot2-4 =16, select Create Curves, select Value for Y axis, select Distance for X axis, select New Plot for Place On, select OK to generate curves, change the header to LINER2, change the label item in Curve Label to 270deg, change the Y axis label to T (°C), and change the X axis label to Distance (mm), use File-Export-Curves, select OK to save the curve data, and use CaptureFrame Area to save the curve plot.

[0143] S18. Extraction of curves from finite element calculations of the circumferential temperature field of the engine cylinder:

[0144] For the first cylinder, the coordinates of the center of the piston's top dead center are (XTO1=(XT2+XT1) / 2=0, 0, Z). max -L=100), the center coordinates of the cylinder piston bottom dead center position circle are (XBO1=(XB2+XB1) / 2=0,0,60), the center coordinates of the cylinder highest temperature point position circle are (XTO1=(XTMP2+XTMP1) / 2=0,0,ZTM1=73.33).

[0145] S181. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. ptop1-1 =17, KPN ptop1-2 =18, KPN ptop1-3 =19, KPN ptop1-4=20 and KPN ptop1-1 =17, retrieve all node numbers and coordinates along the path in the order they appear, as shown in Table 1. Figure 8 As shown, there are 12 nodes in total. The coordinates of the k-th node are KTP. k (XKTP k YKTP k Let k be an integer from 1 to 12, and k be the distance between the kth node and the center of the cylinder piston at top dead center. Then, the angle α between the straight line connecting the kth node and the center of the cylinder piston at top dead center and the positive x-axis is given. top-k For ① when XKTP k ≥0 and YKTP k When ≥0, α top-k =arctan(YKTP) k / XKTP k ) / π*180;②When XKTP k ≥0 and YKTP k When <0, α top-k =[arctan(YKTP)] k / XKTP k )+2π] / π*180;③When XKTP k <0 and YKTP k ≥0 or XKTP k <0 and YKTP k When <0, α top-k =[arctan(YKTP)] k / XKTP k ()+π] / π*180, see Table 2. Select "create curves", select "value" for Y axis, select "distance" for X axis, select "new plot" for "place on", select "ok" to generate curves, change "header" to "LINER1", change the "label" item in "curvelabel" to "TDC", change the "Y axis label" to "T" (°C), and change the "X axis label" to "deg" (°). In "data values", replace the X values ​​with the angle α between the node and the X axis according to the node order. top-k See Table 3. Use file-export-curves, select OK to save the curve data, and use capture frame area to save the curve graph.

[0146] Table 1. Numbers and coordinates of all nodes at the top dead center of the cylinder (arranged in path order).

[0147] Node number coordinate 17 30,0,100 50 25.98,15,100 51 15,25.98,100 18 0,30,100 52 -15,25.98,100 53 -25.98,15,100 19 -30,0,100 54 -25.98,-15,100 55 -15,-25.98,100 20 0,-30,100 56 15,-25.98,100 57 25.98,-15,100 17 30,0,100

[0148] Table 2. Angles between the line connecting the center of the cylinder top dead center and all nodes and the positive X-axis (arranged in path order).

[0149]

[0150]

[0151] Table 3: Replacement of X values ​​in data values ​​(arranged in path order)

[0152] Node number Before replacement After replacement 17 0 0 50 15.53 30 51 31.06 60 18 46.59 90 52 62.12 120 53 77.65 150 19 93.18 180 54 108.71 210 55 124.24 240 20 139.77 270 56 155.3 300 57 170.83 330 17 186.36 360

[0153] S182. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. pbom1-1 =25, KPN pbom1-2 =26, KPN pbom1-3 =27, KPN pbom1-4 =28 and KPN pbom1-1 =25, sorted by path, retrieve the node numbers and coordinates of all nodes on the path. There are 12 nodes in total, and the coordinates of the kth node are KBP. k (XKBP k YKBP k Let k be an integer from 1 to 12. Then, the angle α between the straight line connecting the kth node to the center of the cylinder piston at its bottom dead center and the positive x-axis is given by the line connecting the kth node to the center of the cylinder piston at its bottom dead center. bom-k For ① when XKBP k ≥0 and YKBP k When ≥0, α bom-k =arctan(YKBP) k / XKBP k ) / π*180;②When XKBP k ≥0 and YKBP k When <0, α bom-k =[arctan(YKBP)] k / XKBP k )+2π] / π*180;③When XKBP k <0 and YKBP k ≥0 or XKBP k <0 and YKBP k When <0, α bom-k =[arctan(YKBP)] k / XKBP k)+π] / π*180. Select "Create Curves", choose "Value" for the Y axis and "Distance" for the X axis. Select "New Plot" for "Place On", select "OK" to generate the curve. Change the header to "LINER1", change the label in the curve label to "BDC", change the Y axis label to "T (°)", and change the X axis label to "Deg (°)". In "Data Values", replace the X values ​​with the angle α between the node and the X axis, according to the node order. bom-k Use file-export-curves, select OK to save the curve data, and use capture frame area to save the curve graph;

[0154] S183. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. tmp1-1 =33, KPN tmp1-2 =34, KPN tmp1-3 =35, KPN tmp1-4 =36 and KPN tmp1-1 =33, sorted by path, retrieve the node numbers and coordinates of all nodes on the path. There are 12 nodes in total, and the coordinates of the k-th node are KMP. k (XKMP k YKMP k , 73.33), k is an integer from 1 to 12, then the angle α between the straight line connecting the kth node and the center of the circle where the highest temperature point of the cylinder is located and the positive x-axis is... max-k For ① when XKMP k ≥0 and YKMP k When ≥0, α max-k =arctan(YKMP) k / XKMP k ) / π*180;②When XKMP k ≥0 and YKMP k When <0, α max-k =[arctan(YKMP)] k / XKMP k )+2π] / π*180;③When XKMP k <0 and YKMP k ≥0 or XKMP k <0 and YKMP k When <0, α max-k =[arctan(YKMP)] k / XKMPk )+π] / π*180. Select "Create Curves", choose "Value" for the Y axis and "Distance" for the X axis. Select "New Plot" for "Place On", select "OK" to generate the curve. Change the header to "LINER1", change the label in the curve label to "MAX TEMP", change the Y axis label to "T (°)", and change the X axis label to "Deg (°)". In "Data Values", replace the X values ​​with the angle α between the node and the X axis, according to the node order. max-k Use file-export-curves, select OK to save the curve data, and use capture frame area to save the curve graph.

[0155] For the second cylinder, the coordinates of the center of the piston's top dead center are (XTO2=(XT4+XT3) / 2=70, 0, Z). max -L=100), the center coordinates of the cylinder piston bottom dead center position circle are (XBO2=(XB4+XB3) / 2=70,0,60), the center coordinates of the cylinder highest temperature point position circle are (XTO2=(XTMP4+XTMP3) / 2=70,0,ZTM2=86.67).

[0156] S181. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter K KPN sequentially. ptop2-1 =21, KPN ptop2-2 =22, KPN ptop2-3 =23, KPN ptop2-4 =24 and KPN ptop2-1 =21, retrieve all node numbers and coordinates along the path in the order they appear, as shown in Table 4. Figure 8 As shown, there are 12 nodes in total. The coordinates of the k-th node are KTP. k (XKTP k YKTP k Let k be an integer from 1 to 12, and k be the distance between the kth node and the center of the cylinder piston at top dead center. Then, the angle α between the straight line connecting the kth node and the center of the cylinder piston at top dead center and the positive x-axis is given. top-k For ① when XKTP k -70≥0 and YKTP k When ≥0, α top-k =arctan[YKTP k / (XKTP k -70)] / π*180;②When XKTP k-70≥0 and YKTP k When <0,

[0157] α top-k ={arctan[YKTP] k / (XKTP k -70)]+2π} / π*180; ③When XKTP k -70 < 0 and YKTP k ≥0 or XKTP k -70 < 0 and YKTP k When <0,

[0158] α top-k ={arctan[YKTP] k / (XKTP k -70)]+π} / π*180, see Table 5. Select "create curves", select "value" for Yaxis, select "distance" for Xaxis, select "new plot" for "place on", select "ok" to generate curves, change "header" to "LINER2", change the "label" item in "curve label" to "TDC", change the "Y axis label" to "T(°)", and change the "X axis label" to "deg(°)". In "data values", replace the X values ​​with the angle α between the node and the X-axis according to the node order. top-k See Table 6. Use file-export-curves, select OK to save the curve data, and use capture frame area to save the curve graph.

[0159] Table 4. Numbers and coordinates of all nodes at the top dead center of the cylinder (arranged in path order).

[0160] Node number coordinate 21 100,0,100 58 95.98,15,100 59 85,25.98,100 22 70,30,100 52 55,25.98,100 53 44.02,15,100 19 40,0,100 54 44.02,-15,100 55 55,-25.98,100 20 70,-30,100 56 85,-25.98,100 57 95.98,-15,100 17 100,0,100

[0161] Table 5. Angles between the line connecting the center of the cylinder top dead center and all nodes and the positive X-axis (arranged in path order).

[0162] Node number coordinate <![CDATA[included angle α top > 17 100,0,100 0 50 95.98,15,100 30 51 85,25.98,100 60 18 70,30,100 90 52 55,25.98,100 120 53 44.02,15,100 150 19 40,0,100 180 54 44.02,-15,100 210 55 55,-25.98,100 240 20 70,-30,100 270 56 85,-25.98,100 300 57 95.98,-15,100 330 17 100,0,100 360

[0163] Table 6: X value replacements in data values ​​(arranged in path order)

[0164] Node number Before replacement After replacement 17 0 0 50 15.53 30 51 31.06 60 18 46.59 90 52 62.12 120 53 77.65 150 19 93.18 180 54 108.71 210 55 124.24 240 20 139.77 270 56 155.3 300 57 170.83 330 17 186.36 360

[0165] S182. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. pbom2-1 =29, KPNpbom2-2 =30, KPN pbom2-3 =31, KPN pbom2-4 =32 and KPN pbom2-1 =29, sorted by path, retrieve the node numbers and coordinates of all nodes on the path. There are 12 nodes in total, and the coordinates of the kth node are KBP. k (XKBP k YKBP k Let k be an integer from 1 to 12. Then, the angle α between the straight line connecting the kth node to the center of the cylinder piston at its bottom dead center and the positive x-axis is given by the line connecting the kth node to the center of the cylinder piston at its bottom dead center. bom-k For ① when XKBP k -70≥0 and YKBP k When ≥0,

[0166] α bom-k =arctan[YKBP k / (XKBP k -70)] / π*180;②When XKBP k -70≥0 and YKBP k When <0,

[0167] α bom-k ={arctan[YKBP k / (XKBP k -70)]+2π} / π*180; ③When XKBP k -70 < 0 and YKBP k ≥0 or XKBP k -70 < 0 and YKBP k When <0,

[0168] α bom-k ={arctan[YKBP k / (XKBP k -70)]+π} / π*180. Select "Create Curves," choose "Value" for the Y axis and "Distance" for the X axis. Choose "New Plot" for "Place On," then select "OK" to generate the curve. Change the header to "LINER2," change the label in the curve label to "BDC," change the Y axis label to "T (°)," and the X axis label to "Deg (°)." In the data values, replace the X values ​​with the angle α between the node and the X-axis, according to the node order. bom-k Use file-export-curves, select OK to save the curve data, and use capture frame area to save the curve graph;

[0169] S183. In Hyperview software, access the annotations-measures module, add measuregroups, select node path, and enter KPN sequentially. tmp2-1 =37, KPN tmp2-2 =38, KPN tmp2-3 =39, KPN tmp2-4 =40 and KPN tmp2-1 =37, arranged in path order, retrieve the node numbers and coordinates of all nodes on the path. There are 12 nodes in total, and the coordinates of the k-th node are KMP. k (XKMP k YKMP k , 86.67), k is an integer from 1 to 12, then the angle α between the straight line connecting the kth node and the center of the circle where the highest temperature point of the cylinder is located and the positive x-axis is... max-k For ① when XKMP k -70≥0 and YKMP k When ≥0,

[0170] α max-k =arctan[YKMP k / (XKMP k -70)] / π*180;②When XKMP k -70≥0 and YKMP k When <0, α max-k ={arctan[YKMP k / (XKMP k -70)]+2π} / π*180; ③When XKMP k -70 < 0 and YKMP k ≥0 or XKMP k -70 < 0 and YKMP k When <0, α max-k ={arctan[YKMP k / (XKMP k -70)]+π} / π*180. Select "Create Curves," choose "Value" for the Y axis and "Distance" for the X axis. Choose "New Plot" for "Place On," then select "OK" to generate the curve. Change the header to "LINER2," change the label in the curve label to "MAX TEMP," change the Y axis label to "T (°)," and the X axis label to "deg (°)." In the data values, replace the X values ​​with the angle α between the node and the X-axis, according to the node order. max-kUse file-export-curves, select OK to save the curve data, and use capture frame area to save the curve graph.

[0171] S2. Programming and programming of processes S16-S18 based on TCL scripting language.

[0172] S3. Adjust the geometric coordinate system to satisfy the coordinate system defined in S11: The point where the horizontal plane containing the center line of cylinder bore 1 and the center line of the main journal intersects is defined as the origin of the rectangular coordinate system, such as... Figure 2a , Figure 2b As shown, the X-axis of the rectangular coordinate system is located along the crankshaft from cylinder 1 to cylinder 2, and the Z-axis is located along the main bearing cap to the cylinder block. The Y-axis of the rectangular coordinate system is determined according to the right-hand rule.

[0173] S4. Preprocess the engine cylinder geometry:

[0174] S41. The distance from the top dead center of the engine piston to the top of the cylinder is L = 10mm, and the piston stroke is S = 40mm. Therefore, the distance from the bottom dead center of the engine piston to the top of the cylinder is L + S = 50mm. Using the XY plane as the cutting plane, cut the engine cylinder geometry at the top dead center and bottom dead center positions of the engine piston. Figure 3 As shown;

[0175] S42. Using the XZ and YZ planes as cutting planes respectively, cut the engine cylinder geometry through the engine cylinder geometry center point, as follows: Figure 4 As shown.

[0176] S5. Perform finite element modeling on the engine cylinder bore, placing the finite element meshes on the cylinder bore surface into the same group, named F-LINER: the node distribution of the finite element meshes on the engine cylinder bore surface must satisfy the condition that the node distribution in each XY plane (i.e., the node coordinates Z values ​​are the same) is completely consistent (excluding intermediate nodes), such as... Figure 5a , Figure 5b As shown.

[0177] S6. Perform finite element simulation of engine cylinder block temperature field: This includes defining the material properties of the finite element model, applying loads and boundary conditions to the finite element model, and performing finite element analysis.

[0178] S7. Taking HyperView software as an example, read the finite element simulation calculation result file.

[0179] S8. Call the automatic extraction program for the finite element calculation results of the engine cylinder axial and circumferential temperature fields, which was compiled in S2.

[0180] S81. Use the program to automatically calculate the location of key nodes and obtain the key node number.

[0181] S82. The finite element calculation results curve of the axial temperature field of the engine cylinder is automatically extracted using the program.

[0182] S83. The program automatically extracts the finite element calculation result curve of the circumferential temperature field of the engine cylinder.

[0183] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any modifications, equivalent substitutions, and improvements made by those skilled in the art within the scope of the technology disclosed in the present invention, and within the spirit and principles of the present invention, should be included within the scope of protection of the present invention. Furthermore, all content not described in detail in this specification is prior art known to those skilled in the art.

Claims

1. A method for automatically extracting finite element result curves of engine cylinder temperature field, characterized in that, Includes the following steps: Step 1: Extract standardized rules and define a coordinate system based on the finite element calculation results of the axial and circumferential temperature fields of the engine cylinder. Step 2: Preprocess the engine cylinder geometry; Step 3: Perform finite element modeling on the engine cylinder; Step 4: Perform finite element simulation of the engine block temperature field; Step 5: Use HyperView software to read the finite element simulation calculation result file; Step 6: Calculate the location of key nodes and obtain their numbers; Step six, which involves calculating the locations of key nodes and obtaining their numbers, specifically includes operating on the finite element mesh group F-LINER on the cylinder surface: S1. Obtain the coordinates of the node with the maximum temperature of each cylinder barrel (XTM) j YTM j ZTM j ), j is an integer of 1, 2, 3...n, where n is the number of cylinders; S2. Calculate the location of key nodes at the top of the cylinder and obtain their numbers; select the nodes with the largest Z-value and a Y-value equal to 0, resulting in 2n nodes, where n is the number of cylinders; arrange these 2n nodes in ascending order of their X-values, with the coordinate value of the i-th node being (Xi + Xi). i ,0,Z max (i) is an integer from 1 to 2n; then the top radius of the j-th cylinder is RT. j =(XT) 2j - XT 2j-1 Let KP be the coordinates of the critical node at the top of the j-th cylinder, where KP is the number of cylinders in the range 1, 2, 3, ..., n. topj-1 (XT) 2j ,0,Z max ), KP topj-2 (XT) 2j -RT j RT j Z max ), KP topj-3 (XT) 2j-1 ,0,Z max ), KP topj-4 (XT) 2j -RT j -RT j Z max Based on the coordinates of the critical node at the top of the j-th cylinder, find and obtain the critical node KP at the top of the j-th cylinder. topj-1 KP topj-2 KP topj-3 KP topj-4 The node numbers are KPN in sequence. topj-1 KPN topj-2 KPN topj-3 KPN topj-4 ; S3. Calculate the positions of key nodes at the bottom of the cylinder and obtain their numbers; select the nodes with the smallest Z-value and a Y-value equal to 0, resulting in 2n nodes, where n is the number of cylinders; arrange these 2n nodes in ascending order of their X-values, with the i-th node having a coordinate value of (XB...). i ,0,Z min (i) is an integer from 1 to 2n; then the bottom radius of the j-th cylinder is RB. j =(XB 2j - XB 2j-1 Let KP be the coordinates of the critical node at the bottom of the j-th cylinder, where j is an integer from 1 to n. botj-1 (XB) 2j ,0,Z min ), KP botj-2 (XB) 2j -RB j RB j Z min ), KP botj-3 (XB) 2j-1 ,0,Z min ), KP botj-4 (XB) 2j -RB j -RB j Z min Based on the coordinates of the critical node at the bottom of the j-th cylinder, find and obtain the critical node KP at the bottom of the j-th cylinder. botj-1 KP botj-2 KP botj-3 KP botj-4 The node numbers are KPN in sequence. botj-1 KPN botj-2 KPN botj-3 KPN botj-4 ; S4. Calculate the location of the critical node at the top dead center of the cylinder piston and obtain the critical node number; select the node coordinate value Z as Z. max -L, Y value equals 0, that is, Z=Z max There are 2n nodes with -L,Y=0, where n is the number of cylinders. These 2n nodes are arranged in ascending order of their X-coordinates, and the coordinate of the i-th node is (XTP). i ,0,Z max -L), where i is an integer from 1 to 2, 3 to 2n; then the radius of the piston top dead center position of the j-th cylinder is RTP. j =(XTP) 2j - XTP 2j-1 ) / 2, where j is an integer from 1 to n; the key node coordinates of the top dead center position of the piston in the j-th cylinder are KP. ptopj-1 (XTP) 2j ,0,Z max -L), KP ptopj-2 (XTP) 2j -RTP j RTP j Z max -L), KP ptopj-3 (XTP) 2j-1 ,0,Z max -L), KP ptopj-4 (XTP) 2j -RTP j -RTP j Z max -L); Based on the coordinates of the key node at the top dead center position of the piston in the j-th cylinder, find and obtain the key node KP at the top dead center position of the piston in the j-th cylinder. ptopj-1 KP ptopj-2 KP ptopj-3 KP ptopj-4 The node numbers are KPN in sequence. ptopj-1 KPN ptopj-2 KPN ptopj-3 KPN ptopj-4 ; S5. Calculate the location of the critical node at the bottom dead center of the cylinder piston and obtain the critical node number; select the node coordinate value Z as Z. max -LS, Y value equals 0, i.e. Z=Z max There are 2n nodes with -LS,Y=0, where n is the number of cylinders. These 2n nodes are arranged in ascending order of their X-coordinates, and the coordinate of the i-th node is (XBP). i ,0,Z max -LS), where i is an integer from 1 to 2, 3 to 2n; then the radius of the bottom dead center position of the j-th cylinder piston is RBP. j = (XBP) 2j - XBP 2j-1 Let ) / 2, where j is an integer from 1 to n; then the coordinates of the critical node at the bottom dead center of the piston in the j-th cylinder are KP. pbomj-1 (XBP) 2j ,0,Z max -LS), KP pbomj-2 (XBP) 2j -RBP j RBP j Z max -LS), KP pbomj-3 (XBP) 2j-1 ,0,Z max -LS), KP pbomj-4 (XBP) 2j -RBP j -RBP j Z max -LS); Based on the coordinates of the key node at the bottom dead center of the piston in the j-th cylinder, find and obtain the key node KP at the bottom dead center of the piston in the j-th cylinder. pbomj-1 KP pbomj-2 KP pbomj-3 KP pbomj-4 The node numbers are KPN in sequence. pbomj-1 KPN pbomj-2 KPN pbomj-3 KPN pbomj-4 ; S6. Calculate the location of the critical node at the highest cylinder temperature and obtain the critical node number; for the j-th cylinder, select the node coordinate Z value as ZTM. j The value of Y is equal to 0, that is, Z = ZTM. j There are 2n nodes with Y=0, where n is the number of cylinders and j is an integer from 1 to n. These 2n nodes are arranged in ascending order of their X-coordinates, and the coordinate of the i-th node is (Xi + Xj + Xj). i ,0,ZTM j (i) is an integer from 1 to 2n; then the radius of the highest temperature point in the j-th cylinder is RTMP. j =(XTMP) 2j - XTMP 2j-1 Let ) / 2, where j is an integer from 1 to n; then the coordinates of the key node at the highest temperature point of the j-th cylinder are KP. tmpj-1 (XTMP) 2j ,0,ZTM j ), KP tmpj-2 (XTMP) 2j -RTMP j RTMP j ZTM j ), KP tmpj-3 (XTMP) 2j-1 ,0,ZTM j ), KP tmpj-4 (XTMP) 2j -RTMP j -RTMP j ZTM j Based on the coordinates of the key node at the highest temperature point of the j-th cylinder, find and obtain the key node KP at the highest temperature point of the j-th cylinder. tmpj-1 KP tmpj-2 KP tmpj-3 KP tmpj-4 The node numbers are KPN in sequence. tmpj-1 KPN tmpj-2 KPN tmpj-3 KPN tmpj-4 ; Step 7: Extract the curves from the finite element calculation of the axial temperature field of the engine cylinder. Step 8: Extract the curve of the finite element calculation result of the circumferential temperature field of the engine cylinder.

2. The method for automatically extracting the finite element result curve of the engine cylinder temperature field according to claim 1, characterized in that: Defining the coordinate system as described in step one specifically includes: The point where the center line of cylinder bore 1 intersects the horizontal plane containing the center line of main bearing journal is defined as the origin of the rectangular coordinate system. The direction from cylinder 1 to cylinder 2 along the crankshaft is defined as the X-axis of the rectangular coordinate system, and the direction from the main bearing cap to the cylinder block is defined as the Z-axis of the rectangular coordinate system. The Y-axis of the rectangular coordinate system is determined according to the right-hand rule.

3. The method for automatically extracting the finite element result curve of the engine cylinder temperature field according to claim 1, characterized in that: Step two, which involves preprocessing the engine cylinder geometry, specifically includes the following steps: S1. The distance from the top dead center of the engine piston to the top of the cylinder is L, and the piston stroke is S. Then the distance from the bottom dead center of the engine piston to the top of the cylinder is L+S. Using the XY plane as the cutting plane, the engine cylinder geometry is cut at the top dead center and bottom dead center positions of the engine piston. S2. Using the XZ plane and YZ plane as cutting planes respectively, cut the engine cylinder geometry through the engine cylinder geometry center point.

4. The method for automatically extracting the finite element result curve of the engine cylinder temperature field according to claim 1, characterized in that: Step 3, which describes finite element modeling of the engine cylinder, specifically includes: placing the finite element meshes on the cylinder surface into the same group, named F-LINER: the distribution of the finite element mesh nodes on the engine cylinder surface satisfies that the distribution of nodes in each XY plane is completely consistent, excluding intermediate nodes.

5. The method for automatically extracting the finite element result curve of the engine cylinder temperature field according to claim 1, characterized in that: Step four involves performing finite element simulation of the engine block temperature field, specifically including: defining the material properties of the finite element model, applying loads and boundary conditions to the finite element model, and performing finite element analysis.

6. The method for automatically extracting the finite element result curve of the engine cylinder temperature field according to claim 1, characterized in that: Step seven involves extracting the finite element calculation result curve of the engine cylinder axial temperature field for the j-th cylinder, where j is an integer from 1 to n; this includes the following steps: S1. In Hyperview software, access the annotations-measures module, add measure groups, select node path, and enter KPN sequentially. topj-1 and KPN botj-1 Select "create curves", select "value" for Y axis and "distance" for X axis. Select "new plot" for "place on", select "ok" to generate curves, change "header" to "LINERj", change "label" in "curve label" to "0deg", change "Y axis label" to "T℃", and change "X axis label" to "distance" in mm. Use "file-export-curves", select "ok" to save the curve data, and use "capture frame area" to save the curve plot. S2. In Hyperview software, access the annotations-measures module, add measure groups, select node path, and enter KPN sequentially. topj-2 and KPN botj-2 Select "Create Curves," choose "Value" for the Y-axis and "Distance" for the X-axis. Select "New Plot" for "Place On," then click "OK" to generate the curve. Change the header to "LINERj," change the label in the curve label to "90deg," change the Y-axis label to "T℃," and the X-axis label to "Distance" (unit: mm). Use "File-Export-Curves," click "OK" to save the curve data, and use "Capture Frame Area" to save the graph. S3. In Hyperview software, access the annotations-measures module, add measure groups, select node path, and enter KPN sequentially. topj-3 and KPN botj-3 Select "create curves", select "value" for Y axis and "distance" for X axis. Select "new plot" for "place on", select "ok" to generate curves, change "header" to "LINERj", change "label" in "curve label" to "180deg", change "Y axis label" to "T℃", and change "X axis label" to "distance" in mm. Use "file-export-curves", select "ok" to save the curve data, and use "capture frame area" to save the curve plot. S4. In Hyperview software, access the annotations-measures module, add measure groups, select node path, and enter KPN sequentially. topj-4 and KPN botj-4 Select "Create Curves", choose "Value" for the Y axis and "Distance" for the X axis. Select "New Plot" for "Place On", and click "OK" to generate the curve. Change the header to "LINERj", change the label in the curve label to "270deg", change the Y axis label to "T℃", and change the X axis label to "Distance" in mm. Use "File-Export-Curves", click "OK" to save the curve data, and use "Capture Frame Area" to save the curve plot.

7. The method for automatically extracting the finite element result curve of the engine cylinder temperature field according to claim 1, characterized in that: Step eight describes the extraction of the finite element calculation results curve of the circumferential temperature field of the engine cylinder; for the j-th cylinder, the coordinates of the center of the piston top dead center position are (XTO). j =(XT) 2j + XT 2j-1 ) / 2, 0, Z max -L), the coordinates of the center of the cylinder piston bottom dead center are (XBO) j =(XB 2j + XB 2j-1 ) / 2, 0, Z max -LS), the coordinates of the center of the cylinder's highest temperature point are (XTO). j =(XTMP) 2j + XTMP 2j-1 ) / 2, 0, ZTM j (j is an integer of 1, 2, 3...n), specifically including the following steps: S1. In Hyperview software, access the annotations-measures module, add measure groups, select node path, and enter KPN sequentially. ptopj-1 KPN ptopj-2 KPN ptopj-3 KPN ptopj-4 and KPN ptopj-1 Arrange the nodes in the order of the path, and retrieve the node numbers and coordinates of all nodes on the path. There are m nodes in total, and the coordinates of the k-th node are KTP. k (XKTP) k YKTP k Z max -L), where k is an integer from 1 to m, then the angle between the straight line connecting the kth node and the center of the circle at the top dead center of the cylinder piston and the positive x-axis is α. top-k ; Select "Create Curves," choose "Value" for the Y-axis and "Distance" for the X-axis. Choose "New Plot" for "Place On," and select "OK" to generate the curve. Change the header to "LINERj," and in the curve label section, change the label to "TDC," the Y-axis label to "T℃," and the X-axis label to "Deg." In the data values ​​section, replace the X-values ​​with the angle α between the node and the X-axis, according to the node order. top-k Use file-export-curves, select OK to save the curve data, and use captureframe area to save the curve graph; S2. In Hyperview software, access the annotations-measures module, add measure groups, select node path, and enter KPN sequentially. pbomj-1 KPN pbomj-2 KPN pbomj-3 KPN pbomj-4 and KPN pbomj-1 Arrange the nodes in the order of the path, retrieve the node numbers and coordinates of all nodes on the path. There are m nodes in total, and the coordinates of the k-th node are KBP. k (XKBP) k YKBP k Z max -LS), where k is an integer from 1 to m, then the angle between the straight line connecting the kth node and the center of the circle at the bottom dead center of the cylinder piston and the positive x-axis is α. bom-k ; Select "Create Curves," choose "Value" for the Y-axis and "Distance" for the X-axis. Choose "New Plot" for "Place On," and select "OK" to generate the curve. Change the header to "LINERj," and in the curve label section, change the label to "BDC," the Y-axis label to "T℃," and the X-axis label to "Deg." In the data values ​​section, replace the X-values ​​with the angle α between the node and the X-axis, according to the node order. bom-k Use file-export-curves, select OK to save the curve data, and use captureframe area to save the curve graph; S3. In Hyperview software, access the annotations-measures module, add measure groups, select node path, and enter KPN sequentially. tmpj-1 KPN tmpj-2 KPN tmpj-3 KPN tmpj-4 and KPN tmpj-1 Arrange the nodes in the order of the path, retrieve the node numbers and coordinates of all nodes on the path. There are m nodes in total, and the coordinates of the k-th node are KMP. k (XKMP) k YKMP k ZTM j If k is an integer from 1 to m, then the angle between the straight line connecting the k-th node and the center of the circle at the highest temperature point of the cylinder and the positive x-axis is α. max-k ; Select "Create Curves," choose "Value" for the Y-axis and "Distance" for the X-axis. Choose "New Plot" for "Place On," and select "OK" to generate the curve. Change the header to "LINERj," and in the curve label, change the label to "MAX TEMP." Change the Y-axis label to "T℃," and the X-axis label to "Deg." In the data values, replace the X-values ​​with the angle α between the node and the X-axis, according to the node order. max-k Use file-export-curves, select OK to save the curve data, and use captureframe area to save the curve graph.

8. The method for automatically extracting the finite element result curve of the engine cylinder temperature field according to claim 7, characterized in that: The angle between the x-axis and the positive x-axis is α. top-k ; ①When XKTP k -XTO j ≥0 and YKTP k When ≥0, α top-k =arctan[YKTP k / (XKTP) k -XTO j )] / π*180; ②When XKTP k -XTO j ≥0 and YKTP k When <0, α top-k ={arctan[YKTP k / (XKTP) k -XTO j )]+2π} / π*180; ③ When XKTP k -XTO j <0 and YKTP k ≥0 or XKTP k -XTO j <0 and YKTP k When <0, α top-k ={arctan[YKTP k / (XKTP) k -XTO j )]+π} / π*180; The angle between the x-axis and the positive x-axis is α. bom-k ; ① When XKBP k XBO j ≥0 and YKBP k When ≥0, α bom-k =arctan[YKBP k / (XKBP k XBO j )] / π*180; ②When XKBP k XBO j ≥0 and YKBP k When <0, α bom-k ={arctan[YKBP k / (XKBP k XBO j )]+2π} / π*180; ③ When XKBP k XBO j <0 and YKBP k ≥0 or XKBP k XBO j <0 and YKBP k When <0, α bom-k ={arctan[YKBP k / (XKBP k XBO j )]+π} / π*180; The angle between the x-axis and the positive x-axis is α. max-k : ① When XKMP k -XTO j ≥0 and YKMP k When ≥0, α max-k =arctan[YKMP k / (XKMP k -XTO j )] / π*180; ② When XKMP k -XTO j ≥0 and YKMP k When <0, α max-k ={arctan[YKMP k / (XKMP k -XTO j )]+2π} / π*180; ③ When XKMP k -XTO j <0 and YKMP k ≥0 or XKMP k -XTO j <0 and YKMP k When <0, α max-k ={arctan[YKMP k / (XKMP k -XTO j )]+π} / π*180。 9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program that implements the method as described in any one of claims 1 to 8.

Images