A method and system for constructing three-dimensional temperature sensing lines based on finite element analysis

Abstract

This invention provides a method and system for constructing a three-dimensional temperature sensing wire based on finite element analysis, belonging to the field of data processing technology. The method includes: Step 1, collecting physical parameter data of the temperature sensing wire, including core properties, heating wire temperature coefficient, and preset gap size; Step 2, constructing a virtual three-dimensional geometric model of the temperature sensing wire based on the physical parameter data, including a primary winding body and an intermediate layer structure. This invention accurately constructs and dynamically optimizes the three-dimensional model of the temperature sensing wire through finite element analysis, improving the model's simulation accuracy of actual heat conduction characteristics. It provides a reliable virtual analysis basis for structural design optimization and performance verification of the temperature sensing wire, helping to improve the design efficiency and quality of temperature sensing wire products.

Description

Technical Field

[0001] This invention relates to the field of data processing technology, and in particular to a method and system for constructing three-dimensional temperature sensing lines based on finite element analysis. Background Technology

[0002] As a key component for temperature monitoring in heating blankets, the performance of the temperature sensing wire directly affects the temperature control accuracy and safety of the blanket. Currently, the 3D modeling and thermal analysis of such temperature sensing wires mostly employ computer-aided design (CAD) software and finite element analysis (FEA) technology. The typical process is as follows: a geometric model is constructed based on the macroscopic dimensional parameters and material properties of the temperature sensing wire; then, the model is imported into finite element analysis software for mesh generation and heat conduction simulation calculations.

[0003] However, the layers of the temperature-sensing wire used in heating blankets typically have pre-set non-uniform gaps. The geometry of these gaps and the thermal conductivity of the internal medium can affect the overall heat flow distribution and thermal coupling effect of the temperature-sensing wire. For example, when the temperature-sensing wire is in a folded area of ​​the heating blanket, the gap between its layers may locally shrink, while in a flat area, the gap roughly maintains its initial size. This dynamic change may lead to changes in local thermal resistance. Existing modeling methods still have room for improvement in accurately representing the dynamic changes of such micro-gap and their impact on the heat conduction path. Summary of the Invention

[0004] The technical problem to be solved by the present invention is to provide a three-dimensional construction method and system for temperature sensing wires based on finite element analysis, which improves the mechanical strength and temperature resistance of temperature sensing wires.

[0005] To solve the above-mentioned technical problems, the technical solution of the present invention is as follows:

[0006] Firstly, a method for constructing three-dimensional temperature sensing lines based on finite element analysis, the method comprising:

[0007] Step 1: Collect physical parameter data of the temperature sensing wire, including core properties, heating wire temperature coefficient, and preset gap size;

[0008] Step 2: Based on the physical parameter data, construct a virtual three-dimensional geometric model of the temperature sensing wire, including the primary winding body and the intermediate layer structure;

[0009] Step 3: Set a spatial positioning mark array on the surface of the virtual three-dimensional geometric model. The mark array is evenly distributed around the central axis of symmetry of the blanket body, and the longitudinal spacing between adjacent marks is no more than 0.5 meters. Generate a three-dimensional laying path model based on the mark array, and perform finite element heat conduction mesh division on the three-dimensional laying path model to extract the thermal coupling parameters between nodes.

[0010] Step 4: Calculate the local heat flux density based on the thermal coupling parameters, and dynamically adjust the gap region parameters in the virtual three-dimensional geometric model;

[0011] Step 5: Based on the adjusted gap region parameters, update the virtual three-dimensional geometric model to form the final three-dimensional virtual construction result of the temperature sensing line.

[0012] Secondly, a three-dimensional construction system for temperature sensing lines based on finite element analysis includes:

[0013] The data acquisition module is used to collect physical parameter data such as the properties of the heating wire core, the temperature coefficient of the heating wire, and the preset gap size.

[0014] The modeling module is used to construct a virtual three-dimensional geometric model of the temperature-sensing wire, which includes a primary winding body and an intermediate layer structure, based on physical parameter data.

[0015] The extraction module is used to set a spatial positioning mark array on the surface of a virtual three-dimensional geometric model. The mark array is uniformly distributed around the central axis of symmetry of the blanket, and the longitudinal spacing between adjacent marks is no more than 0.5 meters. Based on the mark array, a three-dimensional laying path model is generated, and the three-dimensional laying path model is divided into finite element heat conduction meshes to extract the thermal coupling parameters between nodes.

[0016] The control module is used to calculate the local heat flux density based on the thermal coupling parameters and dynamically control the parameters of the gap region in the virtual three-dimensional geometric model.

[0017] The update module is used to update the virtual three-dimensional geometric model based on the adjusted gap region parameters, forming the final three-dimensional virtual construction result of the temperature sensing line.

[0018] Thirdly, a computing device includes:

[0019] One or more processors;

[0020] A storage device for storing one or more programs that, when executed by one or more processors, cause the one or more processors to implement the method.

[0021] Fourthly, a computer-readable storage medium storing a program that, when executed by a processor, implements the method.

[0022] The above-described solution of the present invention has at least the following beneficial effects:

[0023] By collecting physical parameters such as preset gap dimensions, the intermediate layer structure is explicitly included when constructing the virtual 3D geometric model, allowing for accurate representation of the gaps between the layers of the sensing wire. This avoids model distortion caused by oversimplification of the gaps and improves the consistency between the 3D model and the actual structure. Local heat flux density is calculated using thermal coupling parameters, and the gap region parameters are dynamically adjusted accordingly. This allows for real-time response to the impact of dynamic changes in the gaps of the sensing wire at different locations on the blanket (such as folds versus flat areas) on heat conduction, ensuring the model accurately reflects the thermal behavior in actual use. A 3D laying path model is generated based on a spatial positioning marker array distributed along the central axis of symmetry of the blanket, and targeted finite element meshing is performed. This ensures the accuracy of heat conduction characteristic analysis while avoiding the drastic increase in computational load caused by indiscriminate fine meshing, improving the efficiency and reliability of thermal simulation results. This provides precise model support for the design optimization and performance evaluation of the sensing wire. Attached Figure Description

[0024] Figure 1 This is a flowchart illustrating a three-dimensional construction method for temperature sensing lines based on finite element analysis, provided by an embodiment of the present invention.

[0025] Figure 2 This is a schematic diagram of a three-dimensional construction system for temperature sensing lines based on finite element analysis, provided by an embodiment of the present invention. Detailed Implementation

[0026] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.

[0027] like Figure 1 As shown, an embodiment of the present invention proposes a three-dimensional construction method for temperature sensing lines based on finite element analysis, the method comprising the following steps:

[0028] Step 1: Collect physical parameter data of the temperature sensing wire, including core properties, heating wire temperature coefficient, and preset gap size;

[0029] Step 2: Based on the physical parameter data, construct a virtual three-dimensional geometric model of the temperature sensing wire, including the primary winding body and the intermediate layer structure;

[0030] Step 3: Set a spatial positioning mark array on the surface of the virtual three-dimensional geometric model. The mark array is evenly distributed around the central axis of symmetry of the blanket body, and the longitudinal spacing between adjacent marks is no more than 0.5 meters. Generate a three-dimensional laying path model based on the mark array, and perform finite element heat conduction mesh division on the three-dimensional laying path model to extract the thermal coupling parameters between nodes.

[0031] Step 4: Calculate the local heat flux density based on the thermal coupling parameters, and dynamically adjust the gap region parameters in the virtual three-dimensional geometric model;

[0032] Step 5: Based on the adjusted gap region parameters, update the virtual three-dimensional geometric model to form the final three-dimensional virtual construction result of the temperature sensing line.

[0033] In this embodiment of the invention, by collecting physical parameters such as preset gap dimensions, the intermediate layer structure is explicitly included when constructing the virtual three-dimensional geometric model. This allows for the accurate representation of the gaps between the layers of the sensing wire, avoiding model distortion caused by oversimplification of the gaps and improving the consistency between the three-dimensional model and the actual structure. By calculating local heat flux density through thermal coupling parameters and dynamically adjusting the gap region parameters accordingly, the model can respond in real time to the impact of dynamic changes in the gaps of the sensing wire at different locations on the blanket (such as folds and flat areas) on heat conduction, enabling the model to accurately reflect the thermal behavior in actual use. A three-dimensional laying path model is generated based on a spatial positioning marker array distributed along the central axis of symmetry of the blanket, and targeted finite element meshing is performed. This ensures the accuracy of heat conduction characteristic analysis while avoiding the problem of a dramatic increase in computational load caused by indiscriminate fine meshing, improving the efficiency and reliability of thermal simulation results. This provides accurate model support for the design optimization and performance evaluation of the sensing wire.

[0034] In a preferred embodiment of the present invention, step 1 involves collecting physical parameter data of the temperature sensing wire, including core properties, heating wire temperature coefficient, and preset gap size; step 2 includes:

[0035] Step 200: Analyze the core attribute parameters to generate the three-dimensional geometry and dimensions of the primary winding structure; analyze the preset gap structure parameters to generate the three-dimensional geometry, dimensions, and material property definitions of the intermediate layer structure.

[0036] Step 201: Spatially assemble the primary winding structure and the intermediate layer structure to form a virtual three-dimensional geometric model.

[0037] In this embodiment of the invention, the basic physical parameters required for constructing the three-dimensional model of the temperature sensing line are obtained, and the specific implementation process is as follows:

[0038] Core property parameter acquisition: Samples are cut from the core and analyzed using an infrared spectrometer. The material type of the core is determined based on the detected spectral characteristics. For example, when a specific nylon characteristic peak is detected, the core material can be identified as nylon, and its density value (e.g., 1.14 g / cm³) and heat resistance temperature range (e.g., -40℃ to 120℃) are recorded simultaneously. Five measurement points are evenly selected along the core's axis, with the first measurement point located 1 cm from one end of the core, and subsequent measurement points selected every 2 cm until 1 cm from the other end. The diameter of each point is measured using a vernier caliper with an accuracy of 0.01 mm. If the five measurements are 2.01 mm, 2.03 mm, 2.02 mm, 2.00 mm, and 2.02 mm respectively, then... Calculate the average value and use it as the reference diameter of the core. Simultaneously, measure the total length of the core using a ruler with a minimum scale of 0.1 mm. If the measured value is 100.2 mm, record this value. Observe the surface of the core using a magnifying glass. If the surface is smooth and without grooves, record it as "smooth surface, no grooves". Also, mark one end of the core 2 cm from the endpoint; this mark is the starting point of the winding. If the core has an irregular shape (such as a hexagonal cross-section), use a 3D scanner to scan the core to obtain cross-sectional contour data, measure the length of each side of the hexagon, and record these key contour dimensions.

[0039] Heating wire temperature coefficient acquisition: A 15cm long sample of heating wire from the same batch as the temperature sensing wire was cut, and its surface insulation layer was carefully removed with a blade. The sample was placed in a temperature control chamber, and the temperature was first set to -10℃. After the temperature stabilized for 30 minutes, its resistance value was measured using a precision resistance meter with an accuracy of 0.001Ω. The measured value was 5.210Ω. The resistance value was then measured every 5℃ increase, and resistance data such as 5.320Ω at 10℃ and 5.380Ω at 15℃ were obtained. Linear fitting was performed on the resistance data corresponding to each recorded temperature point. Specifically, the temperature value was used as the abscissa and the resistance value as the ordinate, and a linear relationship equation R=R0(1+αT) was established between the two, where R is the measured resistance value, R0 is the resistance value at the reference temperature, α is the temperature coefficient, and T is the temperature value. The measurement data were fitted using the least squares method: First, a reference temperature (e.g., 0℃) was selected. The temperature and resistance measurements at each temperature point were substituted into the equations to obtain a set of equations containing R0 and α. Solving the equations yielded the value of α, which is the fitted temperature coefficient. If the fitted temperature coefficient was 0.003Ω / ℃, the fitting error was calculated: For each measurement point, the deviation percentage was calculated as (actual resistance value at that temperature point - R) ÷ actual resistance value at that temperature point × 100%. The arithmetic mean of the deviation percentages for all points was then taken as the fitting error. When the fitting error was 1.5% (meeting the requirement of error rate ≤ 2%), the temperature coefficient was determined to be an effective parameter. The diameter of the heating wire was measured three times at different locations using a microscope. If the measured values ​​were 0.20mm, 0.21mm, and 0.20mm respectively, material analysis determined that the heating wire was a copper-nickel alloy, and the above auxiliary parameters were recorded.

[0040] Preset gap size parameter acquisition: Select an unused complete temperature sensing wire sample, and cut 10 segments along the axial direction, starting from one end, every 10cm, for a total of 10 segments, each 5cm long. Use a slicer to slice each sample, and obtain cross-sectional images of the gap between the winding layers using a microscope. Analyze the cross-sectional images of each sample, and use image measurement software to measure the maximum width, minimum width, and average width of the gap. For example, the maximum width of a certain sample is 0.3mm, the minimum width is 0.1mm, and the average width is 0.2mm. After measuring 10 samples, take the average value as the initial preset gap width. If the average value is 0.18mm, record this value. Observe the distribution pattern of the gap along the axial direction. If the first 2 segments and the last 2 segments are found to have larger gap widths, and the middle 6 segments are more uniform, then mark the first 2 segments as the 0-20cm region, corresponding to a gap width of 0.25mm; the last 2 segments as the 80-100cm region, corresponding to a gap width of 0.23mm; and the middle 6 segments as the 20-80cm region, corresponding to a gap width of 0.15mm. After determining that the material filling the gap is air, its thermal physical properties, such as thermal conductivity (0.026 W / (m·K)) and specific heat capacity (1005 J / (kg·K)), are obtained by querying the material database and recorded.

[0041] Step 200: Three-dimensional geometry generation of the winding structure: Using the core attribute parameters as input, and based on the core's reference diameter (2.016mm) and length (100.2mm), a cylindrical central axis solid is generated in the 3D modeling software. This solid has a diameter of 2.016mm and a length of 100.2mm. The diameter (0.203mm) and winding parameters of the heating wire are analyzed, and the winding method is determined to be a single helix. The winding start point marked on the core is taken as the starting point of the helix. The winding pitch of the temperature sensing wire sample is measured to be 3mm, and the number of winding turns is 30 turns. The guide of the helix is ​​calculated. The lead is 30×3=90mm; in the 3D modeling software, a 3D spiral path with a lead of 90mm is generated based on the core center axis; the diameter of the heating wire (0.203mm) is stretched along the generated spiral path to form a solid heating wire winding structure. This structure is combined with the core center axis to form the complete 3D geometry of the primary winding body; at the same time, the key dimensions of the primary winding body are marked, where the maximum outer diameter is the core diameter plus twice the heating wire diameter (2.016+2×0.203=2.422mm), the total length is 100.2mm, and the number of winding turns is 30 turns.

[0042] 3D geometry and material property generation of the intermediate layer structure: Analyzing the preset gap size parameters, the maximum radius of the single-wound body is 2.422 ÷ 2 = 1.211 mm, the average width of the preset gap is 0.18 mm, and the calculated inner wall radius of the intermediate layer is 1.211 + 0.18 = 1.391 mm; the total length of the temperature sensing wire is 100.2 mm. Based on the gap distribution pattern (wider gaps in the first 20 cm and last 20 cm regions, and uniform gaps in the middle 60 cm region), the inner wall radius of the intermediate layer is adjusted: that is, the inner wall radius of the first 20 cm region is 1.211 + 0.25 = 1.461 mm, the inner wall radius of the last 20 cm region is 1.211 + 0.23 = 1.441 mm, and the inner wall radius of the middle 60 cm region is... The inner wall radius of the region is 1.391 mm. Based on the calculated inner wall radius of each region, assuming the thickness of the outer protective structure is 0.3 mm, the outer wall radii of each region are 1.461 + 0.3 = 1.761 mm, 1.441 + 0.3 = 1.741 mm, and 1.391 + 0.3 = 1.691 mm, respectively, thus forming a ring-shaped three-dimensional solid. For non-uniform gap regions, the geometry is adjusted by local stretching to ensure consistency with the actual gap distribution. The thermal conductivity (0.026 W / (m·K), specific heat capacity (1005 J / (kg·K), and density (1.29 kg / m³) of air are associated with the intermediate layer geometry to establish a mapping relationship between material properties and geometric regions.

[0043] Step 201: Through precise spatial positioning, the primary winding body and the intermediate layer structure are combined into a complete model. The specific implementation process is as follows:

[0044] A three-dimensional spatial coordinate system is established with the central axis of the core as the Z-axis of the global coordinate system and one end of the core as the origin (0, 0, 0). The generated primary winding structure is imported into the above coordinate system. Through movement and rotation operations, its central axis is made to completely coincide with the Z-axis, and its starting end is aligned with the origin (0, 0, 0), ensuring its accurate position in the axial (Z-axis) and radial (X, Y-axis) directions. The intermediate layer structure is imported into the coordinate system, and its position is adjusted with the outer wall of the primary winding as a reference: that is, for the middle 60cm area with uniform gaps, the inner wall of the intermediate layer is kept 0.15mm away from the outer wall of the primary winding; for the non-uniform gap areas of the first 20cm and the last 20cm, the distances are kept at 0.25mm and 0.23mm respectively through local coordinate fine-tuning. The interference check function of the 3D modeling software is used for inspection. If an overlap is found in a certain area, the process is backtracked to step 200 to readjust the geometric dimensions of that area. The total outer diameter of the assembled model is measured. If it is 1.691×2=3.382mm and the total length is 100.2mm, the error is 0.5% compared with the measured value of the actual temperature sensing wire sample (meeting the requirement of error ≤1%). Finally, a complete virtual 3D geometric model containing the primary winding body and the intermediate layer structure is formed. The topological structure of the model (such as the positional relationship between the primary winding body and the intermediate layer) and the dimensional parameters of each part are saved.

[0045] By acquiring parameters in multiple dimensions and with high precision, the comprehensiveness and accuracy of the input data are ensured, avoiding deviations in the structure and performance of the model from the actual temperature sensing line due to missing or incorrect parameters. By generating a three-dimensional structure through structured analytical parameters, the physical parameters are accurately converted into a geometric model, ensuring that the winding shape of the primary winding body and the gap distribution of the intermediate layer are consistent with reality. Spatial assembly based on a unified benchmark ensures the relative positional accuracy of the primary winding body and the intermediate layer, avoiding model distortion caused by component misalignment. The resulting complete virtual model provides a realistic geometric carrier for subsequent steps such as finite element mesh generation and thermal coupling parameter extraction.

[0046] In a preferred embodiment of the present invention, step 3 includes:

[0047] Step 300: On the outer surface of the virtual three-dimensional geometric model, set an array of spatial positioning marks that are uniformly distributed with the central axis of symmetry of the blanket as the reference, wherein the longitudinal spacing between adjacent marks is no more than 0.5 meters;

[0048] Step 301: Based on the spatial coordinates of the spatial positioning marker array, generate a three-dimensional laying path model representing the spatial arrangement of the temperature sensing wires;

[0049] Step 302 involves meshing the 3D laying path model for finite element heat conduction analysis, discretizing the model into a mesh structure composed of nodes and connecting elements. Specifically, this includes:

[0050] Step 3020: Based on the spatial curve features of the three-dimensional laying path model, generate discrete nodes with equal spacing along the path curve.

[0051] Step 3021: Based on the spatial vector relationship between adjacent spatial positioning markers, construct a linear heat conduction unit connecting discrete nodes;

[0052] Step 3022: The linear heat conduction unit is radially expanded into a hexahedral unit group with an equivalent heat conduction cross section. Through the spatial coordinates of the discrete nodes and the topological connection relationship of the hexahedral unit group, a discretized mesh structure containing thermophysical properties is formed.

[0053] Step 303: Extract thermal coupling parameters representing the heat transfer characteristics between adjacent grid nodes from the discretized grid structure. The thermal coupling parameters include inter-node thermal conductivity, inter-node heat capacity, and inter-node thermal resistance.

[0054] In this embodiment of the invention, in 3D modeling software (such as ANSYS, SolidWorks), the "symmetry axis recognition" function is called to automatically identify the central symmetry axis based on the overall geometric features of the virtual model of the heating blanket (such as a symmetrical rectangular outline and a uniformly distributed internal structure). This axis must satisfy the requirement that "the geometric features on both sides are completely symmetrical" (such as the arrangement of the temperature sensing lines on the left and right sides and the distribution of the filling material are completely consistent). After confirmation, the symmetry axis is aligned with the global coordinate system Z-axis established in step 201 (that is, the spatial direction of the symmetry axis coincides with the Z-axis), and the coordinates of the starting point (aligned with the starting end of the temperature sensing line model) and the ending point (aligned with the ending end of the temperature sensing line model) of the axis are recorded.

[0055] Starting from the beginning of the virtual model of the temperature sensing line (the end where the origin (0, 0, 0) is located), use the software's "Point Creation" tool to set marks sequentially along the Z-axis: First, create the first mark on the outer surface of the beginning end, with coordinates (0, 1.691mm, 0) (where "1.691mm" is the radius of the outer wall of the intermediate layer, i.e., the distance from the Z-axis to the outer surface, ensuring the mark is on the outer surface of the model); move 0.3 meters along the Z-axis (the spacing should be selected according to accuracy requirements, ≤0.5 meters), and create the second mark on the outer surface at the corresponding position. Mark the model with coordinates (0, 1.691mm, 300mm). Continue this process, creating a marker every 0.3 meters until the model's entire length is covered (e.g., if the temperature sensing line is 10 meters long, 34 markers are needed: 10 meters ÷ 0.3 meters ≈ 33.3, so use 34 to cover the end). If the distance between the end of the model and the last marker is less than 0.3 meters (e.g., if the end coordinates are 9950mm and the last marker is at 9900mm), then add an extra marker on the outer surface of the end to ensure that the markers cover the entire length of the model without any omissions.

[0056] For the intermediate layer structure of the cylindrical outer surface (outer diameter 3.382mm), each mark must be located on the "genealogy line" of the outer surface (i.e., a straight line parallel to the Z-axis, extending along the side of the cylinder); taking the first mark as an example, the radius of the outer surface of the cylinder (1.691mm) is determined by the software's "measurement" tool, and a point (i.e., genealogy line) is created on the intersection line of the XZ plane passing through the Z-axis and the outer surface, ensuring that the X / Y values ​​of the mark coordinates satisfy "distance to the Z-axis = 1.691mm" (e.g., (1.691mm, 0, 0), (0, 1.691mm, 0) etc.).

[0057] To ensure multi-angle positioning, a circular array of markers needs to be set at each axial position: Centered on the Z-axis, create four markers (forming a circular array) at 90° intervals on the outer surface of each axial position (e.g., Z=0, Z=300mm), with coordinates (1.691mm, 0, z), (0, 1.691mm, z), (-1.691mm, 0, z), and (0, -1.691mm, z) (z is the axial coordinate). After all markers are created, check the three-dimensional coordinate accuracy of each marker using the software's "coordinate verification" function to ensure that the X / Y coordinate error is ≤0.01mm (i.e., the distance deviation to the Z-axis is ≤0.01mm) and the Z-axis coordinate error is ≤0.1mm (i.e., the axial spacing deviation is ≤0.1mm), thus ensuring the accuracy of the marker's spatial positioning.

[0058] Step 301: Based on the coordinates of the spatial positioning markers, fit the actual spatial arrangement path of the temperature sensing wires in the heating blanket. The specific operation is as follows:

[0059] In the 3D modeling software, use the "Coordinate List Export" function to batch extract the 3D coordinates of all spatial positioning marks in step 300. The coordinates of each mark are named as "Mark Number - Axial Position - Ring Position" (e.g., "Mark_1_Z0_0°" represents the first mark with an axial position and a 0° direction). The coordinate values ​​are accurate to 0.001mm (e.g., the coordinates of Mark_1_Z0_0° are (1.691mm, 0.000mm, 0.000mm), and the coordinates of Mark_2_Z300_90° are (0.000mm, 1.691mm, 300.000mm)). Export all coordinate data to an Excel spreadsheet containing columns for "Mark Number, X (mm), Y (mm), Z (mm), Axial Position (Z value), Ring Direction (0° / 90° / 180° / 270°)".

[0060] Group the markers by axial position (e.g., Z=0-300mm, Z=300-600mm, etc.), extract the X / Y coordinates of four ring marks in each group, and calculate the X / Y deviation of different ring marks at the same axial position. Specifically, for four ring marks (0°, 90°, 180°, 270°) at a certain axial position (e.g., Z=300mm), select the X-coordinate with the largest and smallest values ​​from these four marks, and calculate the difference between these two values ​​(i.e., X-coordinate deviation); then select the Y-coordinate with the largest and smallest values ​​from these four marks, and calculate the difference between these two values ​​(i.e., Y-coordinate deviation). If a certain group of marks... If the difference between the maximum and minimum X-coordinate values ​​is ≤0.1mm, and the difference between the maximum and minimum Y-coordinate values ​​is ≤0.1mm (e.g., in the Z=300-600mm segment, the X-coordinate is between 1.691-1.692mm, and the Y-coordinate is between 0.000-0.001mm), it is judged as a "flat area" (the temperature sensing line has no obvious curvature); if the X or Y coordinate deviation of a certain group of marks is >0.1mm (e.g., in the Z=600-900mm segment, the X-coordinate is between 1.691-2.000mm, with a deviation of 0.309mm), it is judged as a "wrinkled area" (the temperature sensing line is bent due to wrinkles in the heating blanket); among them, the flat area is selected using the software's "linear fitting" tool. Generate a straight line from the coordinates of all markers within the same axial segment (e.g., all markers in the Z=300-600mm segment). This straight line must pass through the midpoint marker of the segment and smoothly connect to the adjacent segment paths at both ends. For example, the starting point of the straight line in the Z=300-600mm segment is (1.691mm, 0.000mm, 300mm), and the ending point is (1.691mm, 0.000mm, 600mm), ensuring consistency with the straight line arrangement of the temperature sensing lines in a flat state. For wrinkled areas, use the "cubic spline curve" tool, with the markers within the segment as control points (e.g., the Z=600-900mm segment has three axially positioned markers), and set the "curvature". The "continuity" constraint (the tangent direction of adjacent curve segments is consistent) is used. For example, the coordinates of 600mm are marked as (1.691mm, 0.000mm, 600mm), the coordinates of 750mm (the fold apex) are marked as (2.000mm, 0.500mm, 750mm), and the coordinates of 900mm are marked as (1.700mm, 0.001mm, 900mm). The fitted spline curve needs to smoothly rise from 600mm to 750mm (simulating the fold bulge) and then smoothly descend to 900mm. The maximum curvature of the curve should not exceed the physical bending limit of the temperature sensing line (e.g., radius ≥ 50mm, to avoid excessive bending that could damage the structure).

[0061] After selecting the fitted 3D laying path using the software's "Curvature Analysis" tool, the software will automatically calculate the curvature of each point on the path. First, a target point is selected on the path (e.g., a point at Z=500mm with 3D coordinates (x, y, 500)). Then, two adjacent points are selected before and after this target point (e.g., a point at Z=499mm with coordinates (x1, y1, 499); a point at Z=501mm with coordinates (x2, y2, 501)). The tangent direction vector is then calculated using vector calculations, as follows:

[0062] The tangent direction vector of the preceding adjacent point (the point at Z=499mm) is obtained by subtracting the x1 coordinate of the preceding adjacent point from the x coordinate of the target point, resulting in the x component of the vector; by subtracting the y1 coordinate of the preceding adjacent point from the y coordinate of the target point, resulting in the y component of the vector; and by subtracting the Z coordinate (499mm) of the preceding adjacent point from the Z coordinate (500mm) of the target point, resulting in the Z component of the vector (1mm). The combination of these three components forms the vector from the adjacent point to the target point, i.e., the tangent direction vector of the preceding adjacent point. The tangent direction vector of the following adjacent point (the point at Z=501mm) is obtained by subtracting the x coordinate of the target point from the x2 coordinate of the following adjacent point, resulting in the x component of the vector; by subtracting the y coordinate of the target point from the y2 coordinate of the following adjacent point, resulting in the y component of the vector; and by subtracting the Z coordinate (500mm) of the target point from the Z coordinate (501mm) of the following adjacent point, resulting in the Z component of the vector (1mm). The combination of these three components forms the vector from the target point to the adjacent point, i.e., the tangent direction vector of the following adjacent point. The tangent direction vector of the target point is obtained by first adding the x-components of the vectors obtained from the tangent directions of the preceding and following points, and then adding the y-components of the two vectors. Similarly, the y-components of the two vectors are added together, and the y-components of the two vectors (both 1 mm) are added together, resulting in the y-component (2 mm). The sum of the squares of the sum of the x, y, and Z components is then calculated, and the square root of this sum is taken to obtain the length of the sum vector. Finally, the sum of the x, y, and Z components is divided by the length of the sum vector, and the three values ​​are combined to obtain the tangent direction vector (unit vector) of the target point, reflecting the direction of the curve at the target point.

[0063] When calculating the angle between the tangent vector of the target point and the tangent vector of the preceding adjacent point, and the angle between the tangent vector of the target point and the tangent vector of the following adjacent point, the software uses the following method: For two vectors (such as the tangent direction vector of the target point and the tangent direction vector of the preceding adjacent point), assuming the tangent direction vector of the target point is (0.001, 0.001, 1) and the tangent direction vector of the preceding adjacent point is (0.0008, 0.0009, 1), calculate the dot product of the two vectors: that is, multiply the x-components of the two vectors (0.001 × 0.0008), add the multiplication of the y-components (0.001 × 0.0009), and add the multiplication of the z-components (1 × 1 = 1), the sum is (0.001 × 0.0008). ) + (0.001 × 0.0009) + 1; Calculate the "magnitude" (i.e., the square root of the sum of the squares of each component) of the two vectors respectively. Finally, divide the dot product of the two vectors by the product of the two magnitudes to obtain a cosine value. Using the inverse cosine function formula "angle θ = arccos(dot product of two vectors ÷ (magnitude of the first vector × magnitude of the second vector))", calculate the angle between the target point and the tangent direction vector of the previous adjacent point (e.g., θ1 = 0.5°) and the angle between the target point and the tangent direction vector of the next adjacent point (e.g., θ2 = 0.7°). Then add these two angles and divide by 2 to obtain the average angle. This average angle is used to represent the overall degree of change of the tangent direction at the target point. Measure the arc length (actual path length) between the target point and the previous adjacent point, and between the target point and the next adjacent point. Take the average of the two arc lengths as the reference length. Finally, divide the average angle mentioned above by the reference length to obtain the curvature value of that point (unit: mm). -1 — The larger the curvature value, the more pronounced the curvature of the curve at that point. The software arranges the curvature values ​​of all points in order of their path positions to generate a curvature distribution curve (the horizontal axis represents the position coordinates on the path, such as Z=0mm, Z=100mm, etc., and the vertical axis represents the curvature value at the corresponding position); if there is a sudden change in the curvature value at a certain point in the curvature distribution curve (e.g., the curvature at Z=600mm changes from 0.01mm...), the curvature distribution curve will show a more pronounced curve. -1 It suddenly jumped to 0.1 mm at Z=601 mm. -1 The variation range reached 0.09mm. -1If the abrupt change occurs, the fitting parameters near the abrupt change point need to be adjusted. For example, if the abrupt change point is located between Mark_600_0° (Z=600mm) and Mark_700_0° (Z=700mm), and there is a Mark_650_0° (Z=650mm) in between, the weight of Mark_650_0° can be reduced from 1.0 to 0.8 using the software's "Control Point Weight Adjustment" function (higher weight means the curve is closer to the mark, lower weight means the curve's fit to the mark is weaker), making the curve in that area more gently curved. After adjustment, regenerate the curvature distribution curve and then calculate the rate of change of curvature: select two adjacent points on the curve (e.g., points at Z=550mm and Z=560mm), and read the curvature values ​​of the two points (e.g., 0.03mm each). -1 and 0.05mm -1 ), calculate the difference between the two (0.02mm) -1 ); Measure the actual arc length (e.g., 10mm) of these two points on the path, divide the curvature difference by the arc length to obtain the rate of curvature change, and ensure that the rate of curvature change of all adjacent points is ≤0.02mm. -1 / mm, ultimately making the curvature of the entire path change continuously and gently from one point to the next, which is completely consistent with the shape of the temperature sensing line naturally bending due to the folds when it is actually laid in the heating blanket.

[0064] Step 3020: In the 3D modeling software, call the "Curve Length Measurement" function, select the 3D laying path generated in step 301, and the software automatically calculates the total arc length of the path (i.e., the actual length along the curve). For example, if the path is a combination of a 10-meter-long broken line and a curve, the measurement result is 10002.5mm (including the extra length caused by folds). Record this value as the total path length. Based on the accuracy requirements of the heat conduction analysis (e.g., capturing local hot spots requires high accuracy), set the node spacing. If the analysis focuses on global heat distribution, the spacing can be set to 100mm (total number of nodes = 10002.5mm ÷ 100mm + 1 ≈ 101). If it is necessary to analyze local heat flow in folded areas, reduce the spacing to 50mm in folded sections (e.g., Z = 600-900mm) (100mm remains in non-folded areas) to ensure sufficient local node density. Using the software's "Equal Arc Length Point Creation" tool, discrete nodes are generated along the 3D laying path at set intervals: starting from the beginning of the path (aligned with the beginning of the temperature sensing line, Z=0mm), the coordinates of the first node are the coordinates of the path's starting point (e.g., (1.691mm, 0.000mm, 0.000mm)). Moving 100mm along the path towards the end, the second node is created, and its coordinates are recorded (e.g., (1.691mm, 0.000mm, 100.000mm)). For curved path segments (e.g., folded areas), the software automatically takes points along the curve's arc length to ensure that the arc length of adjacent nodes remains 100mm (even if the distance in the Z-axis direction is less than 100mm). For example, if the Z coordinate difference between two nodes in a folded area is 90mm, the arc length is still 100mm, ensuring that the nodes are evenly distributed on the path. After all nodes are created, they are numbered sequentially, and a coordinate list (containing X / Y / Z values, accurate to 0.01mm) is exported.

[0065] Step 3021: Based on the spatial vector relationship of adjacent spatial positioning markers, connect discrete nodes to form basic heat conduction units. The specific operation is as follows:

[0066] From the list of marker coordinates in step 300, select two adjacent axial markers (such as markers at Z=300mm and Z=600mm), and calculate their spatial vectors: Let marker i (Z=300mm, coordinates (xᵢ, yᵢ, 300mm)) and marker i+1 (Z=600mm, coordinates (xᵢ, yᵢ, 300mm)) be the coordinates of the markers in step 300. i+1 y i+1 Taking 600mm as an example, the vector obtained by using the software's "vector calculation" function is "(x i+1 -xᵢ,y i+1-yᵢ, 300mm) – The direction of this vector (e.g., positive X-axis, negative Y-axis) reflects the overall orientation of the three-dimensional laying path (e.g., tilted to the upper right). For the discrete nodes generated in step 3020, use the software's "Line Element Creation" tool to connect adjacent nodes in numerical order to form linear heat conduction elements. When connecting, ensure that the element direction is consistent with the spatial vector direction of the segment: for example, in the segment marked i to i+1 (vector direction is positive X-axis), that is, connect the 5th discrete node and the 6th discrete node in sequence to form a linear heat conduction element. The heat conduction unit needs to be tilted towards the positive X-axis, with an angle ≤5° with the vector direction (verified using the software's "angle measurement" tool) to ensure that the unit accurately reflects the direction of heat transfer along the three-dimensional laying path. For each linear heat conduction unit, record its geometric parameters: namely, the length (measured by software to measure the spatial distance between adjacent nodes, i.e., arc length, consistent with the node spacing set in step 3020, such as 100mm), the unit's direction vector (consistent with the vector of the line connecting the nodes), and marking whether the unit belongs to a flat area or a wrinkled area (e.g., "wrinkled area -600-900mm").

[0067] Step 3022: Extend the linear heat conduction unit radially to form a three-dimensional unit simulating the actual heat conduction cross section. The specific operation is as follows:

[0068] Based on the total outer diameter of the assembled temperature sensing wire (3.382 mm) in step 201, the radial expansion distance is determined: half of the total outer diameter is 1.691 mm (i.e., the distance from the central axis to the outer surface). Therefore, taking the linear heat conduction unit (central axis) as the reference, the radial expansion is extended 1.691 mm to both the inner and outer sides (in the plane perpendicular to the unit direction) to ensure that the expanded range completely covers the primary winding body (diameter 2.422 mm) and the intermediate layer structure (thickness 0.3 mm). Using the software's "sweep" function, each linear heat conduction unit is radially expanded into a hexahedral unit: For the uniform gap region (the middle 60cm), the cross-section is rectangular (length = equivalent length of the perimeter of the sensing line cross-section, width = 1.691mm), ensuring the cross-sectional area matches the actual heat conduction area (e.g., if the actual cross-sectional area is 8.9mm², the hexahedral cross-sectional area error ≤ 0.1mm²); for the non-uniform regions of the first 20cm (gap 0.25mm) and the last 20cm (gap 0.23mm), the cross-sectional dimensions are adjusted using the software's "local scaling" function—the cross-sectional width of the first 20cm region is increased by 0.1mm compared to the middle region (corresponding to wider gaps), and the last 20cm region is increased by 0.08mm, ensuring the cross-sectional dimensions change synchronously with the gap width. Using the software's "contact setting" function, the connection method between adjacent hexahedral units is defined as "bound contact" (i.e., no relative sliding between units, heat energy can be directly transferred), ensuring a continuous heat energy transfer path. The software's "material library" is invoked to associate the corresponding material properties with the unit: the central area unit (corresponding to the heating wire of the primary winding) is assigned copper-nickel alloy properties (thermal conductivity 40W / (m・K), specific heat capacity 450J / (kg・K)); the outer layer unit (corresponding to the gap between the middle layers) is assigned air properties (thermal conductivity 0.026W / (m・K), specific heat capacity 1005J / (kg・K)). The software's "property verification" function is used to check and ensure that the material properties of each unit are consistent with the corresponding area (primary winding / middle layer).

[0069] Step 303: Obtain the thermal energy transfer characteristic parameters of adjacent nodes from the discretized grid structure. The specific operation is as follows:

[0070] For each adjacent node (e.g., i and i+1), the material properties and geometric dimensions of the hexahedral element connecting them are obtained through the software's "element query" function: the thermal conductivity is read from the element properties (e.g., 0.026 W / (m·K) for an air element); the cross-sectional area is obtained through the software's "cross-sectional measurement" tool (e.g., 0.00001 m²); and the element length, i.e., the node spacing (e.g., 0.1 m). Following the logic of "thermal conductivity = material thermal conductivity × cross-sectional area ÷ element length," the software automatically calculates and records the data.

[0071] The software reads the specific heat capacity (e.g., 1005 J / (kg・K) for air) and density (1.29 kg / m³) from the unit properties; obtains the unit volume (e.g., 0.000001 m³) using the software's "Volume Measurement" tool, and calculates the unit mass (density × volume = 1.29 × 0.000001); automatically calculates and records the thermal resistance (the reciprocal of thermal conductivity) according to the logic of "heat capacity = specific heat capacity × mass" (e.g., 1005 × 1.29 × 0.000001); automatically calculates the thermal resistance (the reciprocal of thermal conductivity) and records the value; and enters the thermal conductivity, heat capacity, and thermal resistance parameters of all adjacent nodes into an Excel spreadsheet in the order of node number, forming a matrix form of "node i - node i+1 - thermal conductivity - heat capacity - thermal resistance".

[0072] A spatial positioning marker array is set up with the central symmetry axis of the blanket as the reference, and the longitudinal spacing between adjacent markers is controlled to accurately locate the spatial position of the temperature sensing lines within the blanket, ensuring that the path model is consistent with the actual layout. A three-dimensional laying path model is generated based on the coordinates of the spatial positioning markers. Through smoothing, the path conforms to the actual bending characteristics of the temperature sensing lines, which can realistically reflect the spatial arrangement of the temperature sensing lines. Mesh generation is performed step by step, generating equally spaced discrete nodes to ensure mesh uniformity. The construction of linear heat conduction elements ensures the accuracy of heat conduction direction, and the expansion to hexahedral element groups improves the accuracy of heat conduction analysis. The final discretized mesh structure can accurately simulate the heat conduction characteristics of the temperature sensing lines.

[0073] In a preferred embodiment of the present invention, step 4 includes:

[0074] Step 400: Based on the inter-node thermal conductivity, inter-node heat capacity, and inter-node thermal resistance parameters, calculate the local heat flux density along the three-dimensional laying path of the sensing wire using finite element heat conduction analysis. Specifically, this includes:

[0075] Step 4010: Based on the inter-node thermal conductivity parameters, inter-node heat capacity parameters, and inter-node thermal resistance parameters, establish the heat conduction path between adjacent grid nodes;

[0076] Step 4011: Determine the direction of heat transfer based on the heat conduction path, and calculate the axial temperature change rate along the three-dimensional laying path.

[0077] Step 4012: Generate the heat flux density values ​​for each local region based on the axial temperature change rate and the inter-node thermal conductivity parameters.

[0078] Step 401: Compare the heat flux density of each local region with the preset heat flux density threshold range to obtain the comparison result; generate gap region parameter adjustment instructions based on the comparison result, wherein when the heat flux density of the corresponding local region is higher than the preset upper limit threshold of heat flux density, an adjustment instruction to increase the gap geometry parameters and decrease the equivalent thermal conductivity of the gap filling material is generated; when the heat flux density of the corresponding local region is lower than the preset lower limit threshold of heat flux density, an adjustment instruction to decrease the gap geometry parameters and increase the equivalent thermal conductivity of the gap filling material is generated.

[0079] Step 402: Execute the gap region parameter adjustment command to dynamically modify the preset gap size parameters of the corresponding region in the virtual three-dimensional geometric model; the preset gap size parameters include gap geometric size parameters and the equivalent thermal conductivity parameters of the gap filling material.

[0080] In this embodiment of the invention, firstly, the "inter-node thermal conductivity, heat capacity, and thermal resistance parameter matrix" (containing the heat transfer characteristic parameters of all adjacent node pairs) generated in step 303 is called, and associated with the coordinate information of the discrete nodes in step 3020, clarifying the spatial positional relationship and heat transfer parameters of each node pair (such as node 1 and node 2, node 2 and node 3, etc.); based on the linear heat conduction unit constructed in step 3021, all adjacent grid nodes with direct connections are determined (i.e., node pairs connected by linear units), excluding nodes that are not directly connected (such as node 1 and node 3, which are not included for now); the "heat transfer capacity weight" of the path is determined based on the inter-node thermal conductivity and thermal resistance parameters. The specific weighting rules are as follows: Node pairs with thermal conductivity ≥ 0.000002 W / Kelvin and thermal resistance ≤ 500000 Kelvin / W have a weight of 0.8-1.0 (the higher the thermal conductivity and the lower the thermal resistance, the closer the weight is to 1.0), indicating that heat is more easily transferred along this path, and therefore they are marked as "preferred heat conduction paths"; Node pairs with thermal conductivity < 0.000002 W / Kelvin or thermal resistance > 500000 Kelvin / W have a weight of 0.3-0.7 (the lower the thermal conductivity and the higher the thermal resistance, the closer the weight is to 0.3), and are marked as "weak heat conduction paths". All adjacent node pairs are connected according to the above weighting relationship to form a heat conduction path network covering the entire three-dimensional laying path. For example, nodes 5 and 6 have a thermal conductivity of 0.000003 W / Kelvin and a thermal resistance of 333333 Kelvin / W, which meets the priority path criteria and has a weight of 0.9, represented by a bold line in the network. On the other hand, a pair of nodes in a certain fold region has a thermal conductivity of 0.000001 W / Kelvin and a thermal resistance of 1,000000 Kelvin / W, which belongs to a weak heat conduction path and has a weight of 0.4, represented by a thin line. This visually distinguishes the heat transfer capabilities of different paths.

[0081] Step 4011: Based on the heat conduction path network constructed in step 4010 and the initial temperature distribution of the sensing wire (e.g., the initial temperature of the sensing wire is set to 30℃ when the starting end is energized and the initial temperature of the ending end is set to the ambient temperature of 25℃), determine the direction of heat transfer for each pair of adjacent nodes. Heat is transferred from the node with higher temperature to the node with lower temperature. For example, if the initial temperature of node 5 is 32℃ and the initial temperature of node 6 is 31℃, then heat is transferred from node 5 to node 6. For adjacent nodes with the same temperature (e.g., both are 30℃), the thermal resistance parameter is used to determine the direction of heat transfer: that is, the side with lower thermal resistance is more likely to receive heat (i.e., heat tends to flow towards the node with lower thermal resistance).

[0082] Select a continuous segment of nodes (e.g., node i, node i+1, node i+2) on the three-dimensional laying path. Use the "Temperature Field Simulation" function of the finite element analysis software to obtain the temperature values ​​of each node under steady-state conditions (e.g., node i is 35℃, node i+1 is 34℃, and node i+2 is 33℃). Calculate the temperature difference between adjacent nodes (e.g., the temperature difference between node i and node i+1 is 1℃), and measure the axial distance between these two nodes along the three-dimensional laying path (i.e., the arc length along the path, e.g., 100 mm). Divide the temperature difference by the axial distance to obtain the axial temperature change rate of this segment of the path (e.g., 1℃ ÷ 100 mm = 0.01℃ / mm), which represents the temperature decrease per millimeter of path length. Repeat the above operation for all continuous node pairs along the entire path to obtain the axial temperature change rate of each local area.

[0083] Step 4012: Obtain the axial temperature change rate of each local area on the three-dimensional laying path calculated in step 4011, and correlate it with the thermal conductivity parameters corresponding to each node pair obtained in step 303 (for example, the axial temperature change rate of a certain fold area is 0.02℃ / mm, and the thermal conductivity of its corresponding node pair is 0.000002W / Kelvin); where heat flux density is the amount of heat passing through a unit area per unit time, and its magnitude is positively correlated with the thermal conductivity of the material (i.e., thermal conductivity parameter) and the rate of temperature change (i.e., axial temperature change rate): the higher the thermal conductivity parameter and the more drastic the axial temperature change, the greater the corresponding heat flux density. The three-dimensional laying path is divided into several local regions along its axial direction (e.g., each region is 500 mm long, or it is divided according to natural morphology such as flat areas and folded areas). For each local region, the thermal conductivity parameters and axial temperature change rate of all node pairs in the region are summarized. The average value of the thermal conductivity parameter (i.e., the average thermal conductivity of the region) and the average value of the axial temperature change rate (i.e., the average axial temperature change rate of the region) are calculated respectively. These two are used as representative parameters of the region (e.g., the average thermal conductivity of a folded area is 0.0000025 W / Kelvin, and the average axial temperature change rate of the region is 0.015℃ / mm). The specific method for obtaining the heat flux density value of the region is as follows: multiply the average thermal conductivity of the region by the average axial temperature change rate of the region to obtain the amount of heat transferred per unit length per unit time; divide the amount of heat transferred per unit length per unit time by the cross-sectional area parameter of the sensing line to obtain the amount of heat passing through a unit area per unit time, which is the heat flux density value of the region. For example, based on the above calculations, the heat flux density value of a certain folded area is "high" (specific value such as 35 watts / square meter), and the heat flux density value of a certain flat area is "medium" (specific value such as 20 watts / square meter). The heat flux density values ​​of each area are recorded according to the area number (e.g., "Area 1 - Heat flux density: 20 watts / square meter" "Area 2 - Heat flux density: 35 watts / square meter").

[0084] Step 401: Based on the safety operating standards of the temperature sensing wire, a preset heat flux density threshold range is set. This range includes an upper threshold and a lower threshold (e.g., the upper threshold is 50 watts / m², and the lower threshold is 10 watts / m²). When the heat flux density exceeds the upper threshold, it is easy to cause overheating in local areas of the temperature sensing wire. When the heat flux density is lower than the lower threshold, it may affect the temperature sensing sensitivity of the temperature sensing wire. The heat flux density values ​​of each local area generated in step 4012 are compared one by one with the preset heat flux density threshold range. If the heat flux density value of a local area is higher than the upper threshold (e.g., 60 watts / m² > 50 watts / m²), the area is marked as an "overheated area". If the heat flux density value of a local area is lower than the lower threshold (e.g., 8 watts / m² < 10 watts / m²), the area is marked as an "insufficient heat transfer area". If the heat flux density value of a local area is between the upper and lower thresholds (e.g., 10 watts / m² ≤ heat flux density value ≤ 50 watts / m²), the area is marked as a "normal area".

[0085] For the "overheated area," to reduce the heat transfer efficiency in this area and avoid localized overheating, the following two gap area parameter adjustment instructions are generated: The first increases the gap geometry in this area (e.g., adjusting the original gap width from 0.2 mm to 0.3 mm) to reduce the contact area between materials, thereby reducing heat transfer efficiency; the second decreases the equivalent thermal conductivity of the gap filling material in this area (e.g., replacing the original air filling with an insulating material with lower thermal conductivity) to further hinder heat transfer. For the "insufficient heat transfer area," to improve the heat transfer efficiency in this area and ensure temperature sensing sensitivity, the following two gap area parameter adjustment instructions are generated: The first decreases the gap geometry in this area (e.g., adjusting the original gap width from 0.3 mm to 0.2 mm) to increase the contact area between materials, thereby enhancing heat transfer efficiency; the second increases the equivalent thermal conductivity of the gap filling material in this area (e.g., replacing the original insulating material with a thermally conductive gel with higher thermal conductivity) to promote heat transfer. For the "normal area", no gap area parameter adjustment command is generated, and the original gap geometry parameters and the equivalent thermal conductivity parameters of the gap filling material are maintained in this area.

[0086] Step 402: Receive the gap region parameter adjustment instruction generated in step 401, specifying the local region to be modified (e.g., "region 3 (Z=1500-2000 mm)"), the parameter adjustment type (including increasing / decreasing gap geometry parameters, increasing / decreasing the equivalent thermal conductivity parameter of gap filling material), and the specific parameter change value (e.g., adjusting the gap width from 0.2 mm to 0.3 mm). In the virtual three-dimensional geometric model, based on the coordinate information of the local region to be modified, locate the gap structure corresponding to that local region (e.g., the gap in region 3 is located between the primary winding body and the intermediate layer). Using the dimension editing function of the 3D modeling software, the gap geometry parameters such as the gap width and thickness of the area are modified according to the specific parameter changes in the gap region parameter adjustment instructions (e.g., from 0.2 mm to 0.3 mm), and the software automatically updates the model's geometry. Using the material property editing function of the 3D modeling software, the equivalent thermal conductivity parameter of the gap filling material in the area is modified according to the specific parameter changes in the gap region parameter adjustment instructions (e.g., from 0.026 W / (m·Kelvin) to 0.01 W / (m·Kelvin)), and the thermal conductivity unit properties of the area are updated simultaneously. After the parameters are modified, the parameter checking function of the 3D modeling software is used to confirm that the gap geometry parameters and the equivalent thermal conductivity parameters of the gap filling material in the area to be modified have been updated according to the gap region parameter adjustment instructions, and that the parameters of the area to be modified transition smoothly with the parameters of the surrounding area (without parameter abrupt changes), to ensure that the virtual 3D geometric model is consistent with the adjusted parameters.

[0087] By establishing heat conduction paths step by step, calculating temperature change rates, and correlating with thermal conductivity, the calculation of local heat flux density is made more consistent with the actual heat transfer characteristics of the sensing wire. By comparing heat flux density with thresholds and generating targeted adjustment instructions, areas that are overheated or underheated can be accurately improved, ensuring uniform heat distribution along the three-dimensional laying path of the sensing wire and avoiding local overheating or insufficient sensitivity. The gap parameters of the virtual model are dynamically modified so that the model can reflect the adjusted heat transfer state in real time, providing an intuitive and operable virtual basis for optimizing the structural design and verifying the safety performance of the sensing wire.

[0088] In a preferred embodiment of the present invention, step 5 includes:

[0089] Step 500: Obtain the dynamically modified preset gap size parameters, including gap geometric size parameters and the equivalent thermal conductivity parameters of the gap filling material;

[0090] Step 501: Based on the gap geometry parameters, reconstruct the three-dimensional geometry of the intermediate layer structure, and assign the equivalent thermal conductivity parameters of the gap filling material to the reconstructed intermediate layer structure to obtain the reconstructed intermediate layer structure.

[0091] Step 502: Spatially assemble the reconstructed intermediate layer structure with the original primary winding structure to obtain the assembled updated virtual three-dimensional geometric model as the final construction result.

[0092] In this embodiment of the invention, adjusted preset gap size parameters are extracted from the virtual three-dimensional geometric model that has been dynamically modified in step 402. These preset gap size parameters include two types:

[0093] Gap geometry parameters: Divide the axial regions of the three-dimensional laying path (e.g., region 1 to region n), and record the specific geometric parameters of the gap in each region, including but not limited to width and thickness (e.g., "region 3 (Z=1500mm-2000mm): gap width 0.3mm, thickness 0.1mm" "region 5 (Z=2500mm-3000mm): gap width 0.15mm, thickness 0.08mm").

[0094] Equivalent thermal conductivity parameters of gap filling materials: For the above axial regions, record the equivalent thermal conductivity of the gap filling materials in each region (e.g., "Region 3: 0.01 W / (m·Kelvin)" "Region 5: 0.03 W / (m·Kelvin)").

[0095] The extracted parameters are verified for completeness to ensure that the parameters of the “overheated area” and “insufficient heat transfer area” marked in step 401 have been updated, and the parameters of the “normal area” remain in their original state. The parameters are summarized in order of axial area to form a “area-gap parameter correspondence table”, clarifying the mapping relationship between each area and gap geometry and material thermal conductivity, and confirming that the parameter units are consistent (geometric size unit is millimeter, thermal conductivity unit is watt / (meter·Kelvin)).

[0096] Step 501: Obtain the original 3D geometric model of the intermediate layer structure constructed in Step 2. This model includes the initial contour, the relative position reference with the primary winding body, and the original gap distribution characteristics. It serves as the basic model for reconstruction. Based on the gap geometric dimension parameters, adjust the geometry of local areas of the intermediate layer: that is, according to the "region-gap parameter correspondence table," locate each target region (e.g., region 3, region 5) in the original intermediate layer model through coordinate matching (e.g., Z-axis coordinate range); for the "overheated region" (e.g., region 3), use the local editing functions of the 3D modeling software (e.g., "offset boundary" "stretch surface") to adjust the inner boundary of the intermediate layer near the primary winding body away from the primary winding body, so that the gap width and thickness reach the target values ​​recorded in Step 500 (e.g., width 0.3m). The gap width and thickness are adjusted to the target values ​​(e.g., 0.15mm width and 0.08mm thickness) while keeping the outer boundary contour of the intermediate layer unchanged. For the "insufficient heat transfer area" (e.g., area 5), ​​the inner boundary of the intermediate layer is adjusted towards the primary winding body using the "shrink boundary" and "offset surface" functions, so that the gap width and thickness are reduced to the target values ​​(e.g., 0.15mm width and 0.08mm thickness). For the "normal area", the original geometric dimensions are kept unchanged, and the boundary between it and the adjacent adjusted area is smoothed (e.g., sharp corners are eliminated using the "rounded corner transition" function) to ensure the continuity of the overall shape of the intermediate layer. The size measurement tool of the software is used to verify whether the adjusted gap width and thickness of each area are consistent with the target parameters. The cross-sectional analysis function is used to check the cross-sectional contour of the intermediate layer to confirm that the geometric shape meets the design requirements.

[0097] In the material properties function of the 3D modeling software, the equivalent thermal conductivity of the gap filling material is associated with each region of the reconstructed intermediate layer. Region 3 is specified as "0.01 W / (m·Kelvin)" and region 5 is specified as "0.03 W / (m·Kelvin)". The original material parameters are kept for the "normal region". The above thermal conductivity parameters are associated with the corresponding geometric regions of the intermediate layer through the attribute binding function.

[0098] Step 502: Obtain the original three-dimensional geometric model of the primary winding body constructed in Step 1. The geometry and spatial position of this model remain unchanged. Use it as the assembly reference part to reconstruct the spatial assembly of the intermediate layer and the primary winding body. The specific process is as follows:

[0099] Using the central axis and end face center of the primary winding body as reference features, corresponding mating references (such as the end face center of the intermediate layer coaxial with the end face center of the primary winding body) are set on the reconstructed intermediate layer structure to ensure spatial alignment between the two. Through the assembly constraint functions of the 3D modeling software (such as "coaxial center" and "face-to-face distance"), the reconstructed intermediate layer structure is spatially assembled with the primary winding body according to a preset relative positional relationship (such as the intermediate layer wrapping around the outside of the primary winding body, with both arranged coaxially). During assembly, the software automatically detects and avoids geometric interference (such as overlap between the intermediate layer and the primary winding body). If interference exists, the process returns to the previous step. 501. Readjust the geometry of the intermediate layer. Use a distance measurement tool to check the actual gap size of each area after assembly (e.g., whether the gap in area 3 is 0.3mm). Use a coaxiality detection tool to confirm the axis deviation between the intermediate layer and the primary winding body (e.g., deviation ≤ 0.01mm). Ensure that the assembly accuracy meets the design requirements. After the assembly is completed and verified to be correct, save the integrated model in a preset three-dimensional format (e.g., STL or STEP format). This model includes the updated intermediate layer structure (including the adjusted geometry and material parameters) and the original primary winding body structure, serving as the final updated virtual three-dimensional geometric model.

[0100] By acquiring dynamically adjusted gap parameters and reconstructing the intermediate layer, the model can accurately reflect the structural features and material properties after gap optimization, ensuring that the model is consistent with actual design requirements. From parameter adjustment to structural reconstruction and then to spatial assembly, a complete design iteration process is formed, so that the optimization results of the intermediate layer are directly reflected in the final model, improving design efficiency. Through benchmark matching and accuracy verification, the relative position of the reconstructed intermediate layer and the primary winding body is ensured to be accurate, avoiding model distortion caused by assembly errors.

[0101] like Figure 2 As shown, embodiments of the present invention also provide a three-dimensional construction system for temperature sensing lines based on finite element analysis, comprising:

[0102] The data acquisition module is used to collect physical parameter data such as the core properties of the temperature sensing wire, the temperature coefficient of the heating wire, and the preset gap size.

[0103] The modeling module is used to construct a virtual three-dimensional geometric model of the temperature-sensing wire, which includes a primary winding body and an intermediate layer structure, based on physical parameter data.

[0104] The extraction module is used to set a spatial positioning mark array on the surface of a virtual three-dimensional geometric model. The mark array is uniformly distributed around the central axis of symmetry of the blanket, and the longitudinal spacing between adjacent marks is no more than 0.5 meters. Based on the mark array, a three-dimensional laying path model is generated, and the three-dimensional laying path model is divided into finite element heat conduction meshes to extract the thermal coupling parameters between nodes.

[0105] The control module is used to calculate the local heat flux density based on the thermal coupling parameters and dynamically control the parameters of the gap region in the virtual three-dimensional geometric model.

[0106] The update module is used to update the virtual three-dimensional geometric model based on the adjusted gap region parameters, forming the final three-dimensional virtual construction result of the temperature sensing line.

[0107] It should be noted that this system is a system corresponding to the above method. All implementation methods in the above method embodiments are applicable to this embodiment and can achieve the same technical effect.

[0108] Embodiments of the present invention also provide a computing device, including: a processor and a memory storing a computer program, wherein the computer program, when executed by the processor, performs the method described above. All implementations in the above method embodiments are applicable to this embodiment and can achieve the same technical effects.

[0109] Embodiments of the present invention also provide a computer-readable storage medium storing instructions that, when executed on a computer, cause the computer to perform the method described above. All implementations in the above method embodiments are applicable to this embodiment and can achieve the same technical effects.

[0110] The above description represents the preferred embodiments of the present invention. It should be noted that those skilled in the art can make various improvements and modifications without departing from the principles of the present invention, and these improvements and modifications should also be considered within the scope of protection of the present invention.

Claims

1. A method for constructing three-dimensional temperature sensing lines based on finite element analysis, characterized in that, The method includes: Step 1: Collect physical parameter data of the temperature sensing wire, including core properties, heating wire temperature coefficient, and preset gap size; Step 2: Based on the physical parameter data, construct a virtual three-dimensional geometric model of the temperature sensing wire, including the primary winding body and the intermediate layer structure; Step 3: Set a spatial positioning mark array on the surface of the virtual three-dimensional geometric model. The mark array is evenly distributed around the central axis of symmetry of the blanket body, and the longitudinal spacing between adjacent marks is no more than 0.5 meters. Generate a three-dimensional laying path model based on the mark array, and perform finite element heat conduction mesh division on the three-dimensional laying path model to extract the thermal coupling parameters between nodes. Step 4: Calculate the local heat flux density based on the thermal coupling parameters, and dynamically adjust the gap region parameters in the virtual three-dimensional geometric model; Step 5: Based on the adjusted gap region parameters, update the virtual three-dimensional geometric model to form the final three-dimensional virtual construction result of the temperature sensing line.

2. The method for constructing three-dimensional temperature sensing lines based on finite element analysis according to claim 1, characterized in that, Step 2 includes: The core properties are analyzed to generate the three-dimensional geometry and dimensions of the primary winding structure; the preset gap structure parameters are analyzed to generate the three-dimensional geometry, dimensions, and material property definitions of the intermediate layer structure. The primary winding structure and the intermediate layer structure are spatially assembled to form a virtual three-dimensional geometric model.

3. The method for constructing three-dimensional temperature sensing lines based on finite element analysis according to claim 2, characterized in that, Step 3 includes: On the outer surface of the virtual three-dimensional geometric model, a spatial positioning mark array is set up with the central symmetry axis of the blanket body as the reference, wherein the longitudinal spacing between adjacent marks is no more than 0.5 meters; Based on the spatial coordinates of the spatial positioning marker array, a three-dimensional laying path model representing the spatial arrangement of the temperature sensing wires is generated. The mesh generation required for finite element heat conduction analysis of the three-dimensional laying path model is to discretize the three-dimensional laying path model into a mesh structure composed of nodes and connecting nodes. From the discretized mesh structure, thermal coupling parameters representing the heat transfer characteristics between adjacent mesh nodes are extracted. These thermal coupling parameters include inter-node thermal conductivity, inter-node heat capacity, and inter-node thermal resistance.

4. The method for constructing three-dimensional temperature sensing lines based on finite element analysis according to claim 3, characterized in that, The mesh generation required for finite element heat conduction analysis of the three-dimensional laying path model involves discretizing the three-dimensional laying path model into a mesh structure composed of nodes and connecting elements, including: Based on the spatial curve features of the three-dimensional laying path model, discrete nodes with equal spacing are generated along the path curve. Based on the spatial vector relationship between adjacent spatial positioning markers, a linear heat conduction unit connecting discrete nodes is constructed; The linear heat conduction element is radially expanded into a hexahedral element group with an equivalent heat conduction cross section. By using the spatial coordinates of the discrete nodes and the topological connection relationship of the hexahedral element group, a discretized mesh structure containing thermophysical properties is formed.

5. The method for constructing three-dimensional temperature sensing lines based on finite element analysis according to claim 4, characterized in that, Step 4 includes: Based on the inter-node thermal conductivity, inter-node heat capacity, and inter-node thermal resistance parameters, the local heat flux density of the sensing wire along the three-dimensional laying path is calculated through finite element heat conduction analysis. The heat flux density of each local region is compared with a preset heat flux density threshold range to obtain the comparison results. Based on the comparison results, adjustment instructions for gap region parameters are generated. Specifically, when the heat flux density of the corresponding local region is higher than the preset upper limit threshold, adjustment instructions are generated to increase the gap geometry parameters and decrease the equivalent thermal conductivity of the gap filling material. When the heat flux density of the corresponding local region is lower than the preset lower limit threshold, adjustment instructions are generated to decrease the gap geometry parameters and increase the equivalent thermal conductivity of the gap filling material. The command to adjust the gap region parameters is executed to dynamically modify the preset gap size parameters of the corresponding region in the virtual three-dimensional geometric model; the preset gap size parameters include gap geometric size parameters and the equivalent thermal conductivity parameters of the gap filling material.

6. The method for constructing three-dimensional temperature sensing lines based on finite element analysis according to claim 5, characterized in that, Based on the inter-node thermal conductivity, inter-node heat capacity, and inter-node thermal resistance parameters, the local heat flux density along the three-dimensional laying path of the sensing wire is calculated using finite element heat conduction analysis, including: Based on the inter-node thermal conductivity parameters, inter-node heat capacity parameters, and inter-node thermal resistance parameters, a heat conduction path between adjacent grid nodes is established. The direction of heat transfer is determined based on the heat conduction path, and the axial temperature change rate along the three-dimensional laying path is calculated. Based on the axial temperature change rate and the inter-node thermal conductivity parameters, the heat flux density values ​​for each local region are generated.

7. The method for constructing three-dimensional temperature sensing lines based on finite element analysis according to claim 6, characterized in that, Step 5 includes: Obtain the dynamically modified preset gap size parameters, including gap geometric size parameters and the equivalent thermal conductivity parameters of the gap filling material; Based on the gap geometry parameters, the three-dimensional geometry of the intermediate layer structure is reconstructed, and the equivalent thermal conductivity parameters of the gap filling material are assigned to the reconstructed intermediate layer structure to obtain the reconstructed intermediate layer structure. The reconstructed intermediate layer structure is spatially assembled with the original primary winding structure to obtain an updated virtual three-dimensional geometric model as the final construction result.

8. A three-dimensional construction system for sensing lines based on finite element analysis, the system implementing the method as described in any one of claims 1 to 7, characterized in that, include: The data acquisition module is used to collect physical parameter data such as the core properties of the temperature sensing wire, the temperature coefficient of the heating wire, and the preset gap size. The modeling module is used to construct a virtual three-dimensional geometric model of the temperature-sensing wire, which includes a primary winding body and an intermediate layer structure, based on physical parameter data. The extraction module is used to set a spatial positioning mark array on the surface of a virtual three-dimensional geometric model. The mark array is uniformly distributed around the central axis of symmetry of the blanket, and the longitudinal spacing between adjacent marks is no more than 0.5 meters. Based on the mark array, a three-dimensional laying path model is generated, and the three-dimensional laying path model is divided into finite element heat conduction meshes to extract the thermal coupling parameters between nodes. The control module is used to calculate the local heat flux density based on the thermal coupling parameters and dynamically control the parameters of the gap region in the virtual three-dimensional geometric model. The update module is used to update the virtual three-dimensional geometric model based on the adjusted gap region parameters, forming the final three-dimensional virtual construction result of the temperature sensing line.

9. A computing device, characterized in that, include: One or more processors; A storage device for storing one or more programs, which, when executed by one or more processors, cause the one or more processors to implement the method as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a program that, when executed by a processor, implements the method as described in any one of claims 1 to 7.

Images