Method, device, equipment and product for predicting interlayer limit load of arc transition area

Abstract

The invention discloses an interlayer limit load prediction method, device, equipment and product of an arc transition area. The method comprises the following steps: carrying out image acquisition on the cross section of a composite metallographic specimen containing arc transition region characteristics to obtain a target metallographic image containing inner and outer contours; identifying each contour point in the image, and dividing the contour points into a mold side contour point set and a product side contour point set according to molding process types; fitting a mold side point set circle, and constructing a plurality of rays according to a preset polar angle interval by taking a circle center as a starting point and a preset reference direction as a zero polar angle; calculating the actual thickness under each polar angle according to the intersection point of each ray on the point sets on the two sides; and in combination with the preset load and the actual thickness, iteratively calculating a target load which enables the minimum value of the safety margin to be zero, and taking the target load as an interlayer pre-estimated limit load. According to the embodiment of the invention, by optimizing the engineering algorithm for estimating the limit load through the interlayer strength test containing the arc transition region characteristics, the interlayer failure risk is accurately quantified, and the accuracy of the interlayer strength test of the composite material is improved.

Description

Technical Field

[0001] This invention relates to the fields of materials science and engineering technology, and in particular to a method, apparatus, equipment and product for predicting the interlayer ultimate load in a circular arc transition zone. Background Technology

[0002] Composite materials, with their significant advantages such as high specific strength, high specific stiffness, and strong structural designability, have become the core material for key load-bearing components in aerospace, high-end equipment, and other fields, and are widely used in the manufacture of important structural components such as beams, ribs, and frames. During the actual service life of composite parts, these components must withstand complex loads for extended periods. However, the inherently low interlaminar strength of composite materials makes them prone to interlaminar failure under load, particularly in the arc transition zone (R-zone) of beam-rib structures. As the connection point between the flange and the web, the arc transition zone bears significant bending loads during load transfer, making it a high-risk area for interlaminar delamination failure. Therefore, during the design phase of beam-rib structures, it is essential to accurately verify the strength of the arc transition zone using engineering algorithms or finite element analysis to ensure the service safety of the components.

[0003] In existing strength tests targeting the interlaminar tensile failure load of the arc transition zone, the theoretical design thickness of the arc transition zone or the measured thickness at a local location is often simplified and used as a uniform thickness parameter for the entire arc transition zone in the calculation. This fails to fully consider the impact of the non-uniform thickness distribution characteristics of the arc transition zone caused by the forming process on its interlaminar stress distribution and strength bearing capacity. This simplified approach easily leads to deviations between the strength estimation results of the arc transition zone and the actual working conditions, making it difficult to accurately reflect the true bearing limit of the arc transition zone. Summary of the Invention

[0004] Based on this, the present invention provides a method, apparatus, equipment and product for predicting the interlayer ultimate load of the arc transition zone, in order to solve the problem that the existing strength verification algorithm does not consider the uneven thickness of the arc transition zone of composite beam rib caused by molding, resulting in large estimation deviation and difficulty in reflecting the true bearing limit.

[0005] In a first aspect, embodiments of the present invention provide a method for predicting the inter-story ultimate load in a circular arc transition zone, comprising:

[0006] Image acquisition is performed on the cross section of the target composite metallographic specimen containing the characteristics of the arc transition region to obtain a complete metallographic image of the target containing the inner and outer contours of the arc transition region.

[0007] In the target metallographic image, each contour point located in the inner and outer contours of the arc transition zone is identified, and according to the molding process type of the composite metallographic specimen, each contour point is divided into a mold side contour point set and a product side contour point set.

[0008] Perform circle fitting on the set of points on the side contour of the mold, and use the center of the fitted circle as the starting point of the ray, and the preset reference direction as the starting direction of the zero polar angle, and construct multiple rays according to the preset polar angle interval.

[0009] Based on the two contour points hit by each ray in the mold side contour point set and the product side contour point set respectively, calculate the actual thickness of the arc transition zone at each polar angle.

[0010] Based on the preset load value and the actual thickness of the arc transition zone at each polar angle, the target load that reduces the minimum safety margin of the arc transition zone to zero is calculated iteratively, and is used as the estimated interlaminar limit load for the interlaminar strength test of the target composite material arc transition zone.

[0011] Secondly, embodiments of the present invention also provide an inter-story ultimate load prediction device for a circular arc transition zone, comprising:

[0012] The image recognition module is used to acquire images of the cross-section of the target composite metallographic specimen containing the features of the arc transition zone, and obtain a complete target metallographic image containing the inner and outer contours of the arc transition zone.

[0013] The edge contour point recognition module is used to identify each contour point located in the inner and outer contours of the arc transition zone in the target metallographic image, and divide the contour points into mold side contour point set and product side contour point set according to the molding process type of the composite metallographic sample.

[0014] The circle fitting module is used to perform circle fitting on the set of points on the side contour of the mold, and to construct multiple rays with the center of the fitted circle as the starting point of the ray and the preset reference direction as the zero polar angle starting direction according to the preset polar angle interval.

[0015] The actual thickness calculation module is used to calculate the actual thickness of the arc transition zone at each polar angle based on the two contour points hit by each ray in the mold side contour point set and the product side contour point set respectively.

[0016] The interlaminar ultimate load prediction module is used to iteratively calculate the target load that reduces the minimum safety margin of the arc transition zone to zero based on the preset load value and the actual thickness of the arc transition zone at each polar angle. This target load serves as the predicted interlaminar ultimate load for the interlaminar strength test of the target composite material arc transition zone.

[0017] Thirdly, embodiments of the present invention also provide an electronic device, the electronic device comprising:

[0018] At least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores a computer program executable by the at least one processor, the computer program being executed by the at least one processor to enable the at least one processor to perform a method for predicting the interlayer ultimate load of a circular arc transition zone according to any embodiment of the present invention.

[0019] Fourthly, embodiments of the present invention also provide a computer-readable storage medium storing computer instructions, which are used to cause a processor to execute and implement the interlayer ultimate load prediction method for a circular arc transition zone as described in any embodiment of the present invention.

[0020] Fifthly, embodiments of the present invention also provide a computer program product, the computer program product including a computer program, which, when executed by a processor, implements a method for predicting the interlayer ultimate load of a circular arc transition zone as described in any embodiment of the present invention.

[0021] This invention, through the complete acquisition of metallographic images of the inner and outer contours of the arc transition zone, avoids the shortcomings of existing technologies in terms of missing local contours, providing comprehensive basic data for subsequent contour recognition and parameter calculation. By combining the molding process to divide the contour point sets on the mold side and the product side, the reference boundary of the arc transition zone and the actual product boundary are clearly defined, solving the boundary confusion problem caused by traditional mixed point sets. Starting from the fitted circle center and constructing rays at polar angle intervals, the thickness calculation of the arc transition zone is realized across all angles, avoiding the randomness of local sampling and fully reflecting the non-uniform thickness characteristics. Based on the iterative calculation of the actual thickness at each polar angle, the mechanical behavior of the initial failure of local weak points in the arc transition zone is accurately matched, solving the deviation problem of traditional load estimation based on theoretical thickness or local thickness.

[0022] It should be understood that the description in this section is not intended to identify key or essential features of the embodiments of the present invention, nor is it intended to limit the scope of the invention. Other features of the invention will become readily apparent from the following description. Attached Figure Description

[0023] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0024] Figure 1 This is a flowchart of a method for predicting the interlayer ultimate load of a circular arc transition zone according to Embodiment 1 of the present invention;

[0025] Figure 2 This is a reference metallographic image of a target that fully includes the inner and outer contours of the arc transition region, applicable to embodiments of the present invention.

[0026] Figure 3 This is an embodiment of the invention applicable to a method based on the identified contour points in the inner and outer contours.

[0027] Data point diagram of the formed circular arc transition area;

[0028] Figure 4 This is a reference diagram showing how a ray hits two contour points respectively, as applicable to an embodiment of the present invention;

[0029] Figure 5 This is a flowchart of another method for predicting the interlayer ultimate load of a circular arc transition zone according to Embodiment 2 of the present invention;

[0030] Figure 6 This is a reference diagram of a polar angle range covering the contour range of a circular arc transition area, applicable to an embodiment of the present invention;

[0031] Figure 7 This is a flowchart of the internal calculation of a standard mechanical model according to Embodiment 2 of the present invention;

[0032] Figure 8 This is a schematic diagram of the structure of an inter-layer ultimate load prediction device for a circular arc transition zone according to Embodiment 3 of the present invention;

[0033] Figure 9 This is a schematic diagram of the structure of an electronic device that implements a method for predicting the interlayer ultimate load of a circular arc transition zone according to an embodiment of the present invention. Detailed Implementation

[0034] To enable those skilled in the art to better understand the present invention, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the scope of protection of the present invention.

[0035] It should be noted that the terms "first," "second," etc., in the specification, claims, and accompanying drawings of this invention are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of the invention described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover a non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0036] Example 1

[0037] Figure 1 This is a flowchart of a method for predicting the failure load of a circular arc transition zone according to Embodiment 1 of the present invention. This embodiment is applicable to situations where the failure strength of the circular arc transition zone of a composite material beam rib is predicted by calculating the ultimate load value of the interlaminar strength test. This method can be executed by a device for predicting the failure load of a circular arc transition zone. This device can be implemented in hardware and / or software and can be configured in an automatic composite material strength testing system. Figure 1 As shown, the method includes:

[0038] S110. Image acquisition is performed on the cross section of the target composite metallographic specimen containing the characteristics of the arc transition zone to obtain a complete target metallographic image containing the inner and outer contours of the arc transition zone; wherein, the metallographic specimen is obtained by metallographic sample preparation process of the interlaminar strength test specimen of the arc transition zone of the target composite material.

[0039] The target composite metallographic specimen refers to the specimen for interlaminar strength testing of the composite material's arc transition zone, which has undergone metallographic preparation processes such as cutting, grinding, and etching, and is used for microstructural observation. Its cross-section includes the arc transition zone to be analyzed. The target metallographic image is an image taken of the specimen cross-section using a metallographic microscope or industrial camera, and must clearly show the inner and outer contours of the arc transition zone. Figure 2 This is a reference image of the target metallographically, showing the complete inner and outer contours of the arc transition zone, as observed using a metallographic fiber microscope. The arc transition zone (R-zone) refers to the arc-shaped transition area at the connection between the flange and web in composite beam-rib structures, such as... Figure 2The rounded corners of the L-shaped structure are critical areas prone to delamination failure under load. Since the rounded transition zone of the composite material is a three-dimensional structure, its thickness is the distance from the inner contour to the outer contour. Missing any contour would render thickness calculations unfounded, while a complete contour can cover the entire angular range of the rounded transition zone from 0° to 90°, avoiding the biases of traditional local sampling measurements.

[0040] S120. In the target metallographic image, identify each contour point located in the inner and outer contours of the arc transition zone, and divide each contour point into a mold side contour point set and a product side contour point set according to the molding process type of the composite metallographic sample.

[0041] Contour points refer to pixels extracted using the algorithm of "grayscale conversion → filtering and noise reduction → Canny edge detection → contour tracking". Background noise points need to be excluded, and only continuous pixels at the inner and outer boundaries of the arc transition area are retained, such as... Figure 3 The diagram shows a data point map of the arc transition zone formed by the contour points in the identified inner and outer contours. The molding process type is the key to determining the mold contact side of the arc transition zone. When the female mold is used for molding, the mold cavity encloses the product, the outer surface of the product is constrained by the mold, and the inner surface is free (therefore, the inner contour is the mold side, and the outer contour is the product side); when the male mold is used for molding, the mold core supports the inner side of the product, the inner surface of the product is constrained by the mold, and the outer surface is free (therefore, the outer contour is the mold side, and the inner contour is the product side).

[0042] Optionally, identifying each contour point located within the inner and outer contours of the arc transition zone may include:

[0043] The target metallographic image is preprocessed by grayscale conversion and filtering to obtain a preprocessed image. Edge pixels in the target metallographic image are obtained by edge recognition of the preprocessed image.

[0044] The edge pixels are searched to obtain at least one closed contour chain, and two target contour chains representing the inner and outer contours of the arc transition area are selected from the contour chains as contour points in the inner and outer contours of the arc transition area.

[0045] Edge detection algorithms, such as the Canny operator, are used to identify pixels in an image where grayscale values ​​change abruptly (i.e., boundary points of contours). Specifically, the contour of a circular transition area exhibits a significant brightness difference from its surroundings (e.g., a sudden change in brightness at the contour). The algorithm calculates the grayscale gradient around each pixel and filters out points whose gradient values ​​exceed a threshold; these points are the edge pixels that constitute the contour. A contour tracking algorithm then connects these discrete edge pixels into a continuous contour chain (i.e., a closed curve with its ends connected). Based on the geometric characteristics of the circular transition area, irrelevant contours are eliminated from multiple contour chains, retaining only the two contour chains that conform to the characteristics of the circular transition area.

[0046] Furthermore, based on the molding process type of the composite metallographic specimen, each contour point is divided into a mold-side contour point set and a product-side contour point set, which may include:

[0047] Obtain the molding process type of the target composite metallographic specimen from the molding process document;

[0048] If the molding process type is female mold molding, then the inner contour point set located in the arc transition area of ​​the target metallographic image is defined as the mold side contour data point set, and the outer contour point set located in the arc transition area of ​​the target metallographic image is defined as the product side contour data point set.

[0049] If the molding process type is positive mold molding, then the outer contour point set located in the arc transition zone of the target metallographic image is defined as the mold side contour data point set, and the inner contour point set located in the arc transition zone of the target metallographic image is defined as the product side contour data point set.

[0050] The molding process document is the authoritative document recording the molding method of composite material specimens, clearly indicating whether the specimen is molded using a female mold or a male mold. The core characteristics of female mold molding are: the mold has a cavity structure, the product is formed within the mold cavity, the inner wall of the mold is in direct contact with the inner surface of the product, and the outer surface of the product is freely formed. Therefore, the inner contour of the arc transition zone in the metallographic image is the boundary in direct contact with the female mold, and its shape is guaranteed by the mold's precision, thus defined as the set of points on the mold's side contour; the outer contour (far) is the boundary where the product is freely formed, thus defined as the set of points on the product's side contour. The core characteristics of male mold molding are: the mold has a core structure, the product is formed by wrapping around the mold core, the outer wall of the mold is in direct contact with the outer surface of the product, and the inner surface of the product is freely formed. Therefore, the outer contour of the arc transition zone in the metallographic image is the boundary in direct contact with the male mold, and its shape is guaranteed by the mold's precision, thus defined as the set of points on the mold's side contour; the inner contour is the boundary where the product is freely formed, thus defined as the set of points on the product's side contour.

[0051] S130. Perform circle fitting on the set of points on the side contour of the mold, and use the center of the fitted circle as the starting point of the ray, and the preset reference direction as the starting direction of the zero polar angle, to construct multiple rays according to the preset polar angle interval.

[0052] Circular fitting refers to using algorithms such as the least squares method to fit the set of points on the mold side profile to obtain an ideal circle that best approximates the distribution of those points. The preset reference direction refers to the angular reference direction of the polar coordinate system, usually defined as the horizontal direction parallel to the sample cross-section (such as the positive x-axis), ensuring the consistency of angle measurements. The zero polar angle starting direction refers to the direction with a polar angle of 0°, i.e., the preset reference direction, which is the starting point of the ray angle. The preset polar angle interval is the angle difference between two adjacent rays; the smaller the interval, the more precise the subsequent thickness calculation. In this embodiment, the ray is a straight line extending from the center of the circle along a specific angle, used to measure the thickness of the arc transition zone in that direction. The reason for choosing to perform circular fitting on the mold side is that the product side profile may deform due to material flow (such as local thickening of the outer profile during female mold forming), and the center of the circle will shift after fitting. However, the mold side profile is guaranteed by the mold itself, and the fitted circle center can accurately reflect the geometric center of the arc transition zone, ensuring that the measured thickness is the true radial thickness.

[0053] S140. Calculate the actual thickness of the arc transition zone at each polar angle based on the two contour points hit by each ray in the mold side contour point set and the product side contour point set.

[0054] Each ray is an infinitely extending straight line. Points within the arc transition zone that are closest to the ray and located within the contour points on either the mold side or the product side need to be selected. These are defined as the first type of intersection (mold side) and the second type of intersection (product side), respectively. Since the ray originates from the center of the circle, this distance is precisely the radial thickness of the arc transition zone at that polar angle, and the radial direction is the primary direction of interlayer shear stress when the arc transition zone bears load. Existing engineering algorithms assume the arc transition zone has a uniform thickness, but the actual thickness after molding may gradually change from 2mm at 0° to 3mm at 90°; different thicknesses correspond to different stresses. Only by calculating the thickness per polar angle can accurate input be provided for subsequent point-by-point stress calculation. If a female mold is used in this embodiment, the mold side contour is the inner contour of the arc transition zone, and the product side contour is the outer contour of the arc transition zone. For example... Figure 4 As shown, the ray originates from the center of the circle and first hits the mold side contour point (the point on the inner contour L1) as the first type of intersection point (point A in the figure); then it hits the product side contour point (the point on the outer contour L2) as the second type of intersection point.

[0055] S150. Based on the preset load value and the actual thickness of the arc transition zone at each polar angle, iteratively calculate the target load that reduces the minimum safety margin of the arc transition zone to zero, which serves as the estimated interlaminar limit load for the interlaminar strength test of the target composite material arc transition zone.

[0056] The preset load value is a reference value for the initial load set at the start of the calculation and serves as the starting point for subsequent iterative calculations. Since the ultimate load of the interlaminar strength test cannot be obtained directly in one step, a value lower than the ultimate load must first be selected, and then gradually adjusted to approach the true ultimate load. The safety margin measures the strength reserve of the arc transition zone (R zone) under the current load, reflecting how much additional load the arc transition zone can withstand before failure. A larger safety margin indicates sufficient strength reserve and greater safety; a smaller safety margin indicates closer proximity to the bearing limit. Because the actual thickness of the arc transition zone is uneven at different polar angle positions, the calculated safety margins at each polar angle position differ. Among all the safety margins at all polar angle positions, the smallest value is the minimum safety margin, corresponding to the weakest point of the arc transition zone.

[0057] The preset load is positively correlated with the actual stress, but the actual stress is also related to the thickness. Therefore, a load with a safety margin of 0 cannot be obtained directly from a single calculation. Thus, an iterative process is needed to approximate the critical load by gradually increasing the preset load to calculate the minimum safety margin. Failure in the circular arc transition zone of composite materials is a progressive failure caused by local weak points, rather than a simultaneous overall failure. Even if the safety margin at other polar angles does not reach 0, as long as the minimum value drops to 0, delamination will occur at the weak point, leading to overall failure. Therefore, using the minimum safety margin as the criterion can accurately capture the critical load at which failure begins.

[0058] This invention, through the complete acquisition of metallographic images of the inner and outer contours of the arc transition zone, avoids the shortcomings of existing technologies in terms of missing local contours, providing comprehensive basic data for subsequent contour recognition and parameter calculation. By combining the molding process to divide the contour point sets on the mold side and the product side, the reference boundary of the arc transition zone and the actual product boundary are clearly defined, solving the boundary confusion problem caused by traditional mixed point sets. Starting from the fitted circle center and constructing rays at polar angle intervals, the thickness calculation of the arc transition zone is realized across all angles, avoiding the randomness of local sampling and fully reflecting the non-uniform thickness characteristics. Based on the iterative calculation of the actual thickness at each polar angle, the mechanical behavior of the initial failure of local weak points in the arc transition zone is accurately matched, solving the deviation problem of traditional load estimation based on theoretical thickness or local thickness.

[0059] Example 2

[0060] Figure 5 This is a flowchart of another method for predicting the inter-story ultimate load in a circular arc transition zone, provided in Embodiment 2 of the present invention. This embodiment is a refinement based on Embodiment 1, specifically as follows: Figure 5 As shown, the method includes:

[0061] S510. Image acquisition is performed on the cross section of the target composite metallographic specimen containing the characteristics of the arc transition zone to obtain a complete target metallographic image containing the inner and outer contours of the arc transition zone; wherein, the metallographic specimen is obtained by metallographic sample preparation process of the interlaminar strength test specimen of the arc transition zone of the target composite material.

[0062] S520. In the target metallographic image, identify each contour point located in the inner and outer contours of the arc transition zone, and divide each contour point into a mold side contour point set and a product side contour point set according to the molding process type of the composite metallographic sample.

[0063] S530. Perform circle fitting on the set of points on the side contour of the mold, and set the center coordinates of the fitted circle as the origin of the polar coordinate system. Define the positive x-axis of the polar coordinate system as the starting direction of the zero polar angle and the positive y-axis as the ending direction of the polar angle 90°, forming a polar angle interval covering the contour range of the arc transition area.

[0064] like Figure 6 As shown, the positive x-axis direction (zero polar angle starting direction) refers to the direction defined horizontally to the right in the polar coordinate system (in the same direction as the horizontal pixel axis of the target metallographic image), and also corresponds to the connection direction between the arc transition zone and the flange in the composite beam-rib structure. The polar angle value of this direction is defined as 0°, which is the starting reference for all ray angles. The positive y-axis direction (polar angle 90° ending direction) refers to the direction defined vertically upward in the polar coordinate system, and also corresponds to the connection direction between the arc transition zone and the web (the web is usually a vertical load-bearing member, and this direction is the end of the load transfer in the arc transition zone). The polar angle value of this direction is defined as 90°, which is the ending reference for the effective analysis range of the arc transition zone.

[0065] S540. Starting from the zero polar angle, multiple rays are generated sequentially along the polar angle increasing direction at preset polar angle intervals, with each ray corresponding to a unique polar angle value.

[0066] Based on the preset polar angle interval, all rays from 0° to 90° are generated sequentially. Each ray starts from the origin of the polar coordinate system and extends to the corresponding polar angle direction beyond the boundary of the target metallographic image. Generating rays at intervals ensures that there is a dedicated ray at each angle position of the arc transition zone, fully capturing the non-uniform distribution characteristics of the thickness of the arc transition zone.

[0067] S550. Based on the two contour points hit by each ray in the mold side contour point set and the product side contour point set, calculate the actual thickness of the arc transition zone at each polar angle.

[0068] Optionally, based on the two contour points hit by each ray in the mold-side contour point set and the product-side contour point set, the actual thickness of the arc transition zone at each polar angle can be calculated, which may include:

[0069] The intersection of each ray with the set of points on the mold side contour is defined as the first type of intersection, and the intersection with the set of points on the product side contour is defined as the second type of intersection.

[0070] Calculate the straight-line distance between the first type of intersection point and the second type of intersection point under each ray;

[0071] The straight-line distances mentioned above are used as the actual thickness of the arc transition zone under the corresponding polar angle for each ray.

[0072] The first type of intersection is the intersection of each ray with the mold side profile point set. The mold side profile point set is the profile directly in contact with the mold; therefore, this intersection point is the boundary point of the arc transition zone closer to the mold. The second type of intersection is the intersection of each ray with the product side profile point set. The product side profile point set is the profile of the product's free-form shape; therefore, this intersection point is the boundary point of the arc transition zone farther from the mold. For each ray, the calculated straight-line distance is directly used as the actual thickness of the arc transition zone at the corresponding polar angle, forming a mapping relationship between polar angle and thickness, thus achieving precise quantification of the non-uniform thickness of the arc transition zone.

[0073] S560. Based on the preset load value and the actual thickness of the arc transition zone at each polar angle, iteratively calculate the target load that reduces the minimum safety margin of the arc transition zone to zero, and use it as the estimated interlaminar limit load for the interlaminar strength test of the target composite material arc transition zone.

[0074] Furthermore, based on the preset load value and the actual thickness of the arc transition zone at each polar angle, the target load that minimizes the safety margin of the arc transition zone to zero is iteratively calculated, including:

[0075] The preset load value is obtained as the current load. The current load and the actual thickness of the arc transition zone at each polar angle are input into the standard mechanical model. The safety margin values ​​corresponding to each polar angle are calculated to form a safety margin distribution set. The current minimum safety margin value is found in the safety margin distribution set.

[0076] Determine whether the current minimum safety margin value meets the margin value critical condition;

[0077] If so, then the current load is determined as the target load;

[0078] If not, the current load is updated using a numerical iterative algorithm, and the updated load value is used as the new current load. The operation of inputting the standard mechanical model to calculate the safety margin value is repeated.

[0079] When the current minimum safety margin value meets the margin value critical condition, the current load corresponding to the current minimum safety margin value will be determined as the target load.

[0080] Reaching the limit state is a critical transition from having strength reserves to having none. However, the relationship between load and safety margin is non-linear, making it impossible to intuitively determine at what load the limit will be reached. An iterative approximation strategy is employed, starting with a preset load lower than the expected load. Through a cycle of "adjusting the preset load → recalculating the safety margin for all polar angles → determining whether the critical point has been reached," the range is gradually narrowed down until the inter-story limit load that makes the minimum safety margin satisfy the critical condition is precisely found. The critical condition is a quantitative standard technically defined as the critical state of failure. When the safety margin is 0, the actual stress equals the allowable stress, which is technically considered as just reaching the failure limit, i.e., satisfying the margin value critical condition.

[0081] Optionally, the current load and the actual thickness of the circular arc transition zone at each polar angle are input into the standard mechanical model to calculate the safety margin values ​​corresponding to each polar angle, including:

[0082] Based on the preset load and the cross-sectional width of the interlaminar strength test specimen in the arc transition zone of the target composite material. Calculate the tensile force flow respectively shear flow and bending moment flow ;in, , , ; This is the maximum actual size of the arc transition zone along the 90° polar angle ending direction; In the direction of zero polar angle, the distance from the center of the fitted circle to the second type of intersection point; and Fixed at 0; To retrieve from the molding process file or to measure manually in advance;

[0083] Obtain the theoretical fillet radius from the molding process document. and the actual thickness at the current polar angle Using the formula Calculate the characteristic radius R between the circular arc transition region and the fitted circle;

[0084] Based on stretch flow Shear flow Bending moment flow And the characteristic radius R, calculate the polar angle. Bending moment flow at the point and polar angle Shear flow at the location ;in, The formula for calculating shear flow is: ;

[0085] Based on the bending moment flow and the shear flow Calculate the polar angles separately Interlaminar shear stress at the location and polar angle Principal stress in the thickness direction at the location The formula for calculating interlayer shear stress is as follows: The formula for calculating the principal stress in the thickness direction is: Based on the interlaminar shear stress and the principal stress in the thickness direction Calculate polar angle Safety margin value at The formula for calculating the safety margin value is as follows: ; To allow maximum shear stress, Both are inherent property parameters of the target composite material, representing the maximum allowable normal stress.

[0086] In this embodiment, we focus on the mechanical scenario of interlaminar failure in the circular arc transition zone of composite materials: in the interlaminar tensile test, the target composite metallographic specimen... Under the influence of the arc transition zone, the interlaminar tensile stress is mainly borne, while the interlaminar shear stress is relatively small. In this embodiment, the arc transition zone is... and Set to 0, so that the entire computational link is completely around the interlayer tensile load. The logical progression of "→inter-story bending moment flow-transverse shear flow→principal stress and inter-story shear stress→inter-story safety margin" allows for a more precise quantification of the risk of inter-story failure, thereby improving the pertinence and accuracy of predictions for this type of failure mode.

[0087] Due to longitudinal force The value is 0, therefore the tensile force flow is zero. The force component per unit width is 0; shear flow. The physical meaning is: to According to cross-sectional width On average, the shear strength per unit width is obtained, ensuring that the shear force matches the width dimension of the circular transition zone; due to the bending moment component... The value is fixed at 0, therefore Only the shear force component passes through the lever arm This generates an additional bending moment, i.e. By integrating the bending moment effect of interlayer shear force, the reference bending moment density of the arc transition zone is obtained; the theoretical fillet radius is obtained from the molding process documents. Combined with the actual thickness at the current polar angle The calculation shows that the circular arc transition region is at the polar angle. The average radius at that point, and the distance from the center of the circle to the product side, is... .

[0088] Bending Moment Flow Polar angle The bending moment strength under all loads is taken as the reference bending moment flow. Based on, superimposed shear flow Due to polar angle The resulting additional bending moment ultimately yields the actual bending moment strength at the current polar angle; shear flow Polar angle The combined shear strength at the point is the original shear flow. The component at the current polar angle ( (This refers to the directional projection coefficient), which yields the actual shear strength borne at the current polar angle; principal stress in the thickness direction. From bending moment flow The resulting thickness-direction normal stress reflects the bending moment, actual thickness, and geometric curvature. , The coupling effect of ) For shear flow The resulting interlaminar shear stress reflects the parabolic distribution characteristics of shear stress in a rectangular cross-section, focusing on the risk of interlaminar shear failure; safety margin value. For comprehensive and This addresses the risk of combined failure, adapting to the "normal stress + shear stress" combined failure mode under interlaminar tensile loading. For ease of understanding... Figure 7 A flowchart of the internal calculations of a standard mechanical model is shown.

[0089] This invention primarily describes the specific calculation process of the safety margin. Through a progressive approach of "load decomposition → geometric coupling → local mechanical quantity synthesis → stress calculation → composite safety margin assessment," the complex geometric non-uniformity of the circular arc transition zone is coupled with the multi-component nature of the load, ultimately achieving precise quantification of the polar angle. The method addresses the issue of low evaluation accuracy caused by neglecting the local geometry and stress coupling in the arc transition zone, which is a safe state.

[0090] Example 3

[0091] Figure 8 This is a schematic diagram of the structure of a device for predicting the interlayer ultimate load in a circular arc transition zone, provided in Embodiment 3 of the present invention. Figure 8 As shown, the device includes:

[0092] The image recognition module 810 is used to acquire images of the cross section of the target composite metallographic specimen containing the features of the arc transition zone, and obtain a complete target metallographic image containing the inner and outer contours of the arc transition zone.

[0093] The edge contour point recognition module 820 is used to identify each contour point located in the inner and outer contours of the arc transition zone in the target metallographic image, and divide the contour points into mold side contour point set and product side contour point set according to the molding process type of the composite metallographic sample.

[0094] The circle fitting module 830 is used to perform circle fitting on the set of points on the side contour of the mold, and to construct multiple rays with the center of the fitted circle as the starting point of the ray and the preset reference direction as the zero polar angle starting direction according to the preset polar angle interval.

[0095] The actual thickness calculation module 840 is used to calculate the actual thickness of the arc transition zone at each polar angle based on the two contour points hit by each ray in the mold side contour point set and the product side contour point set respectively.

[0096] The interlaminar ultimate load prediction module 850 is used to iteratively calculate the target load that reduces the minimum safety margin of the arc transition zone to zero based on the preset load value and the actual thickness of the arc transition zone at each polar angle, and serves as the predicted interlaminar ultimate load for the interlaminar strength test of the target composite material arc transition zone.

[0097] This invention, through the complete acquisition of metallographic images of the inner and outer contours of the arc transition zone, avoids the shortcomings of existing technologies in terms of missing local contours, providing comprehensive basic data for subsequent contour recognition and parameter calculation. By combining the molding process to divide the contour point sets on the mold side and the product side, the reference boundary of the arc transition zone and the actual product boundary are clearly defined, solving the boundary confusion problem caused by traditional mixed point sets. Starting from the fitted circle center and constructing rays at polar angle intervals, the thickness calculation of the arc transition zone is realized across all angles, avoiding the randomness of local sampling and fully reflecting the non-uniform thickness characteristics. Based on the iterative calculation of the actual thickness at each polar angle, the mechanical behavior of the initial failure of local weak points in the arc transition zone is accurately matched, solving the deviation problem of traditional load estimation based on theoretical thickness or local thickness.

[0098] Optionally, based on the above embodiments, the edge contour point recognition module 820 may include:

[0099] An edge recognition unit is used to perform grayscale conversion and filtering and noise reduction preprocessing on the target metallographic image to obtain a preprocessed image, and to obtain edge pixels in the target metallographic image by performing edge recognition on the preprocessed image.

[0100] The contour point filtering unit is used to perform contour search on the edge pixels to obtain at least one closed contour chain, and to filter out two target contour chains representing the inner and outer contours of the arc transition area respectively from the contour chain, as each contour point in the inner and outer contours of the arc transition area.

[0101] Optionally, based on the above embodiments, the edge contour point recognition module 820 may further include:

[0102] A molding process type determination unit is used to obtain the molding process type of the target composite metallographic specimen from the molding process file;

[0103] The first type of contour point definition unit is used to define the inner contour point set located in the arc transition area of ​​the target metallographic image as the mold side contour data point set, and the outer contour point set located in the arc transition area of ​​the target metallographic image as the product side contour data point set if the molding process type is female mold molding.

[0104] The second type of contour point definition unit is used to define the outer contour point set located in the arc transition area of ​​the target metallographic image as the mold side contour data point set if the molding process type is positive mold molding, and to define the inner contour point set located in the arc transition area of ​​the target metallographic image as the product side contour data point set.

[0105] Optionally, based on the above embodiments, the circle fitting module 830 may include:

[0106] The polar angle direction definition unit is used to set the center coordinates of the fitted circle as the origin of the polar coordinate system, define the positive x-axis direction of the polar coordinate system as the zero polar angle starting direction, and the positive y-axis direction as the polar angle 90° ending direction, forming a polar angle interval covering the contour range of the arc transition area.

[0107] The ray generation unit is used to generate multiple rays sequentially from the zero polar angle starting direction along the polar angle increasing direction at preset polar angle intervals, with each ray corresponding to a unique polar angle value.

[0108] Optionally, based on the above embodiments, the actual thickness calculation module 840 may include:

[0109] The intersection point definition unit is used to define the intersection point of each ray with the mold side contour point set as the first type of intersection point, and the intersection point with the product side contour point set as the second type of intersection point;

[0110] The intersection distance calculation unit is used to calculate the straight-line distance between the first type of intersection point and the second type of intersection point under each ray;

[0111] The actual thickness determination unit is used to take the straight line distance as the actual thickness of the arc transition zone under the corresponding polar angle of each ray.

[0112] Optionally, based on the above embodiments, the inter-story ultimate load estimation module 850 may include:

[0113] The safety margin value calculation unit is used to obtain a preset load value as the current load, input the current load and the actual thickness of the arc transition zone at each polar angle into the standard mechanical model, calculate the safety margin values ​​corresponding to each polar angle, form a safety margin distribution set, and find the current minimum safety margin value in the safety margin distribution set.

[0114] The margin value critical condition judgment unit is used to determine whether the current minimum safety margin value meets the margin value critical condition;

[0115] The first destructive load determination unit is used to determine the current load as the target load if the load is such that it is.

[0116] The load value iteration unit is used to, if not, update the current load using a numerical iteration algorithm, take the updated load value as the new current load, and repeatedly execute the operation of inputting the standard mechanical model to calculate the safety margin value.

[0117] The second destructive load determination unit is used to determine the current load corresponding to the current minimum safety margin value as the target load until the current minimum safety margin value meets the margin value critical condition.

[0118] Optionally, based on the above embodiments, the safety margin value calculation unit can also be used to calculate the safety margin value based on the preset load. and the cross-sectional width of the interlaminar strength test specimen in the arc transition zone of the target composite material. Calculate the tensile force flow respectively shear flow and bending moment flow ;in, , , ; This is the maximum actual size of the arc transition zone along the 90° polar angle ending direction; In the direction of zero polar angle, the distance from the center of the fitted circle to the second type of intersection point; and Fixed at 0; To retrieve from the molding process file or to measure manually in advance;

[0119] Obtain the theoretical fillet radius from the molding process document. and the actual thickness at the current polar angle Using the formula Calculate the characteristic radius R between the circular arc transition region and the fitted circle;

[0120] Based on stretch flow Shear flow Bending moment flow And the characteristic radius R, calculate the polar angle. Bending moment flow at the point and polar angle Shear flow at the location ;in, The formula for calculating shear flow is: ;

[0121] Based on the bending moment flow and the shear flow Calculate the polar angles separately Interlaminar shear stress at the location and polar angle Principal stress in the thickness direction at the location The formula for calculating interlayer shear stress is as follows: The formula for calculating the principal stress in the thickness direction is: Based on the interlaminar shear stress and the principal stress in the thickness direction Calculate polar angle Safety margin value at The formula for calculating the safety margin value is as follows: ; To allow maximum shear stress, Both are inherent property parameters of the target composite material, representing the maximum allowable normal stress.

[0122] The inter-layer ultimate load prediction device for a circular arc transition zone provided in this embodiment of the invention can execute the inter-layer ultimate load prediction method for a circular arc transition zone provided in any embodiment of the invention, and has the corresponding functional modules and beneficial effects of the execution method.

[0123] Example 4

[0124] Figure 9 A schematic diagram of an electronic device 10 that can be used to implement embodiments of the present invention is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital assistants, cellular phones, smartphones, wearable devices (e.g., helmets, glasses, watches, etc.), and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely illustrative and are not intended to limit the implementation of the invention described and / or claimed herein.

[0125] like Figure 9As shown, the electronic device 10 includes at least one processor 11 and a memory, such as a read-only memory (ROM) 12 or a random access memory (RAM) 13, communicatively connected to the at least one processor 11. The memory stores computer programs executable by the at least one processor. The processor 11 can perform various appropriate actions and processes based on the computer program stored in the ROM 12 or loaded from storage unit 18 into the RAM 13. The RAM 13 can also store various programs and data required for the operation of the electronic device 10. The processor 11, ROM 12, and RAM 13 are interconnected via a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0126] Multiple components in electronic device 10 are connected to I / O interface 15, including: input unit 16, such as keyboard, mouse, etc.; output unit 17, such as various types of displays, speakers, etc.; storage unit 18, such as disk, optical disk, etc.; and communication unit 19, such as network card, modem, wireless transceiver, etc. Communication unit 19 allows electronic device 10 to exchange information / data with other devices through computer networks such as the Internet and / or various telecommunications networks.

[0127] Processor 11 can be a variety of general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special-purpose artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. Processor 11 performs the various methods and processes described above, such as a method for predicting interlayer limit loads in a circular arc transition zone.

[0128] That is: to acquire images of the cross section of the target composite metallographic specimen containing the features of the arc transition zone, and obtain a complete target metallographic image containing the inner and outer contours of the arc transition zone;

[0129] In the target metallographic image, each contour point located in the inner and outer contours of the arc transition zone is identified, and according to the molding process type of the composite metallographic specimen, each contour point is divided into a mold side contour point set and a product side contour point set.

[0130] Perform circle fitting on the set of points on the side contour of the mold, and use the center of the fitted circle as the starting point of the ray, and the preset reference direction as the starting direction of the zero polar angle, and construct multiple rays according to the preset polar angle interval.

[0131] Based on the two contour points hit by each ray in the mold side contour point set and the product side contour point set respectively, calculate the actual thickness of the arc transition zone at each polar angle.

[0132] Based on the preset load value and the actual thickness of the arc transition zone at each polar angle, the target load that reduces the minimum safety margin of the arc transition zone to zero is calculated iteratively, and is used as the estimated interlaminar limit load for the interlaminar strength test of the target composite material arc transition zone.

[0133] In some embodiments, a method for predicting the inter-layer ultimate load of a circular arc transition zone can be implemented as a computer program tangibly contained in a computer-readable storage medium, such as storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed on electronic device 10 via ROM 12 and / or communication unit 19. When the computer program is loaded into RAM 13 and executed by processor 11, one or more steps of the method for predicting the inter-layer ultimate load of a circular arc transition zone described above can be performed. Alternatively, in other embodiments, processor 11 can be configured to perform a method for predicting the inter-layer ultimate load of a circular arc transition zone by any other suitable means (e.g., by means of firmware).

[0134] Various embodiments of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field-programmable gate arrays (FPGAs), application-specific integrated circuits (ASICs), application-specific standard products (ASSPs), systems-on-a-chip (SoCs), payload-programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments may include implementations in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which may be a dedicated or general-purpose programmable processor, capable of receiving data and instructions from a storage system, at least one input device, and at least one output device, and transmitting data and instructions to the storage system, the at least one input device, and the at least one output device.

[0135] Computer programs used to implement the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, such that when executed by the processor, the computer programs cause the functions / operations specified in the flowcharts and / or block diagrams to be performed. The computer programs may be executed entirely on a machine, partially on a machine, or as a standalone software package, partially on a machine and partially on a remote machine, or entirely on a remote machine or server.

[0136] In the context of this invention, a computer-readable storage medium can be a tangible medium that may contain or store a computer program for use by or in conjunction with an instruction execution system, apparatus, or device. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination thereof. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. More specific examples of machine-readable storage media include electrical connections based on one or more wires, portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fibers, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof.

[0137] To provide interaction with a user, the systems and techniques described herein can be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and pointing device (e.g., a mouse or trackball) through which the user provides input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including sound input, voice input, or tactile input).

[0138] The systems and technologies described herein can be implemented in computing systems that include backend components (e.g., as data servers), or middleware components (e.g., application servers), or frontend components (e.g., user computers with graphical user interfaces or web browsers through which users can interact with implementations of the systems and technologies described herein), or any combination of such backend, middleware, or frontend components. The components of the system can be interconnected via digital data communication of any form or medium (e.g., communication networks). Examples of communication networks include local area networks (LANs), wide area networks (WANs), blockchain networks, and the Internet.

[0139] A computing system can include clients and servers. Clients and servers are generally located far apart and typically interact through communication networks. The client-server relationship is created by computer programs running on the respective computers and having a client-server relationship with each other. The server can be a cloud server, also known as a cloud computing server or cloud host, which is a hosting product within the cloud computing service system to address the shortcomings of traditional physical hosts and VPS services, such as high management difficulty and weak business scalability.

[0140] It should be understood that the various forms of processes shown above can be used, with steps reordered, added, or deleted. For example, the steps described in this invention can be executed in parallel, sequentially, or in different orders, as long as the desired result of the technical solution of this invention can be achieved, and this is not limited herein.

[0141] The specific embodiments described above do not constitute a limitation on the scope of protection of this invention. Those skilled in the art should understand that various modifications, combinations, sub-combinations, and substitutions can be made according to design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of this invention should be included within the scope of protection of this invention.

Claims

1. A method for predicting interlayer limit load of a circular-arc transition region, characterized by, The method comprises the following steps: An image of a cross section of a target composite material metallographic sample containing a circular arc transition region feature is collected to obtain a target metallographic image containing the inner and outer contours of the circular arc transition region; the metallographic sample is obtained by processing a target composite material circular arc transition region interlayer strength test sample by a metallographic sample preparation process; In the target metallographic image, each contour point located in the inner and outer contours of the circular arc transition region is identified, and according to the forming process type of the composite material metallographic sample, the contour points are divided into a mold side contour point set and a product side contour point set; A circle fitting is performed on the mold side contour point set, and a plurality of rays are constructed according to a preset polar angle interval, with the center of the fitted circle as the starting point of the rays and a preset reference direction as the starting direction of the zero polar angle; According to the two contour points hit by each ray in the mold side contour point set and the product side contour point set, the actual thickness of the circular arc transition region at each polar angle is calculated; According to a preset load value and the actual thickness of the circular arc transition region at each polar angle, the target load when the minimum value of the safety margin of the circular arc transition region is reduced to zero is iteratively calculated as the estimated interlayer limit load of the target composite material circular arc transition region interlayer strength test.

2. The method of claim 1, wherein, In the target metallographic image, each contour point located in the inner and outer contours of the circular arc transition region is identified, including: The target metallographic image is preprocessed by grayscale and filter denoising to obtain a preprocessed image, and the edge pixels in the target metallographic image are obtained by edge recognition on the preprocessed image; At least one closed contour chain is obtained by contour finding on the edge pixels, and two target contour chains representing the inner and outer contours of the circular arc transition region are selected from the contour chain as each contour point in the inner and outer contours of the circular arc transition region.

3. The method of claim 1, wherein, According to the forming process type of the composite material metallographic sample, the contour points are divided into a mold side contour point set and a product side contour point set, including: The forming process type of the target composite material metallographic sample is obtained from the forming process file; If the forming process type is negative mold forming, the inner contour point set of the circular arc transition region in the target metallographic image is defined as the mold side contour data point set, and the outer contour point set of the circular arc transition region in the target metallographic image is defined as the product side contour data point set; If the forming process type is positive mold forming, the outer contour point set of the circular arc transition region in the target metallographic image is defined as the mold side contour data point set, and the inner contour point set of the circular arc transition region in the target metallographic image is defined as the product side contour data point set.

4. The method of claim 1, wherein, The center of the fitted circle is taken as the starting point of the rays, and a plurality of rays are constructed according to a preset polar angle interval, with a preset reference direction as the starting direction of the zero polar angle, including: The center coordinates of the fitted circle are taken as the origin of the polar coordinate system, the positive direction of the x-axis of the polar coordinate system is defined as the starting direction of the zero polar angle, and the positive direction of the y-axis is defined as the end direction of the polar angle 90°, forming a polar angle interval covering the contour range of the circular arc transition region; Starting from the zero polar angle starting direction, a plurality of rays are generated in the preset polar angle interval along the increasing direction of the polar angle, and each ray corresponds to a unique polar angle value.

5. The method of claim 4, wherein, According to two profile points hit by each ray in the mold side profile point set and the product side profile point set respectively, the actual thickness of the circular arc transition area under each polar angle is calculated, including: The intersection of each ray with the mold side profile point set is defined as a first type intersection, and the intersection with the product side profile point set is defined as a second type intersection; The straight line distance between the first type intersection and the second type intersection under each ray is calculated; Each straight line distance is taken as the actual thickness of the circular arc transition area corresponding to each polar angle of each ray.

6. The method of claim 1, wherein, According to the preset load value and the actual thickness of the circular arc transition area under each polar angle, the target load is iteratively calculated when the minimum value of the safety margin of the circular arc transition area is reduced to zero, including: The preset load value is taken as the current load, and the current load and the actual thickness of the circular arc transition area under each polar angle are jointly input into the standard mechanical model to calculate the safety margin value corresponding to each polar angle to form a safety margin distribution set, and the current minimum safety margin value is found in the safety margin distribution set; Determine whether the current minimum safety margin value meets the margin value critical condition; If yes, the current load is determined as the target load; If not, the numerical iteration algorithm is used to update the current load, and the updated load value is taken as the new current load, and the operation of inputting into the standard mechanical model to calculate the safety margin value is repeatedly executed; Until the current minimum safety margin value meets the margin value critical condition, the current load corresponding to the current minimum safety margin value is determined as the target load.

7. The method according to any one of claims 1 to 6, characterized in that, The current load and the actual thickness of the circular arc transition area under each polar angle are jointly input into the standard mechanical model to calculate the safety margin value corresponding to each polar angle, including: based on the preset load and the target composite material circular arc transition zone interlayer strength test sample cross section width , respectively calculate the tensile force flow , shear flow and bending moment flow ; wherein, , , ; is the maximum actual size of the circular arc transition zone along the direction of the polar angle 90°; is the distance from the center of the fitting circle to the second type of intersection point in the direction of the zero polar angle; and fixed as 0; is retrieved from the forming process file or pre-manually measured; a theoretical fillet radius from the forming process file and an actual thickness at the current polar angle the characteristic radius R between the arc transition and the fitting circle is calculated using the formula ​ based on the extension flow , the shear flow , the bending moment flow and the characteristic radius R, the bending moment flow at polar angle is calculated and the shear flow at polar angle is calculated ; wherein ; the shear flow calculation formula is: ; Based on the bending moment flow and the shear flow Calculate the polar angles separately Interlaminar shear stress at the location and polar angle Principal stress in the thickness direction at the location The formula for calculating interlayer shear stress is as follows: The formula for calculating the principal stress in the thickness direction is: Based on the interlaminar shear stress and the principal stress in the thickness direction Calculate polar angle Safety margin value at The formula for calculating the safety margin value is as follows: ; To allow maximum shear stress, To allow for maximum normal stress, both are inherent property parameters of the target composite material recorded in the molding process document.

8. A device for predicting the inter-story ultimate load in a circular arc transition zone, characterized in that, The device comprises: An image recognition module is configured to acquire an image of a cross section of a target metallographic sample of a composite material containing a circular arc transition area feature to obtain a target metallographic image containing inner and outer profiles of the circular arc transition area; An edge profile point recognition module is configured to recognize each profile point located in the inner and outer profiles of the circular arc transition area in the target metallographic image, and divide the profile points into a mold side profile point set and a product side profile point set according to a forming process type of the metallographic sample of the composite material; A circle fitting module is configured to perform circle fitting on the mold side profile point set, take a center of the fitted circle as a starting point of a ray, take a preset reference direction as a starting direction of a zero polar angle, and construct a plurality of rays according to a preset polar angle interval; An actual thickness calculation module is configured to calculate the actual thickness of the circular arc transition area under each polar angle according to two profile points hit by each ray in the mold side profile point set and the product side profile point set respectively; An interlayer ultimate load estimation module is configured to iteratively calculate a target load when a minimum value of a safety margin of the circular arc transition area is reduced to zero according to a preset load value and the actual thickness of the circular arc transition area under each polar angle, and take the target load as an estimated interlayer ultimate load of a target composite material circular arc transition area interlayer strength test.

9. An electronic device, comprising: The electronic device comprises: At least one processor; and a memory connected with the at least one processor in communication; wherein the memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor to enable the at least one processor to perform the interlaminar ultimate load prediction method of the circular-arc transition zone according to any one of claims 1-7.

10. A computer program product, characterised in that, The computer program product comprises a computer program which, when executed by a processor, implements the interlaminar ultimate load prediction method of the circular-arc transition zone according to any one of claims 1-7.

Images