A method, device, medium and product for reconstructing a traction-separation curve
The traction-separation curve is reconstructed through the finite element model and gradient descent method, which solves the measurement problems of high cost and large errors in the prior art, and realizes the precise characterization of the interlayer damage mechanism of composite materials.
Patent Information
- Application Number
- CN202410444457.5
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2024-04-12
- Publication Date
- 2025-09-02
- Estimated Expiration
- 2044-04-12
AI Technical Summary
The prior art requires high-cost special equipment when measuring interlayer traction-separation curves of composite materials, and large test errors make it difficult to accurately characterize complex damage mechanisms. Especially in the case of interlayer plasticity and fiber bridges, the simplified curve form leads to poor characterization accuracy of fracture characteristics.
Using the finite element model and gradient descent method, the force-displacement curve of the double cantilever beam specimen is used as input to iteratively calculate the crack expansion amount, reconstruct the traction-separation curve, and combine the energy release rate equilibrium relationship to avoid physical measurement of the crack expansion amount and simplify the measurement process.
It significantly reduces measurement difficulty and cost, improves the accuracy of the traction-separation curve in complex situations, and enables precise characterization of the damage mechanism.
Smart Images

Figure CN118246287B_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the technical field of traction-separation curve analysis, and in particular to a method, device, medium and product for reconstructing a type I traction-separation curve of interlayer cohesion of a fiber reinforced composite material. Background Art
[0002] Composite materials and bonded structures are characterized by lamination, making them susceptible to delamination or debonding under load. The cohesive zone model (CZM) is a commonly used numerical analysis software for simulating delamination and debonding. The CZM requires the definition of a traction-separation curve to describe the relationship between the traction at the interlaminar interface and the crack opening.
[0003] The traction-separation curve can be defined by three key parameters: interlaminar strength, fracture toughness and the shape of the traction-separation curve, e.g. Figure 1 The figure shows the linear pull-separation curve, where the rising edge on the left is the elastic region and the falling edge on the right is the linear pull-separation curve. The elastic region on the left will be ignored in the subsequent analysis.
[0004] In order to obtain the traction-separation curve, it is necessary to obtain it through experimental testing methods. Under the action of type I load, a double cantilever beam tensile test (DCB) can usually be used to test the interlaminar strength and fracture toughness according to the ASTM D5528 / D5528M-21 method.
[0005] To further obtain the shape of the traction-separation curve, it is necessary to measure the crack opening in the experiment and obtain the traction-separation curve through the differential relationship between the J integral value and the crack opening:
[0006]
[0007] Where δ is the crack opening, σ is the interfacial traction, and J is the J-integral value in the cohesive region. δ can be obtained by using testing equipment such as extensometers and DIC, such as Figure 2 Shown are common pull-separation curve forms (including but not limited to) other than the linear pull-separation curve.
[0008] Since testing δ typically requires special testing equipment (such as extensometers and DIC), experimental measurement costs are high, a unified testing standard has yet to be established, and testing errors are generally large. Therefore, most traction-separation curves only use interlaminar strength and fracture toughness as the main material parameters, and the shape of the traction-separation curve is generally simplified by using linear, exponential, or trapezoidal shapes. However, in cases where the damage mechanism is more complex, such as interlaminar plasticity and fiber bridging, the use of simple and common traction-separation curve forms will lead to poor accuracy in characterizing interlaminar fracture characteristics, necessitating actual testing of the shape of the traction-separation curve. Summary of the Invention
[0009] The purpose of the present invention is to provide a method, device, medium and product for reconstructing the traction-separation curve. Only the force-displacement curve of a double cantilever beam specimen is required as input, and there is no need to measure the crack opening in the experiment. This significantly simplifies the measurement difficulty of the traction-separation curve. The reconstructed traction-separation curve can accurately characterize the damage mechanism under complex conditions.
[0010] To achieve the above object, the present invention provides the following solutions:
[0011] In a first aspect, the present invention provides a method for reconstructing a pull-separation curve, comprising:
[0012] Establishing a finite element model of a double cantilever beam specimen, wherein the material of the double cantilever beam specimen is a composite material;
[0013] generating an initial traction-separation curve, wherein the initial traction-separation curve is a correspondence between traction force and crack opening;
[0014] Inputting the initial traction-separation curve and the force-displacement curve into the finite element model to obtain the crack opening of the cohesive region corresponding to each displacement point, wherein the force-displacement curve is obtained based on load data and displacement data between load application points, the load data and the displacement data being obtained by testing the double cantilever beam specimen;
[0015] For each of the displacement points, perform the following steps:
[0016] Calculating the energy release rate at the crack front based on the displacement data and the corresponding load data;
[0017] Calculating the energy release rate of delamination expansion based on the crack opening and the corresponding traction force in the initial traction-separation curve;
[0018] calculating an error function based on the energy release rate of the crack front and the energy release rate of the delamination extension;
[0019] The gradient descent method is used to minimize the error function and obtain the traction force under the current crack opening;
[0020] updating the initial traction-separation curve according to the traction force under the current crack opening to obtain an updated traction-separation curve;
[0021] The updated traction-separation curve is used as a new initial traction-separation curve, and the process returns to step "inputting the initial traction-separation curve and the force-displacement curve into the finite element model to obtain the crack opening in the cohesive region corresponding to each displacement point" until the traction-separation curve approaches a stable state, thereby completing the reconstruction of the traction-separation curve.
[0022] Optionally, after executing the step of “establishing a finite element model of the double cantilever beam specimen”, the reconstruction method further includes:
[0023] The parameters of the finite element model are adjusted so that the force-displacement curve in the linear elastic stage calculated according to the finite element model is consistent with the force-displacement curve in the linear elastic stage obtained by testing the double cantilever beam specimen.
[0024] Optionally, the initial pull-separation curve is a linear pull-separation curve.
[0025] Optionally, the calculation expression of the error function is:
[0026]
[0027] Wherein, C is the error function, L j is the length of the j-th cohesive unit, δ i (j) is the crack opening of the jth cohesive unit in the cohesive region corresponding to the i-th pair of force-displacement data points, x represents the length direction of the cohesive unit, σ i (j) is the traction force calculated based on the traction-separation curve, P i is the load corresponding to the i-th pair of force-displacement data points, U i is the displacement corresponding to the i-th pair of force-displacement data points, E is the longitudinal Young's modulus of the material used in the double cantilever beam specimen, I is the section moment of inertia of the cantilever, and b is the cantilever width.
[0028] Optionally, the calculation expression of the traction force under the current crack opening is:
[0029]
[0030] in, The traction force obtained to minimize the error function is: is the traction force calculated based on the initial traction-separation curve, λ is the learning rate, C is the error function, σ i (j) is the traction force calculated based on the traction-separation curve.
[0031] In a second aspect, the present invention provides a computer device comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of any one of the above-mentioned methods for reconstructing a traction-separation curve.
[0032] In a third aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of any of the above-mentioned methods for reconstructing a traction-separation curve.
[0033] In a fourth aspect, the present invention provides a computer program product, comprising a computer program, which, when executed by a processor, implements the steps of any of the above-mentioned methods for reconstructing a traction-separation curve.
[0034] According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:
[0035] The present invention provides a method, device, medium, and product for reconstructing traction-separation curves. This method uses a finite element model to iteratively update the crack opening within the cohesive region and, combined with the J-integral, calculates the energy release rate of delamination extension. Furthermore, a dual-cantilever beam theory formula is used to calculate the energy release rate at the crack front, establishing an equilibrium relationship based on the energy release rate based on the Dugdale condition. Finally, a gradient descent algorithm is used to reconstruct the traction-separation curve. This method utilizes only DCB force-displacement data as input, eliminating the need for physical measurement of the crack opening displacement at the delamination tip. This significantly simplifies the measurement of the traction-separation curve, and the reconstructed traction-separation curve can accurately characterize damage mechanisms in complex situations. BRIEF DESCRIPTION OF THE DRAWINGS
[0036] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.
[0037] Figure 1 is a schematic diagram of the linear pull-separation curve;
[0038] Figure 2 Schematic diagram of a common nonlinear pull-separation curve;
[0039] Figure 3 A schematic diagram of a process for reconstructing a pull-separation curve provided in Example 1 of the present invention;
[0040] Figure 4 Schematic diagram of a typical DCB experiment in Example 1 of the present invention;
[0041] Figure 5 Schematic diagram of force-displacement curves obtained from a typical DCB test under three different traction-separation curves in Example 1 of the present invention;
[0042] Figure 6 Schematic diagram of the finite element model in Example 1 of the present invention;
[0043] Figure 7 Schematic diagram of extracting crack opening amount based on finite element model in Example 1 of the present invention;
[0044] Figure 8 This is a diagram of the internal structure of a computer device. DETAILED DESCRIPTION
[0045] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.
[0046] The purpose of the present invention is to provide a method, device, medium and product for reconstructing the traction-separation curve. Only the force-displacement curve of a double cantilever beam specimen is required as input, and there is no need to measure the crack opening in the experiment. This significantly simplifies the measurement difficulty of the traction-separation curve. The reconstructed traction-separation curve can accurately characterize the damage mechanism under complex conditions.
[0047] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.
[0048] Example 1
[0049] like Figure 3 As shown, a method for reconstructing a pull-separation curve in this embodiment includes:
[0050] S1: Establishing a finite element model of a double cantilever beam specimen, wherein the material of the double cantilever beam specimen is a composite material.
[0051] S2: generating an initial traction-separation curve, wherein the initial traction-separation curve is a correspondence between traction force and crack opening.
[0052] S3: Input the initial traction-separation curve and the force-displacement curve into the finite element model to obtain the crack opening amount of the cohesive force area corresponding to each displacement point, wherein the force-displacement curve is obtained based on the load data and the displacement data between the load application points, and the load data and the displacement data are obtained by testing the double cantilever beam specimen.
[0053] For each of the displacement points, perform the following steps:
[0054] S4: Calculating the energy release rate of the crack front according to the displacement data and the corresponding load data.
[0055] S5: Calculating the energy release rate of delamination expansion according to the crack opening amount and the corresponding traction force in the initial traction-separation curve.
[0056] S6: Calculating an error function based on the energy release rate of the crack front and the energy release rate of the delamination extension.
[0057] S7 uses the gradient descent method to minimize the error function and obtain the traction force under the current crack opening.
[0058] S8: updating the initial traction-separation curve according to the traction force under the current crack opening amount to obtain an updated traction-separation curve.
[0059] S9: Using the updated traction-separation curve as a new initial traction-separation curve, and returning to S3 until the traction-separation curve approaches a stable state, thereby completing the reconstruction of the traction-separation curve.
[0060] The above-mentioned reconstruction method proposed in this embodiment combines the finite element model and the DCB (double cantilever beam) test to obtain the shape of the traction-separation curve. There is no need to use any special device to measure the crack opening during the DCB test. It is only necessary to measure the loading point displacement and load data of the traditional DCB. The finite element model is used to iteratively calculate the crack opening in the delamination (debonding) area, and the traction-separation curve is obtained by numerical calculation. The basic idea of this embodiment is that the correct traction-separation curve and the crack opening in the delamination (debonding) area need to satisfy the load balance relationship and the energy release rate balance relationship of the DCB structure at the same time. Therefore, it is possible to: 1) establish a crack opening update algorithm (finite element model) based on the load balance relationship; 2) establish a traction-separation curve update algorithm based on the energy release rate balance relationship; through the sequential iterative solution of these two balance relationships, it finally converges to a traction-separation curve that satisfies both balance relationships. The specific calculation process is as follows:
[0061] 1) Carry out DCB test according to ASTM standards, obtain the displacement and load data of the loading point, and then obtain the force-displacement curve. Then calculate the energy release rate of delamination (debonding) under external load according to formula (2), that is, the energy release rate at the crack front mentioned in S4 above:
[0062]
[0063] Where P is the force at the loading point (i.e., the point where the load is applied), U is the displacement of the loading point, E is the longitudinal Young's modulus of the cantilever, b is the cantilever width, and I is the polar moment of inertia of the cantilever.
[0064] like Figure 4 As shown, DCB material is a composite material, and the double cantilever beam is prefabricated in the middle of the cantilever beam with a certain length of layer. Figure 4 The length indicated by a in the figure is the length of the prefabricated layer. E is the longitudinal Young's modulus of the cantilever beam along the beam direction. B is the width of the cantilever beam. Figure 4 It is a two-dimensional cross-sectional view, so b is the width of the cantilever beam perpendicular to the screen direction. I = bh 3 / 12 is the polar moment of inertia of the cantilever beam, h is the thickness of a single cantilever, and 2h is the thickness of the entire cantilever. After the double cantilever beam specimen is clamped by a tensile machine, opposite force loads P are applied to the two cantilevers, and the displacement U between the two load application points is recorded. The force-displacement response curve is as follows: Figure 5 shown.
[0065] In the above DCB test, as the load increases, the precast delamination will further expand, i.e. Figure 4 The a in the equation will become larger. Before the crack propagates, the response of the structure is basically linear. Figure 5 As shown in the left area, the structural response changes when the layers are expanded, as shown in Figure 5 The right area is shown.
[0066] According to the assumptions of fracture mechanics, the energy release rate is the energy released when a crack forms. When a crack forms under an external load, the energy release rate of the crack under the external load can be calculated based on beam theory, i.e., formula (2). At this time, if the stress field near the crack is calculated by J-integral, the energy release rate near the crack can be obtained, i.e., the energy release rate of the delamination expansion mentioned in S7 above. This energy release rate must be equal to the value calculated by formula (2). Therefore, the energy release rate is the basis for constructing formula (3) below in this embodiment.
[0067] 2) See Figure 6 , a finite element model of DCB was established.
[0068] A finite element model of DCB is constructed based on the geometric parameters of DCB, such as length, width and height, as well as the Young's modulus, Poisson's ratio and shear modulus of the material used in DCB.
[0069] The constructed finite element model needs to ensure that the linear elastic stage of the force-displacement curve is consistent with that obtained from the DCB experiment when delamination (debonding) does not occur.
[0070] Because the linear elastic stage occurs when the load is low, delamination does not expand, and the force-displacement (PU) curve at the loading point is essentially linear. Due to differences between the finite element model and the actual physical structure, the force-displacement curve calculated based on the finite element model may not be completely consistent with the experimental test. In this case, the finite element parameters (primarily Young's modulus) need to be adjusted to achieve a basic consistency between the two.
[0071] 3) generating an initial traction-separation curve, which can be generated based on a linear relationship between interlaminar strength and fracture toughness obtained using ASTM standards.
[0072] The cohesive force model is established based on the traction-separation curve, and the force-displacement curve is input into the finite element model. According to the displacement data points of the test data, the crack opening in the corresponding cohesive force area is extracted from the finite element model in sequence, such as Figure 7 As shown in the figure, the small figure in the lower right corner is the extracted crack opening δ, the horizontal axis is the x coordinate in the beam direction, the vertical axis is the δ value, and the discrete expression is an array of δ. i In this case, multiple groups of δ can be obtained i .
[0073] Extract the crack opening in the corresponding cohesive region from the finite element model, including:
[0074] For a pair of force-displacement data (P i -U i ), the crack opening δ in the cohesive region can be obtained from the finite element model force i δ i It represents the crack opening of all cohesive units from the tip of the delamination (debonding) area to the cohesive damage initiation area corresponding to the i-th pair of PU data points. The constitutive model of each cohesive unit is the cohesive model. i The δ of any cohesive unit in i (j) , you can get a corresponding traction value σ i (j) , and the traction force value needs to satisfy the following energy release rate balance relationship:
[0075]
[0076] Wherein, C is the error function, L j is the length of the j-th cohesive unit, δ i (j) is the crack opening of the jth cohesive element in the cohesive region corresponding to the i-th pair of force-displacement data points, x represents the length direction of the cohesive element, that is, the direction of the beam, σ i (j) To convert δ i (j) Substitute the traction force calculated by the traction-separation curve, P i is the load corresponding to the i-th pair of force-displacement data points, U i is the displacement corresponding to the i-th pair of force-displacement data points.
[0077] The DCB experiment records PU data. The subscript i represents the discretized PU data. Therefore, a pair of PU data is expressed as P i -U i The corresponding U can also be extracted in the finite element model. i The crack opening δ in the delamination area is i (Note that this is a vector), the superscript j represents δ i Discrete, where any element δ i (j) According to the traction-separation curve, a σ can be calculated. i (j) , substituted into formula (3). However, formula (3) is valid only when the pull-separation curve is correct; otherwise, the result of formula (3) is not zero. Therefore, the error function of formula (4) is established, and the pull-separation curve is updated by the gradient descent method of formula (5) so that the result of formula (3) is close to 0, thereby achieving the reconstruction of the pull-separation curve.
[0078] 4) Establish the error function shown in formula (4):
[0079]
[0080] 5) Use gradient descent to minimize L, where λ is the learning rate, which is usually a small amount, to obtain the pull-separation curve under the current δ:
[0081]
[0082] Because the crack opening δ in the cohesive region varies for different traction-separation curves, the first update uses the δ obtained from the initial traction-separation curve, and then calculates a new traction-separation curve based on this δ value. When this new traction-separation curve is substituted into the finite element model, the resulting δ will be different from the first. Only through multiple iterations can δ converge to a certain steady-state value, at which point the traction-separation curve will also converge, and the final reconstruction of the traction-separation curve is completed.
[0083] Therefore, the traction-separation curve updated according to formula (5) is substituted into the finite element model to update the crack opening δ in the cohesive region. Repeat steps (4) to (6) and substitute δ into formulas (4) and (5) to update the traction-separation curve. After repeating this iterative process N times, the traction-separation curve will converge to a steady state, which is the final traction-separation curve.
[0084] This embodiment uses a finite element model to numerically calculate the crack opening in the cohesive region. Only the force-displacement curve of the traditional DCB test is required as input, avoiding the need for physical testing of the crack opening in the experiment. The traction-separation curve is reconstructed by iteratively calculating δ using the finite element model. No special equipment such as DIC or extensometer is required to measure δ, significantly simplifying the traction-separation curve testing process, reducing experimental costs, and improving accuracy.
[0085] Example 2
[0086] This embodiment provides a computer device, including: a memory, a processor, and a computer program stored in the memory and executable on the processor. The processor executes the computer program to implement the steps of the method for reconstructing a pull-separation curve in embodiment 1.
[0087] Example 3
[0088] This embodiment provides a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the steps of the method for reconstructing a pull-separation curve in embodiment 1 are implemented.
[0089] Example 4
[0090] This embodiment provides a computer program product, including a computer program. When the computer program is executed by a processor, the steps of the method for reconstructing a pull-separation curve in embodiment 1 are implemented.
[0091] Example 5
[0092] This embodiment proposes a computer device, which may be a database, and its internal structure diagram may be as follows: Figure 8As shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, abbreviated as I / O) and a communication interface. The processor, memory and input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and computer program in the non-volatile storage medium. The database of the computer device is used to store pending transactions. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, a method for reconstructing a traction-separation curve in Example 1 is implemented.
[0093] It should be noted that the object information (including but not limited to object device information, object personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved in the present invention are all information and data authorized by the object or fully authorized by all parties, and the collection, use and processing of relevant data must comply with the relevant laws, regulations and standards of relevant countries and regions.
[0094] Those skilled in the art will appreciate that all or part of the processes in the above-mentioned embodiments can be implemented by guiding the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. In particular, any reference to memory, database or other media used in the embodiments provided by the present invention can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM). The databases involved in the various embodiments provided by the present invention may include at least one of a relational database and a non-relational database. Non-relational databases may include, but are not limited to, distributed databases based on blockchains. The processors involved in the various embodiments provided by the present invention may be, but are not limited to, general-purpose processors, central processing units (CPUs), graphics processing units (GPUs), digital signal processors (DSPs), programmable logic devices (PLDs), data processing logic devices based on quantum computing, and the like.
[0095] The technical features of the above embodiments can be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.
[0096] This document uses specific examples to illustrate the principles and implementation methods of the present invention. The above examples are only intended to help understand the method and core concept of the present invention. At the same time, those skilled in the art will find that the specific implementation methods and application scopes may vary based on the concept of the present invention. In summary, the contents of this specification should not be construed as limiting the present invention.
Claims
1. A method for reconstructing a pull-separation curve, characterized in that: The reconstruction method includes: Establishing a finite element model of a double cantilever beam specimen, wherein the material of the double cantilever beam specimen is a composite material; generating an initial traction-separation curve, wherein the initial traction-separation curve is a correspondence between traction force and crack opening; Inputting the initial traction-separation curve and the force-displacement curve into the finite element model to obtain the crack opening of the cohesive region corresponding to each displacement point, wherein the force-displacement curve is obtained based on load data and displacement data between load application points, and the load data and the displacement data are obtained by testing the double cantilever beam specimen; For each of the displacement points, perform the following steps: Calculating the energy release rate at the crack front based on the displacement data and the corresponding load data; Calculating the energy release rate of delamination expansion based on the crack opening and the corresponding traction force in the initial traction-separation curve; calculating an error function based on the energy release rate of the crack front and the energy release rate of the delamination extension; The gradient descent method is used to minimize the error function and obtain the traction force under the current crack opening; updating the initial traction-separation curve according to the traction force under the current crack opening to obtain an updated traction-separation curve; The updated traction-separation curve is used as a new initial traction-separation curve, and the process returns to step "inputting the initial traction-separation curve and the force-displacement curve into the finite element model to obtain the crack opening in the cohesive region corresponding to each displacement point" until the traction-separation curve reaches a stable state, thereby completing the reconstruction of the traction-separation curve.
2. The method for reconstructing a pull-separation curve according to claim 1, characterized in that: After executing the step of "establishing a finite element model of the double cantilever beam specimen", the reconstruction method further includes: The parameters of the finite element model are adjusted so that the force-displacement curve in the linear elastic stage calculated according to the finite element model is consistent with the force-displacement curve in the linear elastic stage obtained by testing the double cantilever beam specimen.
3. The method for reconstructing a pull-separation curve according to claim 1, characterized in that: The initial pull-separation curve is a linear pull-separation curve.
4. The method for reconstructing a pull-separation curve according to claim 1, characterized in that: The calculation expression of the error function is: Wherein, C is the error function, L j is the length of j cohesive units, is the crack opening of the jth cohesive unit included in the cohesive region corresponding to the i-th pair of force-displacement data points, x represents the length direction of the cohesive unit, is the traction force calculated based on the traction-separation curve, P i is the load corresponding to the i-th pair of force-displacement data points, U i is the displacement corresponding to the i-th pair of force-displacement data points, E is the longitudinal Young's modulus of the material used in the double cantilever beam specimen, I is the section moment of inertia of the cantilever, and b is the cantilever width.
5. The method for reconstructing a pull-separation curve according to claim 1, characterized in that: The calculation expression of the traction force under the current crack opening is: in, The traction force obtained to minimize the error function is: is the traction force calculated based on the initial traction-separation curve, λ is the learning rate, C is the error function, is the traction force calculated based on the traction-separation curve.
6. A computer device comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of the method for reconstructing a pull-separation curve according to any one of claims 1 to 5.
7. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the computer program implements the steps of the method for reconstructing a pull-separation curve according to any one of claims 1 to 5.
8. A computer program product comprising a computer program, characterized in that When the computer program is executed by a processor, the computer program implements the steps of the method for reconstructing a pull-separation curve according to any one of claims 1 to 5.
Citation Information
Patent Citations
Method and device for constructing three-dimensional expansion model of deep-water high-temperature high-pressure stratum fracturing crack
CN115017769A
Method and system for calculating fracture toughness of ultra-high performance concrete reinforced beam
CN117789869A