A method for modeling the performance of a curved layup of a composite laminate layup tape

Abstract

The application belongs to the technical field of performance analysis of a curved lay-up tape lay-up composite laminate, and particularly relates to a performance analysis modeling method of a curved lay-up tape lay-up composite laminate, a composite laminate geometric model is constructed; the fiber lay-up angles of each point in the composite laminate range, the point being located in a lay-up tape, a lay-up tape edge overlapping area or a lay-up tape edge gap area, are calculated; the material properties of each point corresponding to the lay-up tape, the lay-up tape edge overlapping area or the lay-up tape edge gap area, and the fiber lay-up angles are given to the composite laminate geometric model.

Description

Technical Field

[0001] This application belongs to the technical field of performance analysis of composite laminates with laid-up tape curves, and specifically relates to a modeling method for performance analysis of composite laminates with laid-up tape curves. Background Technology

[0002] Composite laminates possess excellent mechanical properties and have wide applications in engineering. Currently, most composite laminates are formed by laying strips side-by-side in straight lines at equal intervals at a certain angle. This does not fully utilize the directional characteristics of the composite strips. Therefore, designs are made to lay composite laminate strips according to a predetermined curved pattern. However, because the curvature of each strip's trajectory is constantly changing, even with equal intervals, the edges of the strips may overlap or gaps may exist. Figure 1 As shown, in order to make the surface of the composite laminate flat, the overlapping parts of the layup strip edges were trimmed in practice.

[0003] Modeling and analyzing the properties of composite laminates can quickly yield their performance and efficiently assist in their design and improvement. However, current modeling of composite laminates with laid-up tape curves is often based on the assumption that the ideal fiber prepreg tape has no width, ignoring the overlap and gaps at the edges of the laid-up tape. The performance of the composite laminates obtained by this technical solution deviates significantly from reality.

[0004] This application is made in view of the aforementioned technical deficiencies.

[0005] It should be noted that the above background information is only used to assist in understanding the inventive concept and technical solution of this invention, and it does not necessarily belong to the prior art of this patent application. In the absence of clear evidence that the above information was disclosed on the filing date of this application, the above background information should not be used to evaluate the novelty and inventiveness of this application. Summary of the Invention

[0006] The purpose of this application is to provide a method for performance analysis and modeling of composite laminates with laid-out curves, in order to overcome or mitigate at least one of the known technical defects.

[0007] The technical solution of this application is:

[0008] A method for performance analysis and modeling of composite laminates with laid-out curves includes:

[0009] Construct a geometric model of the composite laminate;

[0010] The point is located in the kth lay-up tape, the k+1th lay-up tape edge gap region, and the fiber lay-up angle is set as

[0011] The point is located in the kth lay-up tape, the k+1th lay-up tape edge gap region, and the fiber lay-up angle is set as

[0012] According to at least one embodiment of the present application, in the above-mentioned lay-up tape curve lay-up composite laminate performance analysis modeling method, the calculation of the composite laminate range, the point is located in the lay-up tape, the lay-up tape edge overlapping region or the lay-up tape edge gap region, and the fiber lay-up angle, comprises:

[0013] When the composite laminate lay-up tape is not cut:

[0014]

[0015]

[0016]

[0017]

[0018]

[0019] If The point is located in the kth lay-up tape, the k+1th lay-up tape edge gap region, and the fiber lay-up angle is set as

[0020] If The point is located in the kth lay-up tape, the k+1th lay-up tape edge overlapping region, and the fiber lay-up angle is set as

[0021] If The point is located in the kth lay-up tape, and the fiber lay-up angle is

[0022] If The point is located in the k+1th lay-up tape, and the fiber lay-up angle is

[0023] Wherein,

[0024] k is the number of the lay-up tape;

[0025] floor() is the floor function;

[0026] y0(x) is the y-axis coordinate when x is taken as the x-axis coordinate on the 0th lay-up tape center line;

[0027] x, y are coordinates of point P in the range of the composite laminate;

[0028] x0, y0 are coordinates of point on the 0th ply centerline;

[0029] S is the interval between ply centerlines;

[0030] x k , y k are coordinates of point on the kth ply centerline;

[0031] d k is the distance from point P in the range of the composite laminate to the point on the kth ply centerline;

[0032] is the shortest distance from point P in the range of the composite laminate to the point on the kth ply centerline;

[0033] is the x-axis coordinate of the shortest distance from the kth ply centerline to point P in the range of the composite laminate;

[0034] is the fiber placement angle of the kth ply centerline at the x-axis coordinate of the shortest distance from the kth ply centerline to point P in the range of the composite laminate;

[0035] θ(P) is the fiber placement angle of point P in the range of the composite laminate;

[0036] W max is the width of the ply;

[0037] x k+1 , y k+1 are coordinates of point on the k+1th ply centerline;

[0038] d k+1 is the distance from point P in the range of the composite laminate to the point on the k+1th ply centerline;

[0039] is the shortest distance from point P in the range of the composite laminate to the point on the k+1th ply centerline;

[0040] is the x-axis coordinate of the shortest distance from the k+1th ply centerline to point P in the range of the composite laminate;

[0041] is the fiber placement angle of the k+1th ply centerline at the x-axis coordinate of the shortest distance from the k+1th ply centerline to point P in the range of the composite laminate.

[0042] According to at least one embodiment of this application, in the above-described method for performance analysis and modeling of composite laminates with laid-up tape curves, the calculation of the location of each point within the range of the composite laminate, including its location in the laid-up tape, the overlapping area of ​​the laid-up tape edge, or the gap area of ​​the laid-up tape edge, and its fiber lay-up angle, includes:

[0043] When cutting composite laminate layup tapes:

[0044]

[0045]

[0046]

[0047]

[0048]

[0049]

[0050]

[0051] like The point is located in the gap between the edges of the k-th and (k+1)-th laying strips, and its fiber laying angle is set to...

[0052] like exist When the point is located in the overlapping area of ​​the edges of the k-th and (k+1)-th laying strips, its fiber laying angle is set to... exist When the point is located on the kth laying strip, its fiber laying angle is .

[0053] like exist When the point is located in the gap between the edges of the k-th and (k+1)-th laying strips, the fiber laying angle is set to... exist When the point is located on the kth laying strip, its fiber laying angle is .

[0054] like The point is located on the (k+1)th laying strip, and its fiber laying angle is...

[0055] in,

[0056] k is the number of the laying strip;

[0057] floor() rounds down to the nearest integer.

[0058] y0(x) is the y-axis coordinate of the point on the 0th tape centerline when the x-axis coordinate is x;

[0059] x, y are the coordinates of point P in the range of the composite laminate;

[0060] x0, y0 are the coordinates of the point on the 0th tape centerline;

[0061] S is the interval between the tape centerlines;

[0062] x k , y k are the coordinates of the point on the kth tape centerline;

[0063] d k is the distance from point P in the range of the composite laminate to the point on the kth tape centerline;

[0064] is the shortest distance from point P in the range of the composite laminate to the point on the kth tape centerline;

[0065] is the x-axis coordinate of the shortest distance from point P in the range of the composite laminate to the kth tape centerline;

[0066] is the fiber placement angle of the kth tape centerline at the x-axis coordinate of the kth tape centerline;

[0067] θ(P) is the fiber placement angle of point P in the range of the composite laminate;

[0068] W max is the width of the tape;

[0069] x k+1 , y k+1 are the coordinates of the point on the k+1th tape centerline;

[0070] d k+1 is the distance from point P in the range of the composite laminate to the point on the k+1th tape centerline;

[0071] is the shortest distance from point P in the range of the composite laminate to the point on the k+1th tape centerline;

[0072] is the x-axis coordinate of the shortest distance from point P in the range of the composite laminate to the k+1th tape centerline;

[0073] is the fiber placement angle of the k+1th tape centerline at the x-axis coordinate the fiber placement angle of each point in the composite laminated plate;

[0074] d F the critical point of the overlap or gap between the edges of the placement tape;

[0075] W t the width of the fiber prepreg tape in the placement tape;

[0076] FGL is the coverage rate of the placement tape cutting, which is the ratio of the area of the overlapping part of the edges of the placement tape after cutting to the sum of the area of the overlapping part and the area of the gap part;

[0077] the coordinates of the critical point of the overlap or gap between the edges of the kth and k+1th placement tapes;

[0078] the minimum distance from the coordinates of the critical point of the overlap or gap between the edges of the kth and k+1th placement tapes to the center line of the k+1th placement tape.

[0079] According to at least one embodiment of the present application, in the above-mentioned placement tape curve placement composite laminated plate performance analysis modeling method, W max = W t · N t ;

[0080] wherein,

[0081] N t is the number of fiber prepreg tapes in the placement tape.

[0082] According to at least one embodiment of the present application, in the above-mentioned placement tape curve placement composite laminated plate performance analysis modeling method, is obtained by the shortest distance search method.

[0083] The present application has at least the following beneficial technical effects:

[0084] A placement tape curve placement composite laminated plate performance analysis modeling method is provided, which considers the overlap and gap of the placement tape edge part of the placement tape curve placement composite laminated plate, calculates the fiber placement angle of each point in the composite laminated plate, which is located in the placement tape, the placement tape edge overlapping area or the placement tape edge gap area, assigns the material properties corresponding to the placement tape, the placement tape edge overlapping area or the placement tape edge gap area of each point and the fiber placement angle to the composite laminated plate geometric model, completes the modeling of the placement tape curve placement composite laminated plate, and analyzes the performance of the composite laminated plate based on this, which has high accuracy and can effectively assist the design and improvement of the composite laminated plate. BRIEF DESCRIPTION OF DRAWINGS

[0085] Figure 1Figure 1 is a schematic diagram of the composite laminated plate laying tape provided by the embodiment of the present application laid according to the set curve form at equal intervals;

[0086] Figure 2 Figure 2 is a schematic diagram of the process of the composite laminated plate performance analysis modeling method of the laying tape curve laying provided by the embodiment of the present application;

[0087] Figure 3 Figure 3 is a schematic diagram of the 0 laying tape center line provided by the embodiment of the present application;

[0088] Figure 4 Figure 4 is a schematic diagram of the effective width of the laying tape in the y direction provided by the embodiment of the present application;

[0089] Figure 5 Figure 5 is a schematic diagram of the relationship between any point on the composite laminated plate and the laying tape provided by the embodiment of the present application.

[0090] Figure 6 Figure 6 is a schematic diagram of the cutting of the laying tape provided by the embodiment of the present application;

[0091] Figure 7 Figure 7 is a schematic diagram of the critical point of the overlapping or gap part between the edges of the laying tape provided by the embodiment of the present application;

[0092] Figure 8 Figure 8 is a schematic diagram of the laying tape cutting coverage provided by the embodiment of the present application.

[0093] Figure 9 Figure 9 is a schematic diagram of the modeling of the composite laminated plate performance analysis modeling method of the laying tape curve laying when the laying tape is not cut provided by the embodiment of the present application;

[0094] Figure 10 Figure 10 is a schematic diagram of the modeling of the composite laminated plate performance analysis modeling method of the laying tape curve laying when the laying tape is cut provided by the embodiment of the present application. DETAILED DESCRIPTION

[0095] In order to make the technical solutions of the present application and its advantages clearer, the technical solutions of the present application will be further clearly and completely described in detail below, it can be understood that the specific embodiments described here are only part of the embodiments of the present application, and are only used to explain the present application, but not to limit the present application. It should be noted that, in order to facilitate the description, only the parts related to the present application are shown in the drawings, and other related parts can be referred to the general design, and in the case of no conflict, the embodiments in the present application and the technical features in the embodiments can be combined to obtain new embodiments.

[0096] Furthermore, unless otherwise defined, the technical or scientific terms used in this application description shall have the ordinary meaning understood by one of ordinary skill in the art to which this application pertains. The terms "upper," "lower," "left," "right," "center," "vertical," "horizontal," "inner," and "outer," etc., used in this application description to indicate relative direction or positional relationship are used only to indicate relative orientation or positional relationship, and do not imply that the device or component must have a specific orientation, or be constructed and operated in a specific orientation. When the absolute position of the described object changes, its relative positional relationship may also change accordingly, and therefore should not be construed as a limitation on this application. The terms "first," "second," "third," and similar terms used in this application description are used only for descriptive purposes to distinguish different components, and should not be construed as indicating or implying relative importance. The terms "a," "one," or "the," etc., used in this application description should not be construed as an absolute limitation on quantity, but should be construed as indicating the existence of at least one. The terms "including," "comprising," etc., used in this application description mean that the element or object preceding the word covers the element or object listed after the word and its equivalents, without excluding other elements or objects.

[0097] Furthermore, it should be noted that, unless otherwise explicitly specified and limited, terms such as “installation,” “connection,” and “linkage” used in the description of this application should be interpreted broadly. For example, a connection can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection or an indirect connection through an intermediate medium; or it can be a connection within two components. Those skilled in the art can understand its specific meaning in this application according to the specific circumstances.

[0098] The following is in conjunction with the appendix Figures 1 to 10 This application will be described in further detail.

[0099] The layup trajectory of the layup tape is defined using the linear angle layup method, such as... Figure 3 As shown, a rectangular coordinate system is established at the center point of the composite laminate, with the x-axis horizontal and the y-axis vertical. The angle between the fiber layup angle and the x-axis is defined as the layup angle. The angles between the tangent of the curve at the center and boundary of the composite laminate and the positive x-axis are defined as T0 and T1, respectively. This curve passes through the origin and is symmetrical about the origin. The expression for the angle between the curve and the positive x-axis is:

[0100]

[0101] where a is the length of the composite laminate in x-axis direction, and the path curve is marked as <T0|T1>. Obviously, when T0=T1, the variable stiffness composite laminate with the fiber path is converted into a conventional composite laminate, and the fiber path at other positions of the composite laminate can be obtained by translating the reference path curve along the y-axis equidistantly.

[0102] The slope of the reference path curve is:

[0103]

[0104] The trajectory of the reference path curve is:

[0105]

[0106] If the width of a single prepreg tape is Wt, and the number of tapes is Nt, then the maximum width of the prepreg tape in one placement, i.e., the placement tape width, Wmax, is:

[0107] W max = W t · N t ···(4)

[0108] The effective width of the placement tape in the y direction, Weff, is:

[0109]

[0110] where θ is the angle corresponding to a point (x, y) on the center line of the placement tape, as shown in Figure 4 .

[0111] Since the placement tape has a certain width, and the curvature radius of the curve defined by formula (3) is constantly changing, when the reference path curve is translated along the y-axis by a certain distance, a gap or overlap will be generated between the edges of the two placement tapes due to the change in the curvature radius. When the translation distance S of the placement tape changes within a certain range, theoretically, the gap or overlap region generated between the adjacent two placement tapes will also change accordingly, as follows:

[0112] If T0>T1, the translation distance S of the placement tape changes in the interval [W max / cosT1, W max / cosT0], when S=W max / cosT1, the theoretical gap is zero, and the overlap area is maximum; when S=W max / cosT0, the theoretical gap is maximum, and the overlap area is zero.

[0113] Define the center laying strip as laying strip 0, and denote the coordinates of any point on the center line of laying strip 0 as (x0, y0), which satisfies equation (3). The points on the center lines of the other laying strips can be obtained by translating the point (x0, y0) up and down along the y-axis. Let the coordinates of the point on the k-th laying strip corresponding to point (x0, y0) be (x0, yk), and thus:

[0114]

[0115] y k =y0+kS k∈Z···(6)

[0116] like Figure 5 As shown, for any point P(x, y) on the composite laminate, let's assume that point P is located between the center lines of the kth and (k+1)th laying strips. Therefore, k can be obtained from the coordinates of point P using equation (6):

[0117]

[0118] Where floor() is the floor function, we can use it to round down, and thus determine the centerline formula for the k-th and (k+1)-th laying strips:

[0119]

[0120] Let dk and dk+1 be the distances from point P to the center lines of the kth and (k+1)th laying strips, respectively. The calculation formula is as follows:

[0121]

[0122] From equations (8) and (9), the minimum distances from point P to the centerlines of the k-th and (k+1)-th paving strips can be obtained using the shortest distance search numerical method provided later: And the corresponding coordinates are:

[0123] Therefore, it can be determined which of the three cases the location of point P is in: overlapping, gap, or a single layup strip among two layup strips, as well as the possible fiber layup angle θ(P).

[0124] 1. The composite laminate layup strips are not cut:

[0125] like This point is located at the gap defect, so there is no such thing as the fiber layup angle; it can be considered as... The fiber laying angle may not need to be defined.

[0126] like This point is located in the overlapping region, where the thickness increases, and the fiber layup angle is approximately [value missing].

[0127] If The point is located in the k-th lay-up tape, and the fiber lay-up angle is approximately

[0128] If The point is located in the k+1-th lay-up tape, and the fiber lay-up angle is approximately

[0129] II. Cutting of composite laminates lay-up tapes:

[0130] Using non-cutting lay-up tapes is beneficial for structural design and component processing, but it can also lead to a large number of gaps and overlaps. In order to alleviate this situation, the pre-impregnated tape cutting technology is used, which purposefully cuts part of the pre-impregnated tape in the lay-up tape to make the adjacent lay-up tapes more compact, which can significantly reduce the process defects in the structure. The pre-impregnated tape is cut in the direction perpendicular to the fiber, so the sawtooth-shaped lay-up tape boundary is formed, as shown in Figure 6 .

[0131] It is stipulated that only Nt strips of pre-impregnated tapes are cut to form the upper boundary of the entire lay-up tape, as shown in Figure 7 , according to the minimum distance of the point P to the center line of the k-th lay-up tape and the corresponding coordinates , the distance of the critical point Fe or Ge to can be obtained:

[0132]

[0133] where dF is the distance of the critical point Fe or Ge to ;

[0134] FGL is the coverage, which is the ratio of the area of the overlapping part of the cut lay-up tape edge to the sum of the area of the overlapping part and the area of the gap part, as shown in Figure 8 .

[0135] The coordinates of the critical point

[0136]

[0137] The minimum distance of to the center line of the k+1-th lay-up tape is obtained

[0138] If The point is located at the gap defect, and the fiber lay-up angle does not matter, which can be considered as The fiber lay-up angle can also not be defined.

[0139] If This point is located in the overlap defect area without cutting, and it is needed to further determine whether it is in the overlap area when cutting. When , this point is located in the kth laying strip without cutting, and it is needed to further determine whether it is in the gap area when cutting. When , this point is located in the k+1th laying strip, and the fiber laying angle is approximately

[0140] If , this point is located in the kth laying strip without cutting, and it is needed to further determine whether it is in the gap area when cutting. When , this point is located in the gap area after cutting, , the fiber laying angle is approximately

[0141] If , this point is located in the k+1th laying strip, and the fiber laying angle is approximately

[0142] For , the shortest distance can be obtained by searching method, and the numerical method for solving the shortest distance is introduced below taking as an example.

[0143] From equations (8) and (9), we have:

[0144]

[0145] where δ(x k ) is a sign function;

[0146] Let z=x k , The golden section ratio is Φ, then:

[0147]

[0148] Obviously, the problem is transformed into finding the minimum value of f(z) in the interval z∈[x-W max , x+W max ], so the golden search method can be used to solve it, and the detailed process is as follows:

[0149] (1) Initialization

[0150] z1=x-W max , z2=x+W max , z3=z2-Φ(z2-z1), z4=z1+Φ(z2-z1) ··· (14)

[0151] (2) Calculate function value

[0152] f3 = f(z3), f4 = f(z4) ··· (15)

[0153] (3) Update interval

[0154]

[0155] (4) Determine

[0156] If |z2-z1| < TOL, continue (5); if |z2-z1| ≥ TOL, jump to step (2), TOL is a specified tolerance, which can be selected according to requirements, and TOL = 10 can be taken -4 ;

[0157] (5) End, output result

[0158]

[0159] In one specific embodiment, the size of the composite plate is 200mm x 200mm, the lay-up trajectory is <30|60> for example, and the modeling method for analyzing the performance of the composite laminated plate laid by the laid tape curve disclosed in the present application is used for modeling, and the specific steps can be referred to as follows:

[0160] 1) Assuming T0 = 30° and T1 = 60°, according to formula (1), the angle between the reference curve path at each position and the positive direction of the x-axis can be obtained

[0161] 2) Convert T0 and T1 to radians, and the curve trajectory obtained from formula (3) is:

[0162] According to this formula, the value of any point (x, y) can be obtained, so as to draw the initial reference path curve, that is, the initial center line path of the laid tape composed of Nt prepreg tapes, and then translate the reference path curve along the y-axis by a certain distance S, and judge whether the position between the center lines of the adjacent two laid tapes is in the overlapping position, the gap position, or the position of a single laid tape in the two laid tapes;

[0163] 3) According to formula (4) (5), assuming that the width of a single prepreg tape Wt is 6.35mm and the number of tapes Nt is 10, then the width of the laid tape laid at one time Wmax is 63.5mm, and the effective width is which changes with position.

[0164] 4) When the translation distance S takes , only gap defects are included, and when the translation distance S takes When only the overlapping defects are included, at the translation distance of min(S1, S2), there is a point P(10, 50) in the composite laminate range;

[0165] 5) When the composite laminate laying tape is not cut, at the translation distance of S = (S1 + S2) / 2 = 100.2:

[0166] 5.1) The point P(10, 50) is brought into formula (7), and the calculation obtains That is, the P point is in the area between the 0th and 1st laying tapes;

[0167] 5.2) According to formula (9) and the shortest distance search method mentioned in formulas (12)-(17), the minimum distances of the point P to the center lines of the 0th and 1st laying tapes are obtained and and the corresponding coordinates are:

[0168]

[0169] 5.3) Since and (63.5 / 2 is W max / 2), the point is a gap area;

[0170] The above process is programmed through Python, and the 200mm x 200mm range is divided into 1mm2 square units, and the center point coordinates of the units are taken for judgment. After traversing each unit, the laying tape distribution in the entire range can be obtained, as shown in Figure 9 .

[0171] 6) When the composite laminate laying tape is cut, at the translation distance of S = S2 = 73.32, and the coverage rate FGL = 50%, the situation is:

[0172] 6.1) The point P(10, 50) is brought into formula (7), and the calculation obtains That is, the P point is in the area between the 0th and 1st laying tapes;

[0173] 6.2) According to formula (9) and the shortest distance search method mentioned in formulas (12)-(17), the minimum distances of the point P to the center lines of the 0th and 1st laying tapes are obtained and and the corresponding coordinates are:

[0174]

[0175] 6.3) The critical distance dF obtained according to formula (10) is 38.1, and the critical point coordinates are calculated according to formula (11)

[0176] 6.4) According to formula (9) and the shortest distance search method mentioned in formula (12)-(17), the critical point coordinates are obtained The minimum distance to the center line of the first laying tape

[0177] 6.5) Since and (63.5 / 2 is W max / 2), this position is located on the first laying tape, and the fiber angle is 30.88°;

[0178] The above process is programmed by Python, which divides the 200mmx200mm range into 1mm2 square units, and takes the center point coordinates of the units for judgment. After traversing each unit, the laying tape distribution in the entire range can be obtained, such as Figure 10 as shown in the case.

[0179] Each embodiment in the specification is described in a progressive manner, and each embodiment focuses on the difference from other embodiments. The same or similar parts between each embodiment can be referred to each other.

[0180] So far, the technical solution of the present application has been described in combination with the preferred embodiments shown in the drawings. Those skilled in the art should understand that the protection scope of the present application is obviously not limited to these specific embodiments. Without deviating from the principles of the present application, those skilled in the art can make equivalent changes or replacements to related technical features. The technical solutions after these changes or replacements will fall within the protection scope of the present application.

Claims

1. A method of modeling the performance of a curved layup composite laminate layup, characterized by, The method comprises the following steps: constructing a geometric model of the composite laminate; calculating the fiber placement angle of each point in the composite laminate, which is located in the placement tape, the edge overlap area of the placement tape, or the edge gap area of the placement tape; assigning the material properties of each point corresponding to the placement tape, the edge overlap area of the placement tape, or the edge gap area of the placement tape, and the fiber placement angle to the geometric model of the composite laminate; defining the angle between the fiber placement angle of the placement tape and the x-axis as the placement angle, defining the angle between the tangent of the curve at the center and the boundary of the composite laminate and the positive direction of the x-axis as T0 and T1, respectively, the curve passing through the origin and being symmetric about the origin, a being the length of the x-axis of the composite laminate, and defining the center placement tape as the 0th placement tape, and recording the coordinates of any point on the center line of the 0th placement tape as (x0, y0), which satisfies: wherein x and y are the coordinates of point P in the composite laminate; The rest of the points on the center line of the laid strip are obtained by translating the point (x0, y0) up and down along the y-axis. The coordinates of the point corresponding to the point (x0, y0) on the kth laid strip are denoted as (x0, y k ) y k = y0+ kS k∈Z; for any point P(x, y) in the composite laminate, assuming that the point P is located between the center lines of the kth and k+1th placement tapes, k is obtained according to the coordinates of the point P: wherein floor() is the floor function, and y0(x) is the y-axis coordinate when the x-axis coordinate of the center line of the 0th placement tape is x; the placement tape translation distance is S, the center line formulas of the kth and k+1th placement tapes are determined, and the position of the point P is determined to be in which of the three situations of overlap, gap, and a single placement tape in the two placement tapes, and the fiber placement angle θ(P). The distance from the point P to the center line of the kth and k+1th laying tape is d k , d k+1 , which is calculated as The minimum value of the distance from P point to the center line of the kth and k+1th laying tape is obtained by the shortest distance search numerical method, and is respectively: And the corresponding coordinates are:

2. The performance analysis modeling method of the placement tape curve placement composite laminate according to claim 1, wherein the calculation of the fiber placement angle of each point in the composite laminate, which is located in the placement tape, the edge overlap area of the placement tape, or the edge gap area of the placement tape, comprises: wherein θ(P) is the fiber placement angle of point P in the composite laminate. If then the point is located in the edge gap region between the kth layup tape and the k+1th layup tape, and the fiber layup angle is set to If then the point is located in the edge overlap region of the kth layup tape and the k+1th layup tape, and the fiber layup angle is set to If then the point is located on the kth layup tape with a fiber layup angle of If then the point is located on the k+1th layup band with the fiber layup angle 3. The performance analysis modeling method of the placement tape curve placement composite laminate according to claim 1, wherein the calculation of the fiber placement angle of each point in the composite laminate, which is located in the placement tape, the edge overlap area of the placement tape, or the edge gap area of the placement tape, comprises: Dkis the shortest distance from point P in the composite laminate to the centerline of the kth layup tape. xkis the x-axis coordinate of the shortest distance from the kth layup tape centerline to the point P within the composite laminate layup area; the fiber placement angle for the kth layup tape centerline at the x-axis coordinate of xk. when the placement tape of the composite laminate is cut: W max W is the width of the laid tape; x k+1 , y k+1 is the coordinate of the point on the center line of the k+1th laying tape; Dk+1is the shortest distance from point P in the composite laminate to the centerline of the (k+1)th ply; xk+1= the x-axis coordinate of the shortest distance from the centerline of the k+1th layup tape to the point P within the composite laminate layup area; the fiber placement angle of the k+1th layup tape centerline at the x-axis coordinate of xk+1. wherein θ(P) is the fiber placement angle of point P in the composite laminate. FGL is the coverage rate of the cut placement tape, and FGL is the ratio of the area of the edge overlap part to the sum of the area of the overlap part and the area of the gap part after cutting. If then the point is located in the edge gap region between the kth layup tape and the k+1th layup tape, and the fiber layup angle is set to If When , the point is located on the kth laying tape, and the fiber laying angle is set as When , the point is located on the kth laying tape, and the fiber laying angle is set as If When , the point is located at the edge gap region between the kth laying tape and the k+1th laying tape, and the fiber laying angle is set as When , the point is located on the kth laying tape, and the fiber laying angle is If then the point is located on the k+1th layup band with the fiber layup angle 4. The performance analysis modeling method of the placement tape curve placement composite laminate according to claim 3, wherein Dkis the shortest distance from point P in the composite laminate to the centerline of the kth layup tape. xkis the x-axis coordinate of the shortest distance from the kth layup tape centerline to the point P within the composite laminate layup area; the fiber placement angle for the kth layup tape centerline at the x-axis coordinate of where the kth layup tape centerline intersects the x-axis. wherein W max W is the width of the laid tape; Dk+1is the shortest distance from point P in the composite laminate to the centerline of the (k+1)th ply; xk+1= the x-axis coordinate of the shortest distance from the centerline of the k+1th layup tape to the point P within the composite laminate layup area; the fiber placement angle for the k+1th layup tape centerline at the x-axis coordinate of xk+1. d F Critical points where there is an overlap or gap between the edges of the layup tape. W t W is the width of the fiber prepreg tape in the layup tape; ​ Coordinates of critical points for the k, k+1 number of the laid tape edge between the overlapping or gap parts; The minimum distance from the critical point coordinates of the overlapping or gap site between the edges of the kth and k+1th laying tapes to the center line of the k+1th laying tape. ​ W max = W t • N t ; ​ N t The number of fiber prepreg tapes in the layup tape.

Images