An Optimization Design Method for Laminated Glass Structures Based on Near-Field Dynamics Simulation

Abstract

This invention discloses a method for optimizing the design of laminated glass structures based on near-field dynamics simulation, comprising: Step 1, obtaining the expected load on the laminated glass based on its intended installation location; Step 2, initializing the structural dimensions and material parameters of the laminated glass; Step 3, inputting the expected load, structural dimensions, and material parameters of the laminated glass into a near-field dynamics-based simulation model of the laminated glass to obtain the failure conditions of the laminated glass under the expected load; Step 4, determining whether the laminated glass meets safety specifications based on the failure conditions. If it does, obtaining a structural design scheme for the laminated glass under the expected load based on its structural dimensions and material parameters. If it does not meet the specifications, adjusting the structural dimensions and material parameters of the laminated glass within a preset range, and returning to Step 3. This invention can more accurately calculate the response of laminated glass under external loads, resulting in better and more accurate optimization.

Description

Technical Field

[0001] This invention relates to an optimization design method for laminated glass structures based on near-field dynamics simulation, belonging to the technical field of laminated glass structure optimization design. Background Technology

[0002] Laminated glass typically consists of two or more panes of glass bonded together with one or more layers of polymer interlayers. Compared to ordinary glass, the interlayers in laminated glass can absorb the energy of external loads; the polymer interlayers are less prone to breakage than glass and possess a certain degree of strength; in the event of glass breakage, the interlayers can hold the glass together to prevent shards from flying and causing injury. Therefore, laminated glass offers higher safety and strength. These superior properties make laminated glass widely used in various structures. However, the performance of laminated glass is greatly affected by parameters such as the thickness and type of the glass and interlayer. Therefore, optimizing laminated glass structures based on the expected installation location and expected load is both economically and safety-wise significant. The prerequisite for optimizing laminated glass structures is accurately calculating the response of the laminated glass under external loads, especially accurately predicting the failure behavior of laminated glass under external loads.

[0003] The failure modes of laminated glass under load are complex, including glass breakage, debonding at the interface between the laminate and the glass, and failure of the laminate material. Testing the failure behavior of laminated glass experimentally is costly, especially for high-speed impact and explosion tests. With the development of computer technology and the increase in computing speed, numerical methods have become one of the main methods for studying laminated glass failure. Currently, various numerical methods have been applied to the study of laminated glass failure, such as the Finite Element Method (FEM), Finite Discrete Element Method (FDM), Extended Finite Difference Method (FDM), and Extended Finite Element Method (XFEM). However, these methods are based on the framework of continuum mechanics, which presents certain difficulties in simulating discontinuous problems such as failure. Furthermore, mesh-based numerical methods such as the FEM and FDM suffer from significant mesh distortion when solving material failure problems, especially large deformation problems such as high-speed impact and explosion, making it highly dependent on mesh quality and difficult to guarantee the accuracy of experimental results. Summary of the Invention

[0004] To address the limitations and inaccurate results of existing laminated glass structure optimization design methods, this invention proposes a laminated glass structure optimization design method based on near-field dynamics simulation. By simulating the failure of laminated glass using near-field dynamics, the method can more accurately calculate the response of laminated glass under external loads, resulting in better and more accurate optimization.

[0005] To solve the above-mentioned technical problems, the present invention adopts the following technical means:

[0006] This invention proposes an optimization design method for laminated glass structures based on near-field dynamics simulation, comprising the following steps:

[0007] Step 1: Obtain the expected load on the laminated glass based on its intended installation location;

[0008] Step 2: Initialize the structural dimensions and material parameters of the laminated glass;

[0009] Step 3: Input the expected load, structural dimensions and material parameters of the laminated glass into the laminated glass simulation model based on near-field dynamics to obtain the failure condition of the laminated glass under the expected load.

[0010] Step 4: Determine whether the laminated glass meets safety standards based on the failure condition. If it does, obtain the laminated glass structural design scheme under the expected load based on the structural dimensions and material parameters of the laminated glass. If it does not meet the standards, adjust the structural dimensions and material parameters of the laminated glass within the preset range and return to step 3.

[0011] Furthermore, the material parameters include the density, elastic modulus, and Poisson's ratio of the glass in the laminated glass, and the density, elastic modulus, and Poisson's ratio of the interlayer.

[0012] Furthermore, the method for obtaining the failure condition of laminated glass under expected load is as follows:

[0013] Step 301: Input the expected load, structural dimensions and material parameters of the laminated glass into the laminated glass simulation model based on near-field dynamics;

[0014] Step 302: Discretize each part of the laminated glass into near-field dynamic material points and set the calculation time step;

[0015] Step 303: Determine the boundary conditions of the material points based on the expected load;

[0016] Step 304: Search for nearby material points based on the near-field range to obtain a list of near-field ranges;

[0017] Step 305: Obtain the effective range of the interface model based on the interface position of the laminated glass;

[0018] Step 306: Calculate the interaction forces between material points within the same material of the laminated glass based on the near-field dynamics theory and the near-field range list;

[0019] Step 307: Based on the interface model and its range of action, calculate the interfacial forces between material points inside different materials on both sides of the interface of the laminated glass.

[0020] Step 308: Calculate the acceleration of the material point based on the interaction forces and interface forces, and calculate the displacement and velocity of the material point based on the acceleration;

[0021] Step 309: Calculate the damage to the material points using a damage model based on their displacement and velocity.

[0022] Step 310: Determine whether the current time t exceeds the preset simulation duration t. max If t≤t max If the calculation is successful, proceed to the next time step and return to step 306; otherwise, end the calculation. Obtain the overall damage cloud map of the laminated glass based on the material point damage, which represents the failure status of the laminated glass under the expected load.

[0023] Furthermore, the specific operation of step 304 is as follows:

[0024] Randomly select a material point x, and search for material points within its near-field range δ = 3Δx with material point x as the center, generating a list of near-field ranges, where Δx is the discrete spacing.

[0025] Furthermore, in step 305, the interface model's interface scope is:

[0026] δ B = (Δx1 + Δx2) / 2 (1)

[0027] Where, δ B The effective range of the interface model is represented by Δx1 and Δx2, which are the discrete distances on both sides of the laminated glass interface, respectively.

[0028] Furthermore, the condition for considering interfacial interactions between matter point x and matter point x′ is as follows:

[0029] (a) The initial distance between matter point x and matter point x′ is less than the effective range δ of the interface model. B ;

[0030] (b) Material point x and material point x′ are located on opposite sides of the interface of the laminated glass, meaning that material point x and material point x′ have different material properties.

[0031] Furthermore, the acceleration of matter point x under the interaction forces and interfacial forces is:

[0032]

[0033] Where ρ(x) represents the density of the substance point x. H represents the acceleration of point x at time t. x B represents the set of matter points x′ within the near-field range of matter point x. x The interface model represents the set of matter points x′ within the interface range of matter point x. T [x,t] represents the interaction force experienced by a point x within the same material at time t. T [x′,t] represents the interaction force experienced by the material point x′ within the same material at time t, V x′ This represents the volume of the substance point x′. F [x,t] represents the interfacial forces experienced by a point x at time t within the different materials on either side of the interface. F [x′,t] represents the interfacial force experienced by the material point x′ at time t in the different materials on both sides of the interface, and b(x,t) represents the body force density of the material point x.

[0034] Furthermore, in step 308, the displacement and velocity of the material point are calculated using the Verlet explicit integration method, with the following iteration method:

[0035]

[0036] Among them, u n This represents the displacement of the material point at time step n. This represents the velocity of the substance point at time step n. Δt represents the acceleration of the material point at time step n, where Δt is the time step size.

[0037] Furthermore, when material point x interacts with material point x′ within its near-field range, the failure criterion for material point x is that the elongation s is greater than the critical elongation s0. By statistically analyzing the bond breaks of all material points x′ within the near-field range of material point x, the damage to material point x is obtained:

[0038]

[0039]

[0040]

[0041]

[0042] Where y and y′ are the position vectors of matter x and x′, respectively, G0 is the energy release rate of the material to which matter point x belongs, K is the bulk modulus of the material to which matter point x belongs, δ is the near-field range of matter point x, and μ(x,x′) is the discontinuity function of matter point x and other matter points x′ within its near-field range. H represents the damage to material point x at time t. x V represents the set of matter points x′ within the near field. x′ This represents the volume of the substance point x′.

[0043] Furthermore, the failure scenarios of laminated glass under expected loads include glass breakage, debonding of the glass-laminate interface, and lamination failure.

[0044] The following advantages can be obtained by adopting the above technical means:

[0045] This invention proposes a near-field dynamics simulation-based optimization design method for laminated glass structures. Under the expected installation location and load, a near-field dynamics-based simulation model simulates the breakage of laminated glass with different structural dimensions and material parameters, thereby continuously optimizing the laminated glass structure and achieving more accurate and reliable optimization results. This method utilizes near-field dynamics technology, avoiding singularity issues when solving material discontinuity problems, thus more accurately calculating the response of laminated glass under external loads. Furthermore, a near-field dynamics interface model describes the mechanical behavior of the interface between the glass and the interlayer in the laminated glass. Through the interaction between the same and different materials, the stress situation of the laminated glass is better simulated, making the calculation results more accurate. In addition, this invention can represent three failure modes in laminated glass: glass failure, interlayer failure, and interface debonding, reducing costs and preventing failures, thus broadening its application scope and providing results that are more consistent with reality. Attached Figure Description

[0046] Figure 1 This is a flowchart illustrating the steps of an optimized design method for laminated glass structures based on near-field dynamics simulation, as described in this invention.

[0047] Figure 2 This is a simulation flowchart of the laminated glass simulation model based on near-field dynamics of the present invention;

[0048] Figure 3 This is a discrete schematic diagram of the various parts of the laminated glass in this invention;

[0049] Figure 4 This is a schematic diagram illustrating the scope of the interface model in this invention;

[0050] Figure 5 This is a schematic diagram of the crack tensile test in Embodiment 1 of the present invention;

[0051] Figure 6 This is a schematic diagram of the simulation results and experimental results of the crack tensile test in Embodiment 1 of the present invention;

[0052] Figure 7 This is a schematic diagram of the impact test of laminated glass in Embodiment 2 of the present invention;

[0053] Figure 8 This is a simulation result diagram of the laminated glass simulation model in the impact test of laminated glass in Embodiment 2 of the present invention;

[0054] Figure 9 This is a damage diagram from an actual test in the impact test of laminated glass in Embodiment 2 of the present invention;

[0055] Figure 10 This is the velocity-time curve of the impactor in Embodiment 2 of the present invention; Detailed Implementation

[0056] The technical solution of the present invention will be further described below with reference to the accompanying drawings:

[0057] This invention proposes an optimization design method for laminated glass structures based on near-field dynamics simulation, such as... Figure 1 As shown, the specific steps include the following:

[0058] Step 1: Determine the expected load on the laminated glass based on its intended installation location. For example, if the laminated glass is used as architectural glass and is intended to be installed on the exterior wall of a building, its expected load includes wind loads such as wind debris and bird strikes. If the laminated glass is used as aquarium glass and is intended to be installed in the underwater tunnel of an aquarium, its expected load includes water loads. The installation location will result in different loads and different directions of load application.

[0059] Step 2: Initialize the structural dimensions and material parameters of the laminated glass. The structural dimensions include the number of laminates, length and width of the laminated glass, glass thickness, and interlayer thickness. The material parameters include the density, elastic modulus, and Poisson's ratio of the glass and the density, elastic modulus, and Poisson's ratio of the interlayer.

[0060] Step 3: Input the expected load, structural dimensions and material parameters of the laminated glass into the laminated glass simulation model based on near-field dynamics to obtain the failure condition of the laminated glass under the expected load.

[0061] like Figure 2 As shown, the specific steps for step 3 are as follows:

[0062] Step 301: Input the expected load, structural dimensions and material parameters of the laminated glass into the laminated glass simulation model based on near-field dynamics.

[0063] Step 302: Discretize the various parts of the laminated glass (glass and interlayer) into a series of near-field dynamic material points containing material information, such as... Figure 3As shown, the calculation time step is set. During discretization, the model is divided into a series of hexahedral elements, requiring the element division to be as uniform in size as possible. The maximum discretization spacing Δx should not exceed 1 / 100 of the maximum side length of the solid geometric model, and the specific adjustment should be made according to the calculation accuracy requirements.

[0064] Step 303: Determine the boundary conditions of the material point based on the expected load. Determine the location where the external load is applied based on the expected load, mark the location of the external load application, and convert it into corresponding boundary conditions, such as displacement boundary, force boundary, and velocity boundary. Displacement and velocity boundaries can be directly applied to the material point, while force boundaries need to be converted into force density before being applied to the material point.

[0065] Step 304: Search for neighboring material points based on the near-field range to obtain a near-field range list. Since each material point only interacts with other material points within its near-field range, it is necessary to search for neighboring material points and establish a near-field range list. Specifically: randomly select a material point x, and with material point x as the center, search for material points within its near-field range δ = 3Δx to generate a near-field range list.

[0066] Step 305: Obtain the effective range of the interface model based on the interface location of the laminated glass. For example... Figure 4 As shown, in laminated glass, the surface between the glass and the polymer interlayer is the interface, and the effective range of the near-field dynamic interface model is:

[0067] δ B =(Δx1+Δx2) / 2 (8)

[0068] Where, δ B The interface model represents the scope of action of the interface model. Δx1 and Δx2 are the discrete distances on both sides of the laminated glass interface, i.e., the discrete distances between the glass and the interlayer.

[0069] In this embodiment of the invention, the conditions for considering interfacial interactions between material point x and material point x′ are as follows:

[0070] (a) The initial distance between matter point x and matter point x′ is less than the effective range δ of the interface model. B .

[0071] (b) Material point x and material point x′ are located on opposite sides of the interface of the laminated glass, meaning that material point x and material point x′ have different material properties.

[0072] Step 306: Calculate the interaction forces between material points within the same material in the laminated glass, based on the near-field dynamics theory and the near-field range list.

[0073] Step 307: Based on the interface model and its range of action, calculate the interfacial forces between material points inside different materials on both sides of the interface of the laminated glass.

[0074] Step 308: Calculate the acceleration of the material point based on the interaction force and the interface force, and calculate the displacement and velocity of the material point based on the acceleration.

[0075] The acceleration of point x under the interaction forces and interfacial forces is:

[0076]

[0077] Where ρ(x) represents the density of the substance point x. H represents the acceleration of point x at time t. x B represents the set of matter points x′ within the near-field range of matter point x. x The interface model represents the set of matter points x′ within the interface range of matter point x. T [x,t] represents the interaction force experienced by a point x within the same material at time t. T [x′,t] represents the interaction force experienced by the material point x′ within the same material at time t, V x′ This represents the volume of the substance point x′. F [x,t] represents the interfacial forces experienced by a point x at time t within the different materials on either side of the interface. F [x′,t] represents the interfacial force experienced by the material point x′ at time t in the different materials on both sides of the interface, and b(x,t) represents the body force density of the material point x.

[0078] In predicting the failure of laminated glass, a linear elastic model is used for the glass and a viscoelastic model is used for the interlayer.

[0079] The displacement and velocity of a material point are calculated using the Verlet explicit integration method, and the iterative method is as follows:

[0080]

[0081] Among them, u n This represents the displacement of the material point at time step n. This represents the velocity of the substance point at time step n. Δt represents the acceleration of the material point at time step n, where Δt is the time step size.

[0082] Step 309: Calculate the damage to the material points using the damage model based on the displacement and velocity of the material points.

[0083] Based on the displacement and velocity of the material points, a damage model is used to simulate the motion of each material point in the laminated glass under the expected load, and then the damage of the material points is calculated.

[0084] When a material point x interacts with another material point x′ within its near-field range, the failure criterion for material point x is that its elongation s is greater than the critical elongation s0. By statistically analyzing the bond breaks of all material points x′ within the near-field range of material point x, the damage to material point x is obtained:

[0085]

[0086]

[0087]

[0088]

[0089] Where y and y′ are the position vectors of matter x and x′, respectively, G0 is the energy release rate of the material to which matter point x belongs, K is the bulk modulus of the material to which matter point x belongs, δ is the radius of the near-field range of matter point x, and μ(x,x′) is the discontinuity function of matter point x and other matter points x′ within its near-field range. H represents the damage to material point x at time t. x V represents the set of matter points x′ within the near field. x′ This represents the volume of the substance point x′.

[0090] Step 310: Determine whether the current time t exceeds the preset simulation duration t. max If t≤t max If the calculation is successful, proceed to the next time step and return to step 306; otherwise, end the calculation. Obtain the overall damage cloud map of the laminated glass based on the damage at each material point, which represents the failure situation under the expected load.

[0091] Failure scenarios of laminated glass under expected loads include glass breakage, debonding of the glass-laminate interface, and lamination failure.

[0092] Step 4: Determine whether the laminated glass meets safety standards based on the failure condition. If it does, obtain the laminated glass structural design scheme under the expected load based on the structural dimensions and material parameters of the laminated glass. If it does not meet the standards, adjust the structural dimensions and material parameters of the laminated glass within a preset range, which is the range required by the use conditions of the laminated glass, and return to step 3.

[0093] To verify the effectiveness of the method of the present invention, the following two embodiments are given:

[0094] Example 1:

[0095] Example 1 is as follows Figure 5 The tensile test of the laminated glass plate shown is used to verify the applicability of the interface model of the present invention.

[0096] In Example 1, the laminated glass panel is 150mm long and 60mm wide, consisting of two 3mm thick glass panels and a 1.52mm thick PVB interlayer. A pre-installed crack is present between the two glass panels. The glass density is 2700kg / m³. 3 The elastic modulus is 70 GPa, and the Poisson's ratio is 0.20. The density of the PVB interlayer is 1100 kg / m³. 3 The elastic modulus is 0.1 GPa and the Poisson's ratio is 0.49. In the experiment, the bottom of the laminated glass plate was fixed and the top was loaded.

[0097] Analyzing the failure characteristics of laminated glass using the method of this invention:

[0098] (1) Input the above material parameters into the laminated glass simulation model.

[0099] (2) Discretize each part of the laminated glass according to the experimental dimensions. Considering the symmetry of the model, half of it is used for calculation (120×30mm). The glass discretization interval is taken as Δx1=0.60mm, the PVB discretization interval is taken as Δx2=0.15mm, and the time step is 8×10. -8 s, and take the calculation time as t. max =10 -2 s.

[0100] (3) Boundary conditions: The top is subjected to displacement loading at a speed of 1 m / s.

[0101] (4) Search for material points in the near field range of each part. The near field ranges of glass and PVB are δ1=3Δx1=1.80mm and δ2=4Δx2=0.60mm, respectively.

[0102] (5) Scope of application of the interface model δ B =0.42mm, determine the material points that need to be applied to the interface model.

[0103] (6) Solve for the interaction forces of the material points in the laminated glass according to the near-field dynamic constitutive equation and interface model, and solve for the displacement and velocity according to the Verlet explicit integration method.

[0104] (7) According to the near-field dynamic damage model, solve for the damage of material points and interface damage in each part of the laminated glass.

[0105] (8) Repeat (6)-(7) until the calculation time t≥t max The calculation will terminate at that time.

[0106] In Example 1, the deformation and damage (failure) of the laminated glass under this load are as follows: Figure 6As shown, the left side represents the result of a real test with a 10mm loading, and the right side represents the result predicted by the method of this invention. The diagram illustrates the central debonding region. From Figure 6 It is clearly visible that debonding occurred between the PVB and the glass in the middle of the laminated glass plate, indicating failure of the laminated glass interface. The width of the debonded area obtained by this invention is 12.9 mm, which is basically consistent with the actual experimental result of 13.1 mm.

[0107] Example 2:

[0108] Example 2 is as follows Figure 7 The impact test of the laminated glass shown is used to verify that the present invention can more accurately calculate the response of laminated glass under external loads. The laminated glass panel measures 150mm × 100mm and consists of two 2mm thick glass panels and a 1.52mm thick PVB interlayer. The glass density is 2700kg / m³. 3 The elastic modulus is 70 GPa, and the Poisson's ratio is 0.20. The density of the PVB interlayer is 1100 kg / m³. 3 The elastic modulus is 0.1 GPa, and the Poisson's ratio is 0.49. In the experiment, the lower glass layer of the laminated glass was supported by a 2 mm thick rubber pad with a density of 1200 kg / m³. 3 The elastic modulus is 0.042 GPa and the Poisson's ratio is 0.49. The middle of the upper glass plate is impacted by an object with a mass of 2.048 kg and a diameter of 12.7 mm at an impact velocity of 3.83 m / s, and is considered a rigid body.

[0109] Analyzing the failure characteristics of laminated glass using the method of this invention:

[0110] (1) Input the above material parameters into the laminated glass simulation model.

[0111] (2) The laminated glass was discretized according to the test dimensions, with the discretization distance for the glass, PVB, rubber pad, and impactor all set at Δx = 0.50 mm. The time step was 8 × 10⁻⁶. -8 s, and take the calculation time as t. max =3×10 -3 s.

[0112] (3) Boundary conditions: The bottom of the rubber is fixed, that is, the displacement of the bottom surface is 0; the initial velocity of the impacting object is 3.83m / s.

[0113] (4) Search for material points in the near field range of each part. The near field range of the four materials is δ=3Δx=1.50mm.

[0114] (5) Scope of application of the interface model δ B =0.50mm, and determine the material points that need to be applied to the interface model.

[0115] (6) Solve for the interaction forces of the material points in the laminated glass according to the near-field dynamic constitutive equation and interface model, and solve for the displacement and velocity according to the Verlet explicit integration method.

[0116] (7) According to the near-field dynamic damage model, solve for the damage of material points and interface damage in each part of the laminated glass.

[0117] (8) Repeat (6)-(7) until the calculation time t≥t max The calculation will terminate at that time.

[0118] The simulated damage results (failure conditions) of laminated glass using the method of this invention are as follows: Figure 8 As shown, the results of the actual experiment are as follows: Figure 9 As shown in the figure, numerous radial cracks are observed in the center of the upper glass layer (left side), along with several annular cracks near the edge. Compared to the upper glass layer, the radial cracks in the lower glass layer (right side) are more pronounced, but the number of annular cracks is less. Compared to experimental results, the simulation results of this invention capture all the characteristics of the aforementioned laminated glass cracks, demonstrating the applicability and accuracy of this invention in predicting laminated glass failure.

[0119] Figure 10 The figure shows the velocity-time curve of the impactor during the impact process. As can be seen from the figure, the velocity change of the impactor in this invention is basically consistent with the experimental results. Before 1.0 ms, the velocity of the impactor in both the experiment and the results of this invention decreases rapidly. After 1.0 ms, the decreasing trend slows down. It is evident that the glass layers of the laminated glass have failed by about 1.0 ms, and the PVB layer mainly plays a role in the subsequent impact process.

[0120] Compared with existing technologies, the method of this invention utilizes peridynamics, a nonlocal particle-based approach that uses integral equations instead of differential equations for modeling. This avoids singularity issues when solving material discontinuities, thus more accurately calculating the response of laminated glass under external loads. Furthermore, the model includes a description of material damage, naturally and spontaneously simulating crack initiation and propagation under external loads. This invention describes the mechanical behavior of the interface between the glass and the interlayer in laminated glass through a peridynamic interface model. By observing the interactions between the same and different materials, it better simulates the stress conditions of laminated glass, making the calculation results more accurate and reliable. In addition, this invention can represent three failure modes in laminated glass: glass failure, interlayer failure, and interface debonding, resulting in a wider range of applications and more realistic results.

[0121] The above description is only a preferred embodiment of the present invention. It should be noted that for those skilled in the art, several improvements and modifications can be made without departing from the technical principles of the present invention, and these improvements and modifications should also be considered within the scope of protection of the present invention.

Claims

1. A method for optimizing the design of laminated glass structures based on near-field dynamics simulation, characterized in that, Includes the following steps: Step 1: Obtain the expected load on the laminated glass based on its intended installation location; Step 2: Initialize the structural dimensions and material parameters of the laminated glass; Step 3: Input the expected load, structural dimensions and material parameters of the laminated glass into the laminated glass simulation model based on near-field dynamics to obtain the failure condition of the laminated glass under the expected load. Step 4: Determine whether the laminated glass meets the safety specifications based on the failure condition. If it does, obtain the laminated glass structural design scheme under the expected load based on the structural dimensions and material parameters of the laminated glass. If it does not meet the specifications, adjust the structural dimensions and material parameters of the laminated glass within the preset range and return to step 3. Step 3 includes the following specific operations: Step 301: Input the expected load, structural dimensions and material parameters of the laminated glass into the laminated glass simulation model based on near-field dynamics; Step 302: Discretize each part of the laminated glass into near-field dynamic material points and set the calculation time step; Step 303: Determine the boundary conditions of the material points based on the expected load; Step 304: Search for nearby material points based on the near-field range to obtain a list of near-field ranges; Step 305: Obtain the effective range of the interface model based on the interface position of the laminated glass; Step 306: Calculate the interaction forces between material points within the same material of the laminated glass based on the near-field dynamics theory and the near-field range list; Step 307: Based on the interface model and its range of action, calculate the interfacial forces between material points inside different materials on both sides of the interface of the laminated glass. Step 308: Calculate the acceleration of the material point based on the interaction forces and interface forces, and calculate the displacement and velocity of the material point based on the acceleration; Step 309: Calculate the damage to the material points using the damage model based on their displacement and velocity. Step 310: Determine whether the current time t exceeds the preset simulation duration. ,like If the calculation is successful, proceed to the next time step and return to step 306; otherwise, end the calculation. Based on the damage at each material point, obtain the overall damage cloud map of the laminated glass, which represents the failure status of the laminated glass under the expected load.

2. The method for optimizing the design of laminated glass structures based on near-field dynamics simulation according to claim 1, characterized in that, The material parameters include the density, elastic modulus, and Poisson's ratio of the glass in the laminated glass, and the density, elastic modulus, and Poisson's ratio of the interlayer.

3. The method for optimizing the design of laminated glass structures based on near-field dynamics simulation according to claim 1, characterized in that, The specific operation of step 304 is as follows: Randomly select a matter point x, and search its near-field range with matter point x as the center. For the matter points within the range, generate a list of near-field ranges, where... The interval is discrete.

4. The method for optimizing the design of laminated glass structures based on near-field dynamics simulation according to claim 1, characterized in that, In step 305, the scope of the interface model is: ； in, Indicates the scope of the interface model. and These represent the discrete distances on both sides of the laminated glass interface.

5. The method for optimizing the design of laminated glass structures based on near-field dynamics simulation according to claim 4, characterized in that, Matter point x and matter point The conditions for considering interface interactions are: (a) Matter point With matter point The initial distance is less than the effective range of the interface model. ; (b) Matter Points and matter points On both sides of the interface of the laminated glass, i.e., material points and matter points The materials have different properties.

6. The method for optimizing the design of laminated glass structures based on near-field dynamics simulation according to claim 1, characterized in that, The acceleration of point x under the interaction forces and interfacial forces is: ； in, This represents the density of a point x. This represents the acceleration of point x at time t. This represents the matter point within the near-field range of matter point x. The set, The interface model representing the material point x represents the material points within the interface range. The set, This represents the interaction force experienced by a point x within the same material at time t. Represents the material point at time t The interaction forces that are borne within the same material Representing a point of matter volume, This represents the interfacial force experienced by a point x at time t within the different materials on either side of the interface. Represents the matter point at time t The interfacial forces experienced by the different materials on both sides of the interface. Let x represent the volumetric density of a point x.

7. The method for optimizing the design of laminated glass structures based on near-field dynamics simulation according to claim 1, characterized in that, In step 308, the displacement and velocity of the material point are calculated using the Verlet explicit integration method, and the iteration method is as follows: ； in, This represents the displacement of the material point at time step n. This represents the velocity of the substance point at time step n. This represents the acceleration of the substance point at time step n. For time step.

8. The method for optimizing the design of laminated glass structures based on near-field dynamics simulation according to claim 1, characterized in that, When matter point x is related to matter points within its near field During interaction, the failure criterion for material point x is that the elongation s is greater than the critical elongation. For all matter points within the near field range of matter point x The damage at material point x is obtained by statistically analyzing the broken bonds: ； ； ； ； in, and They are the material point x and The position vector, G0 is the energy release rate of the material to which the material point x belongs, and K is the bulk modulus of the material to which the material point x belongs. Let x be the near-field range of the matter point. Let matter point x be the other matter points within its near-field range. discontinuous functions, This represents the damage to the material point x at time t. Represents a point of matter within the near field. The set, Representing a point of matter The volume.

9. The method for optimizing the design of laminated glass structures based on near-field dynamics simulation according to claim 1, characterized in that, Failure scenarios of laminated glass under expected loads include glass breakage, debonding of the glass-laminate interface, and lamination failure.

Images