A CAE mechanical simulation method and device based on a BIM information model and a medium

Abstract

The application discloses a CAE mechanical simulation method and device based on a BIM information model and a medium, relates to the technical field of engineering simulation, and comprises the following steps: reading the BIM information model, extracting component basic mechanical correlation information of a main bearing component, generating a component semantic table, a structure correlation table and a cross-section transition table; determining a main force transmission path based on the component semantic table, combining the structure correlation table to connect local structure areas in series, generating a structure buffer propagation chain and a chain segment order table; calculating a chain segment buffer release ratio of each chain segment based on the structure buffer propagation chain, the chain segment order table and the cross-section transition table, and calibrating a chain segment type and a corresponding analysis expression mode through the chain segment buffer release ratio; constructing a CAE analysis model based on the corresponding analysis expression mode of each chain segment, performing mechanical solving, extracting a local high-response area, and combining a propagation chain form to determine a real danger area. The application improves the force flow recognition accuracy, the local mechanical behavior expression capability and the danger area determination reliability.

Description

Technical Field

[0001] This invention relates to the field of engineering simulation technology, and in particular to a CAE mechanical simulation method, equipment and medium based on BIM information model. Background Technology

[0002] With the development of digital design and structural analysis technologies for building engineering, CAE mechanical simulation based on BIM information models has become an important technical approach in the field of engineering simulation. Conventional methods typically involve first extracting the geometry, materials, connections, and boundary information of components from the BIM model, and then establishing corresponding finite element or other CAE analysis models to solve for the displacement, stress, and strain responses of the structure under load. This assists in completing structural stress analysis, identifying local weak points, and design verification, and has been widely applied to the engineering analysis of building structures, prefabricated components, and complex nodes.

[0003] However, under complex load-bearing systems, the above-mentioned conventional methods usually still have two problems: On the one hand, although the local structure, cross-sectional changes and component connection relationships in the BIM model can be extracted, CAE modeling often lacks an orderly organization for the main force transmission process, resulting in insufficient utilization of local mechanical correlation information; on the other hand, the simulation expression methods for different structural sections are usually relatively uniform, making it difficult to fully reflect the differentiated mechanical characteristics such as local diffusion, transition and convergence, thus affecting the refinement and reliability of CAE mechanical simulation results. Summary of the Invention

[0004] In view of the aforementioned existing problems, the present invention is proposed.

[0005] Therefore, this invention provides a CAE mechanical simulation method based on BIM information model to solve the problems of insufficient utilization of local mechanical correlation information and insufficient simulation expression capability of chain segment differences in existing technologies.

[0006] To solve the above-mentioned technical problems, the present invention provides the following technical solution: In a first aspect, the present invention provides a CAE mechanical simulation method based on a BIM information model, comprising: reading the BIM information model, extracting the component foundation mechanical correlation information of the main load-bearing component, and generating a component semantic table, a structural correlation table, and a section transition table; determining the main force transmission path based on the component semantic table, and combining the structural correlation table to connect local structural areas in sequence along the main force transmission path to generate a structural buffer propagation chain and a chain segment sequence table; calculating the chain segment mitigation ratio of each chain segment based on the structural buffer propagation chain, the chain segment sequence table, and the section transition table, and calibrating the chain segment type and corresponding analysis expression method through the chain segment mitigation ratio; constructing a CAE analysis model based on the corresponding analysis expression method of each chain segment, performing mechanical solution, extracting local high-response areas, and determining the propagation chain morphology according to the chain segment type and arrangement order; matching the local high-response areas, the main force transmission path, the structural buffer propagation chain, and the propagation chain morphology, and determining the actual danger zone according to different scenarios.

[0007] As a preferred embodiment of the CAE mechanical simulation method based on BIM information model described in this invention, the generation of component semantic table, structural association table, and cross-section transition table includes: reading the BIM information model file corresponding to the target project; extracting the component category, spatial location, material properties, cross-section properties, connection relationships, support relationships, and load attribution relationships of the main load-bearing components to generate a component semantic table; identifying the local structural areas corresponding to the main load-bearing components and generating a structural association table based on the association relationship between the local structural areas and the main load-bearing components; identifying cross-section change sections along the main direction of the main load-bearing components and recording the start and end positions of the cross-section transition and the corresponding change parameters to generate a cross-section transition table.

[0008] As a preferred embodiment of the CAE mechanical simulation method based on BIM information model described in this invention, the buffer propagation chain includes: determining candidate load-bearing components according to load attribution relationships, and determining candidate support components according to support relationships; starting from the candidate load-bearing components, selecting the next level of main force-transmitting components step by step according to the component connection relationships until reaching the candidate support components, thereby generating a main force transmission path; obtaining the effective local structural areas of each main load-bearing component on the main force transmission path, and sorting them according to the order along the main force transmission path; and sequentially connecting the effective local structural areas that maintain continuous buffer association to generate a structural buffer propagation chain.

[0009] As a preferred embodiment of the CAE mechanical simulation method based on BIM information model described in this invention, the chain segment sequence table includes: arranging the effective local structural areas in ascending order of their projection starting positions on the main force transmission path; comparing the interval distance between two adjacent effective local structural areas; when the interval distance is not greater than the characteristic size of the component's local cross-section, determining that the two adjacent effective local structural areas maintain a continuous buffer connection on the main force transmission path and connecting them in series; when the interval distance is greater than the characteristic size of the component's local cross-section, determining that the buffer connection between two adjacent effective local structural areas is interrupted, and using the interruption position as the termination position of a structural buffer propagation chain; defining the main direction segment between two adjacent effective local structural areas in the same structural buffer propagation chain as a chain segment; recording the sequence number of each chain segment to generate a chain segment sequence table.

[0010] As a preferred embodiment of the CAE mechanical simulation method based on BIM information model described in this invention, the step of calibrating the segment type and corresponding analysis expression method by the segment release ratio includes: for each segment, calculating the structural unfolding scale, material continuity scale, geometric abrupt change scale, and stress contraction scale; calculating the segment release ratio of the segment based on the structural unfolding scale, material continuity scale, geometric abrupt change scale, and stress contraction scale; when the segment release ratio is not greater than a first release threshold, the current segment is calibrated as a convergent segment; when the segment release ratio is greater than the first release threshold but not greater than a second release threshold, the current segment is calibrated as a transition segment; when the segment release ratio is greater than the second release threshold, the current segment is calibrated as a diffusion segment; and according to the segment type, the connection area and transition area of ​​the corresponding segment are calibrated as standard connection expression, local transition expression, or local diffusion expression, respectively.

[0011] As a preferred embodiment of the CAE mechanical simulation method based on BIM information model described in this invention, the calculation of the structural unfolding scale, material continuity scale, geometric abrupt change scale, and stress contraction scale includes: determining the structural unfolding scale based on the correspondence between the local structural regions at both ends of the chain segment and the effective force transmission width in the middle of the chain segment; determining the material continuity scale based on the continuity of material parameters in the local structural regions at both ends of the chain segment; determining the geometric abrupt change scale based on the degree of cross-sectional transition characterized by the thickness change, effective force transmission width change, and cross-sectional centroid offset within the chain segment; and determining the stress contraction scale based on the change in effective force transmission width between the inlet and outlet cross-sections of the chain segment.

[0012] As a preferred embodiment of the CAE mechanical simulation method based on BIM information model described in this invention, the determination of the propagation chain morphology includes: constructing a continuous analysis skeleton according to the component connection sequence on the main force transmission path; constructing a CAE analysis model for the corresponding chain segment range using standard connection expression, local transition expression, or local diffusion expression based on the analysis expression method corresponding to each chain segment; performing mechanical solution and extracting local high-response regions; and determining the propagation chain morphology as a peak sensitive chain, a fluctuating transition chain, or a stable diffusion chain based on the arrangement order of chain segment types in the same structural buffer propagation chain.

[0013] As a preferred embodiment of the CAE mechanical simulation method based on BIM information model described in this invention, the determination of the true hazard zone includes: projecting the local high-response area along the main force transmission path direction and matching its position with each structural buffer propagation chain to determine the target structural buffer propagation chain and chain segment; when the target structural buffer propagation chain is a peak-sensitive chain, determining the true hazard zone or pseudo-hazard zone based on whether the local high-response area is confined to the convergence segment and combined with local CAE analysis; when the target structural buffer propagation chain is a fluctuating transition chain, determining the true hazard zone or propagation offset type pseudo-hazard zone based on the consistency between the peak position of the area and the fluctuating boundary position; when the target structural buffer propagation chain is a stable diffusion chain, determining the true hazard zone or diffusion discontinuity type pseudo-hazard zone based on the continuous distribution of the local high-response area within the diffusion segment.

[0014] In a second aspect, the present invention provides a computer device including a memory and a processor, wherein the memory stores a computer program, wherein when the computer program is executed by the processor, it implements any step of the CAE mechanical simulation method based on BIM information model as described in the first aspect of the present invention.

[0015] Thirdly, the present invention provides a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements any step of the CAE mechanical simulation method based on BIM information model as described in the first aspect of the present invention.

[0016] The beneficial effects of this invention are as follows: by generating a structural buffer propagation chain and a chain segment sequence table, an orderly expression of the actual force transmission process and the local structural buffering effect is realized, which improves the accuracy of force flow identification; by calculating the chain segment release ratio and calibrating the chain segment type and corresponding analysis expression method, differentiated modeling of the local mechanical behavior of different chain segments is realized, which improves the accuracy of CAE analysis. Attached Figure Description

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

[0018] Figure 1 This is a flowchart of a CAE mechanical simulation method based on a BIM information model.

[0019] Figure 2 This is a flowchart of the generation of the main force transmission path and the construction of a buffer propagation chain.

[0020] Figure 3 This is a flowchart for calculating the chain segment sustained-release ratio and labeling the chain segment type.

[0021] Figure 4 This is a flowchart for identifying the propagation chain morphology and determining the actual danger zone. Detailed Implementation

[0022] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings.

[0023] Many specific details are set forth in the following description in order to provide a full understanding of the invention. However, the invention may also be practiced in other ways different from those described herein, and those skilled in the art can make similar extensions without departing from the spirit of the invention. Therefore, the invention is not limited to the specific embodiments disclosed below.

[0024] Secondly, the term "one embodiment" or "embodiment" as used herein refers to a specific feature, structure, or characteristic that may be included in at least one implementation of the present invention. The phrase "in one embodiment" appearing in different places in this specification does not necessarily refer to the same embodiment, nor is it a single or selective embodiment that is mutually exclusive with other embodiments.

[0025] Reference Figures 1-4 This is one embodiment of the present invention, which provides a CAE mechanical simulation method based on a BIM information model, including the following steps: S1. Read the BIM information model, extract the component foundation mechanical association information of the main load-bearing components, and generate the component semantic table, structural association table and cross-section transition table.

[0026] Read the BIM information model file corresponding to the target project.

[0027] The BIM information model file includes at least component object data, material object data, section object data, connection relationship data, support relationship data, load attribution data, and opening and local structural data.

[0028] After importing the BIM information model file into the parsing environment, the project reference coordinate system, floor reference elevation, and component instance number are read. The spatial coordinates of all components are then standardized so that all components are expressed in the same engineering coordinate system.

[0029] According to the component instance number, each component's category identifier, geometric outline, material name, material parameters, cross-sectional dimensions, connection location, support location, load location, and local structural location are read one by one to form the original component record set.

[0030] Furthermore, the primary load-bearing components are screened on the original component record set.

[0031] Specifically, components categorized as beams, columns, supports, slabs, walls, foundations, and load-bearing connectors are considered as candidate load-bearing components.

[0032] For each candidate load-bearing component, read the load-bearing attribute marker, upstream connection object, and downstream connection object. If the candidate load-bearing component has at least one of the following: force transfer relationship, support relationship, and load attribution relationship, it is determined to be the main load-bearing component.

[0033] Assign a unique component code to the determined main load-bearing component and write it into the main load-bearing component set.

[0034] Furthermore, the structural mechanical semantic information of each main load-bearing component is extracted one by one.

[0035] For linear main load-bearing components, after reading the starting point coordinates, ending point coordinates, and cross-sectional outer contour, the initial axis is determined by connecting the starting point coordinates and the ending point coordinates. Then, the initial axis is checked for eccentricity by combining the center of the cross-sectional outer contour to obtain the component axis. The length parameter is established along the component axis direction with the starting point of the component axis as the origin.

[0036] For planar main load-bearing components, after reading the outer boundary, thickness direction, and mid-surface position, the outer boundary is closed and the mid-surface contour is extracted. The in-plane main extension direction is determined based on the degree of expansion of the mid-surface contour in two orthogonal directions, and the normal direction is established in combination with the thickness direction.

[0037] For foundation and load-bearing connectors, after reading the outer contour, contact boundary and thickness dimensions, first identify the upstream contact boundary and downstream contact boundary, then establish a local reference direction based on the direction from the center of the upstream contact boundary to the center of the downstream contact boundary, and establish a thickness reference direction in combination with the thickness dimensions.

[0038] The component code, component category, spatial location, material properties, cross-sectional properties, connection relationships, support relationships, and load attribution relationships of each main load-bearing component are organized to generate a component semantic table.

[0039] Furthermore, local structural regions are identified.

[0040] Using the connection position, support transition position, opening adjacent position, and cross-sectional change position of each main load-bearing component as the center, local structural data is read within a search range (e.g., 450mm) around the component.

[0041] The local structural data includes at least the node expansion zone, the local thickening zone, the stiffening structural zone, the transition zone at the tunnel entrance edge, the support transition layer, and the connecting transition section.

[0042] Project each local structural area onto the main direction of the corresponding main load-bearing component to obtain the projection start point, projection end point, and projection coverage length of the local structural area.

[0043] The shortest distance from the center of the local structural zone to the axis of the main load-bearing member is calculated as the deviation distance. The material type and elastic modulus difference between the local structural zone and the main load-bearing member are read. The local structural zone code, main load-bearing member code, associated position, deviation distance and material difference are recorded accordingly to form a candidate associated record of the local structural zone.

[0044] To filter out local structural regions with practical impact on mechanical simulation from the candidate correlation records of local structural regions, the structural correlation degree between the local structural regions and the main load-bearing members is calculated, and the expression is as follows: ; in, Indicated in the main load-bearing component Length position Local structural region With main load-bearing components The degree of structural correlation between them Indicates axial coverage. Indicates the cross-sectional proximity ratio. Indicates the center deviation rate. This represents the stiffness mutation rate.

[0045] It should be noted that the axial coverage rate is calculated by the ratio of the projected coverage length of the local structural zone in the main direction of the main load-bearing member to the length of the current local reference section; the section proximity rate is calculated by the ratio of the effective projected width of the local structural zone in the current section to the effective force transmission width of the section; the center deviation rate is calculated by the ratio of the shortest distance from the center of the local structural zone to the reference axis of the main load-bearing member to the effective force transmission width of the section; and the stiffness mutation rate is calculated by the ratio of the absolute value of the difference between the axial stiffness of the local structural zone and the axial stiffness of the main load-bearing member to their sum.

[0046] It should be noted that the length of the local reference section is determined by taking the projection position of the center of the local structural area onto the main direction of the main load-bearing member as the center, and then merging the lengths of the local cross-section feature dimensions on both sides along the main direction; the effective force transmission width is obtained by extracting the width of the solid outline of the main load-bearing member participating in force transmission within the current cross-section, and deducting the remaining width after deducting holes, gaps, and non-load-bearing edges; the effective projection width is determined by projecting the local structural area into the current cross-section along a direction perpendicular to the main direction of the main load-bearing member, and reading the projection span within the effective force transmission range of the main load-bearing member.

[0047] when When the threshold for determining the association is not less than the threshold for determining the local structure region, the local structure region is determined. With main load-bearing components If a valid construction association exists, the corresponding component code, local construction area code, association location, and association type will be written into the construction association table.

[0048] It should be noted that the threshold for determining the structural association is determined by statistically analyzing the distribution boundary between the valid and invalid structural association samples that have been manually confirmed in historical finite element modeling projects, and taking the lowest structural association degree within the distribution boundary that can distinguish between valid and invalid structural associations. The preferred value range is [0.45, 0.65].

[0049] Furthermore, a cross-sectional transition table is generated.

[0050] Specifically, the cross-sections are discretized along the main direction of each main load-bearing component at fixed step lengths (e.g., 50 mm) to obtain multiple discrete cross-sections.

[0051] At each discrete cross-section, the effective force transmission width, effective thickness, and the position of the cross-section centroid are read.

[0052] Compare the effective force transmission width difference, effective thickness difference, and centroid offset of two adjacent discrete sections. When any change exceeds the corresponding transition judgment condition, mark the adjacent section as a section transition section, and record the section transition start position, section transition end position, thickness change, width change, and centroid offset to generate a section transition table.

[0053] It should be noted that the thickness change determination criterion is determined by statistically analyzing the thickness difference sequence of adjacent discrete sections of similar main load-bearing components within the geometrically stable section, and taking the maximum thickness difference that appears continuously and stably in the thickness difference sequence (e.g., the corresponding change amount in at least 3 consecutive adjacent discrete sections); the width change determination criterion is determined by statistically analyzing the effective force transmission width difference sequence of adjacent discrete sections of similar main load-bearing components within the geometrically stable section, and taking the maximum width difference that appears continuously and stably in the effective force transmission width difference sequence; the centroid offset determination criterion is determined by statistically analyzing the centroid offset distance sequence of adjacent discrete sections of similar main load-bearing components within the geometrically stable section, and taking the maximum offset distance that appears continuously and stably in the centroid offset distance sequence.

[0054] S2. Based on the component semantic table, determine the main force transmission path, and combine it with the structural association table to connect the local structural areas in sequence along the main force transmission path to generate the structural buffer propagation chain and the chain segment sequence table.

[0055] The main load-bearing components with load attribution relationships are identified as candidate components for load action.

[0056] The main load-bearing components with a support relationship are identified as candidate support components.

[0057] For each candidate component subject to load, read the load assignment location and use the component corresponding to the load assignment location as the current starting component.

[0058] Furthermore, starting from the current starting component, all adjacent main load-bearing components that have a direct connection with the current starting component are read as downstream candidate components.

[0059] For each downstream candidate component, extract the output direction of the current starting component pointing to the connection position, the input direction of the downstream candidate component extending from the connection position, the approximation distance of the downstream candidate component to the nearest supporting candidate component, and the center offset of the contact boundary between the current starting component and the downstream candidate component at the connection position.

[0060] The main force continuity determination value is calculated based on the output direction, input direction, approximation distance, and contact boundary center offset. The expression is: ; in, Indicates the current starting component With downstream candidate components The continuous determination value of the main transmitted force between them Indicates the directional matching rate. Indicates the support approximation ratio. Indicates the contact deviation rate. This represents the shunting deflection rate.

[0061] It should be noted that the directional matching rate is determined by the absolute value of the cosine of the angle between the output direction of the current starting component and the input direction of the downstream candidate component; the support approximation rate is determined by reading the path distance from the center of the current starting component to the center of the nearest support candidate component, and the path distance from the center of the downstream candidate component to the center of the same support candidate component, respectively, and calculating the ratio of the difference between these two distances to the path distance from the center of the current starting component to the center of the nearest support candidate component. When the path distance from the center of the downstream candidate component to the center of the same support candidate component is less than the path distance from the center of the current starting component to the center of the nearest support candidate component, the support approximation rate increases with the increase of the distance difference; the contact deviation rate is determined by the ratio of the distance between the center of the contact boundary of the current starting component and the downstream candidate component at the connection position to the effective contact span of the current connection position; the shunting deflection rate is determined by reading the output direction of the current starting component pointing to all downstream candidate components, calculating the angular deviation between each output direction and the average output direction, and taking the sum of the maximum and minimum values ​​of the angular deviation. The ratio is determined such that when there is only one downstream candidate component, the shunting deflection rate is 0.

[0062] It should be noted that the effective contact span is determined by reading the boundary length or contact width of the contact points at the current connection location.

[0063] When there are multiple downstream candidate components, the downstream candidate component with the largest continuous main force transmission value is selected as the next-level main force transmission component of the current starting component.

[0064] When the maximum main force transmission value corresponds to two or more downstream candidate components, the component with the shortest approach distance to the support candidate component is selected as the next level main force transmission component.

[0065] Connect the current starting component and the next selected main force transmission components in series according to the connection sequence until the support candidate component is reached, and obtain the main force transmission path from the load application position to the support position.

[0066] The main force transmission path is repeatedly obtained for all candidate components subjected to loads. When there are multiple main force transmission paths, the main force transmission path with the largest corresponding load range is retained as the target main force transmission path for the current load application location.

[0067] It should be noted that the load attribution range is determined by reading the load attribution boundary of the candidate component currently subjected to load in the BIM information model, and then projecting the load attribution boundary onto the force-bearing surface where the starting component of the main force transmission path is located.

[0068] After obtaining the target main force transmission path, read all effective local structural areas corresponding to each main load-bearing component on the construction association table and the target main force transmission path.

[0069] Extract the projection start point, projection end point, and associated location type of each effective local structural region in the main direction of the target main force transmission path.

[0070] All valid local structural regions are sorted in ascending order of their projection starting positions to obtain a sequence of local structural regions.

[0071] Furthermore, for two adjacent effective local structural regions in the sequence of local structural regions, the interval distance between the projection endpoint of the previous effective local structural region and the projection starting point of the next effective local structural region is compared.

[0072] When the interval distance is not greater than the characteristic size of the local section of the component, it is determined that two adjacent effective local structural zones maintain a continuous buffer relationship on the main force transmission path.

[0073] When the distance between two adjacent effective local structural zones is greater than the characteristic size of the local section of the component, the buffer association between them is determined to be interrupted, and the interruption position is taken as the termination position of a structural buffer propagation chain.

[0074] It should be noted that the local section feature dimensions are determined by reading the effective force transmission width and effective thickness of the main load-bearing member section at the current associated position on the main force transmission path, and taking the minimum value of the two.

[0075] All effective local structural regions that maintain continuous buffer associations are sequentially connected in series to generate a structural buffer propagation chain.

[0076] Furthermore, the main directional segment between two adjacent effective local structural regions in the construction buffer propagation chain is defined as a chain segment. The upstream local structural region code, downstream local structural region code, chain segment start position, chain segment end position, and chain segment sequence number of each chain segment are recorded to generate a chain segment sequence table.

[0077] Among them, the main force transmission path is used to characterize the main transmission sequence of loads between the main load-bearing components; the structural buffer propagation chain is used to characterize the distribution of local structural zones around the main force transmission path that continuously participate in force diffusion and stiffness transition; and the chain segment sequence table is used to characterize the sequential relationship of each chain segment in the structural buffer propagation chain.

[0078] S3. Calculate the chain segment release ratio of each chain segment based on the constructed buffer propagation chain, chain segment sequence table, and cross-sectional transition table, and calibrate the chain segment type and corresponding analysis expression method through the chain segment release ratio.

[0079] Furthermore, the structural unfolding scale, material continuity scale, geometric abrupt change scale, and stress contraction scale of the current chain segment are calculated.

[0080] Specifically, the structural development scale is determined by reading the effective projected widths of the upstream and downstream local structural regions within their corresponding cross sections and taking the ratio of the maximum value of the two to the effective force transmission width at the middle position of the current chain segment.

[0081] The material continuity scale is determined by reading the elastic modulus of the upstream local structural region material and the elastic modulus of the downstream local structural region material, and by the ratio of their geometric mean to their arithmetic mean.

[0082] By reading the thickness change, effective force transmission width change, and cross-sectional centroid offset of the current chain segment, and normalizing them with the effective thickness, effective force transmission width, and effective force transmission width of the upstream position of the current chain segment, the maximum value of the three is taken to determine the geometric change scale.

[0083] The force contraction scale is determined by reading the effective force transmission width of the inlet section and the effective force transmission width of the outlet section of the current chain segment. When the effective force transmission width of the inlet section is greater than the effective force transmission width of the outlet section, the ratio of the effective force transmission width of the inlet section to the effective force transmission width of the outlet section is used as the force contraction scale. When the effective force transmission width of the inlet section is not greater than the effective force transmission width of the outlet section, the force contraction scale is 1.

[0084] It should be noted that the inlet section is the discrete section corresponding to the starting position of the chain segment, and the outlet section is the discrete section corresponding to the ending position of the chain segment. When there are multiple section transition records within the chain segment, the thickness change, the effective force transmission width change, and the section centroid offset are taken as the maximum values ​​in the corresponding records.

[0085] Furthermore, the chain segment release ratio is calculated based on the structural unfolding scale, material continuity scale, geometric abrupt change scale, and stress contraction scale, and the expression is: ; in, Indicates the first The first in the construction of the buffer propagation chain The segment-to-segment release ratio of each chain segment, Indicates the first The first in the construction of the buffer propagation chain The scale of the construction and expansion of each chain segment. Indicates the first The first in the construction of the buffer propagation chain The material continuity scale of each chain segment Indicates the first The first in the construction of the buffer propagation chain The geometric aberration scale of each chain segment Indicates the first The first in the construction of the buffer propagation chain The stress contraction scale of each chain segment.

[0086] It should be noted that the chain segment release ratio is used to characterize the ability of the current chain segment to support the diffusion of local response along the main force transmission path. Among them, the structural expansion scale and the material continuity scale jointly characterize the promoting effect of local structural regions on the expansion and continuous transmission of force flow, while the geometric abrupt change scale and the stress contraction scale jointly characterize the enhancing effect of geometric discontinuities and stress convergence within the chain segment on the concentration of local response. When the chain segment has a large degree of structural expansion and good material continuity, and at the same time has small geometric abrupt changes and weak stress contraction, the chain segment release ratio increases, thus indicating that the chain segment is more conducive to weakening the local high response along the main force transmission path.

[0087] Furthermore, a chain segment type label is performed for each chain segment.

[0088] When the chain segment release ratio is not greater than the first release threshold, the current chain segment is marked as the closure chain segment.

[0089] When the sustained-release ratio of a segment is greater than the first sustained-release threshold but not greater than the second sustained-release threshold, the current segment is designated as a transition segment.

[0090] When the sustained-release ratio of the chain segment is greater than the second sustained-release threshold, the current chain segment is designated as a diffusion chain segment.

[0091] It should be noted that the first sustained-release threshold is determined by reading the chain segment sustained-release ratio sequence corresponding to samples whose local high response range remains within a single chain segment in historical finite element modeling projects, and taking the largest chain segment sustained-release ratio in the chain segment sustained-release ratio sequence that still limits the local high response to a single chain segment, preferably in the range of [0.80, 1.10]. The second sustained-release threshold is determined by reading the chain segment sustained-release ratio sequence corresponding to samples whose local high response continuously crosses two adjacent chain segments in historical finite element modeling projects, and taking the smallest chain segment sustained-release ratio in the chain segment sustained-release ratio sequence that first causes the local high response to continuously cross two adjacent chain segments, preferably in the range of [1.20, 1.80].

[0092] Furthermore, the expression method of synchronous calibration chain segment analysis is also discussed.

[0093] When the chain segment type is a convergence chain segment, the connection region and transition region corresponding to the current chain segment are marked as standard connection expressions.

[0094] When the segment type is a transition segment, the connection region and transition region corresponding to the current segment are marked as local transition expressions.

[0095] When the segment type is a diffuse segment, the connection region and transition region corresponding to the current segment are marked as locally diffused expression.

[0096] The chain segment sequence number, chain segment sustained-release ratio, chain segment type, and chain segment analysis expression method are recorded accordingly, generating a chain segment sustained-release ratio table, a chain segment type table, and a chain segment analysis expression method table.

[0097] S4. Construct a CAE analysis model based on the corresponding analysis expression of each chain segment, perform mechanical solution, extract local high response regions, and determine the propagation chain morphology according to the chain segment type and arrangement order.

[0098] According to the connection sequence of components on the main force transmission path, extract the component code, component category, spatial location, cross-sectional properties and material properties of each main load-bearing component in sequence.

[0099] Based on the component category, extract the component axis of the linear main load-bearing component and establish a linear analysis component; extract the mid-surface of the planar main load-bearing component and establish a planar analysis component; and establish solid analysis components for the foundation and load-bearing connectors according to the solid contour.

[0100] Analysis nodes are established based on adjacent connection positions on the main force transmission path, and the main load-bearing components are connected into a continuous analysis skeleton using the analysis nodes as boundaries.

[0101] Furthermore, according to the chain segment sequence number in the chain segment sequence table, the chain segment analysis expression method corresponding to each chain segment is read one by one.

[0102] When the chain segment analysis is expressed as a standard connection expression, the original connection stiffness and original cross-sectional properties of the corresponding connection and transition regions of the chain segment remain unchanged.

[0103] When the chain segment analysis is expressed as a local transition expression, multiple transition sections are inserted at fixed intervals (e.g., 20mm when the chain segment length is 120mm) along the main force transmission path between the chain segment start position and the chain segment end position. Based on the thickness change and effective force transmission width change in the section transition table, the effective thickness and effective force transmission width of each transition section are linearly interpolated according to the distance ratio from the chain segment start position to the chain segment end position. The local stiffness of each transition section is calculated based on the interpolated section area or section moment of inertia and then gradually assigned.

[0104] It should be noted that the transition section is a virtual analysis section located within the chain segment and perpendicular to the direction of the main force transmission path.

[0105] When the chain segment analysis is expressed as a local diffusion expression, the effective local force transmission range is extended between the chain segment start position and the chain segment end position, and the effective projection width of the upstream local structural region and the downstream local structural region is written into the corresponding analysis section, forming a local force transmission expression that extends continuously along the main force transmission path direction.

[0106] Furthermore, the continuous analysis skeleton is discretized into a mesh.

[0107] For linear analysis components, linear discrete segments are divided along the component axis with a fixed length step.

[0108] For planar analysis components, rectangular discrete blocks are divided along the main in-plane extension direction and the secondary in-plane extension direction.

[0109] For solid analysis components, the solid discrete volume is divided along three orthogonal directions.

[0110] When a discrete position falls within the connection zone, support transition zone, cross-sectional change zone, or the range of a chain segment expressed as a local transition or local diffusion expression, the discrete scale of the corresponding position is reduced to half of the ordinary discrete scale to improve the accuracy of local response recognition.

[0111] It should be noted that the ordinary discrete scale is determined by reading the local cross-sectional characteristic dimensions of the current main load-bearing component, preferably 0.5 to 1 times the local cross-sectional characteristic dimensions; the connection area, support transition area and cross-sectional change area are directly determined by the corresponding positions in the main force transmission path and cross-sectional transition table.

[0112] Furthermore, the support relationships and load attribution relationships in the component semantic table are read and written into the CAE analysis model.

[0113] Specifically, for positions with support relationships, the support type of the corresponding support boundary is read. When the support type is a fixed support, the translational and rotational degrees of freedom of the corresponding analysis node in three orthogonal directions are constrained. When the support type is a sliding support, the translational degrees of freedom perpendicular to the support surface are constrained, while the translational and rotational degrees of freedom along the support surface are retained.

[0114] For locations with load attribution relationships, first read the load magnitude, load direction, and load range corresponding to the load attribution boundary.

[0115] When the load assignment boundary corresponds to a linear analysis member, the load is written into the corresponding analysis section as a line load.

[0116] When a load is assigned to a surface analysis component corresponding to the boundary, the load is written into the corresponding analysis surface as a surface load.

[0117] When a load is assigned to a boundary element corresponding to a solid analysis component, the load is written to the corresponding solid boundary in the manner of boundary pressure.

[0118] The discretized mesh data is written into the CAE analysis model according to the correspondence between nodes and elements. Material parameters are assigned to the corresponding elements, boundary constraints are assigned to the corresponding analysis nodes, and load conditions are assigned to the corresponding analysis sections, analysis surfaces, or solid boundaries. Based on the node equilibrium relationship, mechanical solution is performed to obtain the displacement response values ​​of all discrete nodes and the stress response values ​​and strain response values ​​of the center positions of all discrete elements.

[0119] The mechanical solution based on nodal equilibrium relationships involves writing the discretized mesh data, material parameters, boundary constraints, and load conditions into the CAE analysis model according to the correspondence between nodes, elements, and boundary objects. The element stiffness matrix is ​​calculated for each discrete element, and the overall stiffness matrix is ​​assembled according to the node connection relationship. The overall stiffness matrix and load array are corrected by combining the boundary constraints. The nodal equilibrium relationship between the overall stiffness matrix, the nodal displacement array, and the load array is established, and the nodal equilibrium relationship is solved to obtain the displacement of each discrete node. Based on the displacement of each discrete node, the strain response value and stress response value at the center position of all discrete elements are calculated.

[0120] After obtaining all response values, the local response concentration ratio is calculated for each discrete location, expressed as: ; in, Indicates the first The local response concentration ratio at discrete locations Indicates the first The current response value at each discrete location. Indicates the first The median of the response values ​​within the neighborhood of a discrete location.

[0121] It should be noted that the current response value is preferentially taken as the equivalent stress value; when the equivalent stress value cannot be directly read, the absolute value of the principal strain or the displacement gradient value is taken; the neighborhood is determined by selecting a fixed number (e.g., 2) of consecutive adjacent discrete positions upstream and downstream along the main force transmission path with the current discrete position as the center.

[0122] Furthermore, when the local response concentration ratio is not less than the local high response determination threshold, the current discrete location is marked as a high response location.

[0123] Perform connectivity merging along the main force transmission path for all high-response locations to obtain local high-response regions.

[0124] It should be noted that the local high response judgment threshold is determined by statistically analyzing the local response concentration ratio sequence of the discrete positions corresponding to the geometrically stable sections in the completed finite element modeling projects in the past, and taking the maximum local response concentration ratio corresponding to the local response concentration ratio sequence in the continuous stable response section (e.g., maintaining a stable response state for 3 consecutive times). The preferred value range is [1.20, 1.80].

[0125] Furthermore, the propagation chain morphology is determined based on the chain segment type table and the chain segment sequence table.

[0126] When there are consecutive convergence segments in the same structural buffer propagation chain, and the coverage of the main direction corresponding to the consecutive convergence segments reaches the peak sensitivity determination condition, the structural buffer propagation chain is marked as a peak sensitive chain.

[0127] When the converging segment and the diffusion segment alternate along the chain segment sequence in the same structural buffer propagation chain, the structural buffer propagation chain is marked as a wave transition chain.

[0128] When there are continuous diffusion chain segments in the same structural buffer propagation chain, and the main direction coverage of the continuous diffusion chain segments meets the diffusion determination criteria, the structural buffer propagation chain is marked as a stable diffusion chain.

[0129] Record the corresponding propagation chain code, propagation chain morphology category, and chain segment sequence number to generate a propagation chain morphology table.

[0130] It should be noted that the peak sensitivity determination criterion is determined by statistically analyzing the sequence of continuous convergence chain segment lengths corresponding to peak sensitive chain samples in the completed finite element modeling project, and taking the smallest segment length in the segment length sequence that still allows the local high response to remain within the control range of the convergence chain segment; the diffusion determination criterion is determined by statistically analyzing the sequence of continuous diffusion chain segment lengths corresponding to stable diffusion chain samples in the completed finite element modeling project, and taking the smallest segment length in the segment length sequence that allows the local high response to continuously expand along the propagation chain direction.

[0131] Further, extract the location markers of key chain segments.

[0132] When the propagation chain is a spike-sensitive chain, the starting position of the first convergence segment and the ending position of the last convergence segment are read as the position markers of the convergence segment.

[0133] When the propagation chain is a wave transition chain, the boundary position between the adjacent convergence chain segment and the diffusion chain segment is read and used as the wave boundary position marker.

[0134] When the propagation chain is a stable diffusion chain, the start and end positions of the continuous diffusion chain segments are read as the location markers of the diffusion segments.

[0135] Record the propagation chain code, critical segment category, critical segment start position, and critical segment end position accordingly to generate critical segment position markers.

[0136] S5. Match the local high-response area, the main force transmission path, the construction of buffer propagation chain, and the propagation chain morphology to determine the real danger zone according to different situations.

[0137] For each local high-response region, extract the outer boundary contour, the region center position, and the region peak position. Then, project the outer boundary contour along the main force transmission path direction onto the main force transmission path to obtain the region's main direction coverage area and the corresponding projection start and end points.

[0138] Position matching is performed based on the overlap length between the projection start point and projection end point and the range of each structural buffer propagation chain segment. The structural buffer propagation chain with the largest overlap length is determined as the target structural buffer propagation chain, and the range of chain segments in the target structural buffer propagation chain that continuously overlap with the projection start point and projection end point is determined as the target chain segment.

[0139] Furthermore, the propagation chain morphology category of the target construction buffer propagation chain in the propagation chain morphology table is read, and the start position and end position of the corresponding key chain segment are read from the key chain segment position marker.

[0140] When the propagation chain morphology type is a spike-sensitive chain, it is further determined whether the current local high-response region only covers the chain segment range corresponding to the position marker of the convergence section.

[0141] When a local high-response region only covers the chain segment range corresponding to the location marker of the convergence section, the current local high-response region is recorded as a candidate for a restricted danger zone.

[0142] Read the analysis expression method of the candidate chain segment corresponding to the restricted danger area in the chain segment analysis expression method table.

[0143] When the analysis expression method of the candidate chain segment corresponding to the restricted hazard area is standard connection expression or local transition expression, the analysis expression method of the candidate chain segment corresponding to the restricted hazard area is switched to local diffusion expression, and a local CAE analysis sub-model is re-established only for the chain segment range where the candidate restricted hazard area is located.

[0144] The switched local diffusion expression is written into the local CAE analysis sub-model corresponding to the chain segment range where the candidate confined danger zone is located. The element stiffness matrix of the local CAE analysis sub-model is recalculated and assembled to form a local global stiffness matrix. The mechanical solution is re-executed in combination with the original boundary constraints and load conditions to obtain the displacement response values ​​of local discrete nodes and the stress response values ​​and strain response values ​​of the center position of local discrete elements. Based on the re-solved response values, high response positions are extracted and connected and merged to obtain the adjusted local high response region.

[0145] Furthermore, by comparing the main direction coverage length of the local high-response area before and after adjustment, the ratio of the coverage length of the candidate restricted hazard area in the main force transmission path direction after adjustment to the coverage length of the candidate restricted hazard area in the main force transmission path direction before adjustment is taken as the area shrinkage ratio.

[0146] When the area shrinkage ratio is not greater than the area shrinkage judgment threshold, the current restricted danger zone candidate is judged as a false danger zone.

[0147] When the area shrinkage ratio is greater than the area shrinkage judgment threshold, the current restricted hazard zone candidate is judged as a real hazard zone.

[0148] It should be noted that the regional shrinkage judgment threshold is determined by statistically analyzing the regional shrinkage ratio sequence of pseudo-hazardous area samples before and after local diffusion expression verification in historical finite element modeling projects, and taking the largest regional shrinkage ratio in the regional shrinkage ratio sequence that can maintain the stable shrinkage characteristics of the pseudo-hazardous area. The preferred value range is [0.60, 0.85].

[0149] Furthermore, when the propagation chain morphology is a wave transition chain, the wave boundary position marker is read, and the regional peak position of the current local high response region is extracted. The main direction distance between the regional peak position and the wave boundary position marker is compared, and the change in the chain segment release ratio within the corresponding chain segment of the current local high response region is used to determine whether the peak position is consistent with the wave boundary position.

[0150] Specifically, the distance difference between the peak position and the boundary of the fluctuation is extracted, and the distance difference is normalized by the coverage length of the current local high response area in the direction of the main force transmission path to obtain the position deviation rate.

[0151] Read the chain segment release ratio of adjacent chain segments on both sides of the wave boundary and compare the direction of change of the two.

[0152] When the peak position of the region falls within the neighborhood of the wave boundary, and the sustained release ratio of the chain segments on both sides of the peak position increases on one side and decreases on the other side, it is determined that the current local high response region is consistent with the wave boundary, and thus it is identified as a real danger zone.

[0153] When the peak position of the region is separated from the neighborhood of the fluctuation boundary, or when the slow release ratio of the chain segments on both sides of the peak position does not change alternately, it is identified as a propagation-offset pseudo-hazard zone.

[0154] It should be noted that the neighborhood of the wave boundary is determined by reading the effective force transmission width of the cross section corresponding to the wave boundary and expanding outwards to both sides of the main force transmission path based on the effective force transmission width.

[0155] Furthermore, when the propagation chain morphology is a stable diffusion chain, the location marker of the diffusion segment is read, and it is determined whether the current local high response region is continuously distributed along the location marker of the diffusion segment.

[0156] Specifically, the overlap length of the current local high response area within the diffusion segment and the total coverage length of the current local high response area are calculated respectively to obtain the diffusion overlap rate.

[0157] When the diffusion overlap rate is not less than the diffusion judgment threshold, and there is no interruption point along the main force transmission path in the current local high response area, the current local high response area is directly judged as a real danger zone.

[0158] When the diffusion overlap rate is less than the diffusion judgment threshold, or when there is an interruption in the current local high response area along the main force transmission path, the current local high response area is judged as a diffusion discontinuity type pseudo-hazard area.

[0159] It should be noted that the diffusion overlap rate is determined by the ratio of the coverage length of the current local high response region falling within the diffusion segment location mark range to the total coverage length of the current local high response region; the diffusion judgment threshold is determined by the diffusion overlap rate sequence corresponding to the stable diffusion chain samples in the historical finite element modeling projects, and the minimum diffusion overlap rate that can maintain the stability of the continuous diffusion response identification result in the diffusion overlap rate sequence is selected, with the preferred value range being [0.65, 0.90].

[0160] Furthermore, all real danger zones are classified into danger levels.

[0161] Specifically, the peak stress response value, peak strain response value, and area coverage length within each actual danger zone are read. The allowable stress of the material in the main load-bearing component where the actual danger zone is located is read. The ratio of the peak stress response value to the allowable stress of the material is calculated. The ratio is compared with the threshold of the danger level range. When the ratio falls into the lower danger level range, it is determined to be a low danger level. When the ratio falls into the middle danger level range, it is determined to be a medium danger level. When the ratio falls into the higher danger level range, it is determined to be a high danger level.

[0162] It should be noted that the allowable stress of a material refers to the upper limit of the stress that can be used in the structural design for the material corresponding to the main load-bearing component where the actual danger zone is located.

[0163] It should be noted that the hazard level intervals are obtained by reading the stress utilization rate sequences corresponding to low-hazard, medium-hazard, and high-hazard samples from historical finite element modeling projects, and by statistically analyzing the concentrated distribution intervals of stress utilization rates for the three types of samples. The concentrated distribution intervals corresponding to low-hazard samples are determined as low-hazard level intervals, the concentrated distribution intervals corresponding to medium-hazard samples are determined as medium-hazard level intervals, and the concentrated distribution intervals corresponding to high-hazard samples are determined as high-hazard level intervals.

[0164] This embodiment also provides a computer device applicable to the CAE mechanical simulation method based on BIM information model, including: a memory and a processor; the memory is used to store computer-executable instructions, and the processor is used to execute the computer-executable instructions to realize the CAE mechanical simulation method based on BIM information model as proposed in the above embodiment.

[0165] The computer device can be a terminal, comprising a processor, memory, communication interface, display screen, and input devices connected via a system bus. The processor provides computing and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system and computer programs. The internal memory provides an environment for the operation of the operating system and computer programs stored in the non-volatile storage media. The communication interface is used for wired or wireless communication with external terminals; wireless communication can be achieved through Wi-Fi, carrier networks, NFC (Near Field Communication), or other technologies. The display screen can be an LCD screen or an e-ink screen. The input devices can be a touch layer covering the display screen, buttons, a trackball, or a touchpad on the computer device's casing, or an external keyboard, touchpad, or mouse.

[0166] This embodiment also provides a storage medium storing a computer program. When executed by a processor, the program implements the CAE mechanical simulation method based on the BIM information model as proposed in the above embodiments. The storage medium can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as Static Random Access Memory (SRAM), Electrically Erasable Programmable Read-Only Memory (EEPROM), Erasable Programmable Read Only Memory (EPROM), Programmable Red-Only Memory (PROM), Read-Only Memory (ROM), magnetic storage, flash memory, magnetic disk, or optical disk.

[0167] In summary, this invention achieves an ordered expression of the actual force transmission process and the local structural buffering effect by generating a structural buffer propagation chain and a chain segment sequence table, thereby improving the accuracy of force flow identification. By calculating the chain segment release ratio and calibrating the chain segment type and corresponding analysis expression method, it achieves differentiated modeling of the local mechanical behavior of different chain segments, thereby improving the accuracy of CAE analysis.

[0168] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit it. Although the present invention has been described in detail with reference to preferred embodiments, those skilled in the art should understand that modifications or equivalent substitutions can be made to the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, and all such modifications or substitutions should be covered within the scope of the claims of the present invention.

Claims

1. A CAE mechanical simulation method based on BIM information model, characterized in that, include: Read the BIM information model, extract the component foundation mechanical association information of the main load-bearing components, and generate component semantic table, structural association table and section transition table; The main force transmission path is determined based on the component semantic table, and local structural regions are connected in sequence along the main force transmission path in combination with the structural association table to generate structural buffer propagation chain and chain segment sequence table. The chain segment release ratio of each chain segment is calculated based on the constructed buffer propagation chain, chain segment sequence table and cross-sectional transition table, and the chain segment type and corresponding analysis expression method are identified by the chain segment release ratio. A CAE analysis model is constructed based on the corresponding analysis expression of each chain segment, mechanical solution is performed, local high response regions are extracted, and the propagation chain morphology is determined according to the chain segment type and arrangement order. By matching local high-response areas, main force transmission paths, construction of buffer propagation chains, and propagation chain morphology, the actual danger zone is determined according to different scenarios.

2. The CAE mechanical simulation method based on BIM information model as described in claim 1, characterized in that, The generated component semantic table, construction association table, and section transition table include: Read the BIM information model file corresponding to the target project; Extract the component category, spatial location, material properties, cross-sectional properties, connection relationships, support relationships, and load assignment relationships of the main load-bearing components to generate a component semantic table; Identify the local structural regions corresponding to the main load-bearing components, and generate a structural association table based on the relationship between the local structural regions and the main load-bearing components; Identify the cross-sectional change sections along the main direction of the main load-bearing member, record the start and end positions of the cross-sectional transition and the corresponding change parameters, and generate a cross-sectional transition table.

3. The CAE mechanical simulation method based on BIM information model as described in claim 1 or 2, characterized in that, The buffer propagation chain includes: Candidate members for load action are determined based on load attribution relationships, and candidate members for support are determined based on support relationships. Starting with the candidate component under load, the next level of main force transmission component is selected step by step according to the component connection relationship until the support candidate component is reached, thus generating the main force transmission path. Obtain the effective local structural areas of each main load-bearing component along the main force transmission path and sort them according to the order along the main force transmission path; By sequentially connecting effective local structural regions that maintain continuous buffer associations, a structural buffer propagation chain is generated.

4. The CAE mechanical simulation method based on BIM information model as described in claim 3, characterized in that, The chain segment sequence list includes: Arranged in ascending order of the starting position of the effective local structural zone's projection on the main force transmission path; Compare the interval distance between two adjacent effective local structural regions; When the interval is not greater than the characteristic dimension of the local section of the component, it is determined that two adjacent effective local structural zones maintain a continuous buffer connection on the main force transmission path and are connected in series; When the interval distance is greater than the local cross-sectional characteristic size of the component, the buffer association between two adjacent effective local structural zones is determined to be interrupted, and the interruption position is taken as the termination position of a structural buffer propagation chain. The main directional segment between two adjacent effective local structural regions in the same structural buffer propagation chain is defined as a chain segment. Record the sequential number of each chain segment and generate a chain segment sequence table.

5. The CAE mechanical simulation method based on BIM information model as described in claim 1 or 4, characterized in that, The method of calibrating segment types and corresponding analytical expressions by segment sustained-release ratio includes: For each chain segment, calculate the structural unfolding scale, material continuity scale, geometric abrupt change scale, and stress contraction scale; The chain segment release ratio is calculated based on the structural unfolding scale, material continuity scale, geometric abrupt change scale, and stress contraction scale. When the chain segment release ratio is not greater than the first release threshold, the current chain segment is marked as the closure chain segment; When the sustained-release ratio of the chain segment is greater than the first sustained-release threshold but not greater than the second sustained-release threshold, the current chain segment is marked as a transitional chain segment. When the sustained-release ratio of the chain segment is greater than the second sustained-release threshold, the current chain segment is marked as a diffused chain segment; Based on the segment type, the connection region and transition region of the corresponding segment are labeled as standard connection expression, local transition expression, or local diffusion expression.

6. The CAE mechanical simulation method based on BIM information model as described in claim 5, characterized in that, The calculated structural unfolding scale, material continuity scale, geometric abrupt change scale, and stress contraction scale include: The structural development scale is determined based on the correspondence between the local structural regions at both ends of the chain segment and the effective force transmission width in the middle of the chain segment. The material continuity scale is determined based on the continuity of material parameters in the local structural regions at both ends of the chain segment; The geometric abrupt change scale is determined based on the degree of cross-sectional transition characterized by changes in chain segment thickness, effective force transmission width, and cross-sectional centroid offset. The force contraction scale is determined based on the change in the effective force transmission width between the inlet and outlet sections of the chain segment.

7. The CAE mechanical simulation method based on BIM information model as described in claim 3, characterized in that, The determination of the propagation chain form includes: Construct a continuous analysis skeleton according to the connection sequence of components on the main force transmission path; Based on the analysis and expression methods corresponding to each chain segment, a CAE analysis model is constructed for the corresponding chain segment range using standard connection expression, local transition expression, or local diffusion expression. Perform mechanical solutions and extract locally high-response regions; Based on the order of the chain segment types in the same structural buffer propagation chain, the propagation chain morphology is determined as a spike-sensitive chain, a fluctuating transition chain, or a stable diffusion chain.

8. The CAE mechanical simulation method based on BIM information model as described in claim 7, characterized in that, The determination of the actual danger zone includes: Project the local high-response region along the main force transmission path and match its position with each structural buffer propagation chain to determine the target structural buffer propagation chain and chain segment. When the target constructs a buffer propagation chain that is a spike-sensitive chain, the real danger zone or the false danger zone is determined based on whether the local high-response region is limited by the convergence segment and in combination with local CAE analysis. When the target constructs a buffer propagation chain that is a wave transition chain, the true danger zone or the propagation deviation type pseudo danger zone is determined based on the consistency between the peak position of the region and the boundary position of the wave. When the target constructs a buffer propagation chain that is a stable diffusion chain, the real danger zone or the pseudo-danger zone with diffusion discontinuity is determined based on the continuous distribution of local high-response areas within the diffusion segment.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the CAE mechanical simulation method based on BIM information model as described in any one of claims 1 to 8.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the steps of the CAE mechanical simulation method based on the BIM information model as described in any one of claims 1 to 8.

Images