A parameterized finite element simulation calculation method and storage medium for improving design efficiency of point-supported glass curtain walls

Abstract

The present application relates to the technical field of building curtain wall design, and more particularly to a parameterized finite element simulation calculation method of point support glass curtain wall for improving design efficiency and a storage medium, comprising: defining an input parameter set and setting it as a global variable; obtaining a segmented geometric model; giving the segmented geometric model unit and material properties; performing mesh division on the segmented geometric model to generate finite element meshes, so that corresponding finite element nodes exist at geometric key point positions; identifying the finite element nodes and applying displacement constraints, while applying uniform pressure loads; obtaining mechanical response data, when any variable is changed, repeating modeling calculation based on the changed parameters; extracting stress distribution, deformation and support reaction data from the finite element solution calculation results. The present application realizes parameterized automatic simulation, accurately captures stress concentration; optimizes the use of computing resources while accurately simulating the large deformation behavior of glass brittle materials, and improves design efficiency and structural safety.

Description

Technical Field

[0001] This invention relates to the field of building curtain wall design technology, and in particular to a parametric finite element simulation calculation method and storage medium for point-supported glass curtain walls that improves design efficiency. Background Technology

[0002] With the rapid development of building curtain wall technology, point-supported glass curtain walls have been widely used in modern building engineering due to their good transparency and aesthetic appeal. Against this backdrop, to ensure the safety and reliability of the curtain wall structure, stress analysis of the point-supported glass panels has become a crucial step in the design process. Traditional techniques typically employ two calculation schemes: simplified mechanical models based on analytical methods and manual modeling and calculation methods based on general-purpose finite element software.

[0003] The simplified mechanical model based on analytical methods simplifies the point-supported glass panel into a mechanical model. Based on the theory of elastic thin plates, stress and deformation are approximated by looking up tables or empirical formulas. This method divides the glass panel into several calculation regions, calculates the bending moment and deflection of each region separately, and then obtains the overall response through the superposition principle. However, analytical methods rely on many simplification assumptions and cannot fully consider the refined influence of differences in the location, number, and opening size of point supports on the overall mechanical properties. Furthermore, they are difficult to accurately reflect the true stress distribution and deformation characteristics of the glass under actual stress conditions, especially in stress concentration areas near the point supports.

[0004] The scheme based on the manual modeling and calculation method of general-purpose finite element software uses general-purpose finite element analysis software (such as SAP2000, ANSYS Workbench, etc.) for numerical simulation. This scheme requires manually creating the geometric model of the glass panel in the software interface, defining material properties, manually creating constraint points according to the point support location, manually generating the mesh, applying load conditions, and finally executing the solution calculation and manually extracting the results. However, in practical applications, this scheme involves a cumbersome and inefficient modeling process. Each calculation of a glass model with different dimensions or support positions requires a complete re-execution of the modeling process, including geometric drawing, mesh generation, and boundary condition setting. This results in a large amount of repetitive work and low computational efficiency. Furthermore, this scheme has a high operational threshold, relying heavily on professional experience. Operators need in-depth knowledge of finite element theory and software operation experience, making it difficult for non-professionals to master and use, thus limiting its widespread application in engineering design. In addition, the consistency of the calculation results is difficult to guarantee: manual operation is prone to modeling differences due to human factors, such as mesh generation quality and constraint application location, affecting the accuracy and comparability of the calculation results. Moreover, the scheme is difficult to compare and iterate: when comparing multiple schemes or performing parameter sensitivity analysis, it is necessary to repeatedly modify the model manually and recalculate, making rapid iterative optimization difficult.

[0005] In practical engineering applications, the size, thickness, and support point layout of point-supported glass vary greatly, requiring frequent adjustments during the design phase. However, traditional methods, once the calculation scheme is determined, are difficult to adjust flexibly. Repeated modeling leads to extended design cycles and increased costs, resulting in significant shortcomings in terms of economy and practicality. It is evident that existing technologies in the stress analysis of point-supported glass curtain walls either fail to accurately capture stress concentrations near the support points due to excessive simplification assumptions, creating potential structural safety hazards; or they rely on manual operation, which is not only inefficient and involves a large amount of repetitive work, but also results in inconsistent modeling quality due to varying operator skill levels, with constraint application positions not matching design drawings. In severe cases, this can lead to calculation distortions and even curtain wall engineering accidents. Summary of the Invention

[0006] The purpose of this invention is to address the shortcomings of existing technologies by providing a parametric finite element simulation calculation method for point-supported glass curtain walls that improves design efficiency, including:

[0007] S1: Define a set of input parameters including the geometric dimensions of the glass panel, the coordinates of the point support positions, the glass thickness, and the external load, and set the set of input parameters as a global variable;

[0008] S2: Generate an initial geometric model of the glass panel based on the geometric dimensions, generate geometric key points on the initial geometric model with the point support position coordinates as a reference, and obtain a segmented geometric model defined by the geometric key points through Boolean segmentation;

[0009] S3: Based on the glass thickness in the global variables, define the shell element type, material properties, and cross-sectional thickness, and assign element and material properties to the segmented geometric model;

[0010] S4: Mesh the segmented geometric model to generate a finite element mesh, and ensure that there are corresponding finite element nodes at the locations of the geometric key points;

[0011] S5: Automatically identify the finite element nodes corresponding to the geometric key point positions, apply displacement constraints to the identified finite element nodes, and simultaneously apply a uniformly distributed pressure load defined by the global variable to the glass panel surface.

[0012] S6: Perform finite element calculation to obtain the mechanical response data of the glass panel. When any variable in the input parameter set changes, repeat steps S2 to S6 based on the changed parameters.

[0013] S7: Extract stress distribution, deformation and support reaction data from the finite element solution results, and generate corresponding contour maps and numerical reports.

[0014] Preferably, in step S1, the input parameter set is defined, including:

[0015] Physical units can be explicitly set using the unit system definition function of finite element software.

[0016] The variable definition function defines the glass panel width and height as geometric dimension variables, the coordinates of each support point in the panel plane as coordinate parameter variables, the glass thickness as thickness parameter variables, and the uniformly distributed pressure load value acting on the glass surface as external load parameter variables.

[0017] Preferably, in step S2, generating geometric key points includes:

[0018] The geometric modeling kernel of the finite element software is invoked to create geometric key points using the coordinates of the point support position as parameters;

[0019] Create an auxiliary dividing line that passes through each geometric key point and is parallel to the coordinate axis;

[0020] The initial geometric model is cut using the auxiliary dividing lines, so that each geometric key point becomes a common vertex or edge point of adjacent sub-faces, thus becoming a topological node of the segmented geometric model.

[0021] Preferably, in step S3, assigning element and material properties to the segmented geometric model includes:

[0022] By defining the element type, the element type of the segmented geometric model is defined as a shell element suitable for thin shell structures;

[0023] The material property definition function allows you to define material property parameters, including the elastic modulus, Poisson's ratio, and density of glass materials.

[0024] The glass thickness parameter in the global variable is assigned to the cross-sectional thickness attribute of the shell unit by using the cross-sectional attribute definition function.

[0025] Preferably, in step S4, the segmented geometric model is meshed, including:

[0026] A quadrilateral mesh is generated using a mapping mesh generation algorithm, and the geometric key points are used as vertices or edge points of the geometric model, which will inevitably generate finite element nodes at the corresponding positions.

[0027] The generated finite element mesh is subjected to quality checks, which include at least the evaluation of element warpage, aspect ratio, and interior angle.

[0028] If the inspection results do not meet the preset quality threshold, the mesh generation parameters will be automatically adjusted and the mesh generation will be re-executed until a finite element mesh that meets the quality requirements is generated.

[0029] Preferably, in step S5, applying displacement constraints and uniformly distributed pressure loads includes:

[0030] Obtain the coordinates of the geometric key points of each support point, set the numerical tolerance based on the coordinates, and search for nodes with matching coordinates in the finite element node set;

[0031] Verify the number of nodes found. If it is unique, then the node is the finite element node at the corresponding support point location. Apply displacement constraints to the identified node.

[0032] Select the geometric surface of the glass panel using the surface selection function, and apply the uniformly distributed pressure load to the geometric surface using the surface load application function.

[0033] Preferably, in step S6, performing the finite element solution calculation includes:

[0034] The analysis type is specified by the analysis type definition function, and the analysis type includes at least static analysis;

[0035] Activate the geometric nonlinear large deformation effect through the geometric nonlinearity setting function;

[0036] Perform the solution calculation to obtain the mechanical response results, including the stress, deformation, and support reaction of the glass panel.

[0037] Preferably, in step S6, the finite element solution calculation is performed, including:

[0038] In response to a change in any variable value in the input parameter set, a threshold judgment is made on the amount of change. If the change exceeds the preset engineering tolerance threshold, the current finite element model is automatically cleared, and steps S2 to S6 are re-executed based on the changed parameters. If the change is within the tolerance threshold range, the current finite element model remains unchanged.

[0039] Preferably, in step S7, stress distribution, deformation, and support reaction data are extracted from the finite element solution results, and corresponding contour maps and numerical reports are generated, including:

[0040] By using the node-based cloud map display function, stress component parameters can be specified to generate and output stress distribution cloud maps of the glass panel.

[0041] By using the node-based cloud map display function, specifying displacement component parameters, the deformation result cloud map of the glass panel is generated and output;

[0042] Using the reaction force solution list function, filter the support point nodes with applied displacement constraints, calculate and list the support reaction force data.

[0043] Based on the same concept, the present invention also provides a storage medium storing computer-readable instructions, which, when executed by one or more processors, cause the one or more processors to perform the steps of a parametric finite element simulation calculation method for a point-supported glass curtain wall that improves design efficiency, as described in any one embodiment.

[0044] Compared with the prior art, the beneficial effects of the present invention are:

[0045] (1) This invention defines a set of input parameters, including the geometric dimensions of the glass panel, the coordinates of the point support position, the glass thickness and the external load, and sets them as global variables. Based on the set of input parameters, it automatically performs the entire process of geometric modeling, mesh generation, constraint application, solution calculation and result output. This allows for the rapid solution of different specifications of models by only modifying the parameter variables. It solidifies the repetitive labor of traditional manual modeling into an automated process, thereby greatly improving the design efficiency of point-supported glass curtain walls. At the same time, it eliminates modeling errors caused by differences in manual operation and reduces the safety hazards of curtain wall structure design from the source.

[0046] (2) This invention generates geometric key points based on the coordinates of the support position, obtains the segmented geometric model defined by the geometric key points by Boolean segmentation, makes each geometric key point a common vertex or edge point of adjacent sub-faces, and ensures that there are corresponding finite element nodes at the position of the geometric key points during mesh generation, so as to achieve a strict one-to-one correspondence between the support point position and the finite element node, avoids the constraint application error caused by the need for node coupling or multi-point constraint approximation processing due to the support point falling inside the element in the traditional method, thereby accurately capturing the stress concentration near the support point and ensuring the reliability of the safety assessment of the glass curtain wall structure.

[0047] (3) This invention activates the geometric nonlinear large deformation effect to adapt to the real mechanical response of the glass panel under large deflection conditions, and performs threshold judgment on the change when any variable in the input parameter set changes. Only when the change exceeds the preset engineering tolerance threshold is the full process recalculation triggered, avoiding invalid calculation due to small parameter changes. It can accurately simulate the brittle failure characteristics and large deformation behavior of glass materials, and optimize the efficiency of computing resources. Thus, while ensuring the accuracy of structural analysis, it supports rapid comparison and iterative optimization of multiple schemes, and has significant advantages of high efficiency, high reliability and high safety. Attached Figure Description

[0048] Various other advantages and benefits will become apparent to those skilled in the art upon reading the following detailed description of preferred embodiments. The accompanying drawings are for illustrative purposes only and are not intended to limit the invention.

[0049] Figure 1This is a flowchart of a parametric finite element simulation calculation method for point-supported glass curtain walls that improves design efficiency according to the present invention.

[0050] Figure 2 This is another flowchart of a parametric finite element simulation calculation method for point-supported glass curtain walls that improves design efficiency according to the present invention.

[0051] Figure 3 This is a typical geometric model diagram of point-supported glass;

[0052] Figure 4 This is a distribution diagram of the stress calculation results for the point-supported glass of the present invention;

[0053] Figure 5 This is a diagram showing the calculation results of the deformation of the point-supported glass in this invention;

[0054] Figure 6 The diagram shows the calculation results of the point-supported glass constraint reaction force according to the present invention. Detailed Implementation

[0055] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention. Obviously, the described embodiments are only some, not all, of the embodiments described in this application. All other embodiments obtained by those skilled in the art based on the embodiments in this application without creative effort are within the scope of protection of this application.

[0056] Those skilled in the art will understand that, unless otherwise stated, the singular forms “a” and “an” used herein, and “the”, may also include the plural forms. It should be further understood that the term “comprising” as used in this specification means the presence of the stated features, integers, steps, operations, elements, and / or components, but does not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and / or groups thereof.

[0057] First Embodiment

[0058] Please see Figure 1 and Figure 2 As shown, this embodiment provides a parametric finite element simulation calculation method for point-supported glass curtain walls to improve design efficiency, including the following steps:

[0059] S1: Define a set of input parameters including the geometric dimensions of the glass panel, the coordinates of the point support positions, the glass thickness, and the external loads, and set the set of input parameters as a global variable.

[0060] Preferably, in step S1, the input parameter set is defined, including:

[0061] The physical units are explicitly set using the unit system definition function of the finite element software. Specifically, in this embodiment, the ANSYS analysis environment is started, the current database is cleared, and a new work file name and analysis title are set. The physical unit system used by the analysis environment is explicitly defined using the APDL / UNITS command or an equivalent method.

[0062] The glass panel width and height are defined as geometric dimension variables using the variable definition function (variable assignment statements in ANSYS Parametric Design Language (APDL)). The coordinates of each support point within the panel plane are defined as coordinate parameter variables, the glass thickness is defined as a thickness parameter variable, and the uniformly distributed pressure load acting on the glass surface is defined as an external load parameter variable. Specifically, in this embodiment, the / CLEAR command clears all model and result data from the current database, creating a completely new analysis environment. The / FILNAME command sets the workfile name to "PointSupportedGlass," the / TITLE command sets the analysis title to "Analysis of PointSupported Glass Panel," and comments explicitly specify that the physical unit system used in the analysis environment is N-mm-Kg (Newton-millimeters-kilogram).

[0063] The code is as follows:

[0064] FINISH! End the current module.

[0065] / CLEAR,START ! Clears the database and starts a new analysis.

[0066] / FILNAME,BL ! Sets the working file name to BL

[0067] / TITLE,BL / UNINT:N,mm,Kg !Sets the title and specifies the unit system (Newton, millimeter, kilogram).

[0068] More preferably, the point support position coordinate parameters are as follows:

[0069] x1=85 !x coordinate (left constraint point)

[0070] x2 = 1115 !x coordinate (right constraint point)

[0071] x3 = 1200 !x coordinate (right side of the glass)

[0072] y1=150 !y coordinate (lower constraint point)

[0073] y2=1000 !y coordinate (upper constraint point)

[0074] y3 = 1150 !y coordinate (top of the glass)

[0075] th_bl=6 ! Calculate the thickness of the glass

[0076] press_bl=1.0e-3 ! Load, unit MPa.

[0077] S2: Generate an initial geometric model of the glass panel based on the geometric dimensions. Generate geometric key points on the initial geometric model based on the coordinates of the point support positions. Obtain a segmented geometric model defined by the geometric key points through Boolean segmentation. Specifically, in this embodiment, the modeling command flow of the ANSYS parametric design language APDL is called to automatically construct the initial geometric model.

[0078] Preferably, in step S2, generating geometric key points includes:

[0079] The geometric modeling kernel of the finite element software is invoked to create geometric key points using the coordinates of the point support positions as parameters. Specifically, in this embodiment, based on the geometric dimension parameters, the key point creation command K of the ANSYS parametric design language APDL is invoked to create key points representing the four corner points of the glass panel.

[0080] Even better, the K command creates keypoints at the four corners of the glass panel, as shown in the following code snippet:

[0081] K,1,0,0

[0082] K,2,x1,Ky(1)

[0083] K,3,x2,Ky(1)

[0084] K,4,x3,Ky(1) .

[0085] Using each geometric key point as the center, an auxiliary dividing line is created that passes through the point and is parallel to the coordinate axis. Specifically, in this embodiment, the key points are connected by the line creation command L to generate the boundary line of the glass panel, and the geometric plane representing the original glass panel is generated from the boundary line by the surface creation command A.

[0086] Even better, the L command connects key points to generate lines, as shown in the following code snippet:

[0087] L,1,2

[0088] L,2,3

[0089] L,3,4

[0090] LGEN,2,1,3,1, ,y1,,,0

[0091] LGEN,2,1,3,1, ,y2,,,0

[0092] LGEN,2,1,3,1, ,y3,,,0 !Copy lines to generate the boundaries of the glass plate (LGEN command)

[0093] L,1,5

[0094] L,5,9

[0095] L,9,13

[0096] LGEN,2,13,15,1,X1,,,,0

[0097] LGEN,2,13,15,1,X2,,,,0

[0098] LGEN,2,13,15,1,X3,,,,0

[0099] Nummrg,kp !Merge keypoints (Nummrg,kp)

[0100] NUMCMP,all !Compression number (NUMCMP,all).

[0101] The initial geometric model is cut using auxiliary dividing lines, making each geometric key point a common vertex or edge point of adjacent sub-faces, thus becoming a topological node of the divided geometric model. Specifically, in this embodiment, command A is called to generate a geometric plane representing the original glass panel from the region enclosed by these four boundary lines. Part of the code is shown below:

[0102] A,1,4,16,13 ! Create face (A command).

[0103] Based on the coordinate parameters of the point support positions, the K command is used to create the positioning key points corresponding to each point support on the geometric plane. Then, the ASBL command is used to perform Boolean segmentation on the geometric plane with the positioning key points as the segmentation reference. Finally, a finite element geometric model with geometric segmentation points at the support points is formed, which can be used for subsequent mesh generation and constraint application.

[0104] More preferably, the ASBL command is invoked to perform Boolean segmentation on the original geometric surface using the cutting lines generated through the aforementioned support points. This ensures that the support point positions are precisely preserved as independent geometric entities in the geometric model, ultimately generating a model like... Figure 3 The finite element geometric model of the glass panel, which includes the point support location features, is shown below, with some code examples:

[0105] ASBL,1,ALL ! Create a glass plate area by cutting the surface through ASBL.

[0106] Traditional mapping mesh generation methods often result in support points falling inside or on the edges of elements when generating meshes, failing to automatically generate corresponding finite element nodes. This leads to the need for approximation through node coupling or multi-point constraints when applying displacement constraints. This not only introduces constraint application errors but also makes it impossible to accurately capture stress concentrations near support points. For glass curtain walls made of brittle materials, this error directly leads to an underestimation of peak stress, posing a structural safety hazard.

[0107] This embodiment uses a forced Boolean partitioning operation to cut the original geometric surface at the support points, making each support point a common vertex or edge point of adjacent geometric sub-surfaces. Therefore, during subsequent mapping and mesh generation, these geometric key points will inevitably generate corresponding finite element nodes. This design ensures a strict one-to-one correspondence between the support point positions and the finite element nodes. Displacement constraints can be directly applied to the nodes without any approximation, eliminating positional errors in constraint application at the geometric topology level. This guarantees the accuracy of stress concentration calculations near the support points and lays a reliable mechanical foundation for the safety assessment of glass curtain wall structures.

[0108] S3: Based on the glass thickness in the global variables, define the shell element type, material properties and cross-sectional thickness, and assign element and material properties to the segmented geometric model.

[0109] Preferably, in step S3, assigning element and material properties to the segmented geometric model includes:

[0110] The element type definition function defines the element type of the segmented geometric model as a shell element suitable for thin shell structures. Specifically, in this embodiment, the element type definition command ET of the ANSYS parametric design language APDL defines the element type of the geometric model as a SHELL181 shell element. The key options of the SHELL181 shell element are set according to the mechanical analysis requirements of the glass panel through the element key option setting command KEYOPT. The key options include, but are not limited to, element coordinate system, large deformation effect switch, and integration scheme.

[0111] The material property definition function defines material property parameters, including the elastic modulus, Poisson's ratio, and density of glass materials. Specifically, in this embodiment, the material property parameters are defined using the material property definition command MP.

[0112] The glass thickness parameter in the global variable is assigned to the section thickness attribute of the shell element through the section attribute definition function. Specifically, in this embodiment, the section type corresponding to the shell element is created by the section type definition command SECTYPE and the section number is specified. The glass thickness parameter is assigned to the thickness attribute of the shell section by the section parameter assignment command SECDATA.

[0113] More preferably, the ET command is used to define the element type as SHELL181. This element is suitable for thin-shell structure analysis and can accurately simulate the bending and membrane stress behavior of glass panels.

[0114] ET,1,181 !Use SHELL181 unit (for thin shell structures).

[0115] Even better, the KEYOPT command can be used to set key element options to activate layer data storage, facilitating post-processing for viewing stress distribution.

[0116] KEYOPT,1,8,2 !Set cell options (storage layer data).

[0117] More preferably, the MP command is used to define the linear elastic properties of the glass material. Taking ordinary float glass as an example, the elastic modulus, Poisson's ratio, and density are defined as follows:

[0118] MP,ex,1,72000 !Elastic modulus 72 GPa (unit N / mm²)

[0119] MP,prxy,1,0.2 !Poisson's ratio 0.2

[0120] MP,dens,1,2.56e-9 !Density 2.56e-9 kg / mm³.

[0121] The SECTYPE command is used to create the section type corresponding to the shell element and specify the section number. The SECDATA command is used to assign the defined glass thickness parameter variable to the section thickness, thereby realizing the dynamic association between the element section properties and key design parameters.

[0122] SECDATA, th_bl,1,0.0,3 ! Defines the shell cross-sectional thickness.

[0123] S4: Mesh the segmented geometric model to generate a finite element mesh, and ensure that there are corresponding finite element nodes at the locations of geometric key points.

[0124] Preferably, in step S4, the segmented geometric model is meshed, including:

[0125] A quadrilateral mesh is generated using a mapping meshing algorithm, making geometric key points the vertices or edge points of the geometric model. This inevitably generates finite element nodes at the corresponding locations. Specifically, in this embodiment, the ANSYS APDL element size setting command ESIZE is called to set the global element size according to the geometric dimensions and analysis accuracy requirements (setting the global element size to 50mm as the basic parameter for controlling mesh density: ESIZE,50,0 ! Element size 50mm). Smart meshing is disabled using the SMRTSIZE command. The APDL mesh shape specification command MSHAPE is called to specify the mesh shape as quadrilateral. The mapping meshing algorithm is selected using the MSHKEY command to obtain a regular, high-precision structured mesh: MSHAPE,0,2d ! Quadrilateral mesh; MSHKEY,1 ! Mapped mesh. The APDL face meshing command AMESH is called to perform automatic meshing on the geometric model, generating a computational mesh: AMESH,all ! Mesh all faces.

[0126] The generated finite element mesh is subjected to quality checks, which include at least the evaluation of element warpage, aspect ratio, and interior angle. Specifically, in this embodiment, the APDL mesh quality check command CHECK is called to perform quality checks on the generated computational mesh.

[0127] If the inspection results do not meet the preset quality threshold, the mesh generation parameters will be automatically adjusted and the mesh generation will be re-executed until a finite element mesh that meets the quality requirements is generated.

[0128] S5: Automatically identify finite element nodes corresponding to the geometric key point locations, apply displacement constraints to the identified finite element nodes, and simultaneously apply a uniformly distributed pressure load defined by global variables to the glass panel surface.

[0129] Preferably, in step S5, applying displacement constraints and uniformly distributed pressure loads includes:

[0130] The coordinates of the geometric key points of each support point are obtained, and the numerical tolerance is set based on these coordinates. The nodes with matching coordinates are searched in the finite element node set. Specifically, in this embodiment, the node selection function of APDL is called to locate the node corresponding to the support point position in the computational grid.

[0131] Verify the number of nodes found. If it is unique, then the node is the finite element node at the corresponding support point position. Apply displacement constraints to the identified nodes. Specifically, in this embodiment, displacement constraints are applied to the nodes through the key point displacement constraint command DK or the node displacement constraint command D.

[0132] The geometric surface of the glass panel is selected by the surface selection function, and a uniform pressure load is applied to the geometric surface by the surface load application function. Specifically, in this embodiment, the geometric surface representing the entire glass panel is selected by the surface selection command ASEL, and a uniform pressure load defined by external load parameters is applied to the geometric surface by the surface load application command SFA. DK,10,UX,,,,UY,UZ,,,ROTZ !Fix key point 10 (UX,UY,UZ,ROTZ);

[0133] DK,11,UX,,,,UY,UZ,,,ROTZ !Fixed keypoint 11;

[0134] DK,6,UX,,,,,UZ,,,ROTZ !Constrain key point 6 (UX,UZ,ROTZ);

[0135] DK,7,UX,,,,,UZ,,,ROTZ !Constrain key point 7;

[0136] sf,all,pres,press_bl !Apply pressure load;

[0137] GPLOT;

[0138] !ACEL,0,9.8e3,0, !Specifies gravitational acceleration; the elevation should be specified by y.

[0139] In this step, the program automatically identifies finite element nodes that strictly coincide with the support point and applies completely fixed displacement constraints (constraining the three translational degrees of freedom and the rotational degree of freedom about the surface normal). This approach fully considers the characteristics of brittle fracture of glass: the constraint position is consistent with the actual stress state of the support point, avoiding false stress concentration or constraint stiffness distortion caused by approximate constraint treatment.

[0140] S6: Perform finite element analysis to obtain the mechanical response data of the glass panel. When any variable in the input parameter set changes, repeat steps S2 to S6 based on the changed parameters.

[0141] Preferably, in step S6, the finite element solution calculation is performed, including:

[0142] The analysis type is specified through the analysis type definition function. The analysis type includes at least static analysis. Specifically, in this embodiment, the analysis type is specified through the APDL analysis type definition command ANTYPE.

[0143] The geometric nonlinearity setting function is used to activate the geometric nonlinearity large deformation effect. Specifically, in this embodiment, the geometric nonlinearity setting command NLGEOM is used to activate the geometric nonlinearity large deformation effect, NLGEOM,1 ! Considering large deformation.

[0144] When glass panels are subjected to wind loads or concentrated loads, even if the stress level has not reached the ultimate strength, the deflection can be several times the thickness. In this case, linear analysis will significantly overestimate stiffness and underestimate deformation, leading to a dangerous design. This invention automatically adapts to the large deformation characteristics of glass curtain walls through a built-in geometric nonlinear switch, enabling the calculation results to truly reflect the mechanical behavior of the structure under ultimate conditions, providing a reliable basis for safety assessment.

[0145] The solution calculation is performed to obtain the mechanical response results, including the stress, deformation, and support reaction of the glass panel. Specifically, in this embodiment, the solution calculation is performed using the SOLVE command to obtain the mechanical response results, including the stress, deformation, and support reaction of the glass panel. ASEL,all ! Select all faces; SOLVE ! Solve; exit the solver module using the FINISH command; fini ! End the solver.

[0146] Preferably, in step S6, the finite element solution calculation is performed, including:

[0147] In response to any change in the value of any variable in the input parameter set, a threshold judgment is made on the amount of change. If the change exceeds the preset engineering tolerance threshold, the current finite element model is automatically cleared, and steps S2 to S6 are re-executed based on the changed parameters. If the change is within the tolerance threshold range, the current finite element model is kept unchanged or only a local update is performed to avoid unnecessary full-process recalculation triggered by minor parameter changes (such as modeling rounding errors), thereby optimizing the efficiency of computing resource utilization and keeping the values ​​of unchanged parameters unchanged.

[0148] This solution also incorporates a parameter change monitoring mechanism. When a user modifies any input parameter, the system automatically calculates the rate of change and compares it with preset engineering tolerance thresholds (e.g., geometric dimension changes exceeding 1mm, thickness changes exceeding 0.1mm, load changes exceeding 1%). A complete recalculation is triggered only when the change exceeds the threshold; minor changes due to unit conversion or rounding are automatically ignored, and existing calculation results are retained. This mechanism effectively avoids a large amount of invalid calculations caused by parameter fine-tuning, significantly improving the utilization efficiency of computing resources in multi-solution comparison scenarios.

[0149] S7: Extract stress distribution, deformation and support reaction data from the finite element solution results, and generate corresponding contour maps and numerical reports.

[0150] Preferably, in step S7, stress distribution, deformation, and support reaction data are extracted from the finite element solution results, and corresponding contour maps and numerical reports are generated, including:

[0151] By using the node de-clouding display function, stress component parameters are specified respectively, and stress distribution cloud maps of the glass panel are generated and output. Specifically, in this embodiment, the ANSYS general post-processor module is entered and the result file is loaded. The stress component parameters are specified respectively through the node de-clouding display command PLNSOL, and stress distribution cloud maps of the glass panel are generated and output.

[0152] By using the node-based cloud map display function, the displacement component parameters are specified, and the deformation result cloud map of the glass panel is generated and output. Specifically, in this embodiment, the displacement component parameters are specified using the PLNSOL command, and the deformation result cloud map of the glass panel is generated and output.

[0153] The support reaction force data is calculated and output in a list by using the reaction force solution list function to filter the support point nodes with applied displacement constraints. Specifically, in this embodiment, the support point nodes with applied displacement constraints are filtered by using the reaction force solution list command PRRSOL, and the support reaction force data is calculated and output in a list.

[0154] More preferably, the / POST1 command is invoked to enter the general post-processor module. The program automatically reads the result file generated after step S6 is completed, as shown below:

[0155] / POST1

[0156] / EFACET,1 .

[0157] Please see Figure 4 As shown, the PLNSOL command is invoked, specifying the display result as the first principal stress, to plot the stress distribution on the model using a continuous contour plot, as shown below:

[0158] PLNSOL, S,1, 0,1.0.

[0159] Please see Figure 5 As shown, the PLNSOL command is invoked, specifying the display result as the total displacement (U, SUM), to generate and output the deformation pattern and displacement contour map of the glass panel under load, as shown below:

[0160] !Show deformed image

[0161] / POST1

[0162] / EFACET,1

[0163] PLNSOL, U,SUM, 0,1.0.

[0164] Please see Figure 6 As shown, the PRRSOL command is invoked, and the reaction forces acting on the support points with applied displacement constraints are selected as output items. The reaction force values ​​for each support are calculated and output in list form.

[0165] Display support reaction force

[0166] / POST1

[0167] PRRSOL,FZ! This indicates the support reaction force.

[0168] Second Embodiment

[0169] Based on the same concept, this embodiment also provides a computer device, including a memory and a processor. The memory stores computer-readable instructions, which, when executed by the processor, cause the processor to perform the steps of a parametric finite element simulation calculation method for a point-supported glass curtain wall that improves design efficiency, as described in the embodiment.

[0170] Based on the same concept, the present invention also provides a storage medium storing computer-readable instructions, which, when executed by one or more processors, cause the one or more processors to perform the steps of a parametric finite element simulation calculation method for a point-supported glass curtain wall that improves design efficiency, as described in any one embodiment.

[0171] It is understandable that the parametric finite element simulation calculation method for point-supported glass curtain walls mentioned above, which improves design efficiency, can be stored in a computer-readable storage medium if implemented as software functional modules and sold or used as independent products. Based on this understanding, the technical solution of this invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer server or a network device, etc.) to execute all or part of the steps of the methods in the various embodiments of this invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0172] Computer-readable storage media may include data signals propagated in baseband or as part of a carrier wave, carrying readable program code. Such propagated data signals may take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. A readable storage medium may also be any readable medium other than a readable storage medium that can transmit, propagate, or transfer a program for use by or in connection with an instruction execution system, apparatus, or device. The program code contained on the readable storage medium may be transmitted using any suitable medium, including but not limited to wireless, wired, optical fiber, RF, etc., or any suitable combination thereof.

[0173] The above description is merely a preferred embodiment of the present invention. The scope of protection of the present invention is not limited to the above embodiments. All technical solutions falling within the scope of the present invention's concept are within the scope of protection of the present invention. It should be noted that for those skilled in the art, any improvements and modifications made without departing from the principles of the present invention should also be considered within the scope of protection of the present invention.

Claims

1. A parametric finite element simulation calculation method for point-supported glass curtain walls to improve design efficiency, characterized in that, Includes the following steps: S1: Define a set of input parameters including the geometric dimensions of the glass panel, the coordinates of the point support positions, the glass thickness, and the external load, and set the set of input parameters as a global variable; S2: Generate an initial geometric model of the glass panel based on the geometric dimensions, generate geometric key points on the initial geometric model with the point support position coordinates as a reference, and obtain a segmented geometric model defined by the geometric key points through Boolean segmentation; S3: Based on the glass thickness in the global variables, define the shell element type, material properties, and cross-sectional thickness, and assign element and material properties to the segmented geometric model; S4: Mesh the segmented geometric model to generate a finite element mesh, and ensure that there are corresponding finite element nodes at the locations of the geometric key points; S5: Automatically identify the finite element nodes corresponding to the geometric key point positions, apply displacement constraints to the identified finite element nodes, and simultaneously apply a uniformly distributed pressure load defined by the global variable to the glass panel surface. S6: Perform finite element calculation to obtain the mechanical response data of the glass panel. When any variable in the input parameter set changes, repeat steps S2 to S6 based on the changed parameters. S7: Extract stress distribution, deformation and support reaction data from the finite element solution results, and generate corresponding contour maps and numerical reports.

2. The parametric finite element simulation calculation method for point-supported glass curtain walls to improve design efficiency according to claim 1, characterized in that, In step S1, the input parameter set is defined, including: Physical units can be explicitly set using the unit system definition function of finite element software. The variable definition function defines the glass panel width and height as geometric dimension variables, the coordinates of each support point in the panel plane as coordinate parameter variables, the glass thickness as thickness parameter variables, and the uniformly distributed pressure load value acting on the glass surface as external load parameter variables.

3. The parametric finite element simulation calculation method for point-supported glass curtain walls to improve design efficiency according to claim 1, characterized in that, In step S2, geometric key points are generated, including: The geometric modeling kernel of the finite element software is invoked to create geometric key points using the coordinates of the point support position as parameters; Create an auxiliary dividing line that passes through each geometric key point and is parallel to the coordinate axis; The initial geometric model is cut using the auxiliary dividing lines, so that each geometric key point becomes a common vertex or edge point of adjacent sub-faces, thus becoming a topological node of the segmented geometric model.

4. The parametric finite element simulation calculation method for point-supported glass curtain walls to improve design efficiency according to claim 1, characterized in that, In step S3, the segmented geometric model is assigned element and material properties, including: By defining the element type, the element type of the segmented geometric model is defined as a shell element suitable for thin shell structures; The material property definition function allows you to define material property parameters, including the elastic modulus, Poisson's ratio, and density of glass materials. The glass thickness parameter in the global variable is assigned to the cross-sectional thickness attribute of the shell unit by using the cross-sectional attribute definition function.

5. The parametric finite element simulation calculation method for point-supported glass curtain walls to improve design efficiency according to claim 1, characterized in that, In step S4, the segmented geometric model is meshed, including: A quadrilateral mesh is generated using a mapping mesh generation algorithm, and the geometric key points are used as vertices or edge points of the geometric model, which will inevitably generate finite element nodes at the corresponding positions. The generated finite element mesh is subjected to quality checks, which include at least the evaluation of element warpage, aspect ratio, and interior angle. If the inspection results do not meet the preset quality threshold, the mesh generation parameters will be automatically adjusted and the mesh generation will be re-executed until a finite element mesh that meets the quality requirements is generated.

6. The parametric finite element simulation calculation method for point-supported glass curtain walls to improve design efficiency according to claim 1, characterized in that, In step S5, displacement constraints and uniformly distributed pressure loads are applied, including: Obtain the coordinates of the geometric key points of each support point, set the numerical tolerance based on the coordinates, and search for nodes with matching coordinates in the finite element node set; Verify the number of nodes found. If it is unique, then the node is the finite element node at the corresponding support point location. Apply displacement constraints to the identified node. Select the geometric surface of the glass panel using the surface selection function, and apply the uniformly distributed pressure load to the geometric surface using the surface load application function.

7. The parametric finite element simulation calculation method for point-supported glass curtain walls to improve design efficiency according to claim 1, characterized in that, In step S6, performing the finite element solution calculation includes: The analysis type is specified by the analysis type definition function, and the analysis type includes at least static analysis; Activate the geometric nonlinear large deformation effect through the geometric nonlinearity setting function; Perform the solution calculation to obtain the mechanical response results, including the stress, deformation, and support reaction of the glass panel.

8. The parametric finite element simulation calculation method for point-supported glass curtain walls to improve design efficiency according to claim 1, characterized in that, In step S6, the finite element method (FEM) calculation is performed, including: In response to a change in any variable value in the input parameter set, a threshold judgment is made on the amount of change. If the change exceeds the preset engineering tolerance threshold, the current finite element model is automatically cleared, and steps S2 to S6 are re-executed based on the changed parameters. If the change is within the tolerance threshold range, the current finite element model remains unchanged.

9. The parametric finite element simulation calculation method for point-supported glass curtain walls to improve design efficiency according to claim 1, characterized in that, In step S7, stress distribution, deformation, and support reaction data are extracted from the finite element solution results, and corresponding contour plots and numerical reports are generated, including: By using the node-based cloud map display function, stress component parameters can be specified to generate and output stress distribution cloud maps of the glass panel. By using the node-based cloud map display function, specifying displacement component parameters, the deformation result cloud map of the glass panel is generated and output; Using the reaction force solution list function, filter the support point nodes with applied displacement constraints, calculate and list the support reaction force data.

10. A storage medium storing computer-readable instructions, characterized in that, When the computer-readable instructions are executed by one or more processors, the one or more processors perform the steps of a parametric finite element simulation calculation method for point-supported glass curtain walls that improves design efficiency, as described in any one of claims 1 to 9.

Images