Braid 3D simulation system based on finite element analysis

Abstract

The invention discloses a woven tape 3D simulation system based on finite element analysis, and relates to the technical field of woven tape simulation, multidirectional stress on a grid area can be accurately calculated by applying different types of force, comprehensive stress is constructed by using an equivalent stress theory, and the system can identify whether a stress field is smooth or not based on the comprehensive stress, so that the woven tape 3D simulation system can be used for simulating the woven tape. And a stress concentration area and a potential instability risk can be quickly found. The structure evaluation module monitors the relative slippage condition among yarns in the multilayer structure of the braid by regulating and controlling different vibration load conditions, calculates the interlayer shear stress, and fits and outputs a structure stability score in combination with an effect evaluation model, and the design can dynamically reflect the stability performance of the braid under the complex working condition, so that the stability of the braid under the complex working condition is improved. And a scientific basis is provided for optimizing design parameters (such as interweaving point density, yarn material and the like) of the braid.

Description

Technical Field

[0001] The present invention relates to the technical field of ribbon simulation, and in particular to a ribbon 3D simulation system based on finite element analysis. Background Art

[0002] As a key product widely used in the textile industry, webbing is used in many industries, including clothing, industrial lifting, automotive safety, aerospace, etc. With the continuous development of advanced manufacturing technology, the precise analysis and optimization design of the mechanical properties of webbing has gradually become a research hotspot. Especially in the development of high-performance webbing, it is particularly important to accurately simulate and evaluate the mechanical response and stability of the webbing structure.

[0003] At present, in the practice of webbing design and performance evaluation, traditional methods mainly rely on experimental tests and experience, which have problems such as long test cycles, high costs, and poor repeatability. At the same time, experimental methods are difficult to intuitively present the stress distribution and interlayer slip behavior of the internal structure of the webbing, especially under vibration loads and multi-axial stress conditions. Traditional methods are difficult to effectively capture details such as stress concentration and interlayer shear stress.

[0004] The main reasons for the above-mentioned status quo and shortcomings are the complexity of the webbing interlaced structure, the mechanical coupling between the multi-layer composite yarns, and the nonlinear behavior of yarn slip under dynamic conditions. These factors make it difficult for traditional methods to fully capture the true mechanical response of the webbing under complex working conditions. Summary of the invention

[0005] In view of the deficiencies in the prior art, the present invention provides a 3D ribbon simulation system based on finite element analysis, which solves the problems in the above-mentioned background technology.

[0006] To achieve the above objectives, the present invention is implemented through the following technical solutions: a 3D simulation system for webbing based on finite element analysis, including a model building module, a discrete division module, a finite element analysis module, a structure evaluation module and a feedback module;

[0007] The model building module is used to generate a webbing geometric model using a three-dimensional modeling software in advance, and import the webbing geometric model into a finite element simulation software, and input different working conditions and related material properties of the webbing into the finite element simulation software;

[0008] The discrete division module is used to identify each interlacing point from the webbing geometric model in the finite element simulation software, and perform a finite element mesh division operation based on each interlacing point to generate a plurality of groups of finite element mesh regions;

[0009] The finite element analysis module is used to start the simulation program according to different working conditions set in the model construction module, apply force to the webbing, analyze the deformation of each finite element mesh area in the webbing under different working conditions, and obtain multi-directional stress based on the deformation, and use the equivalent stress theory to construct the comprehensive stress Zyz of each finite element mesh area. Based on the comprehensive stress Zyz of each finite element mesh area, identify whether the stress field between each finite element mesh area in the webbing is smooth;

[0010] The structural evaluation module is used to adjust different vibration loads according to different working conditions set in the model building module, monitor the relative slip between yarns in the multi-layer structure of the webbing, so as to obtain relevant webbing state data, and construct the interlayer shear stress Jyz under different vibration load conditions based on the relevant webbing state data, and combine the trained effect evaluation model, and after linear normalization, fit the output structure stability score Jwpf;

[0011] The feedback module is used to estimate the stability performance of the internal structure of the webbing based on the structural stability score Jwpf value, and feed it back to the background center to formulate an optimization plan.

[0012] Preferably, the model building module includes a modeling unit and a data input unit;

[0013] The modeling unit is used to obtain the relevant geometric feature data of the ribbon through various detection instruments in advance, and simulate the geometric model of the ribbon in proportion to the ribbon according to the relevant geometric feature data and the three-dimensional modeling software, wherein the relevant geometric feature data includes the width, thickness, yarn arrangement, yarn color, yarn thickness, weaving pattern, spatial distribution, interlayer thickness of the ribbon , the coordinates of the interlacing points of the warp and weft yarns , the paths and interweaving density of warp and weft yarns; 3D modeling software includes SolidWorks and AutoCAD;

[0014] The data input unit is used to save the geometric model of the webbing in the STEP file format, and import the geometric model of the webbing into the finite element simulation software, and input different working conditions and related material properties of the webbing into the finite element simulation software, wherein the related material properties of the webbing include elastic modulus , Poisson's ratio , yarn tension and shear modulus ; Different working conditions include load type and load size, among which load types include longitudinal tension, lateral tension, torsional load and vibration load.

[0015] Preferably, the discrete division module includes a yarn geometric distribution analysis unit and a grid division unit;

[0016] The yarn geometric distribution analysis unit is used to identify each interlacing point from the webbing geometric model, and determine the warp space coordinates around each interlacing point in the webbing in combination with relevant geometric feature data. and the spatial position of the weft yarn , based on the spatial coordinates of the warp yarns around each interlacing point in the ribbon and the spatial position of the weft yarn , respectively obtain the curvature of the warp and weft yarns around each interlacing point and mark them as the warp curvature and weft curvature , which is obtained by the following formula: ;

[0017] In the formula, represents the velocity vector of the warp yarn, represents the acceleration vector of the warp yarn, and t represents the position point on the yarn; is the cross product modulus of the velocity vector and the acceleration vector of the warp, is the magnitude of the velocity vector of the warp yarn; ;

[0018] In the formula, represents the velocity vector of the weft yarn, represents the acceleration vector of the weft yarn, and t represents the position point on the yarn; is the cross product modulus of the velocity vector and the acceleration vector of the weft yarn; is the magnitude of the velocity vector of the weft yarn.

[0019] Preferably, the grid division unit is used to calculate the warp curvature obtained by the yarn geometric distribution analysis unit according to the relevant material properties. and weft curvature , get the discrete spacing in each direction around each interlacing point, including the discrete spacing on the warp yarn and discrete spacing on the weft yarns The specific acquisition method is as follows: ; ;

[0020] In the formula, Indicates the pressure exerted by adjacent warp yarns; Indicates the pressure exerted by adjacent weft yarns; It indicates the average pressure exerted by adjacent weft yarns at each interlacing point; It indicates the average pressure exerted by adjacent warp yarns at each interlacing point; represents the warp curvature; represents the weft yarn curvature; represents a constant; represents the adjustment coefficient;

[0021] According to the coordinates of the intersection points , along the warp direction, according to the tangent vector of the warp and discrete spacing on the warp , generate each node on the warp yarn, along the weft yarn direction, according to the curve tangent vector of the weft yarn and discrete spacing on the weft yarns , generate each node on the weft yarn, and obtain the coordinates of each node on the warp yarn in the following way And the coordinates b of each node on the weft yarn: ; ;

[0022] In the formula, represents the coordinates of the i-th node on the warp around the corresponding interweaving point, represents the jth node coordinates on the weft yarn around the corresponding interweaving point, and j represents the node number;

[0023] By setting the coordinates of each node on the warp And the coordinates b of each node on the weft yarn are combined to form several groups of finite element mesh areas.

[0024] Preferably, the finite element analysis module includes a multi-directional stress analysis unit and a macro state unit;

[0025] The multi-directional stress analysis unit is used to start the simulation program according to different working conditions set in the model building module, and respectively perform longitudinal tension application conditions, transverse tension application conditions and rotational force application conditions on the webbing, and analyze the deformation of each finite element mesh area in the webbing under different working conditions, so as to obtain the width value Kdz of each finite element mesh area under the longitudinal tension application condition, the thickness value Hdz of each finite element mesh area under the transverse tension application condition, and the width value Kzz of each finite element mesh area under the rotational force application condition, which is specifically obtained by the following formula: ; ; ;

[0026] In the formula, represents the initial width of the finite element mesh region, represents Poisson's ratio, represents the initial length value of the finite element mesh region, represents the initial thickness value of the finite element mesh area, represents the lateral tensile force, represents the elastic modulus, represents the original cross-sectional area of ​​the finite element mesh region, represents the effective radius of the cross section, It indicates the angle when the yarn is bent. Indicates the length change value It represents the tension applied to the yarn;

[0027] According to the width value Kdz of each finite element mesh area under the longitudinal tension application condition, the thickness value Hdz of each finite element mesh area under the transverse tension application condition, and the width value Kzz of each finite element mesh area under the rotational force application condition, the longitudinal stress Zy, transverse stress Hy and shear stress Qy are obtained respectively, which are specifically obtained by the following formulas: ; ; ;

[0028] In the formula, represents the tension applied to the yarn, represents the lateral tensile force, is the torque acting on the cross section, represents the distance from the point on the cross section to the neutral axis, represents the polar moment of inertia of the cross section;

[0029] Construct the comprehensive stress Zyz of each finite element mesh area. The specific acquisition method is as follows: ;

[0030] Preferably, the macro state unit is used to analyze the continuity of the stress gradient distribution according to the comprehensive stress Zyz of each finite element mesh area to construct a continuity coefficient Lxxs, and based on the continuity coefficient Lxxs, identify whether the stress field between each finite element mesh area in the webbing is smooth, and the continuity coefficient Lxxs is obtained by the following formula: ;

[0031] In the formula, represents the comprehensive stress in the current finite element mesh area, represents the combined stress in the adjacent finite element mesh region, Represents the distance between two finite element mesh regions.

[0032] Preferably, based on the value of the continuity coefficient Lxxs, the continuity coefficient Lxxs is compared with a preset continuity threshold to identify whether the stress field between the finite element mesh regions in the webbing is smooth, and the specific contents are as follows:

[0033] If the continuity coefficient Lxxs exceeds the preset continuity threshold, it is identified that the stress field between the finite element mesh areas in the webbing is not smooth, and the instability risk instruction will be triggered;

[0034] If the continuity coefficient Lxxs does not exceed the preset continuity threshold, it is recognized that there is no unevenness in the stress field between the finite element mesh areas in the webbing, and the instability risk instruction will not be triggered temporarily.

[0035] Preferably, the structural assessment module includes a vibration load condition analysis unit and a comprehensive analysis unit;

[0036] The vibration load condition analysis unit is used to receive the instability risk instruction, and according to the different working conditions set in the model building module, adjust different vibration loads, monitor the relative slip between the yarns in the multi-layer structure of the webbing, so as to obtain relevant webbing state data, and the relevant webbing state data includes the interlayer thickness And the interlayer slip displacement difference under different vibration loads f , based on the relevant webbing state data, the interlayer shear stress Jyz under different vibration load conditions is constructed, which is obtained by the following formula: ;

[0037] In the formula, represents the shear modulus of the ribbon, Indicates the interlayer thickness, It represents the interlayer slip displacement difference under different vibration loads f.

[0038] Preferably, the comprehensive analysis unit is used to construct an effect evaluation model of the webbing in the simulation process by using deep learning technology, by inputting the interlayer shear stress Jyz and the continuity coefficient Lxxs under different vibration load conditions into the effect evaluation model, and after linear normalization processing, fitting and outputting the knot stability score Jwpf, the knot stability score Jwpf is obtained by the following formula: ;

[0039] In the formula, and are weight values, is the correction constant, where 0< <1,0< <1, and The specific value is set by the user according to the situation.

[0040] Preferably, the feedback module is used to pre-set an evaluation threshold U, and compare and analyze the evaluation threshold U with the structural stability score Jwpf to estimate the stability performance of the internal structure of the webbing, and the specific contents are as follows:

[0041] If the structural stability score Jwpf exceeds the assessment threshold U, the stability performance of the current internal structure of the webbing is judged to be in an abnormal state. At this time, the judgment result is fed back to the background center in real time, and the density of the yarn interlacing points is increased to 1.1 times the original. The damping effect is introduced in the vibration analysis to reduce the impact of vibration on interlayer slip by increasing the internal friction of the webbing;

[0042] If the structural stability score Jwpf does not exceed the evaluation threshold U, it is judged that the stability performance of the current internal structure of the webbing is not in an abnormal state. At this time, the judgment result is fed back to the background center in real time, and the density of the yarn interlacing points is reduced to 1.1 times the original. The structural stability score Jwpf is analyzed again to see whether it does not exceed the evaluation threshold U, so as to continuously analyze the changing trend of the structural stability score Jwpf.

[0043] The present invention provides a 3D simulation system for a webbing based on finite element analysis, which has the following beneficial effects:

[0044] (1) The system automatically identifies the interweaving points in the webbing geometric model through the discrete partitioning module, and generates high-quality finite element mesh areas based on the interweaving points, so that the simulation analysis can fully capture the mechanical properties of the webbing. Combined with the finite element analysis module, by applying different types of forces (such as longitudinal tension, lateral tension, rotational force, etc.), it can accurately calculate the multi-directional stress on the mesh area, and use the equivalent stress theory to construct the comprehensive stress. Based on the comprehensive stress, the system can identify whether the stress field is smooth, helping to quickly find stress concentration areas and potential instability risks. The structural evaluation module monitors the relative slip between yarns in the multi-layer structure of the webbing by adjusting different vibration load conditions, calculates the interlayer shear stress, and combines the effect evaluation model to fit and output the structural stability score. This design can dynamically reflect the stability performance of the webbing under complex working conditions and provide a scientific basis for optimizing the design parameters of the webbing (such as interweaving point density, yarn material, etc.). The feedback module uses the structural stability score to evaluate the stability performance of the internal structure of the webbing in real time to determine whether the design needs to be optimized. When the structural stability score exceeds the threshold, the system can feedback the instability risk information and generate an optimization plan (such as increasing the density of interlacing points or adjusting vibration damping). This mechanism ensures the safety and reliability of the webbing design, while improving R&D efficiency and product performance. In short, the system further realizes the closed loop of the entire process from modeling to optimization, and provides efficient analysis tools for complex structures and working conditions of webbings, which not only improves simulation accuracy and design efficiency, but also provides strong technical support for the safety performance and reliability of webbings.

[0045] (2) By introducing the curvature calculation formula, the curvature of the warp and weft yarns around each interlacing point is obtained respectively, and the local bending strength and trajectory change rate of the yarn are quantified. In the formula, the local curvature of the curve is reflected by the cross product modulus of the yarn velocity vector and the acceleration vector; the tangent change rate of the curve at the current point is reflected by the modulus of the velocity vector. This method can not only accurately describe the geometric characteristics of the yarn in the complex interlacing structure, but also provide a scientific basis for the optimization of the discrete spacing. The accurate calculation of the yarn geometric characteristics (such as curvature) enables the discrete partitioning module to adaptively adjust the discrete spacing according to the local bending strength and change rate of the curve, thereby generating a high-quality finite element mesh area around the interlacing point. The fine mesh partitioning can more realistically reflect the mechanical behavior of the yarn and the stress distribution of the interlacing point, significantly improving the accuracy and reliability of the simulation calculation. Through automated interlacing point recognition and yarn curvature analysis, complex ribbon geometric models can be quickly processed without manual intervention, providing an efficient solution for multi-batch and multi-specification ribbon simulation analysis, significantly improving the application efficiency of the system. In summary, the present invention enhances the characterization capability of geometric characteristics and physical behaviors in the finite element model of the webbing through precise yarn geometric distribution analysis, and provides strong technical support for improving simulation accuracy and subsequent performance analysis.

[0046] (3) By combining the curvature information of the warp and weft yarns and the interaction pressure between the yarns, the discrete spacing in the warp and weft directions is dynamically adjusted, so that the discrete spacing automatically shrinks in the area where the yarn is more curved, and appropriately expands in the straight area, which improves the adaptability and flexibility of the grid division to the actual ribbon structure. By generating each node in the warp and weft directions according to the discrete spacing and the curve tangent vector, and combining the warp and weft node coordinates, a high-precision finite element grid area is formed. The node coordinate calculation method follows the yarn geometry, so that the generated grid can carefully reflect the local characteristics around the interweaving point, providing a fine calculation basis for subsequent mechanical simulation and performance evaluation. Material properties (such as yarn tension, yarn bending angle) and control parameters are introduced into the discrete spacing formula. By dynamically adjusting the discrete spacing, it can effectively adapt to ribbon structures with different material properties and weaving methods. Under different loads and working conditions, the system can flexibly adjust the grid density to ensure the precise matching of grid division and mechanical behavior.

[0047] (4) The system can reflect the stress state of the webbing in different force directions in many aspects through the multi-directional stress analysis method, which enhances the reliability of the evaluation results. By constructing comprehensive stress and conducting a holistic analysis of the stress field in each finite element mesh area, it helps to accurately locate the stress concentration area, identify possible failure risk points, and provide data support for the structural design optimization of the webbing.

[0048] (5) By calculating the gradient distribution of the comprehensive stress and constructing the continuity coefficient, the stress field changes between adjacent finite element mesh areas can be accurately reflected. When the stress gradient is discontinuous, the abnormal area can be quickly detected and marked, triggering the instability risk instruction, ensuring that the potential failure area of ​​the webbing structure is identified in time. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] Figure 1 This is a block diagram of a 3D simulation system for a webbing based on finite element analysis according to the present invention. DETAILED DESCRIPTION

[0050] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0051] Example 1

[0052] See also Figure 1, the present invention provides a 3D simulation system for a webbing based on finite element analysis, including a model building module, a discrete division module, a finite element analysis module, a structure evaluation module and a feedback module;

[0053] The model building module is used to generate a webbing geometric model using a three-dimensional modeling software in advance, and import the webbing geometric model into a finite element simulation software, and input different working conditions and related material properties of the webbing into the finite element simulation software;

[0054] The discrete division module is used to identify each interlacing point from the webbing geometric model in the finite element simulation software, and perform a finite element mesh division operation based on each interlacing point to generate a plurality of groups of finite element mesh regions;

[0055] The finite element analysis module is used to start the simulation program according to different working conditions set in the model construction module, apply force to the webbing, analyze the deformation of each finite element mesh area in the webbing under different working conditions, and obtain multi-directional stress based on the deformation, and use the equivalent stress theory to construct the comprehensive stress Zyz of each finite element mesh area. Based on the comprehensive stress Zyz of each finite element mesh area, identify whether the stress field between each finite element mesh area in the webbing is smooth;

[0056] The structural evaluation module is used to adjust different vibration loads according to different working conditions set in the model building module, monitor the relative slip between yarns in the multi-layer structure of the webbing, so as to obtain relevant webbing state data, and construct the interlayer shear stress Jyz under different vibration load conditions based on the relevant webbing state data, and combine the trained effect evaluation model, and after linear normalization, fit the output structure stability score Jwpf;

[0057] The feedback module is used to estimate the stability performance of the internal structure of the webbing based on the structural stability score Jwpf value, and feed it back to the background center to formulate an optimization plan.

[0058] During the operation of this system, the system uses the model building module to generate the geometric model of the webbing in advance using the 3D modeling software and imports it into the finite element simulation software. The relevant material properties of the webbing (such as elastic modulus, Poisson's ratio, etc.) and different working conditions (such as longitudinal tension, lateral tension, vibration load, etc.) are input into the simulation software, so that the system can efficiently build the webbing model and carry out accurate simulation, providing a reliable basis for subsequent performance analysis. The discrete partitioning module automatically identifies the interweaving points in the webbing geometric model and divides the finite element mesh based on this, thereby generating several groups of finite element mesh areas. The finite element analysis module combines the set working conditions, applies external force to the webbing, analyzes its deformation under different conditions, and calculates multi-directional stress values. Using the equivalent stress theory, the comprehensive stress of each grid area is constructed, and by analyzing the stress gradient distribution between grids, it is identified whether the internal stress field of the webbing is smooth. This precise analysis provides strong support for the identification of stress concentration areas and the judgment of potential instability problems. The structural evaluation module further expands the functions of the system. It can monitor the relative slip between yarns in the multi-layer structure of the webbing by regulating the vibration load, and calculate the interlayer shear stress based on these state data. Combined with the effect evaluation model constructed by deep learning, the system can output the structural stability score after linear normalization. This score provides a quantitative basis for the stability of the webbing structure under different working conditions, helping users to grasp the mechanical properties of the webbing in real time. The feedback module evaluates the stability of the internal structure of the webbing according to the structural stability score. If the structural stability score is abnormal, the system can promptly feedback to the background and start the optimization process, such as adjusting the density of the yarn interlacing points or introducing the damping effect, so as to reduce the interlayer slip and improve the impact of vibration on stability. This closed-loop feedback mechanism maximizes the intelligence of the system and ensures the reliability and safety of the webbing design. Through modular design, this system is compatible with different finite element simulation software (such as ANSYS and ABAQUS), and can be widely used in the webbing performance optimization needs in many fields such as clothing, safety, industrial lifting and aerospace. Compared with traditional methods, the system significantly shortens the design iteration cycle and reduces experimental testing costs. At the same time, it improves the performance and reliability of webbing products through high-precision simulation, bringing significant economic benefits and social value to industrial applications.

[0059] Example 2

[0060] Please refer to Figure 1 ,Specifically: the model building module includes a modeling unit and a data input unit;

[0061] The modeling unit is used to obtain the relevant geometric feature data of the ribbon through various detection instruments in advance, and simulate the geometric model of the ribbon in proportion to the ribbon according to the relevant geometric feature data and the three-dimensional modeling software, wherein the relevant geometric feature data includes the width, thickness, yarn arrangement, yarn color, yarn thickness, weaving pattern, spatial distribution, interlayer thickness of the ribbon , the coordinates of the interlacing points of the warp and weft yarns , the paths and interweaving density of warp and weft yarns; 3D modeling software includes SolidWorks and AutoCAD;

[0062] The data input unit is used to save the geometric model of the webbing in the STEP file format, ensuring that the geometric model of the webbing can be recognized and correctly parsed by the corresponding simulation software to support further operations of the finite element analysis, realize seamless data flow between different software platforms, and avoid model information loss or simulation failure due to file incompatibility; and import the geometric model of the webbing into the finite element simulation software, and input different working conditions and related material properties of the webbing into the finite element simulation software, wherein the related material properties of the webbing include elastic modulus , Poisson's ratio , yarn tension and shear modulus ; Different working conditions include load types and load sizes, among which load types include longitudinal tension, lateral tension, torsional load and vibration load; finite element simulation software includes ANSYS and ABAQUS.

[0063] In this embodiment, the modeling unit uses a detection instrument to obtain the geometric feature data of the webbing, such as width, thickness, yarn arrangement, yarn color and thickness, weaving pattern, spatial distribution, coordinates of warp and weft interlacing points, path and interlacing density, and combines the three-dimensional modeling software to generate a three-dimensional geometric model that is proportional to the actual webbing. This method can accurately reproduce the real geometric characteristics of the webbing, especially the microscopic characteristics of the complex interlacing structure and yarn distribution, and provide high-precision input data for subsequent finite element analysis to ensure the accuracy of the simulation results. The data input unit saves the generated geometric model in a standardized STEP file format to ensure that the geometric model can be correctly identified and parsed by mainstream finite element simulation software such as ANSYS and ABAQUS. Seamless data transfer is achieved between different software platforms, avoiding model information loss or simulation failure due to incompatible file formats, improving the reliability and efficiency of data transmission, and providing support for cross-platform collaboration and result verification. Input the relevant material properties of the webbing (such as elastic modulus, Poisson's ratio, yarn tension and shear modulus) and various working conditions (such as longitudinal tension, lateral tension, torsional load and vibration load) into the simulation software, so that the model can adapt to different mechanical analysis scenarios. By simulating and analyzing the performance of the webbing under multiple working conditions, the behavior of the webbing under static and dynamic loading conditions can be comprehensively evaluated, providing a scientific basis for design optimization and failure analysis. From geometric feature data acquisition to 3D modeling, and then to the import of simulation software and the input of material properties, this module realizes a highly automated process, effectively reducing manual intervention and possible errors, and improving the efficiency of simulation modeling. The systematic workflow can quickly respond to different analysis needs and is suitable for the simulation needs of multiple batches and multiple specifications of webbing. By simulating the mechanical properties of the webbing under different working conditions, the strength, stability and instability risk of the webbing can be evaluated. In summary, the present invention realizes a seamless connection of the entire process from geometric model generation to simulation data input. Through precise data acquisition and modeling, efficient data transmission and simulation analysis of multiple working conditions, the accuracy, efficiency and applicability of webbing simulation are significantly improved, providing strong technical support for the design and evaluation of high-performance webbing.

[0064] Example 3

[0065] Please refer to Figure 1 ,Specifically: the discrete division module includes a yarn geometric distribution analysis unit and a grid division unit;

[0066] The yarn geometric distribution analysis unit is used to identify each interlacing point from the webbing geometric model, and determine the warp space coordinates around each interlacing point in the webbing in combination with relevant geometric feature data. and the spatial position of the weft yarn , based on the spatial coordinates of the warp yarns around each interlacing point in the ribbon and the spatial position of the weft yarn , respectively obtain the curvature of the warp and weft yarns around each interlacing point and mark them as the warp curvature and weft curvature , which is obtained by the following formula: ;

[0067] In the formula, represents the velocity vector of the warp yarn (first derivative), represents the acceleration vector of the warp yarn (second-order derivative), and t represents the position point on the yarn; is the cross product modulus of the velocity vector and the acceleration vector of the warp yarn, representing the local bending strength of the curve, is the modulus of the velocity vector of the warp, indicating the rate of change (magnitude) of the curve in the tangent direction of the current point; ;

[0068] In the formula, represents the velocity vector of the weft yarn (first derivative), represents the acceleration vector (second-order derivative) of the weft yarn, and t represents the position point on the yarn. Corresponding to the physical scene, t is a one-dimensional parameter of the yarn trajectory, indicating a position of the yarn on the spatial path. is the cross product modulus of the velocity vector and the acceleration vector of the weft yarn, representing the local bending strength of the curve, is the modulus of the velocity vector of the weft yarn, indicating the rate of change (magnitude) of the curve in the tangent direction of the current point;

[0069] Velocity Vector Represents the velocity vector, that is, the first-order derivative of the curve with respect to the position point t on the yarn, describing the tangent direction and rate of change of the curve in space. The acquisition method is for the spatial coordinates of the warp yarns around each interlacing point in the ribbon. and the spatial position of the weft yarn , can be calculated using the finite difference method: ;

[0070] The analytical method uses calculus to find the first-order derivative, where Indicates the position point on the yarn The yarn space coordinates, Indicates the position difference on the yarn;

[0071] Acceleration vector represents the acceleration vector, that is, the second derivative of the curve at the position point t on the yarn, and represents the rate of change of the velocity vector: ;

[0072] Acquisition method: For discrete coordinate points, use the difference method to find the second-order derivative:

[0073] Among them, yarn includes warp yarn and weft yarn;

[0074] In this embodiment, the yarn geometric distribution analysis unit can accurately identify the interweaving points of the warp and weft yarns from the webbing geometric model, and extract the coordinates and paths of the warp and weft yarns in space in combination with relevant geometric feature data. The curvature of the warp and weft yarns is calculated separately using the curvature formula, thereby accurately reflecting the local bending strength and trajectory change rate of the yarn. Through the calculation of the curvature, the deformation characteristics of the yarn near each interweaving point can be quantified, providing a scientific basis for finite element mesh division. The present invention uses the cross-product modulus of the velocity vector and the acceleration vector and the tangent direction change rate to describe the yarn curvature, which can fully capture the local bending characteristics and deformation behavior of the yarn at the interweaving point. This multi-dimensional data description method improves the accuracy of discrete division, makes the mesh division more refined, and can more realistically simulate the complexity of the webbing interweaving structure, providing a guarantee for the accuracy of subsequent finite element analysis. By utilizing the calculation results of the curvature of the warp and weft yarns, the mesh division unit can adaptively adjust the discrete density near each interweaving point. For areas with larger curvature, a finer grid is divided; for areas with smaller curvature, a coarser grid is divided. This curvature-based discrete partitioning strategy can optimize the number of grids, reduce computational complexity, and improve simulation efficiency while ensuring simulation accuracy. The accurately calculated curvatures of warp and weft yarns provide reliable basic data for performance analysis under dynamic conditions. The curvature results can not only help identify stress concentration areas on the structure, but also assist in evaluating deformation modes under different loading conditions, thereby providing a scientific basis for webbing design and optimization. In summary, the discrete partitioning module not only improves the modeling capabilities of complex interwoven structures by accurately extracting yarn geometric properties and adaptively partitioning the grid, but also greatly improves the accuracy and efficiency of simulation analysis, providing strong technical support for the research and development of high-performance webbing.

[0075] Example 4

[0076] Please refer to Figure 1 Specifically: the grid division unit is used to calculate the warp curvature obtained by the yarn geometric distribution analysis unit according to the relevant material properties. and weft curvature , get the discrete spacing in each direction around each interlacing point, including the discrete spacing on the warp yarn and discrete spacing on the weft yarns The specific acquisition method is as follows: ; ;

[0077] In the formula, Indicates the pressure exerted by adjacent warp yarns; Indicates the pressure exerted by adjacent weft yarns; It indicates the average pressure exerted by adjacent weft yarns at each interlacing point; It indicates the average pressure exerted by adjacent warp yarns at each interlacing point; represents the warp curvature; represents the weft yarn curvature; Represents a constant, used to control the base value of discrete density; represents the adjustment coefficient, which is used to amplify the effect of pressure on the discrete spacing; represents the initial discrete spacing on the warp; represents the initial discrete spacing on the weft yarn; where the pressure P exerted by adjacent yarns is obtained by the following formula: in, Indicates the yarn tension, It indicates the angle when the yarn is bent;

[0078] According to the coordinates of the intersection points , along the warp direction, according to the tangent vector of the warp and discrete spacing on the warp , generate each node on the warp yarn, along the weft yarn direction, according to the curve tangent vector of the weft yarn and discrete spacing on the weft yarns , generate each node on the weft yarn, and obtain the coordinates of each node on the warp yarn in the following way And the coordinates b of each node on the weft yarn: ; ;

[0079] In the formula, represents the coordinates of the i-th node on the warp around the corresponding interweaving point, represents the jth node coordinates on the weft yarn around the corresponding interweaving point, and j represents the node number, which is an integer and increases from 0 until the entire length between two interlacing points is covered;

[0080] The curve tangent vector of the warp yarn and the tangent vector of the weft yarn Refers to the unit tangent vector in the yarn direction, describing the extension direction of the warp yarn, which can be obtained by the following formula: ;

[0081] By setting the coordinates of each node on the warp and the coordinates b of each node on the weft yarn to form several groups of finite element mesh regions;

[0082] In this embodiment, according to the warp and weft curvatures, as well as the pressure applied by adjacent yarns, the discrete spacing in the warp and weft directions is adaptively adjusted, the influence of yarn curvature and pressure is combined through a formula, the discrete spacing is dynamically adjusted, and the high curvature area and high pressure area have a finer grid division. This optimization enables the grid division to more accurately characterize the complex mechanical behavior around the interweaving point, and provides a high-resolution finite element basis for simulation calculation. The pressure applied by the adjacent yarns is calculated by using the yarn tension and the angle when the yarn is bent, and the force transmission characteristics of the yarn interweaving point are reflected by the pressure. Combining the pressure parameter with the discrete spacing ensures that the grid division can fully consider the mechanical differences of the ribbon in different areas, and provide accurate data support for subsequent simulation analysis. Through the discrete spacing and curve tangent vectors in the warp and weft directions, accurate node coordinates are generated, and the node coordinates of the warp and weft yarns are combined to form a finite element grid area around the interweaving point. This method reflects higher accuracy and adaptability in the grid division of the interweaving point, and can more realistically simulate the geometric characteristics and stress state inside the ribbon. Dynamic adjustment of discrete spacing and high-precision node coordinate generation make the finite element mesh area division more reasonable, avoid unnecessary over-dense or over-sparse meshes, and significantly improve the accuracy and efficiency of simulation calculations. Especially in the complex interweaving point area, this method can better capture the stress distribution and deformation details, laying a solid foundation for subsequent simulation analysis and structural evaluation. In summary, the present invention successfully achieves efficient mapping between the geometric characteristics and mechanical properties of the webbing through optimized mesh division unit design, providing a strong guarantee for the accuracy and reliability of the webbing 3D simulation system.

[0083] Example 5

[0084] Please refer to Figure 1 ,Specifically: the finite element analysis module includes a multi-directional stress analysis unit and a macro state unit;

[0085] The multi-directional stress analysis unit is used to start the simulation program according to different working conditions set in the model building module, and respectively perform longitudinal tension application conditions, transverse tension application conditions and rotational force application conditions on the webbing, and analyze the deformation of each finite element mesh area in the webbing under different working conditions, so as to obtain the width value Kdz of each finite element mesh area under the longitudinal tension application condition, the thickness value Hdz of each finite element mesh area under the transverse tension application condition, and the width value Kzz of each finite element mesh area under the rotational force application condition, which is specifically obtained by the following formula: ; ; ;

[0086] In the formula, represents the initial width of the finite element mesh region, represents Poisson's ratio, represents the initial length value of the finite element mesh region, represents the initial thickness value of the finite element mesh area, represents the lateral tensile force, represents the elastic modulus, represents the original cross-sectional area of ​​the finite element mesh region, represents the effective radius of the cross section, It indicates the angle when the yarn is bent. Indicates the length change value, which is obtained by the following formula: Indicates the tension applied to the yarn; effective radius of the cross section It is a parameter used to describe the equivalent geometric characteristics of the cross section of the ribbon. In the actual ribbon, the cross section may be an irregular shape (such as an ellipse or a rectangle), and the effective radius of the cross section is equivalent to a circular cross section to simplify the calculation. The effective radius of the cross section is the local width of the webbing and thickness To calculate, it is obtained by the following formula: in, is the constant in the area formula of a circle, used to convert a rectangle into a circle; It represents shear strain, which reflects the degree of deformation of the material under shear stress;

[0087] It should be noted that the reason why the longitudinal tension affects the width is that the dominant longitudinal stress and the Poisson effect cause the transverse dimension (width) to shrink.

[0088] The reason why transverse tension affects thickness is that the dominant transverse stress, Poisson effect, causes compression in the thickness direction of the yarn.

[0089] The reason why the rotational force affects the width is that the dominant shear stress, slip and shear deformation mainly affect the weft yarn spacing, resulting in width changes.

[0090] The reason why the effect is directional is that the weaving structure: the warp and weft directions respond differently under stress conditions. Mechanical transmission path: the dominant stress of different working conditions is transmitted through the direction of the yarn, making the deformation of width or thickness the main manifestation. Deformation sensitivity: the mechanical sensitivity of width and thickness in different directions is different, resulting in deformation mainly concentrated in a certain direction.

[0091] According to the width value Kdz of each finite element mesh area under the longitudinal tension application condition, the thickness value Hdz of each finite element mesh area under the transverse tension application condition, and the width value Kzz of each finite element mesh area under the rotational force application condition, the longitudinal stress Zy, transverse stress Hy and shear stress Qy are obtained respectively, which are specifically obtained by the following formulas: ; ; ;

[0092] In the formula, represents the tension applied to the yarn, represents the lateral tensile force, is the torque acting on the cross section, represents the distance from the point on the cross section to the neutral axis, represents the polar moment of inertia of the cross section;

[0093] The above torque acting on the cross section It is a physical quantity that describes the rotation effect of a force on an axis. The torque acting on a cross section is the cumulative effect of shear force on the cross section caused by external loads or internal torques. Based on finite element analysis, the load conditions are input and the torque distribution of yarn contact friction is calculated;

[0094] Distance from a point on the cross section to the neutral axis It refers to the area of ​​the cross section that is not deformed under the torsional load. The distance from a point to the neutral axis refers to the shortest distance between any point on the cross section and the neutral axis. This parameter affects the shear stress borne at the point. The farther the distance, the greater the shear stress borne. Acquisition method: 3D scanner: obtain the geometric shape and size of the ribbon cross section; Image processing and CAD modeling: extract the cross-sectional shape through image recognition software, calculate the neutral axis position using modeling software, and mark the distance from the point to the neutral axis;

[0095] Polar moment of inertia of cross section It is a geometric characteristic parameter that describes the cross-sectional shape's resistance to torsional deformation. It depends on the cross-sectional shape and its size. Acquisition method: 3D scanner and image recognition software: obtain the cross-sectional shape and use the software to calculate the polar moment of inertia.

[0096] Among them, shear strain The relationship between the shear stress Qy is: in, It represents the shear modulus of the webbing (unit: Pa), which indicates the ability of the webbing to resist shear deformation;

[0097] Construct the comprehensive stress Zyz of each finite element mesh area. The specific acquisition method is as follows: ;

[0098] In this embodiment, by simulating three typical working conditions of longitudinal tension, transverse tension and rotational force, the width change and thickness change of the webbing under different load conditions are accurately analyzed, which provides comprehensive data support for understanding the structural deformation and mechanical behavior of the webbing in a complex use environment, making the evaluation result more reliable. Using the equivalent stress theory, combined with the width, thickness and cross-sectional changes under different working conditions, the longitudinal stress, transverse stress and shear stress are calculated respectively. The coupling relationship between shear stress and shear strain further refines the analysis of shear performance, thereby accurately reflecting the comprehensive stress distribution of the webbing under complex load conditions, providing an important basis for subsequent optimization design. By equating the irregular cross section to a circular cross section and introducing the effective radius formula, the calculation difficulty of the complex structure is significantly simplified, while ensuring the calculation accuracy. The introduction of the cross-section effective radius parameter solves the interference of the complex cross-sectional shape on the calculation of mechanical properties and improves the ability of the simulation model to restore the real situation. By constructing the comprehensive stress Zyz of the finite element mesh area, the multi-directional stress influence is systematically integrated. The analysis of the comprehensive stress field makes the identification of stress concentration areas more accurate, which helps to evaluate the stability of the webbing structure and prevent potential instability risks. In summary, the present invention systematically evaluates the deformation and stress distribution of the webbing under complex stress conditions through the multi-directional stress analysis and comprehensive stress construction of the simulation module, greatly improves the accuracy of mechanical property analysis, and provides strong technical support for the design optimization, material selection and performance evaluation of the webbing.

[0099] Example 6

[0100] Please refer to Figure 1 Specifically: the macro state unit is used to analyze the continuity of the stress gradient distribution according to the comprehensive stress Zyz of each finite element mesh area to construct the continuity coefficient Lxxs. Based on the continuity coefficient Lxxs, it is identified whether the stress field between each finite element mesh area in the webbing is smooth. The continuity coefficient Lxxs is obtained by the following formula: ;

[0101] In the formula, represents the comprehensive stress in the current finite element mesh area, represents the combined stress in the adjacent finite element mesh region, Represents the distance between two finite element mesh regions.

[0102] Based on the value of the continuity coefficient Lxxs, the continuity coefficient Lxxs is compared with a preset continuity threshold to identify whether the stress field between the finite element mesh regions in the webbing is smooth. The specific contents are as follows:

[0103] If the continuity coefficient Lxxs exceeds the preset continuity threshold, it is identified that the stress field between the finite element mesh areas in the webbing is not smooth, and the instability risk instruction will be triggered;

[0104] If the continuity coefficient Lxxs does not exceed the preset continuity threshold, it is recognized that there is no unevenness in the stress field between the finite element mesh areas in the webbing, and the instability risk instruction will not be triggered temporarily.

[0105] In this embodiment, by defining the continuity coefficient Lxxs, combined with the comprehensive stress difference and distance between the current grid area and the adjacent area, the continuity degree of the stress field is quantified. This coefficient effectively reveals the changing law of the stress field inside the ribbon, and provides a scientific basis for identifying potential non-smooth areas. By comparing the continuity coefficient Lxxs with the set continuity threshold, it is accurately judged whether the stress field has non-smoothness. When the continuity coefficient Lxxs exceeds the threshold, the system automatically triggers the instability risk instruction, prompting potential structural instability problems. This real-time monitoring mechanism can actively discover design defects during the simulation process, avoid mechanical failure of the ribbon in actual use, and improve the safety of the product. In the case of a smooth stress field, by analyzing the changing trend of the continuity coefficient Lxxs, it can help designers further optimize the stress distribution of each area in the ribbon, reduce local stress concentration, and enhance the overall structural stability of the ribbon. This provides data support for the design optimization of the ribbon, which is particularly important under complex working conditions. The present invention uses the automated analysis function of the simulation software and takes the continuity coefficient Lxxs as the core parameter to realize the intelligent evaluation of the stress field, avoiding the subjective bias and accuracy problems that may be caused by manual judgment, and improving the analysis efficiency and result reliability of the system. After identifying the non-smooth area of ​​the stress field, the instability risk instruction is triggered in time to provide a reference basis for subsequent design improvement and optimization. This mechanism can greatly improve the credibility of the simulation results and ensure the stability and reliability of the ribbon design in practical applications.

[0106] Example 7

[0107] Please refer to Figure 1 ,Specifically: the structural evaluation module includes a vibration load condition analysis unit and a comprehensive analysis unit;

[0108] The vibration load condition analysis unit is used to receive the instability risk instruction, and according to the different working conditions set in the model building module, adjust different vibration loads, monitor the relative slip between the yarns in the multi-layer structure of the webbing, so as to obtain relevant webbing state data, and the relevant webbing state data includes the interlayer thickness And the interlayer slip displacement difference under different vibration loads f , based on the relevant webbing state data, the interlayer shear stress Jyz under different vibration load conditions is constructed, which is obtained by the following formula: ;

[0109] In the formula, represents the shear modulus of the ribbon, Indicates the interlayer thickness, It represents the interlayer slip displacement difference under different vibration loads f.

[0110] The comprehensive analysis unit is used to construct an effect evaluation model of the webbing in the simulation process by using deep learning technology. The interlayer shear stress Jyz and the continuity coefficient Lxxs under different vibration load conditions are input into the effect evaluation model, and after linear normalization, the knot stability score Jwpf is fitted and output. The knot stability score Jwpf is obtained by the following formula: ;

[0111] In the formula, and are weight values, is the correction constant, where 0< <1,0< <1, and The specific value is set by the user according to the situation.

[0112] The feedback module is used to pre-set the evaluation threshold U, and compare and analyze the evaluation threshold U with the structural stability score Jwpf to estimate the stability performance of the internal structure of the webbing. The specific contents are as follows:

[0113] If the structural stability score Jwpf exceeds the evaluation threshold U, the stability performance of the current internal structure of the ribbon is judged to be in an abnormal state. At this time, the judgment result is fed back to the background center in real time, and the density of the yarn interlacing points is increased to 1.1 times the original (increasing the intersection points of the warp and weft yarns through simulation), and the damping effect is introduced in the vibration analysis to reduce the influence of vibration on interlayer slip by increasing the internal friction of the ribbon;

[0114] If the structural stability score Jwpf does not exceed the evaluation threshold U, it is judged that the stability performance of the current internal structure of the webbing is not in an abnormal state. At this time, the judgment result is fed back to the background center in real time, and the density of the yarn interlacing points is reduced to 1.1 times the original (reducing the intersection points of the warp and weft yarns through simulation), and the structural stability score Jwpf is analyzed again to see whether it does not exceed the evaluation threshold U, so as to continuously analyze the changing trend of the structural stability score Jwpf.

[0115] In this embodiment, the vibration load condition analysis unit accurately obtains the interlayer thickness and slip displacement difference by monitoring the relative slip of the yarn in the multilayer structure, and further calculates the interlayer shear stress Jyz. This real-time data collection and analysis method can fully reflect the structural changes under different vibration load conditions and help identify potential instability risks. The comprehensive analysis unit constructs an effect evaluation model through deep learning technology, comprehensively analyzes the interlayer shear stress Jyz and the continuity coefficient Lxxs, and outputs the structural stability score Jwpf through linear normalization processing. This multi-dimensional evaluation method effectively improves the accuracy of the webbing performance analysis and provides a scientific basis for subsequent optimization design and performance improvement. The feedback module dynamically adjusts the interweaving point density of the webbing according to the comparison result of the stability score Jwpf and the preset threshold U. When an abnormal state is found, the vibration resistance of the structure is enhanced by increasing the density of the warp and weft interweaving points and introducing the damping effect; under normal conditions, the material utilization efficiency is optimized by appropriately reducing the interweaving point density. This adaptive optimization mechanism balances the contradiction between structural stability and material cost. Through real-time monitoring of instability risks, data feedback and dynamic adjustment, the present invention forms a complete closed-loop control process to ensure the stable performance of the webbing under complex working conditions. At the same time, the continuous analysis of the score change trend further verifies the effectiveness of the optimization measures. In summary, the present invention significantly improves the stability and reliability of the webbing under vibration loads through the multifunctional design of the structural evaluation module, combined with deep learning and adaptive optimization technology, and provides important technical support for the design, manufacture and practical application of the webbing.

[0116] Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present invention, and that the scope of the present invention is defined by the appended claims and their equivalents.

Claims

1. A 3D simulation system for webbing based on finite element analysis, characterized in that: It includes model building module, discrete partitioning module, finite element analysis module, structural evaluation module and feedback module; The model building module is used to generate a webbing geometric model using a three-dimensional modeling software in advance, and import the webbing geometric model into a finite element simulation software, and input different working conditions and related material properties of the webbing into the finite element simulation software; The discrete division module is used to identify each interlacing point from the webbing geometric model in the finite element simulation software, and perform a finite element mesh division operation based on each interlacing point to generate a plurality of groups of finite element mesh regions; The finite element analysis module is used to start the simulation program according to different working conditions set in the model construction module, apply force to the webbing, analyze the deformation of each finite element mesh area in the webbing under different working conditions, and obtain multi-directional stress based on the deformation, and use the equivalent stress theory to construct the comprehensive stress Zyz of each finite element mesh area. Based on the comprehensive stress Zyz of each finite element mesh area, identify whether the stress field between each finite element mesh area in the webbing is smooth; The structural evaluation module is used to adjust different vibration loads according to different working conditions set in the model building module, monitor the relative slip between yarns in the multi-layer structure of the webbing, so as to obtain relevant webbing state data, and construct the interlayer shear stress Jyz under different vibration load conditions based on the relevant webbing state data, and combine the trained effect evaluation model, and after linear normalization, fit the output structure stability score Jwpf; The feedback module is used to estimate the stability performance of the internal structure of the webbing based on the structural stability score Jwpf value, and feed it back to the background center to formulate an optimization plan.

2. A 3D simulation system for webbing based on finite element analysis according to claim 1, characterized in that: The model building module includes a modeling unit and a data input unit; The modeling unit is used to obtain the relevant geometric feature data of the ribbon through various detection instruments in advance, and simulate the geometric model of the ribbon in proportion to the ribbon according to the relevant geometric feature data and the three-dimensional modeling software, wherein the relevant geometric feature data includes the width, thickness, yarn arrangement, yarn color, yarn thickness, weaving pattern, spatial distribution, interlayer thickness of the ribbon , the coordinates of the interlacing points of the warp and weft yarns , the paths and interweaving density of warp and weft yarns; 3D modeling software includes SolidWorks and AutoCAD; The data input unit is used to save the geometric model of the webbing in the STEP file format, and import the geometric model of the webbing into the finite element simulation software, and input different working conditions and related material properties of the webbing into the finite element simulation software, wherein the related material properties of the webbing include elastic modulus , Poisson's ratio , yarn tension and shear modulus ; Different working conditions include load type and load size, among which load types include longitudinal tension, lateral tension, torsional load and vibration load.

3. The 3D simulation system for webbing based on finite element analysis according to claim 2, characterized in that: The discrete division module includes a yarn geometric distribution analysis unit and a grid division unit; The yarn geometric distribution analysis unit is used to identify each interlacing point from the webbing geometric model, and determine the warp space coordinates around each interlacing point in the webbing in combination with relevant geometric feature data. and the spatial position of the weft yarn , based on the spatial coordinates of the warp yarns around each interlacing point in the ribbon and the spatial position of the weft yarn , respectively obtain the curvature of the warp and weft yarns around each interlacing point and mark them as the warp curvature and weft curvature , which is obtained by the following formula: ； In the formula, represents the velocity vector of the warp yarn, represents the acceleration vector of the warp yarn, and t represents the position point on the yarn; is the cross product modulus of the velocity vector and the acceleration vector of the warp, is the magnitude of the velocity vector of the warp yarn; ； In the formula, represents the velocity vector of the weft yarn, represents the acceleration vector of the weft yarn, and t represents the position point on the yarn; is the cross product modulus of the velocity vector and the acceleration vector of the weft yarn; is the magnitude of the velocity vector of the weft yarn.

4. The 3D simulation system for webbing based on finite element analysis according to claim 3, characterized in that: The mesh division unit is used to calculate the warp curvature obtained by the yarn geometric distribution analysis unit according to the relevant material properties. and weft curvature , get the discrete spacing in each direction around each interlacing point, including the discrete spacing on the warp yarn and discrete spacing on the weft yarns The specific acquisition method is as follows: ； ； In the formula, Indicates the pressure exerted by adjacent warp yarns; Indicates the pressure exerted by adjacent weft yarns; It indicates the average pressure exerted by adjacent weft yarns at each interlacing point; It represents the average pressure exerted by adjacent warp yarns at each interlacing point; represents the warp curvature; represents the weft yarn curvature; represents a constant; represents the adjustment coefficient; According to the coordinates of the intersection points , along the warp direction, according to the tangent vector of the warp and discrete spacing on the warp , generate each node on the warp yarn, along the weft yarn direction, according to the curve tangent vector of the weft yarn and discrete spacing on the weft yarns , generate each node on the weft yarn, and obtain the coordinates of each node on the warp yarn in the following way And the coordinates b of each node on the weft yarn: ； ； In the formula, represents the coordinates of the i-th node on the warp around the corresponding interweaving point, represents the jth node coordinates on the weft yarn around the corresponding interweaving point, and j represents the node number; By setting the coordinates of each node on the warp And the coordinates b of each node on the weft yarn are combined to form several groups of finite element mesh areas.

5. The 3D simulation system for webbing based on finite element analysis according to claim 4, characterized in that: The finite element analysis module includes a multi-directional stress analysis unit and a macro state unit; The multi-directional stress analysis unit is used to start the simulation program according to different working conditions set in the model building module, and respectively perform longitudinal tension application conditions, transverse tension application conditions and rotational force application conditions on the webbing, and analyze the deformation of each finite element mesh area in the webbing under different working conditions, so as to obtain the width value Kdz of each finite element mesh area under the longitudinal tension application condition, the thickness value Hdz of each finite element mesh area under the transverse tension application condition, and the width value Kzz of each finite element mesh area under the rotational force application condition, which is specifically obtained by the following formula: ； ； ； In the formula, represents the initial width of the finite element mesh region, represents Poisson's ratio, represents the initial length value of the finite element mesh region, represents the initial thickness value of the finite element mesh area, represents the lateral tensile force, represents the elastic modulus, represents the original cross-sectional area of ​​the finite element mesh region, represents the effective radius of the cross section, It indicates the angle when the yarn is bent. Indicates the length change value It represents the tension applied to the yarn; According to the width value Kdz of each finite element mesh area under the longitudinal tension application condition, the thickness value Hdz of each finite element mesh area under the transverse tension application condition, and the width value Kzz of each finite element mesh area under the rotational force application condition, the longitudinal stress Zy, transverse stress Hy and shear stress Qy are obtained respectively, which are specifically obtained by the following formulas: ； ； ； In the formula, represents the tension applied to the yarn, represents the lateral tensile force, is the torque acting on the cross section, represents the distance from the point on the cross section to the neutral axis, represents the polar moment of inertia of the cross section; Construct the comprehensive stress Zyz of each finite element mesh area. The specific acquisition method is as follows: 。 6. The 3D simulation system for webbing based on finite element analysis according to claim 5, characterized in that: The macro state unit is used to analyze the continuity of the stress gradient distribution according to the comprehensive stress Zyz of each finite element mesh area to construct the continuity coefficient Lxxs. Based on the continuity coefficient Lxxs, it is identified whether the stress field between each finite element mesh area in the webbing is smooth. The continuity coefficient Lxxs is obtained by the following formula: ； In the formula, represents the comprehensive stress in the current finite element mesh area, represents the combined stress in the adjacent finite element mesh region, Represents the distance between two finite element mesh regions.

7. The 3D simulation system for webbing based on finite element analysis according to claim 6, characterized in that: Based on the value of the continuity coefficient Lxxs, the continuity coefficient Lxxs is compared with a preset continuity threshold to identify whether the stress field between the finite element mesh regions in the webbing is smooth. The specific contents are as follows: If the continuity coefficient Lxxs exceeds the preset continuity threshold, it is identified that the stress field between the finite element mesh areas in the webbing is not smooth, and the instability risk instruction will be triggered; If the continuity coefficient Lxxs does not exceed the preset continuity threshold, it is recognized that there is no unevenness in the stress field between the finite element mesh areas in the webbing, and the instability risk instruction will not be triggered temporarily.

8. The 3D simulation system for webbing based on finite element analysis according to claim 7, characterized in that: The structural assessment module includes a vibration load condition analysis unit and a comprehensive analysis unit; The vibration load condition analysis unit is used to receive the instability risk instruction, and according to the different working conditions set in the model building module, adjust different vibration loads, monitor the relative slip between the yarns in the multi-layer structure of the webbing, so as to obtain relevant webbing state data, and the relevant webbing state data includes the interlayer thickness And the interlayer slip displacement difference under different vibration loads f , based on the relevant webbing state data, the interlayer shear stress Jyz under different vibration load conditions is constructed, which is obtained by the following formula: ； In the formula, represents the shear modulus of the ribbon, Indicates the interlayer thickness, It represents the interlayer slip displacement difference under different vibration loads f.

9. The 3D simulation system for webbing based on finite element analysis according to claim 8, characterized in that: The comprehensive analysis unit is used to construct an effect evaluation model of the webbing in the simulation process by using deep learning technology. The interlayer shear stress Jyz and the continuity coefficient Lxxs under different vibration load conditions are input into the effect evaluation model, and after linear normalization, the knot stability score Jwpf is fitted and output. The knot stability score Jwpf is obtained by the following formula: ； In the formula, and are weight values, is the correction constant, where 0< <1,0< <1, and The specific value is set by the user according to the situation.

10. The 3D simulation system for webbing based on finite element analysis according to claim 1, characterized in that: The feedback module is used to pre-set the evaluation threshold U, and compare and analyze the evaluation threshold U with the structural stability score Jwpf to estimate the stability performance of the internal structure of the webbing. The specific contents are as follows: If the structural stability score Jwpf exceeds the assessment threshold U, the stability performance of the current internal structure of the webbing is judged to be in an abnormal state. At this time, the judgment result is fed back to the background center in real time, and the density of the yarn interlacing points is increased to 1.1 times the original. The damping effect is introduced in the vibration analysis to reduce the impact of vibration on interlayer slip by increasing the internal friction of the webbing; If the structural stability score Jwpf does not exceed the evaluation threshold U, it is judged that the stability performance of the current internal structure of the webbing is not in an abnormal state. At this time, the judgment result is fed back to the background center in real time, and the density of the yarn interlacing points is reduced to 1.1 times the original. The structural stability score Jwpf is analyzed again to see whether it does not exceed the evaluation threshold U, so as to continuously analyze the changing trend of the structural stability score Jwpf.

Images