A seismic analysis method for reinforced concrete frame structures based on Revit-OpenSEES

Abstract

A seismic analysis method for reinforced concrete frame structures based on Revit‑OpenSEES belongs to the field of civil engineering technology. Technical solution: Convert the finite element model from Revit's BIM model to OpenSEES, and perform seismic analysis settings and output settings; group, organize and extract the analysis data output by OpenSEES, and display and output them in the form of graphs, tables or documents; perform secondary analysis of seismic damage to the structure based on the concrete and steel strain data output by OpenSEES; automatically store key structural seismic damage information in the Revit structural model; display the seismic damage status through a three-dimensional view, or perform virtual roaming through VR software. Beneficial effects: The present invention realizes the integration of BIM commonly used software and finite element analysis programs, realizes seismic analysis and damage evaluation of reinforced concrete frame structures directly based on BIM models, and provides a solution for structural seismic analysis to participate in the collaborative work of BIM.

Description

Technical Field

[0001] The invention belongs to the technical field of civil engineering and relates to a seismic analysis method of a reinforced concrete frame structure based on Revit-OpenSEES. Background Art

[0002] At present, BIM technology is developing rapidly in the field of civil engineering at home and abroad, but the functions of various mainstream BIM software are limited. Generally, data conversion between various professional building software is required through the IFC standard to achieve BIM-based data sharing and collaborative work.

[0003] With the development of urbanization, earthquakes pose a more and more serious potential threat to economic development and people's lives, so the continuous development and improvement of building seismic theories and methods is particularly important. At present, the research and application of BIM-based seismic design and analysis of building structures are still in their infancy. The model conversion between BIM and structural seismic analysis software is still limited to the conversion of geometric information and simple material information, while vertical and horizontal loads, analysis units and other information still need to be set in the structural software. In addition, most research work has not yet achieved the integration of BIM and seismic analysis functions. How to integrate the design, analysis and evaluation functions of building structures into the BIM platform to achieve shared data and collaborative analysis directly based on BIM is a problem that needs to be solved in further research. Summary of the invention

[0004] In order to solve the technical problem of BIM integration of seismic analysis technology and realize seismic analysis and evaluation directly based on BIM building structure model, the present invention provides a seismic analysis method of reinforced concrete frame structure based on Revit-OpenSEES, which can realize the direct setting of material properties, component sections, analysis units, load conditions, analysis and output in Revit software, and call OpenSEES program to perform elastic-plastic seismic analysis of reinforced concrete frame structure; directly output the analysis results of OpenSEES program in the form of graphics, tables and data files in Revit software; display the seismic damage results of reinforced concrete frame structure in Revit software through three-dimensional model view and virtual roaming; and automatically store the key seismic damage information obtained by analysis in Revit software, and view the damage information corresponding to different earthquake levels through the structural component attribute table.

[0005] The technical solution is as follows:

[0006] A seismic analysis method for reinforced concrete frame structures based on Revit-OpenSEES, the steps are as follows:

[0007] S1. In the pre-analysis processing module, the finite element model is converted from Revit's BIM model to OpenSEES. After the model conversion is completed, seismic analysis and output settings are performed;

[0008] S2. In the post-analysis processing module, the analysis data output by OpenSEES is grouped, sorted and extracted, and displayed and output in the form of graphs, tables or documents;

[0009] S3. In the post-analysis processing module, a secondary seismic damage analysis is performed on the structure based on the concrete and steel strain data output by OpenSEES to determine the seismic damage level of each structural component that cannot be directly obtained by OpenSEES.

[0010] S4. In the post-analysis processing module, key structural earthquake damage information is automatically stored in the Revit structural model as BIM data for collaborative analysis and structural analysis and information management throughout the building life cycle;

[0011] S5. In the earthquake damage visualization display module, the earthquake damage status is displayed through a three-dimensional view, or virtual roaming is performed through VR software.

[0012] Further, the specific steps of step S1 are as follows:

[0013] Node conversion: extract the analysis nodes and corresponding coordinates in the Revit structural model, provide numbers through the numbering class, obtain the center coordinates of the entire floor and the center coordinates of each local floor separated by the beam as the dividing line;

[0014] Constraint setting: Filter the nodes that need to be subject to bottom constraints, and use the center points of each floor slab that have been obtained to define the floor slab constraints of each floor;

[0015] Material setting: Develop and set functions in Revit according to the material mechanical information required for OpenSEES analysis. When setting materials, read the material type in the Revit structural model, and then match the set mechanical parameters to the corresponding materials in the model according to the material type;

[0016] Section setting: For the steel fiber of the section, the coordinates and material properties of the steel bar are read in the Revit model and converted into the corresponding TCL command stream; for concrete fiber, it is necessary to first divide the fiber mesh in the section, and generate concrete fiber by selecting the number of fibers in two orthogonal directions of the section; for the component section, the data extraction algorithm developed by Revit secondary development technology is used to automatically extract the material information of the steel bar, the relative position information of the steel bar, the material information of the concrete, and the number of concrete meshes;

[0017] Unit settings: The setting interface of beam and column units is separated to set different unit types for beams and columns;

[0018] Vertical load setting: The vertical dead loads of components such as beams, columns and walls are automatically calculated by reading the geometric data and material type information of the Revit model; the dead and live loads of floors and roofs are automatically transferred to the supporting beams in the form of line loads according to the load transfer rules of two-way slabs and one-way slabs, and the vertical loads of each floor are determined according to the combination rules of the representative values ​​of gravity loads; the representative values ​​of gravity loads of each floor are automatically calculated and redistributed to each node according to the load-bearing area allocated to the node;

[0019] Output and analysis settings: including static output settings and interface output settings. Static output settings are to preset node displacement and unit shear force output command streams in the background code to provide basic data for subsequent analysis. Interface output settings are for users to select the target record type, set the output format, select the save path, determine whether to open the time series, select the node, and select the output data type in the dialog box of the front-end interface. The program will automatically generate the corresponding TCL command stream. According to the analysis rules and parameter setting requirements of the OpenSEES program, modal analysis, gravity analysis, dynamic time history analysis and Pushover analysis are performed respectively.

[0020] Analysis call: While outputting the TCL command stream, directly call OpenSEES on the Revit platform to perform finite element analysis.

[0021] Furthermore, in the section setting, the data extraction algorithm developed by Revit secondary development technology is used to automatically extract the required information. The specific process is as follows:

[0022] (1) Use the element collector FilteredElementCollector, the category filter ElementCategoryFilter and the function WherePasses to filter the longitudinal reinforcement;

[0023] (2) Screening of longitudinal reinforcements belonging to the same component;

[0024] (3) Record the coordinates of each longitudinal steel bar relative to the cross section. The Y and Z coordinates of the steel bar are recorded in the X direction of the component cross section, and the X and Z coordinates of the steel bar are recorded in the Y direction.

[0025] (4) Obtain the diameter of the steel bar and use the circle area calculation formula to calculate the cross-sectional area of ​​each longitudinal steel bar;

[0026] (5) Write a concrete fiber generation function CreateConcreteFiber(), pass the cross-sectional dimensions, the number of fibers in the X direction of the cross section, and the number of fibers in the Y direction of the cross section as parameters into the function, and calculate the concrete fiber coordinates corresponding to the cross section;

[0027] (6) The unique attributes such as Familysymbol ID, steel bar material ID, steel bar size, and concrete material ID are combined into a matching identification code for the section to prevent repeated acquisition of the same section of different components.

[0028] Further, step S2 includes displaying and outputting the vertex displacement time history curve, and the implementation steps are as follows:

[0029] (1) Use the Windows File class to read the OpenSEES analysis node displacement file and filter the time history data of the top target node;

[0030] (2) Store the data in the string[] string array, use the string splitting function to cut out the data and convert it into double type;

[0031] (3) Store time data in List <double>In the times generic collection, the displacement data is stored in List <double>disp generic collection;

[0032] (4) Use the Chart control in the winform control library to draw a curve, and use times and disp as coordinate points to fit the curve.

[0033] Further, step S2 includes displaying and outputting the base shear force time history curve, and the implementation steps are as follows:

[0034] (1) Use the Windows File class to read the OpenSEES analysis unit shear data file and filter the time history data of the shear force of the bottom unit. The shear force data of each unit is 13 columns. The first column is time, the second to seventh columns are the x, y, z direction shear force and Mx, My, My bending moments at one end of the unit, and the eighth to thirteenth columns are the x, y, z direction shear force and Mx, My, Mz bending moments at the other end.

[0035] (2) Store the data in the string[] string array, use the string splitting function to cut out the data and convert it into double type;

[0036] (3) Store time data in List <double>In the times generic collection, the displacement data is stored in List <double>force in generic collection;

[0037] (4) Use the Chart control in the Winform control library to draw a curve, and use times and force as coordinate points to fit the curve.

[0038] Further, step S2 includes displaying and outputting Pushover analysis data, and the implementation steps are as follows:

[0039] (1) Use the Windows File class to read the OpenSEES pushover analysis vertex displacement file, where the first column is the load multiple and the second column is the displacement data, and filter out the above two columns of data;

[0040] (2) Store the data in the string[] string array, use the string splitting function to cut out the data and convert it into double type;

[0041] (3) Store the displacement data in List <double>In the disp generic collection, the load multiple data is stored in List <double>loads generic collection;

[0042] (4) Use the Chart control in the winform control library to draw a curve, and use disp and loads as coordinate points to draw the curve.

[0043] Further, in step S3, the steps of implementing earthquake damage level analysis are as follows:

[0044] (1) Data extraction: Use the File class to loop through the strain data and store it into the string[] string array;

[0045] (2) Data screening: Use the string segmentation function to segment the target data and convert it into double type;

[0046] (3) Take the maximum value: calculate the maximum value of the absolute value of the strain time history data;

[0047] (4) Determine the damage level: By using the damage level judgment function based on the judgment criteria in the table below, the damage interval where the maximum value is located can be determined to obtain the damage condition of the component;

[0048]

[0049] (5) Damage data display: Use Winform's DataGridView control to display damage data in the earthquake damage visualization display module.

[0050] Further, in step S4, the steps for storing through Revit project parameters are as follows:

[0051] (1) Set the shared parameter file address, parameter group name, parameter type and other conditions to create a shared parameter file;

[0052] (2) Read the shared parameter file, bind the shared parameters to the model components, and add project parameters to the specified family;

[0053] (3) Finally, assign values ​​to the project parameters and store the analysis data in the corresponding model components. The earthquake damage result data of any component can be viewed in the property bar interface in Revit.

[0054] Further, in step S5, the steps of displaying the earthquake damage status through the three-dimensional view are as follows:

[0055] (1) Filtering family instances: Use the element collector class, element category filter, and Wherepasses() function to filter the family instance FamilyInstance;

[0056] (2) Matching components: Match the damage analysis results with the component family instances in Revit;

[0057] (3) Find family instance attributes: filter the attributes that control component surface coloring;

[0058] (4) Modify color: Modify the component color according to the damage result to show the damage status of the component.

[0059] Further, in step S5, the color-annotated Revit model is loaded into Revit roaming, 3DMax roaming, or Fuzor roaming to perform virtual roaming inside the building structure model, so as to observe the damage level and distribution of each component in the building structure and view the earthquake damage information of the corresponding components.

[0060] The beneficial effects of the present invention are:

[0061] The seismic analysis method of reinforced concrete frame structure based on Revit-OpenSEES described in the present invention has the following beneficial effects:

[0062] The present invention uses Revit software as a platform, uses C# language, and develops a structural seismic elastoplastic analysis and evaluation program RTO that integrates Revit and OpenSEES through a Revit API interface.

[0063] After the operator opens the Revit building structure model in the Revit software, the RTO program plug-in is started. In the program, material mechanical properties, analysis unit selection, component section settings, structural constraint settings, vertical load settings, various analysis settings and output settings are set for the model, and a TCL command stream is automatically generated; OpenSEES is called in the RTO program to read the TCL command stream for seismic finite element analysis and output analysis data. The operator performs post-processing analysis in the program according to the analysis output settings, including the display and output of various curves, tables and document data. In addition, the visual display of the seismic damage state of the structure can be achieved by automatically adjusting the color of the Revit model, and the key information of seismic damage can be stored in the properties of the Revit component model as part of the BIM information for BIM collaborative analysis and the management of the entire life cycle of the building.

[0064] The present invention realizes the integration of common BIM software and finite element analysis program, realizes seismic analysis and damage evaluation of reinforced concrete frame structure directly based on BIM model, and provides a solution for structural seismic analysis to participate in the collaborative work of BIM. In addition, the present invention also provides pre- and post-processing functions for OpenSEES, which has certain application value in the fields of BIM application and building structure seismic resistance. BRIEF DESCRIPTION OF THE DRAWINGS

[0065] Figure 1 This is the overall design idea diagram of the RTO program;

[0066] Figure 2 This is a schematic diagram of the three-layer architecture of RTO;

[0067] Figure 3 The overall logic diagram extracted for BIM data;

[0068] Figure 4 This is a schematic diagram of the analytical nodes in Revit;

[0069] Figure 5 This is a schematic diagram of the entire floor slab and the separated partial floor slab;

[0070] Figure 6 Revit model diagram for reinforcement of beam components;

[0071] Figure 7 Revit model diagram for reinforcement of column section;

[0072] Figure 8 Flowchart of algorithm for automatic application of vertical loads;

[0073] Fig. 9 Store flow charts for Revit project parameters;

[0074] Fig.10 Reference map for Fuzor walkthrough;

[0075] Fig.11 This is the overall architecture diagram of the RTO program;

[0076] Fig.12 This is the architecture diagram of the pre-analysis processing module;

[0077] Fig.13 This is the architecture diagram of the post-processing module for analysis;

[0078] Fig.14 It is the front-end operation interface of the RTO system;

[0079] Fig.15 Provides a constraint setting interface for the RTO system;

[0080] Fig.16 Set up the interface for the RTO concrete material;

[0081] Fig.17 Create interfaces for concrete fiber mesh;

[0082] Fig.18 Select the interface for the unit type;

[0083] Fig.19 Set up interfaces for vertical loads;

[0084] Fig. 20 Set the record type interface for output;

[0085] Fig.21 Select the interface for a single node;

[0086] Fig. 22 Select the interface for the scope node;

[0087] Fig.23 It is the node data type interface;

[0088] Fig.24 It is the unit data type interface;

[0089] Fig.25 Set the preview window interface for output;

[0090] Fig.26 Set up the interface for gravity analysis;

[0091] Fig. 27 Set up the interface for Pushover analysis;

[0092] Fig.28 Set up the interface for dynamic time history analysis;

[0093] Fig.29 It is the top displacement time history curve interface;

[0094] Fig.30 It is the Pushover capacity curve interface in top displacement-base shear format;

[0095] Fig.31 This is a comparison diagram of top displacement time history curves;

[0096] Fig.32 Comparison diagram of the maximum inter-story displacement angle distribution curve

[0097] Fig.33 This is a comparison chart of Pushover capability curves;

[0098] Fig.34 Provides interface for damage analysis and output data tables;

[0099] Fig.35 A Revit model diagram that shows earthquake damage;

[0100] Fig.36 A walkthrough diagram of the BIM model with earthquake damage annotated;

[0101] Fig.37 An interface for seismic damage information stored in Revit beam component properties;

[0102] Fig.38 Interface for seismic damage information stored in Revit column component properties. DETAILED DESCRIPTION

[0103] The following will refer to the attached Figure 1-38 The seismic analysis method of reinforced concrete frame structure based on Revit-OpenSEES described in the present invention is described in more detail.

[0104] Example 1

[0105] 1. Programming ideas

[0106] The overall design idea of ​​the RTO program is as follows Figure 1 As shown, it specifically includes: (1) pre-analysis processing module, (2) post-analysis processing module, and (3) earthquake damage visualization display module.

[0107] 2 Program Architecture

[0108] The RTO program is developed based on a three-tier development framework and is appropriately adjusted according to the application scenarios and functions of the software. RTOForms is the presentation layer; RTOBLL is the business logic layer; RTODAL is the data access layer; RTODataBase is the data class that replaces the database; RTOCommon is the common command layer; and RTOModel is used to define custom classes. The relationship between the layers is as follows: Figure 2 shown.

[0109] 3. Technical points of pre-processing module

[0110] The RTO program pre-processing module includes the model data conversion and seismic analysis setting functions between Revit and OpenSEES. The conversion from Revit's BIM model to OpenSEES's finite element model can be understood as the BIM data extraction process. Its overall logical process is as follows: Figure 3 After completing the model conversion, in order to call the OpenSEES program for finite element analysis, it is also necessary to develop seismic analysis settings and output settings functions in the RTO program.

[0111] 3.1 Node Conversion

[0112] (1) Analysis node extraction

[0113] In OpenSEES modeling, it is necessary to establish nodes for each beam, column, wall and other components. The command flow for defining a 3D node in OpenSEES is: node number x coordinate y coordinate z coordinate. The Revit structural model has corresponding analysis nodes, such as Figure 4 As shown. In Revit software, the structural analysis model can be enabled, and the connection between components will generate analysis nodes (ReferencePoint). The analysis nodes and corresponding coordinates in the Revit structural model can be obtained through the element collector (FilteredElementCollector), category filter (ElementCategoryFilter) and WherePasses function.

[0114] (2) Numbering

[0115] The three-dimensional coordinates can be directly obtained from the spatial coordinates of the analysis node. Each node needs to be assigned a unique number, so a numbering class UniqueTag is designed, which provides a method GetTag() to provide the number. Each Node can be generated through a loop statement, and the Insert function of the DAL layer can be called in the BLL layer of the three-layer architecture to store the Node in the properties of the DataBase data class.

[0116] (3) Floor slab nodes

[0117] The floor in OpenSEES is defined by the center node of the floor. The floor in the Revit model provides the coordinate Origin of the center node. However, in actual research, it is found that not only the center coordinate of the entire floor is needed, but also the center coordinate of each local floor separated by the beam as the dividing line, such as Figure 5 shown.

[0118] 3.2 Constraint Settings

[0119] The TCL constraint command flow of OpenSEES is: fix node number xyz Mx My Mz. There are two kinds of constraints to consider, namely the bottom constraint and the floor constraint.

[0120] (1) Bottom constraint

[0121] Calculate the minimum Z coordinate MinZ among all nodes, and use MinZ as a condition to filter out nodes with the same Z coordinate as MinZ. It should be noted here that since the coordinate data is stored in double type, there is a precision problem, so "==" cannot be used directly to compare the same, and AB<0.00001 (the precision can be controlled by yourself) is needed to determine whether the coordinates are the same. The filtered nodes are the nodes that need to apply bottom constraints.

[0122] (2) Floor constraints

[0123] The floor constraints need to be applied in conjunction with the rigid plate assumption. The floor constraints of each layer are defined using the center points of the floor plates that have been obtained. In the rigid plate assumption, the main node is the rigid center of the floor plate, and the subsidiary nodes are other nodes on the same layer.

[0124] 3.3 Material settings

[0125] Revit software only provides simple information such as material type, but does not include the material mechanics information required for finite element analysis, so it cannot be directly converted into the material model data required for OpenSEES analysis. To solve this problem, a material setting function was developed in the front end of the RTO program, that is, the setting function was developed in Revit according to the material mechanics information required for OpenSEES analysis, such as the compressive yield strength, peak strain, residual strength after yield and other parameters required for concrete materials.

[0126] During material installation, the RTO program can read the material type in the Revit structural model, then match the set mechanical parameters to the corresponding material in the model according to the material type, and form a TCL command stream about the material properties.

[0127] 3.4 Section settings

[0128] Considering that the fiber model of OpenSEES can well simulate the mechanical properties of the plastic hinge area of ​​reinforced concrete frame structural components, a setting function for the fiber section is developed in the RTO program. For the steel fiber of the section, the program can automatically read the coordinates and material properties of the steel bar in the Revit model and convert them into the corresponding TCL command stream; for concrete fiber, it is necessary to first divide the fiber mesh in the section, so a concrete fiber setting function is developed in the RTO program to generate concrete fibers by selecting the number of fibers in two orthogonal directions of the section.

[0129] For the cross section of the component, it is necessary to obtain the material information of the steel bars, the relative position information of the steel bars, the material information of the concrete, and the number of concrete grids. The data extraction algorithm developed using Revit secondary development technology can automatically extract the required information. The specific process is as follows:

[0130] (1) Use the element collector FilteredElementCollector, the category filter ElementCategoryFilter and the function WherePasses to filter the longitudinal reinforcement;

[0131] (2) Screening of longitudinal reinforcements belonging to the same component;

[0132] (3) Record the coordinates of each longitudinal steel bar relative to the cross section. The Y and Z coordinates of the steel bar are recorded in the X direction of the component cross section, and the X and Z coordinates of the steel bar are recorded in the Y direction.

[0133] (4) Obtain the diameter of the steel bar and calculate the cross-sectional area of ​​each longitudinal steel bar using the circle area calculation formula;

[0134] (5) Write a concrete fiber generation function CreateConcreteFiber(), pass the cross-sectional dimensions, the number of fibers in the X direction of the cross section, and the number of fibers in the Y direction of the cross section as parameters into the function, and calculate the concrete fiber coordinates corresponding to the cross section;

[0135] The unique attributes such as Familysymbol ID, steel bar material ID, steel bar size, concrete material ID, etc. are combined into a matching identification code for the section to prevent repeated acquisition of the same section for different components.

[0136] 3.5 Unit Settings

[0137] According to the commonly used analysis units in OpenSEES, the RTO program provides four unit setting functions. Taking the dispBeamColumn unit as an example, its TCL command flow is: element dispBeamColumn$eleTag$iNode$jNode$numIntgrPts$secTag$transfTag.

[0138] The meanings of the parameters in the above command flow are shown in Table 1:

[0139] Table 1 Parameters defined by the dispBeamColumn unit

[0140]

[0141] Among them, the number, the two end node number, the section number and the coordinate transformation number can be automatically associated through the program, and the user only needs to set the unit integration points. In addition, in order to meet different setting requirements, the setting interface of the beam and column units is separated, and different unit types can be set for beams and columns.

[0142] 3.6 Vertical load setting

[0143] For the vertical dead loads of components such as beams, columns and walls, the RTO program can automatically calculate them by reading the geometric data and material type information of the Revit model. For the dead and live loads of floors and roofs, they must be automatically transferred to the supporting beams in the form of line loads according to the load transfer rules of two-way slabs and one-way slabs, respectively, and the vertical loads of each floor must be determined according to the combination rules of the representative values ​​of gravity loads. However, in the OpenSEES program, the structural vertical loads are dispersed at each node of the model. Therefore, the RTO program will further automatically calculate the representative value of the gravity load for each floor, and redistribute it to each node according to the load-bearing area allocated to the node. The vertical load generation algorithm based on Revit model information is as follows: Figure 8 shown.

[0144] The principle of point load calculation is similar to that of mass source calculation. On the basis of line load, the line load on the beam is multiplied by the length to obtain the concentrated load, and the concentrated load is evenly distributed to the nodes at both ends. The same applies to the column.

[0145] 3.7 Output and Analysis Settings

[0146] The output settings of the RTO program include static output settings and interface output settings.

[0147] (1) Static output setting is to preset the output command streams such as node displacement and unit shear force in the background code to provide basic data for subsequent analysis. This part is a fixed output and cannot be viewed on the front end. The setting process is completed in the background, which will increase the calculation time a little, but will not increase the calculation difficulty.

[0148] (2) Interface output settings refer to the user selecting the target record type, setting the output format, selecting the save path, determining whether to enable the time series, selecting the node, selecting the output data type, etc. in the dialog box of the front-end interface. The program will automatically generate the corresponding command flow.

[0149] According to the analysis rules and parameter setting requirements of the OpenSEES program, the analysis setting interfaces for modal analysis, gravity analysis, dynamic time history analysis and Pushover analysis were developed respectively.

[0150] 3.8 Analysis call

[0151] The OpenSEES calling function has been developed in RTO. While RTO outputs the TCL command stream, OpenSEES can be directly called on the Revit platform for finite element analysis.

[0152] 4 Technical points of post-processing module

[0153] In the RTO system, the post-processing module can complete the following operations:

[0154] (1) Group, organize and extract the analysis data output by OpenSEES, and display and output them in the form of graphs, tables or documents.

[0155] (2) A secondary seismic damage analysis of the structure is performed based on the concrete and steel strain data output by OpenSEES to determine the seismic damage levels of each structural component that cannot be directly obtained by OpenSEES.

[0156] (3) Automatically store key structural earthquake damage information in the Revit structural model as BIM data for collaborative analysis and structural analysis and information management throughout the building life cycle.

[0157] 4.1 Chart and document output

[0158] Through C# programming, the analysis data of OpenSEES is automatically extracted and output in the form of graphics, tables or documents after sorting and processing. Taking three functions as examples, the ideas and implementation of the development of the chart output function of the post-processing module are explained.

[0159] (1) Vertex displacement time history curve

[0160] 1) Use the Windows File class to read the OpenSEES analysis node displacement file and filter the time history data of the top target node.

[0161] 2) Store the data in the string[] string array, use the string splitting function to cut out the data and convert it into double type.

[0162] 3) Store time data in List <double>In the times generic collection, the displacement data is stored in List <double>disp generic collection.

[0163] 4) Use the Chart control in the winform control library to draw a curve, and use times and disp as coordinate points to fit the curve.

[0164] (2) Base shear force time history curve

[0165] 1) Use the Windows File class to read the OpenSEES analysis unit shear data file and filter the time-history data of the shear force of the bottom unit. It should be noted here that the shear force data of a unit is 13 columns. The first column is time, columns 2-7 are the x, y, z direction shear force and Mx, My, My bending moments at one end of the unit, and columns 8-13 are the x, y, z direction shear force and Mx, My, Mz bending moments at the other end.

[0166] 2) Store the data in the string[] string array, use the string splitting function to cut out the data and convert it into double type.

[0167] 3) Store time data in List <double>In the times generic collection, the displacement data is stored in List <double>force generic collection.

[0168] 4) Use the Chart control in the Winform control library to draw a curve, and use times and force as coordinate points to fit the curve.

[0169] (3) Pushover data

[0170] 1) Use the Windows File class to read the OpenSEES pushover analysis vertex displacement file, where the first column is the load multiple and the second column is the displacement data. Filter out these two columns of data.

[0171] 2) Store the data in the string[] string array, use the string splitting function to cut out the data and convert it into double type.

[0172] 3) Store the displacement data in List <double>In the disp generic collection, the load multiple data is stored in List <double>loads generic collection.

[0173] 4) Use the Chart control in the winform control library to draw the curve, and use disp and loads as coordinate points to draw the curve.

[0174] 4.2 Earthquake damage analysis

[0175] According to the "Building Seismic Resilience Evaluation Standard", the seismic damage state of reinforced concrete frame structural components can be judged by the material strain of the positive section, and the material strain judgment standards corresponding to each damage level are shown in Table 2. OpenSEES can output the strain time history data of the fiber of the component under the action of earthquake. Substituting the maximum strain values ​​of the longitudinal steel bars and concrete into Table 2, the seismic damage level of each structural component can be obtained.

[0176] Based on the above principles, an earthquake damage analysis algorithm is written to extract the strain data of the steel fiber and concrete fiber of all components, and the damage level is determined according to the strain value. The algorithm implementation process is as follows:

[0177] 1) Data extraction. Use the File class to loop through the strain data and store it into the string[] string array;

[0178] 2) Data screening. Use the string segmentation function to segment the target data and convert it into double type;

[0179] 3) Take the maximum value. Calculate the maximum absolute value of the strain time history data.

[0180] 4) Determine the damage level. By using the damage level judgment function to determine the damage interval where the maximum value is located, the damage condition of the component can be obtained.

[0181] 5) Damage data display: Use Winform's DataGridView control to display damage data in the RTO post-processing interface.

[0182] Table 2 Strain-based seismic damage criteria for reinforced concrete frame structural components

[0183]

[0184] Note: ε1 is the principal tensile strain of steel bars; ε2 is the principal tensile strain of concrete; ε p and ε cu are the peak strain and ultimate strain of concrete under uniaxial compression; ε y is the yield strain of the steel bar.

[0185] 4.3 Damage Information Storage

[0186] After obtaining the earthquake damage level of each structural component, this information can be automatically stored in the Revit component model and used as BIM data for structural collaborative analysis and integrated management of the entire life cycle. The RTO program uses the Revit project parameter storage method, which can use Revit secondary development technology to achieve rapid batch storage. Users can directly view, add and modify data in the Revit property bar interface.

[0187] When storing through Revit project parameters, first set the shared parameter file address, parameter group name, parameter type and other conditions to create a shared parameter file. Then, read the shared parameter file, bind the shared parameters to the model components, and add project parameters to the specified family. Finally, assign values ​​to the project parameters and store the analysis data in the corresponding model components. You can view the seismic evaluation result data of any component in the property bar interface in Revit. The implementation process of Revit project parameter storage is as follows: Fig. 9 shown.

[0188] 5 Technical points for damage visualization

[0189] In the BIM-based display of structural seismic damage status, in addition to using the powerful three-dimensional viewing function of BIM software, the BIM model can also be imported into VR software for virtual roaming to give full play to its visualization advantages.

[0190] 5.1 3D View Display

[0191] The corresponding relationship between earthquake damage level and Revit model color is shown in Table 3. Revit secondary development technology is used to automatically change the color of structural components according to their damage level, and the damage status and distribution characteristics of the structure under earthquake action are displayed in a three-dimensional view. The specific algorithm is as follows:

[0192] (1) Filter family instances. Use the element collector class, element category filter, and Wherepasses() function to filter family instances FamilyInstance;

[0193] (2) Matching components. Match the damage analysis results with the component family instances in Revit;

[0194] (3) Find the family instance attributes. Filter the attributes that control the surface coloring of the component;

[0195] Modify color. Modify the component color according to the damage result to show the damage status of the component.

[0196] Table 3 Revit model marking colors corresponding to earthquake damage levels

[0197] Damage status level Level 0 Level 1 Level 2 Level 3 Level 4 Model annotation color Shell color yellow orange color red Dark Red

[0198] 5.2 Virtual Tour Display

[0199] The color-coded Revit model can be loaded into virtual tour software, such as Revit Tour, 3DMax Tour, and Fuzor Tour. Users can conduct virtual tours inside the building structure model, so that they can more clearly observe the damage level and distribution of each component in the building structure, and can view the earthquake damage information of the corresponding components. Fig.10 shown.

[0200] Example 2

[0201] This embodiment describes the specific interface of the RTO system running in the Revit software.

[0202] Fig.14 This is the interface of the RTO system running in Revit software. The left side is the function switching area, the upper right side is the parameter setting area, and the lower right side is the information prompt area, which provides information prompts for user operations.

[0203] 1 System pre-processing operations

[0204] 1.1 Model node conversion

[0205] Start the RTO system in the program that opens the Revit structural model. The RTO system will automatically identify the nodes in the Revit structural model in the background and convert them into node parameters that can be read by OpenSEES.

[0206] 1.2 Constraint Settings

[0207] In the function switching area of ​​the front-end interface, click Constraint Settings. In the parameter setting area, you can add foundation bottom constraints and select floor assumptions, such as Fig.15 As shown, the floor is assumed to correspond to the rigid floor unit model and membrane unit model in OpenSEES, and the two models are connected to the load-bearing system using different constraints.

[0208] 1.3 Material Settings

[0209] The RTO system can automatically read the material type in the Revit model, and set the mechanical properties of concrete and steel materials in the front-end interface, such as Fig.16 shown.

[0210] 1.4 Section settings

[0211] The RTO system uses the fiber model of OpenSEES to simulate the mechanical properties of the cross section of reinforced concrete components, so the cross section here refers to the fiber cross section. The RTO system can automatically read the coordinates and material information of the steel bars in the Revit model and convert them into a cross-section steel fiber command stream that can be read by OpenSEES. For the cross-section concrete fiber setting, it is necessary to perform fiber meshing in the cross-section setting of the front-end interface and convert it into a cross-section concrete fiber command stream based on the concrete material information. Fig.17 Interface for cross-section concrete fiber meshing.

[0212] 1.5 Unit Settings

[0213] Four commonly used OpenSEES units can be selected in the RTO program to simulate the mechanical properties of beams and columns.

[0214] The unit setting function allows you to select the analysis unit type and assign it to the selected structural component. Fig.18 This is the setting interface for the column unit type.

[0215] 1.6 Vertical load setting

[0216] Enter the representative value of the gravity load in the load condition parameter setting interface, and the program will call the load distribution algorithm to apply the gravity load to the beam or column through a line load or a point load and pass it downward. Fig.19 Set up the interface for vertical loads.

[0217] 1.7 Output Settings

[0218] The output interface can be used to set the type, method, range and storage of output information required for various analyses, such as Fig. 20 As shown, it includes the selection of record type, output format, save address, range, data type, etc. and the automatically generated TCL command stream window.

[0219] (1) Record type selection

[0220] Record types include nodes, elements, sections, and fibers. You can choose to output a certain type of data result according to your needs. For example, if you select the node type, the node reaction data will be output. Since different record types require different parameters to be set, the interface will automatically switch when switching record types.

[0221] (2) Output format

[0222] Two output format settings are provided. The file type is the default type, which records the output as a .out type file, and the xml type records the output as a txt type file.

[0223] (3) Range type

[0224] In order to reduce the amount of calculation, you can set the output range in the output settings, for example, select several key nodes of the structure through a single node or a range node, and output the reaction data of these nodes, such as Fig.21 and Fig. 22 shown.

[0225] (4) Data type

[0226] Each record type corresponds to a different data type. For example, if the record type is node, you can choose to output data types such as absolute displacement, velocity, acceleration, etc. When switching the record type, the data type will automatically match the record type. Fig.23 and Fig.24 The data types correspond to the node type and unit type respectively.

[0227] (5) Output setting command stream preview

[0228] The output command stream preview window is provided on the right side of the output setting interface. The defined output settings are automatically displayed in the list, so that errors can be found and corrected in time. Fig.25 shown.

[0229] 1.8 Structural Analysis Settings

[0230] (1) Gravity analysis settings

[0231] In the gravity analysis interface, you need to set the integrator, matrix bandwidth processing, convergence criterion, iteration algorithm, maximum iteration step and analysis type. Fig.26 shown.

[0232] (2) Pushover analysis settings

[0233] In the Pushover analysis settings, in addition to the analysis settings such as integrator, matrix bandwidth processing, and convergence criteria, the program also provides the function of setting the pushover force direction and loading mode. The loading modes include triangular distribution, uniform distribution, parabolic distribution, and vibration mode distribution. For example, Fig. 27 shown.

[0234] (3) Dynamic time history analysis settings

[0235] The setting interface of dynamic time history analysis is as follows: Fig.28 It is necessary to link the seismic wave recording file, set the amplitude modulation coefficient of the seismic wave, the time interval of the seismic wave data, and set the structural damping ratio and other dynamic parameter settings.

[0236] 2 System post-processing operations

[0237] 2.1 Chart and data output

[0238] (1) Chart output

[0239] The RTO system post-processing can output the corresponding time history analysis curves, including: displacement time history of each node, shear time history of each floor, hysteresis energy dissipation time history of each floor, etc. Pushover capacity curves in top displacement-base shear format and top displacement-load multiple format can be output respectively, where the base shear data is calculated by multiplying the sum of the lateral pushover force by the load multiple. In addition, the RTO system post-processing can also output modal analysis data tables. Fig.29 and Fig.30 The top displacement time history curve and Pushover capacity curve output by the program are given respectively.

[0240] (2) Document data output

[0241] The post-processing of the RTO system can also output analysis data in the form of documents, making it easier for users to conduct in-depth analysis of the data. Fig.31 , Fig.32 and Fig.33 An analysis chart drawn using other software based on the data output by the system.

[0242] 2.2 Damage secondary analysis, output and display

[0243] OpenSEES can output the fiber strain data of components under earthquake action, and the seismic damage level of each component can be evaluated through secondary damage analysis. Fig.34 It is the display interface for damage analysis and damage identification data.

[0244] In order to clearly display the structural earthquake damage status in the Revit model, each model component in Revit is marked according to the color corresponding to each damage level in Table 3. Fig.35 and Fig.36 The earthquake damage conditions are displayed through Revit 3D view and virtual tour respectively.

[0245] 2.3 Earthquake damage data annotation

[0246] Through elastic-plastic analysis, the strain data of the longitudinal reinforcement and concrete in the structural components are obtained. The earthquake damage level of each component is determined according to Table 2, and the damage information can be automatically added to the property column of each component in the Revit model, such as Fig.37 and Fig.38 shown.

[0247] The above description is only a preferred specific implementation manner of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the technical field can make equivalent replacements or changes according to the technical solutions and concepts of the present invention within the technical scope disclosed by the present invention, which should be covered by the protection scope of the present invention.< / double> < / double> < / double> < / double> < / double> < / double> < / double> < / double> < / double> < / double> < / double> < / double>

Claims

1. A seismic analysis method for reinforced concrete frame structure based on Revit-OpenSEES, characterized in that: Here are the steps: S1. In the pre-analysis processing module, the finite element model is converted from the Revit BIM model to OpenSEES. After the model conversion is completed, the seismic analysis settings and output settings are performed; The specific steps of step S1 are as follows: Node conversion: extract the analysis nodes and corresponding coordinates in the Revit structural model, provide numbers through the numbering class, obtain the center coordinates of the entire floor and the center coordinates of each local floor separated by the beam as the dividing line; Constraint setting: Filter the nodes that need to be subject to bottom constraints, and use the center points of each floor slab that have been obtained to define the floor slab constraints of each floor; Material setting: Develop and set functions in Revit according to the material mechanical information required for OpenSEES analysis. When setting materials, read the material type in the Revit structural model, and then match the set mechanical parameters to the corresponding materials in the model according to the material type; Section setting: For the steel fiber of the section, the coordinates and material properties of the steel bar are read in the Revit structural model and converted into the corresponding TCL command stream; for concrete fiber, it is necessary to first divide the fiber mesh in the section, and generate concrete fiber by selecting the number of fibers in two orthogonal directions of the section; for the component section, the data extraction algorithm developed by Revit secondary development technology is used to automatically extract the material information of the steel bar, the relative position information of the steel bar, the material information of the concrete, and the number of concrete meshes; Unit settings: The setting interface of beam and column units is separated to set different unit types for beams and columns; Vertical load setting: The vertical dead loads of beams, columns and wall components are automatically calculated by reading the geometric data and material type information of the Revit structural model; the dead and live loads of floors and roofs are automatically transferred to the supporting beams in the form of line loads according to the load transfer rules of two-way slabs and one-way slabs, and the vertical loads of each floor are determined according to the combination rules of the representative values ​​of gravity loads; The representative value of gravity load of each floor will be automatically calculated and distributed to each node according to the load-bearing area allocated to the node; Output and analysis settings: including static output settings and interface output settings. Static output settings are to preset node displacement and unit shear force output command streams in the background code to provide basic data for subsequent analysis. Interface output settings are for users to select the target record type, set the output format, select the save path, determine whether to open the time series, select the node, and select the output data type in the dialog box of the front-end interface. The program will automatically generate the corresponding TCL command stream. According to the analysis rules and parameter setting requirements of the OpenSEES program, modal analysis, gravity analysis, dynamic time history analysis and Pushover analysis are performed respectively. Analysis call: While outputting the TCL command stream, directly call OpenSEES on the Revit platform to perform finite element analysis; S2. In the post-analysis processing module, the analysis data output by OpenSEES is grouped, sorted and extracted, and displayed and output in the form of graphs, tables or documents; S3. In the post-analysis processing module, a secondary seismic damage analysis is performed on the structure based on the concrete and steel strain data output by OpenSEES to determine the seismic damage level of each structural component that cannot be directly obtained by OpenSEES. In step S3, the steps for implementing earthquake damage level analysis are as follows: (1) Data extraction: Use the File class to loop through the strain data and store it into the string[] string array; (2) Data screening: Use the string segmentation function to segment the target data and convert it into double type; (3) Take the maximum value: calculate the maximum value of the absolute value of the strain time history data; (4) Determine the damage level: By using the damage level judgment function based on the judgment criteria in the table below, the damage interval where the maximum value is located can be determined to obtain the damage condition of the component; Where ε1 is the principal tensile strain of the steel bar; ε2 is the principal tensile strain of the concrete; p and ε cu are the peak strain and ultimate strain of concrete under uniaxial compression; ε y is the yield strain of the steel bar; (5) Damage data display: Use Winform's DataGridView control to display damage data in the earthquake damage visualization display module; S4. In the post-analysis processing module, key structural earthquake damage information is automatically stored in the Revit structural model as BIM data for collaborative analysis and structural analysis and information management throughout the building life cycle; S5. In the earthquake damage visualization display module, the earthquake damage status is displayed through a three-dimensional view, or virtual roaming is performed through VR software.

2. The seismic analysis method for reinforced concrete frame structure based on Revit-OpenSEES as claimed in claim 1, characterized in that: In the section setting, the data extraction algorithm developed by Revit secondary development technology is used to automatically extract the required information. The specific process is as follows: (1) Use the element collector FilteredElementCollector, the category filter ElementCategoryFilter and the function WherePasses to filter the longitudinal reinforcement; (2) Screening of longitudinal reinforcements belonging to the same component; (3) Record the coordinates of each longitudinal steel bar relative to the cross section. The Y and Z coordinates of the steel bar are recorded in the X direction of the component cross section, and the X and Z coordinates of the steel bar are recorded in the Y direction. (4) Obtain the diameter of the steel bar and use the circle area calculation formula to calculate the cross-sectional area of ​​each longitudinal steel bar; (5) Write a concrete fiber generation function CreateConcreteFiber(), pass the cross-sectional dimensions, the number of fibers in the X direction of the cross section, and the number of fibers in the Y direction of the cross section as parameters into the function, and calculate the concrete fiber coordinates corresponding to the cross section; (6) The Familysymbol ID, steel bar material ID, steel bar size, and concrete material ID are combined into a matching identification code for the section to prevent repeated acquisition of the same section of different components.

3. The seismic analysis method for reinforced concrete frame structure based on Revit-OpenSEES as claimed in claim 1, characterized in that: Step S2 includes displaying and outputting the vertex displacement time history curve, and the implementation steps are as follows: (1) Use the Windows File class to read the OpenSEES analysis node displacement file and filter the time history data of the top target node; (2) Store the data in the string[] string array, use the string splitting function to cut out the data and convert it into double type; (3) Store time data in List <double>In the times generic collection, the displacement data is stored in List <double> disp generic collection;< / double> < / double> (4) Use the Chart control in the winform control library to draw a curve, and use times and disp as coordinate points to fit the curve.

4. The seismic analysis method for reinforced concrete frame structure based on Revit-OpenSEES as claimed in claim 1, characterized in that: Step S2 includes displaying and outputting the base shear time history curve, and the implementation steps are as follows: (1) Use the Windows File class to read the OpenSEES analysis unit shear data file and filter the time history data of the shear force of the bottom unit. The shear force data of each unit is 13 columns. The first column is time, the second to seventh columns are the x, y, z direction shear force and Mx, My, My bending moments at one end of the unit, and the eighth to thirteenth columns are the x, y, z direction shear force and Mx, My, Mz bending moments at the other end. (2) Store the data in the string[] string array, use the string splitting function to cut out the data and convert it into double type; (3) Store time data in List <double>In the times generic collection, the displacement data is stored in List <double> force in generic collection;< / double> < / double> (4) Use the Chart control in the Winform control library to draw a curve, and use times and force as coordinate points to draw the curve.

5. The seismic analysis method for reinforced concrete frame structure based on Revit-OpenSEES as claimed in claim 1, characterized in that: Step S2 includes the display and output of Pushover analysis data, and the implementation steps are as follows: (1) Use the Windows File class to read the OpenSEES pushover analysis vertex displacement file, where the first column is the load multiple and the second column is the displacement data, and filter out the above two columns of data; (2) Store the data in the string[] string array, use the string splitting function to cut out the data and convert it into double type; (3) Store the displacement data in List <double>In the disp generic collection, the load multiple data is stored in List <double> loads generic collection;< / double> < / double> (4) Use the Chart control in the winform control library to draw a curve, and use disp and loads as coordinate points to fit the curve.

6. The seismic analysis method for reinforced concrete frame structure based on Revit-OpenSEES as claimed in claim 1, characterized in that: In step S4, the steps for storing through Revit project parameters are as follows: (1) Set the shared parameter file address, parameter group name, and parameter type conditions to create a shared parameter file; (2) Read the shared parameter file, bind the shared parameters to the model components, and add project parameters to the specified family; (3) Finally, assign values ​​to the project parameters, store the analysis data in the corresponding model components, and view the earthquake damage result data of any component in the property bar interface in Revit.

7. The seismic analysis method for reinforced concrete frame structure based on Revit-OpenSEES as claimed in claim 1, characterized in that: In step S5, the steps for displaying the earthquake damage status through the three-dimensional view are as follows: (1) Filtering family instances: Use the element collector class, element category filter, and Wherepasses() function to filter the family instance FamilyInstance; (2) Matching components: Match the damage analysis results with the component family instances in Revit; (3) Find family instance attributes: filter the attributes that control component surface coloring; (4) Modify color: Modify the component color according to the damage result to show the damage status of the component.

8. The seismic analysis method for reinforced concrete frame structure based on Revit-OpenSEES as claimed in claim 1, characterized in that: In step S5, the color-annotated Revit model is loaded into Revit roaming, 3DMax roaming, or Fuzor roaming to perform virtual roaming inside the building structure model, so as to observe the damage level and distribution of each component in the building structure and view the earthquake damage information of the corresponding components.

Images