A simulation method and system for elliptical fiber winding line design

Abstract

The present invention discloses a simulation method and system for the design of elliptical fiber winding lines, relating to the technical fields of composite fiber winding and computer-aided design simulation. The method comprises: establishing a core mold model based on core mold parameters; obtaining a transition section slip coefficient and a barrel winding angle based on a preset transition section slip coefficient range and the length of the winding transition section, or based on a preset barrel winding angle and a preset transition section slip coefficient range; calculating winding equal division points based on the core mold parameters, a preset yarn width, and the barrel winding angle; obtaining the center angle of the entire round trip of the equal division point based on the winding angle, the barrel semi-axis, and the transition section slip coefficient; determining the center angle of the line based on a preset deviation angle and the center angle of the entire round trip of the equal division point; optimizing the transition section slip coefficient based on the two center angles; obtaining a doffing point based on the optimized slip coefficient; and performing simulation based on the doffing point and the core mold model. The present invention can improve the efficiency of fiber winding line design and production.

Description

Technical Field

[0001] The present invention relates to the technical field of composite material filament winding and computer-aided design simulation, and in particular to a simulation method and system for linear design of elliptical body filament winding. Background Art

[0002] Filament winding technology is an advanced manufacturing technology that uses winding machine equipment to control the relative movement between the nozzle and the core mold, and according to the set winding parameters, the resin-impregnated fiber yarn is evenly and stably wound on the surface of the core mold according to a certain pattern to achieve the strength performance requirements, and then the overall curing and molding is carried out.

[0003] The development of fiber winding technology is inseparable from the support of software. The use of winding simulation design software can greatly improve product R&D efficiency, save design and development costs, simplify complex winding theories, and facilitate research and analysis by technical personnel. However, the research and development of winding simulation software only covers axisymmetric rotational body winding products with circular cross-sections (such as pressure vessels, bottles, cones, spheres, straight tubes). The line design for non-axisymmetric body winding products cannot be effectively and evenly distributed. The line design method is relatively limited, and it is not easy to find a suitable tangent line for production. It is necessary to repeatedly iterate and change the design parameters to calculate the line. In the simulation, the yarn width cannot be effectively displayed, and it is impossible to determine whether the line is evenly distributed. Actual production tests are often required, resulting in low line design and production efficiency. Summary of the Invention

[0004] The purpose of the present invention is to provide a simulation method and system for elliptical body fiber winding line design, which can improve the efficiency of fiber winding line design and production.

[0005] To achieve the above object, the present invention provides the following solutions:

[0006] A simulation method for designing a fiber winding line for an elliptical body, comprising:

[0007] Obtaining core mold parameters, and establishing a core mold model using OpenGL drawing functions according to the core mold parameters; the core mold parameters include the left transition section length, the barrel section length, the right transition section length, the barrel section major semi-axis and the barrel section minor semi-axis of the core mold;

[0008] Obtaining the barrel winding angle, the left transition section slip coefficient, and the right transition section slip coefficient according to a preset transition section slip coefficient range and a preset winding transition section length or according to a preset barrel winding angle and a preset transition section slip coefficient range;

[0009] Calculating winding equal points according to the core mold parameters, the preset yarn width and the barrel winding angle;

[0010] Obtaining the center angle of each winding equally divided point throughout the entire round trip according to the coordinates of each winding equally divided point, the barrel winding angle, the semi-major axis of the barrel segment, the semi-minor axis of the barrel segment, the left transition section slip coefficient, and the right transition section slip coefficient;

[0011] Determine the center angle of the line type based on the preset deviation angle and the center angle of each winding equally divided point throughout the entire round trip;

[0012] Optimizing the left transition section slip coefficient and the right transition section slip coefficient according to the entire round-trip center angle of each winding equally divided point and the center angle of the line to obtain an optimized left transition section slip coefficient and an optimized right transition section slip coefficient;

[0013] Obtaining the doffing points on the entire back-and-forth fiber path at each equally divided winding point according to the optimized left transition section slip coefficient and the optimized right transition section slip coefficient;

[0014] Linear simulation is performed based on the doffing points on the entire back-and-forth fiber path at each winding equally divided point and the core mold model.

[0015] Optionally, the calculating of winding equal division points according to the core mold parameters, the preset yarn width and the barrel winding angle specifically includes:

[0016] Calculating the number of cycles according to the core mold parameters, the preset yarn width and the barrel winding angle;

[0017] The coordinates of the winding equal division points are determined according to the number of cycles.

[0018] Optionally, determining the center angle of the line shape according to the preset deviation angle and the center angle of each winding equally divided point in the entire round trip specifically includes:

[0019] Determine the uniform distribution line type parameters based on the preset deviation angle and the center angle of each winding equal-division point throughout the entire round trip;

[0020] The center angle of the line type is determined according to the uniformly distributed line type parameters.

[0021] Optionally, the calculating the number of cycles according to the core mold parameters, the preset yarn width and the barrel winding angle specifically includes:

[0022] The perimeter of the barrel section is obtained according to the semi-major axis and the semi-minor axis of the barrel section;

[0023] The number of cycles is calculated based on the cross-sectional perimeter of the barrel, the preset yarn width and the barrel winding angle.

[0024] Optionally, the linear simulation is performed based on the doffing points on the entire back-and-forth fiber path at each winding equally divided point and the core mold model, specifically including:

[0025] If yarn width simulation is not performed, linear simulation is performed based on the doffing points on the entire back-and-forth fiber path at each winding equally divided point and the core mold model;

[0026] If yarn width simulation is performed, the unit normal vector of the doffing point on the fiber path of each winding equal-division point is calculated based on the center angle of the entire round trip of each winding equal-division point and the parallel circle radius of the core mold section;

[0027] The unit tangent vector of the doffing point on the fiber path of each winding equal-division point is calculated according to the parallel circle radius of the core mold section, the center angle of the entire round trip of each winding equal-division point, the major semi-axis and the minor semi-axis of the barrel section;

[0028] Calculate the unit offset vector of the doffing point on the entire round-trip fiber path at each equally divided winding point according to the unit normal vector of the doffing point on the entire round-trip fiber path at each equally divided winding point and the unit tangent vector of the doffing point on the entire round-trip fiber path at each equally divided winding point;

[0029] Obtaining an offset doffing point corresponding to each winding equally divided point according to the preset yarn width and the unit offset vector of the doffing point on the entire back-and-forth fiber path at each winding equally divided point;

[0030] Linear simulation is performed based on the offset doffing points corresponding to each winding equally divided point, the doffing points on the entire back-and-forth fiber path of each winding equally divided point, and the core mold model.

[0031] A simulation system for elliptical body fiber winding line design, comprising:

[0032] An acquisition module is used to obtain core mold parameters and establish a core mold model using OpenGL drawing functions according to the core mold parameters; the core mold parameters include the left transition section length, the barrel section length, the right transition section length, the major semi-axis of the barrel section, and the minor semi-axis of the barrel section;

[0033] A slip coefficient and barrel winding angle calculation module is used to obtain the barrel winding angle, the left transition section slip coefficient and the right transition section slip coefficient according to a preset transition section slip coefficient range and a preset winding transition section length or according to a preset barrel winding angle and a preset transition section slip coefficient range;

[0034] A winding equal division point calculation module is used to calculate the winding equal division points according to the core mold parameters, the preset yarn width and the barrel winding angle;

[0035] A winding equally divided point center angle calculation module is used to obtain the center angle of each winding equally divided point throughout the entire round trip based on the coordinates of each winding equally divided point, the barrel winding angle, the semi-major axis of the barrel segment, the semi-minor axis of the barrel segment, the left transition section slip coefficient, and the right transition section slip coefficient;

[0036] A linear center angle calculation module is used to determine the linear center angle based on a preset deviation angle and the center angle of each winding equally divided point throughout the entire round trip;

[0037] an optimization module, configured to optimize the left transition section slip coefficient and the right transition section slip coefficient according to the entire round-trip center angle of each winding equally divided point and the center angle of the line, to obtain an optimized left transition section slip coefficient and an optimized right transition section slip coefficient;

[0038] a doffing point calculation module, configured to obtain the doffing points on the entire round-trip fiber path at each winding equally divided point according to the optimized left transition section slip coefficient and the optimized right transition section slip coefficient;

[0039] The simulation module is used to perform linear simulation based on the yarn drop points on the entire back-and-forth fiber path at each winding equally divided point and the core mold model.

[0040] Optionally, the winding equal division point calculation module specifically includes:

[0041] A cycle number calculation unit, used to calculate the cycle number according to the core mold parameters, the preset yarn width and the barrel winding angle;

[0042] The winding equally divided point calculation unit is used to determine the coordinates of the winding equally divided point according to the number of cycles.

[0043] Optionally, the linear center angle calculation module specifically includes:

[0044] A uniformly distributed linear parameter calculation unit is used to determine the uniformly distributed linear parameters according to a preset deviation angle and the center angle of each winding equally divided point throughout the entire round trip;

[0045] The line center angle calculation unit is used to determine the line center angle according to the uniformly distributed line parameters.

[0046] Optionally, the cycle number calculation unit specifically includes:

[0047] The barrel section perimeter calculation subunit is used to obtain the barrel section perimeter based on the barrel section's semi-major axis and the barrel section's semi-minor axis;

[0048] The cycle number calculation subunit is used to calculate the cycle number according to the perimeter of the barrel cross section, the preset yarn width and the barrel winding angle.

[0049] Optionally, the simulation module specifically includes:

[0050] A first simulation unit is configured to perform line simulation based on the doffing points on the entire back-and-forth fiber path at each winding equally divided point and the core model if yarn width simulation is not performed;

[0051] The unit normal vector calculation unit of the doffing point is used to calculate the unit normal vector of the doffing point on the entire fiber path of each winding equally divided point according to the center angle of the entire round trip of each winding equally divided point and the parallel circle radius of the core mold cross section when performing yarn width simulation;

[0052] The unit tangent vector calculation unit of the doffing point is used to calculate the unit tangent vector of the doffing point on the entire back-and-forth fiber path of each winding equally divided point based on the parallel circle radius of the core mold cross section, the center angle of the entire back-and-forth of each winding equally divided point, the major semi-axis and the minor semi-axis of the barrel segment;

[0053] a unit offset vector calculation unit for a doffing point, for calculating a unit offset vector of a doffing point on the entire round-trip fiber path at each equally divided winding point based on a unit normal vector of the doffing point on the entire round-trip fiber path at each equally divided winding point and a unit tangent vector of the doffing point on the entire round-trip fiber path at each equally divided winding point;

[0054] A calculation unit, configured to obtain an offset doffing point corresponding to each equally divided winding point according to the preset yarn width and a unit offset vector of the doffing point on the entire back-and-forth fiber path at each equally divided winding point;

[0055] The second simulation unit is used to perform linear simulation according to the offset doffing points corresponding to each winding equally divided point, the doffing points on the entire back-and-forth fiber path of each winding equally divided point, and the core mold model.

[0056] According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects: the present invention establishes a core mold model according to the core mold parameters; obtains the transition section slip coefficient and the barrel winding angle according to the preset transition section slip coefficient range and the winding transition section length or according to the preset barrel winding angle and the preset transition section slip coefficient range; calculates the winding equal division point according to the core mold parameters, the preset yarn width and the barrel winding angle; obtains the center angle of the equal division point for the entire round trip according to the winding angle, the barrel semi-axis and the transition section slip coefficient; determines the center angle of the line type according to the preset deviation angle and the center angle of the equal division point for the entire round trip; optimizes the transition section slip coefficient according to the two center angles; obtains the yarn dropping point according to the optimized slip coefficient; and performs simulation according to the yarn dropping point and the core mold model, thereby improving the efficiency of fiber winding line type design and production. BRIEF DESCRIPTION OF THE DRAWINGS

[0057] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0058] Figure 1A flow chart of the simulation method for elliptical body fiber winding line design provided by the present invention;

[0059] Figure 2 Schematic diagram of the core mold model;

[0060] Figure 3 This is a simulation diagram of the fiber-free yarn wide line using method one;

[0061] Figure 4 This is a simulation diagram showing the yarn width line shape using method 1;

[0062] Figure 5 This is the simulation diagram of the fiber without yarn width line using method 2;

[0063] Figure 6 This is a simulation diagram showing the yarn width line shape using method 2. DETAILED DESCRIPTION

[0064] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. 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 making creative efforts are within the scope of protection of the present invention.

[0065] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0066] like Figure 1 As shown, an embodiment of the present invention provides a simulation method for designing an elliptical body fiber winding line, comprising:

[0067] Obtain core mold parameters, and use OpenGL drawing functions to build a core mold model based on the core mold parameters. The core mold model is as follows: Figure 2 As shown; the core mold parameters include the left transition section length, the barrel section length, the right transition section length, the long semi-axis of the barrel section and the short semi-axis of the barrel section.

[0068] The barrel winding angle, the left transition section slip coefficient and the right transition section slip coefficient are obtained according to the preset transition section slip coefficient range and the preset winding transition section length or according to the preset barrel winding angle and the preset transition section slip coefficient range, wherein the left and right transition sections adopt non-geodesic winding and the barrel adopts geodesic winding.

[0069] The winding equal division points are calculated according to the core mold parameters, the preset yarn width and the barrel winding angle.

[0070] According to the coordinates of each of the winding equally divided points, the barrel winding angle, the long semi-axis of the barrel section, the short semi-axis of the barrel section, the slip coefficient of the left transition section and the slip coefficient of the right transition section, the center angle of each of the winding equally divided points in the entire round trip is obtained. The center angle of the entire round trip of the winding equally divided point refers to the total center angle of the fiber completing a round trip winding starting from the winding equally divided point. The coordinates of the winding equally divided point are the coordinates of the first yarn dropping point among all the yarn dropping points that make up the entire round trip of the winding equally divided point.

[0071] The center angle of the line is determined based on the preset deviation angle and the center angle of the entire round trip of each winding equally divided point.

[0072] The left transition section slip coefficient and the right transition section slip coefficient are optimized according to the entire round center angle of each winding equally divided point and the linear center angle to obtain the optimized left transition section slip coefficient and the optimized right transition section slip coefficient.

[0073] The doffing points on the entire round-trip fiber path at each equally divided winding point are obtained according to the optimized left transition section slip coefficient and the optimized right transition section slip coefficient.

[0074] Linear simulation is performed based on the doffing points on the entire back-and-forth fiber path at each winding equally divided point and the core mold model.

[0075] In practical applications, the barrel winding angle, the left transition section slip coefficient, and the right transition section slip coefficient are obtained according to a preset transition section slip coefficient range and a preset winding transition section length or according to a preset barrel winding angle and a preset transition section slip coefficient range, specifically including:

[0076] Method 1: Preset the transition section slip coefficient range and the winding transition section length, and perform linear calculation starting from the end of the transition section to determine whether there is an intersection between the winding angles at the left and right equators in this case. If not, optimize the transition section slip coefficient (left transition section slip coefficient and right transition section slip coefficient) or core mold parameters, expand the transition section slip coefficient range or modify the left and right transition section lengths until there is an intersection. Then, the barrel winding angle, left transition section slip coefficient and right transition section slip coefficient under the current situation are obtained. The judgment process is as follows: Given the transition section slip coefficient range and the guaranteed winding transition section length, the barrel winding angle range can be obtained by solving the equations of the following center angle and winding angle. According to the two barrel winding angle ranges obtained from the left and right transition sections, determine whether there is a right intersection. Whether a winding angle within the allowable error range can be obtained at the left and right equators depends on the friction coefficient range and the transition section length. The winding angle error is 10 -4 .

[0077] Method 2: Preset the barrel winding angle and transition section slip coefficient range, perform linear calculation, and judge whether the left transition section length and the right transition section length meet the requirements to make the winding angle transition from the barrel winding angle to 90°. If not, increase the barrel winding angle or the transition section slip coefficient range or increase the barrel winding angle and transition section slip coefficient range at the same time until the requirements are met. The barrel winding angle, left transition section slip coefficient and right transition section slip coefficient in the current situation are obtained. The judgment process is based on the following center angle and winding angle joint calculation equation. The left and right transition section lengths are known, which is the equation Z value. The slip coefficient and the barrel winding angle are known. The transition length required for the winding angle to transition to 89.99° can be calculated. If the length is less than the transition section length in the core mold parameters, it meets the requirements. If it is greater than the transition section length in the core mold parameters, it does not meet the requirements.

[0078] In practical applications, the calculation of the winding equal division points according to the core mold parameters, the preset yarn width and the barrel winding angle specifically includes:

[0079] The number of cycles is calculated according to the core mold parameters, the preset yarn width and the barrel winding angle.

[0080] The coordinates of the winding equal division points are determined according to the number of cycles. After the number of cycles is determined, the cross section of the barrel is divided into equal parts by the number of cycles to obtain the coordinates of the winding equal division points.

[0081] In practical applications, the center rotation angle of each winding equally divided point is obtained according to the coordinates of each winding equally divided point, the barrel winding angle, the semi-major axis of the barrel segment, the semi-minor axis of the barrel segment, the left transition section slip coefficient and the right transition section slip coefficient, specifically:

[0082] Step 1: Solve the center angle calculation equation and the winding angle together to obtain the center angle of the doffing point. The solution equation is:

[0083]

[0084]

[0085] Among them: α is the winding angle of the doffing point; θ is the central rotation angle of the doffing point; Z is the axial coordinate of the doffing point. The axial coordinate refers to the coordinate position of the core mold with the central axis as the axis. This value depends on the position of the origin coordinate system defined by the user. Generally, the coordinate system is taken at the axis center where the left transition section and the barrel section intersect; μ is the slip coefficient including the slip coefficient of the left transition section and the slip coefficient of the right transition section; a is the major semi-axis of the ellipse; b is the minor semi-axis of the ellipse.

[0086] The fiber path of the entire back and forth winding equal-division point is composed of several yarn dropping points. The specific number of yarn dropping points is determined manually. The number of yarn dropping points is related to the linear accuracy and will not cause deviation in the center angle of the entire back and forth winding equal-division point. The division of the yarn dropping points is divided by the axial coordinate z value. The number of yarn dropping points included in the entire back and forth of each winding equal-division point is the same, and the axial coordinates of the yarn dropping points are also the same. The only change is the center angle value calculated based on each winding equal-division point. Substituting the left transition section slip coefficient into the above formula, several yarn dropping points of the left transition section can be obtained. The center angle of the last yarn dropping point is the center angle of the entire left transition section. Similarly, the center angle of the last yarn dropping point of the barrel section is the center angle of the entire barrel section, and the center angle of the last yarn dropping point of the right transition section is the center angle of the entire right transition section. Although only the left slip coefficient and the right slip coefficient can be substituted into the formula each time, the solution of the differential equation requires given initial conditions, which are determined according to the cross-sectional position of each winding equal-division point.

[0087] Step 2. The center angle of the entire winding equal-division point is the center angle of the last dropping yarn point among all the dropping yarn points that make up the entire winding equal-division point in the left transition section when the line goes out, the center angle of the last dropping yarn point among all the dropping yarn points that make up the entire winding equal-division point in the barrel section, and the center angle of the last dropping yarn point among all the dropping yarn points that make up the entire winding equal-division point in the right transition section, and the center angle of the last dropping yarn point among all the dropping yarn points that make up the entire winding equal-division point in the right transition section when the line returns, the center angle of the last dropping yarn point among all the dropping yarn points that make up the entire winding equal-division point in the barrel section, and the center angle of the last dropping yarn point among all the dropping yarn points that make up the entire winding equal-division point in the left transition section.

[0088] In practical applications, the barrel winding angle refers to the winding angle of the yarn point in the barrel section. Since the slip coefficient of the barrel section is 0, the winding angles of all the doffing points in the entire barrel section are equal. Since the left and right transition sections have slip coefficients, the winding angles of the left and right transition sections are constantly changing. Assuming that the winding angle of the left end doffing point of the left transition section is 89.99°, and the right end doffing point coincides with the left end doffing point of the barrel section, the winding angle of the right end doffing point of the left transition section is the barrel winding angle determined according to method one or method two. Similarly, the winding angle of the left end doffing point of the right transition section is the barrel winding angle determined according to method one or method two.

[0089] Based on the left and right transition slip coefficients and the barrel winding angle determined by Method 1, the center angle of the yarn point in the left transition section is calculated as follows: Starting from the left end of the left transition section and moving toward the right end, the length of the left transition section is assumed to be 100 mm. The number of doffing points is manually determined. For example, if 11 doffing points are set, the axial z coordinates of these 11 doffing points are 0 mm, 10 mm, 20 mm, and so on. The initial conditions are set, such as z = 0 mm at the left end of the left transition section, the initial doffing angle (i.e., the winding angle at the winding equalizer) is α = 89.99°, and the center angle of the doffing point (i.e., the center angle at the winding equalizer) is 0°. Substituting this into the differential equation, the computer automatically iterates the calculation. The iterative increment is also manually set, and the size of the increment is related to the calculation accuracy. The calculation is completed when z = 100 mm. The data for the 11 doffing points, including the z coordinate value, winding angle, and center angle of the doffing point, are then output based on the set number of doffing points.

[0090] In practical applications, for the left transition slip coefficient, right transition slip coefficient and barrel winding angle determined by method 2, the center turning angle of the yarn point in the left transition section is calculated as follows: starting from the right end of the left transition section and calculating towards the left end of the left transition section, the number of yarn dropping points is determined manually, for example, 11 yarn dropping points are set. The initial conditions are determined, for example, the initial point z = 0mm at the right end of the left transition section, the initial yarn dropping point winding angle is the barrel winding angle α = 70° obtained by method 2, and the center turning angle of the yarn dropping point = 0°. Substituting into the differential equation, the computer can automatically iterate the calculation. The iterative increment step is also set manually. The increment step size is related to the calculation accuracy. When the winding angle α is calculated to be greater than or equal to 89.99°, the calculation is completed, and the axial coordinate z value of the left transition section is recorded. If the z value at this time is 100mm, the coordinates of the 11 yarn dropping points are 0mm, 10mm, 20mm...100mm. Finally, the data of 11 doffing points are output according to the set number of doffing points, including the z coordinate value, winding angle and center angle of the doffing point.

[0091] The calculation process of the center angle of the yarn drop point of the barrel section is the same as above and will not be repeated here.

[0092] In practical applications, the method of determining the center angle of the line according to the preset deviation angle and the center angle of each winding equally divided point in the entire round trip specifically includes:

[0093] The uniform distribution line type parameters are determined based on the preset deviation angle and the center angle of each winding equally divided point throughout the entire round trip.

[0094] The center angle of the line type is determined according to the uniformly distributed line type parameters.

[0095] In practical applications, the linear simulation is performed based on the doffing points on the entire back-and-forth fiber path at each winding equally divided point and the core model, specifically including:

[0096] If yarn width simulation is not performed, line type simulation is performed by drawing line types with the help of OpenGL according to the yarn drop points on the entire back-and-forth fiber path at each winding equally divided point and the core mold model.

[0097] If yarn width simulation is performed, the unit normal vector of the yarn drop point on the entire round-trip fiber path of each winding equally divided point is calculated based on the center turning angle of the entire round-trip of each winding equally divided point and the parallel circle radius of the core mold section.

[0098] The unit tangent vector of the doffing point on the entire round-trip fiber path of each winding equally divided point is calculated according to the parallel circle radius of the core mold section, the center turning angle of the entire round-trip of each winding equally divided point, the major semi-axis and the minor semi-axis of the barrel segment.

[0099] The unit offset vector of the yarn dropping point on the entire round-trip fiber path at each winding equally divided point is calculated based on the unit normal vector of the yarn dropping point on the entire round-trip fiber path at each winding equally divided point and the unit tangent vector of the yarn dropping point on the entire round-trip fiber path at each winding equally divided point.

[0100] The offset doffing point corresponding to each winding equally divided point is obtained according to the preset yarn width and the unit offset vector of the doffing point on the entire back-and-forth fiber path of each winding equally divided point.

[0101] Linear simulation is performed based on the offset doffing points corresponding to each winding equally divided point, the doffing points on the entire back-and-forth fiber path of each winding equally divided point, and the core mold model.

[0102] In practical applications, the calculation of the number of cycles according to the core mold parameters, the preset yarn width and the barrel winding angle specifically includes:

[0103] The circumference of the barrel section is obtained based on the long semi-axis and the short semi-axis of the barrel section.

[0104] The number of cycles is calculated based on the cross-sectional perimeter of the barrel, the preset yarn width and the barrel winding angle.

[0105] In practical applications, the uniform distribution line parameters are determined based on the preset deviation angle and the center angle of each winding dividing point throughout the entire round trip, including:

[0106] According to the continued fraction theory, the uniform distribution of the line type is optimized and the uniform distribution line type parameters are recorded. The calculation steps are as follows:

[0107] Step A1: Increase the number of achievable line types and appropriately expand the range of the number of loops to [N min ,N max ], given the deviation angle Δθ.

[0108] Step A2: The maximum center angle θ of the center angles of the winding points is calculated. maxThe maximum and minimum number of mandrel turns per winding are calculated with the deviation angle Δθ. The calculation formula is as follows:

[0109] r cmin =(θ max -Δθ) / 360°

[0110] r cmax =(θ max +Δθ) / 360°

[0111] Step A3: By [r cmin ,r cmax ] and [N min ,N max ]Determine the value range of span d[d min ,d max ], the calculation formula is:

[0112] r c =W+d / N

[0113] Where: W is the integer number of turns of the mandrel during one round trip.

[0114] Step A4: Screen based on the condition that d and N are mutually prime, eliminate the combinations that do not meet the requirements, and apply the continued fraction Euclidean algorithm to d / N to obtain a series of winding line parameters such as the number of tangent points, series, total number of loops, and the center angle of the cart's round trip.

[0115] Step A5: Determine uniform fill parameters based on the winding line parameters. Uniform fill parameters include, but are not limited to, winding line parameters. Uniform fill parameters include maximum center angle, winding angle, left transition section slip coefficient, right transition section slip coefficient, number of tangent points, number of loops, coverage, etc. These parameters generally refer to functions required for line design that can be obtained.

[0116] In practical applications, the center angle of the line type is determined according to the uniformly distributed line type parameters, specifically including:

[0117] The center angle of the uniformly filled line is obtained based on the parameters of the uniformly filled line. The transition slip coefficient is optimized so that the center angle of the winding bisection point is equal to the center angle of the selected line. The transition slip coefficient is iteratively optimized using the bisection method within the range [x-0.1, x+0.1]. When the error between the center angle of the winding bisection point corresponding to the optimized transition slip coefficient and the center angle of the uniformly filled line is within 0.0001, the two are considered equal. x refers to the input transition slip coefficient.

[0118] In practical applications, the doffing points on the entire back-and-forth fiber path at each winding equally divided point are obtained according to the optimized left transition section slip coefficient and the optimized right transition section slip coefficient, specifically including:

[0119] The winding angle, center angle and axial coordinate of the doffing point on the entire back-and-forth fiber path of all winding equally divided points are calculated based on the optimized transition section slip coefficient.

[0120] The coordinates of all doffing points are calculated based on the winding angles, center angles and axial coordinates of the doffing points on the entire back-and-forth fiber path of all winding equally divided points. All doffing points can be called in subsequent steps to obtain the winding trajectory line shape.

[0121] In practical applications, the number of cycles is calculated based on the cross-sectional perimeter of the barrel, the preset yarn width, and the barrel winding angle, specifically including:

[0122] The number of loops is calculated according to the formula: N = ceil(C*cosα / B), where N is the number of loops, C is the circumference of the barrel section, and B is the preset yarn width. ceil(.) returns the smallest integer greater than .

[0123] In practical applications, the unit normal vector of the doffing point on the fiber path of each winding equal point is calculated based on the center angle of the entire round trip of each winding equal point and the parallel circle radius of the core mold section. Specifically include:

[0124] According to the formula:

[0125]

[0126] calculate,

[0127] in: is the core surface equation, θ is the central angle of the doffing point; is the first-order derivative with respect to θ; is the first-order derivative with respect to z; r' is the first-order derivative of the parallel circle radius; r refers to the equivalent circle radius at a certain point on the core mold cross section.

[0128] In practical applications, the unit tangent vector of the doffing point on the fiber path of each winding equal-division point is calculated based on the parallel circle radius of the core mold section, the center angle of the entire round trip of each winding equal-division point, the major semi-axis and the minor semi-axis of the barrel segment. Specifically include:

[0129] According to the formula:

[0130]

[0131] in, calculate.

[0132] In practical applications, the unit offset vector of the doffing point on the entire round-trip fiber path at each equally divided winding point is calculated based on the unit normal vector of the doffing point on the entire round-trip fiber path at each equally divided winding point and the unit tangent vector of the doffing point on the entire round-trip fiber path at each equally divided winding point, specifically including:

[0133] Unit normal vector of the doffing point on the entire back-and-forth fiber path at each winding equal division point The unit tangent vector of the doffing point on the entire round-trip fiber path at each winding dividing point The cross product calculation obtains the unit offset vector of the yarn drop point on the entire back and forth fiber path at each winding equal division point

[0134] In practical applications, the offset doffing point corresponding to each winding equally divided point is obtained according to the preset yarn width and the unit offset vector of the doffing point on the entire back-and-forth fiber path of each winding equally divided point, specifically including:

[0135] One doffing point gets two equidistant offset doffing points in the positive direction of the offset vector, and two equidistant offset doffing points in the negative direction of the offset vector. The distance between the outermost offset doffing point in the positive direction and the outermost offset doffing point in the negative direction is the yarn width. The four offset doffing points plus the original doffing points make a total of five doffing points. Assuming that a winding trajectory has 20 doffing points, plus the corresponding offset doffing points, there are a total of 100 doffing points.

[0136] In practical applications, linear simulation is performed based on the offset doffing points corresponding to each winding equal-division point, the doffing points on the entire back-and-forth fiber path of each winding equal-division point, and the core model, specifically including:

[0137] The coordinate components of the yarn dropping points are read in a counterclockwise direction with three adjacent yarn dropping points as triangle vertices. Once all the coordinate components are read, a fiber line can be formed into a fiber surface sheet, thereby completing the fiber yarn width simulation.

[0138] In accordance with the above method, an embodiment of the present invention provides a simulation system for designing a filament winding line shape of an elliptical body, comprising:

[0139] The acquisition module is used to obtain the core mold parameters and establish a core mold model based on the core mold parameters using OpenGL drawing functions; the core mold parameters include the left transition section length, barrel section length, right transition section length, barrel section major semi-axis and barrel section minor semi-axis of the core mold.

[0140] The slip coefficient and barrel winding angle calculation module is used to obtain the barrel winding angle, the left transition section slip coefficient and the right transition section slip coefficient according to the preset transition section slip coefficient range and the preset winding transition section length or according to the preset barrel winding angle and the preset transition section slip coefficient range.

[0141] The winding equal division point calculation module is used to calculate the winding equal division point according to the core mold parameters, the preset yarn width and the barrel winding angle.

[0142] The winding equally divided point center angle calculation module is used to obtain the center angle of each winding equally divided point in the entire round trip based on the coordinates of each winding equally divided point, the barrel winding angle, the major semi-axis of the barrel segment, the minor semi-axis of the barrel segment, the left transition section slip coefficient and the right transition section slip coefficient.

[0143] The linear center angle calculation module is used to determine the linear center angle according to the preset deviation angle and the center angle of each winding equally divided point throughout the entire round trip.

[0144] The optimization module is used to optimize the left transition section slip coefficient and the right transition section slip coefficient according to the entire round-trip center angle of each winding equally divided point and the center angle of the line to obtain the optimized left transition section slip coefficient and the optimized right transition section slip coefficient.

[0145] The yarn drop point calculation module is used to obtain the yarn drop point on the entire round-trip fiber path of each winding equal division point according to the optimized left transition section slip coefficient and the optimized right transition section slip coefficient.

[0146] The simulation module is used to perform linear simulation based on the yarn drop points on the entire back-and-forth fiber path at each winding equally divided point and the core mold model.

[0147] In practical applications, the winding equal division point calculation module specifically includes:

[0148] The cycle number calculation unit is used to calculate the cycle number according to the core mold parameters, the preset yarn width and the barrel winding angle.

[0149] The winding equally divided point calculation unit is used to determine the coordinates of the winding equally divided point according to the number of cycles.

[0150] In practical applications, the linear center angle calculation module specifically includes:

[0151] The uniformly distributed linear parameter calculation unit is used to determine the uniformly distributed linear parameter according to the preset deviation angle and the center angle of each winding equally divided point.

[0152] The line center angle calculation unit is used to determine the line center angle according to the uniformly distributed line parameters.

[0153] In practical applications, the cycle number calculation unit specifically includes:

[0154] The barrel section perimeter calculation subunit is used to obtain the barrel section perimeter based on the long semi-axis and the short semi-axis of the barrel section.

[0155] The cycle number calculation subunit is used to calculate the cycle number according to the perimeter of the barrel cross section, the preset yarn width and the barrel winding angle.

[0156] In practical applications, the simulation module specifically includes:

[0157] The first simulation unit is used to perform line simulation according to the yarn dropping points on the entire back-and-forth fiber path at each winding equally divided point and the core mold model if yarn width simulation is not performed.

[0158] The unit normal vector calculation unit of the yarn drop point is used to calculate the unit normal vector of the yarn drop point on the entire round-trip fiber path of each winding equal-division point based on the center turning angle of the entire round-trip of each winding equal-division point and the parallel circle radius of the core mold section if yarn width simulation is performed.

[0159] The unit tangent vector calculation unit of the yarn dropping point is used to calculate the unit tangent vector of the yarn dropping point on the entire back-and-forth fiber path of each winding equal-division point based on the parallel circle radius of the core mold section, the center turning angle of the entire back-and-forth of each winding equal-division point, the major semi-axis and the minor semi-axis of the barrel segment.

[0160] The unit offset vector calculation unit of the yarn dropping point is used to calculate the unit offset vector of the yarn dropping point on the entire round-trip fiber path of each winding equally divided point based on the unit normal vector of the yarn dropping point on the entire round-trip fiber path of each winding equally divided point and the unit tangent vector of the yarn dropping point on the entire round-trip fiber path of each winding equally divided point.

[0161] The calculation unit is used to obtain the offset doffing point corresponding to each winding equally divided point according to the preset yarn width and the unit offset vector of the doffing point on the entire back and forth fiber path of each winding equally divided point.

[0162] The second simulation unit is used to perform linear simulation according to the offset doffing points corresponding to each winding equally divided point, the doffing points on the entire back-and-forth fiber path of each winding equally divided point, and the core mold model.

[0163] The present invention further provides a C++-based filament winding line design simulation system based on the above-mentioned elliptical body filament winding line design simulation method. The system includes the following modules:

[0164] 1. Core model parameter modeling module, used to receive user input parameters to establish a model.

[0165] 2. Material parameter input module, used to receive user input material parameters to obtain fiber yarn width and single layer thickness.

[0166] 3. Line type design module, used to output a list of evenly distributed line types and obtain the corresponding line type fiber paths.

[0167] 4. Linear simulation module, used for linear static simulation and dynamic simulation.

[0168] 5. G code generation module, which outputs G code that matches the winding machine tool for actual production.

[0169] The present invention provides a more specific embodiment to illustrate the above method:

[0170] The cross section of the core mold in this embodiment is an elliptical cross section, and the core mold with an elliptical cross section is non-axisymmetric. The length of the barrel section in this embodiment is L=600mm, the long semi-axis a=80mm, the short semi-axis b=40mm, the left and right heads are extended heads, the length of the left and right transition sections is L1,2=120mm, and the left and right transition winding lengths are determined to be 120mm using method 1, the slip coefficient μ=0.12, the yarn width B=10mm, and the barrel winding angle is calculated to be 50.5°, and the center angle is 2470.9°. The schematic diagram of the core mold is as follows: Figure 2 As shown, the linear yarn-free width simulation diagram is as follows Figure 3 As shown, the linear yarn width simulation diagram is as follows Figure 4 shown.

[0171] The present invention also provides a more specific embodiment to illustrate the above method:

[0172] The cross section of the core mold in this embodiment is an elliptical cross section, and the elliptical cross section core mold is non-axisymmetric. The length of the barrel section in this embodiment is L = 600mm, the long semi-axis a of the barrel section is 80mm, the short semi-axis b is 40mm, the left and right heads are extended heads, the length of the left and right transition sections is L1,2 = 180mm, the barrel winding angle is determined to be 55° using method 2, the slip coefficient μ = 0.12, the yarn width B = 10mm, and the calculated left and right transition winding lengths are 102.4mm, and the center angle is 2688.3°. The linear yarn width simulation diagram is as follows Figure 5 As shown, the linear yarn width simulation diagram is as follows Figure 6 shown.

[0173] The beneficial effects of the present invention are:

[0174] 1. Method 1 of the present invention calculates the winding angle from the left and right head polar holes to the cylinder body by giving the core mold parameters and slip coefficient. By comparing whether there is an intersection between the winding angles at the left and right equators, it can be quickly determined whether a linear shape can be generated, whether the given core mold parameters and slip coefficient are appropriate, and whether measures need to be taken to increase the friction coefficient to generate a linear shape, thereby effectively improving linear design and production efficiency.

[0175] 2. The present invention uses the continued fraction theory to optimize the uniform distribution of line types. Adding the deviation angle parameter and expanding the range of the number of cycles can effectively increase the number of line types that meet the requirements, making it easier to find suitable tangent line types for production.

[0176] 3. The present invention develops a fiber winding linear design simulation system through C++ programming, quickly models the system by inputting core mold parameters, quickly designs the linear parameters by inputting linear parameters, and quickly outputs the G code that matches the winding machine tool through the G code generation module for actual production, thereby effectively improving production efficiency. At the same time, in the simulation, the offset vector of the yarn dropping point is obtained by calculating its normal vector and tangent vector, and a number of offset yarn dropping points are obtained according to the offset vector. The triangular unit is constructed in the counterclockwise direction to realize the display of the fiber yarn width. The simulation can quickly determine whether the fibers are evenly distributed and whether the winding requirements are met.

[0177] 4. The present invention introduces fiber yarn width display in the linear design simulation, which can intuitively show the fiber coverage and effectively improve the fiber winding linear design and production efficiency.

[0178] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on the differences from other embodiments. Reference can be made to the common and similar parts between the various embodiments. For the systems disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the description is relatively simple, and the relevant parts can be referred to the method description.

[0179] This document uses specific examples to illustrate the principles and implementation methods of the present invention. The above examples are only intended to help understand the method and core concept of the present invention. At the same time, those skilled in the art will find that the specific implementation methods and application scopes may vary based on the concept of the present invention. In summary, the contents of this specification should not be construed as limiting the present invention.

Claims

1. A simulation method for designing the winding line of an elliptical body, characterized in that: include: Obtaining core mold parameters, and establishing a core mold model using OpenGL drawing functions according to the core mold parameters; the core mold parameters include the left transition section length, the barrel section length, the right transition section length, the barrel section major semi-axis and the barrel section minor semi-axis of the core mold; Obtaining the barrel winding angle, the left transition section slip coefficient, and the right transition section slip coefficient according to a preset transition section slip coefficient range and a preset winding transition section length or according to a preset barrel winding angle and a preset transition section slip coefficient range; Calculating winding equal points according to the core mold parameters, the preset yarn width and the barrel winding angle; Obtaining the center angle of each winding equally divided point throughout the entire round trip according to the coordinates of each winding equally divided point, the barrel winding angle, the semi-major axis of the barrel segment, the semi-minor axis of the barrel segment, the left transition section slip coefficient, and the right transition section slip coefficient; Determine the center angle of the line type based on the preset deviation angle and the center angle of each winding equally divided point throughout the entire round trip; Optimizing the left transition section slip coefficient and the right transition section slip coefficient according to the entire round-trip center angle of each winding equally divided point and the center angle of the line to obtain an optimized left transition section slip coefficient and an optimized right transition section slip coefficient; Obtaining the doffing points on the entire back-and-forth fiber path at each equally divided winding point according to the optimized left transition section slip coefficient and the optimized right transition section slip coefficient; The linear simulation is performed based on the doffing points on the entire back-and-forth fiber path at each winding equally divided point and the core model, specifically including: If yarn width simulation is not performed, linear simulation is performed based on the doffing points on the entire back-and-forth fiber path at each winding equally divided point and the core mold model; If yarn width simulation is performed, the unit normal vector of the doffing point on the fiber path of each winding equal-division point is calculated based on the center angle of the entire round trip of each winding equal-division point and the parallel circle radius of the core mold section; The unit tangent vector of the doffing point on the fiber path of each winding equal-division point is calculated according to the parallel circle radius of the core mold section, the center angle of the entire round trip of each winding equal-division point, the major semi-axis and the minor semi-axis of the barrel section; Calculate the unit offset vector of the doffing point on the entire round-trip fiber path at each equally divided winding point according to the unit normal vector of the doffing point on the entire round-trip fiber path at each equally divided winding point and the unit tangent vector of the doffing point on the entire round-trip fiber path at each equally divided winding point; Obtaining an offset doffing point corresponding to each winding equally divided point according to the preset yarn width and the unit offset vector of the doffing point on the entire back-and-forth fiber path at each winding equally divided point; Linear simulation is performed based on the offset doffing points corresponding to each winding equally divided point, the doffing points on the entire back-and-forth fiber path of each winding equally divided point, and the core mold model.

2. The simulation method for designing the elliptical fiber winding line according to claim 1, characterized in that: The calculating of the winding equal division points according to the core mold parameters, the preset yarn width and the barrel winding angle specifically includes: Calculating the number of cycles according to the core mold parameters, the preset yarn width and the barrel winding angle; The coordinates of the winding equal division points are determined according to the number of cycles.

3. The simulation method for designing the elliptical fiber winding line according to claim 1, characterized in that: The method of determining the center angle of the line according to the preset deviation angle and the center angle of each winding equally divided point in the entire round trip specifically includes: Determine the uniform distribution line type parameters based on the preset deviation angle and the center angle of each winding equal-division point throughout the entire round trip; The center angle of the line type is determined according to the uniformly distributed line type parameters.

4. The simulation method for designing the elliptical fiber winding line according to claim 2, characterized in that: The calculating of the number of cycles according to the core mold parameters, the preset yarn width and the barrel winding angle specifically includes: The perimeter of the barrel section is obtained according to the semi-major axis and the semi-minor axis of the barrel section; The number of cycles is calculated based on the cross-sectional perimeter of the barrel, the preset yarn width and the barrel winding angle.

5. A simulation system for elliptical fiber winding line design, characterized in that: include: An acquisition module is used to obtain core mold parameters and establish a core mold model using OpenGL drawing functions according to the core mold parameters; the core mold parameters include the left transition section length, the barrel section length, the right transition section length, the major semi-axis of the barrel section, and the minor semi-axis of the barrel section; A slip coefficient and barrel winding angle calculation module is used to obtain the barrel winding angle, the left transition section slip coefficient and the right transition section slip coefficient according to a preset transition section slip coefficient range and a preset winding transition section length or according to a preset barrel winding angle and a preset transition section slip coefficient range; A winding equal division point calculation module is used to calculate the winding equal division points according to the core mold parameters, the preset yarn width and the barrel winding angle; A winding equally divided point center angle calculation module is used to obtain the center angle of each winding equally divided point throughout the entire round trip based on the coordinates of each winding equally divided point, the barrel winding angle, the semi-major axis of the barrel segment, the semi-minor axis of the barrel segment, the left transition section slip coefficient, and the right transition section slip coefficient; A linear center angle calculation module is used to determine the linear center angle based on a preset deviation angle and the center angle of each winding equally divided point throughout the entire round trip; an optimization module, configured to optimize the left transition section slip coefficient and the right transition section slip coefficient according to the entire round-trip center angle of each winding equally divided point and the center angle of the line, to obtain an optimized left transition section slip coefficient and an optimized right transition section slip coefficient; a doffing point calculation module, configured to obtain the doffing points on the entire round-trip fiber path at each winding equally divided point according to the optimized left transition section slip coefficient and the optimized right transition section slip coefficient; A simulation module, for performing line simulation based on the doffing points on the entire back-and-forth fiber path at each winding equally divided point and the core mold model; The simulation module specifically includes: A first simulation unit is configured to perform line simulation based on the doffing points on the entire back-and-forth fiber path at each winding equally divided point and the core model if yarn width simulation is not performed; The unit normal vector calculation unit of the doffing point is used to calculate the unit normal vector of the doffing point on the entire fiber path of each winding equally divided point according to the center angle of the entire round trip of each winding equally divided point and the parallel circle radius of the core mold cross section when performing yarn width simulation; The unit tangent vector calculation unit of the doffing point is used to calculate the unit tangent vector of the doffing point on the entire back-and-forth fiber path of each winding equally divided point based on the parallel circle radius of the core mold cross section, the center angle of the entire back-and-forth of each winding equally divided point, the major semi-axis and the minor semi-axis of the barrel segment; a unit offset vector calculation unit for a doffing point, for calculating a unit offset vector of a doffing point on the entire round-trip fiber path at each equally divided winding point based on a unit normal vector of the doffing point on the entire round-trip fiber path at each equally divided winding point and a unit tangent vector of the doffing point on the entire round-trip fiber path at each equally divided winding point; A calculation unit, configured to obtain an offset doffing point corresponding to each equally divided winding point according to the preset yarn width and a unit offset vector of the doffing point on the entire back-and-forth fiber path at each equally divided winding point; The second simulation unit is used to perform linear simulation according to the offset doffing points corresponding to each winding equally divided point, the doffing points on the entire back-and-forth fiber path of each winding equally divided point, and the core mold model.

6. The simulation system for elliptical body fiber winding line design according to claim 5, characterized in that: The winding equal division point calculation module specifically includes: A cycle number calculation unit, used to calculate the cycle number according to the core mold parameters, the preset yarn width and the barrel winding angle; The winding equally divided point calculation unit is used to determine the coordinates of the winding equally divided point according to the number of cycles.

7. The simulation system for elliptical body fiber winding line design according to claim 5, characterized in that: The linear center angle calculation module specifically includes: A uniformly distributed linear parameter calculation unit is used to determine the uniformly distributed linear parameters according to a preset deviation angle and the center angle of each winding equally divided point throughout the entire round trip; The line center angle calculation unit is used to determine the line center angle according to the uniformly distributed line parameters.

8. The simulation system for elliptical body fiber winding line design according to claim 6, characterized in that: The cycle number calculation unit specifically includes: The barrel section perimeter calculation subunit is used to obtain the barrel section perimeter based on the barrel section's semi-major axis and the barrel section's semi-minor axis; The cycle number calculation subunit is used to calculate the cycle number according to the perimeter of the barrel cross section, the preset yarn width and the barrel winding angle.

Images