CN116090078B - BIM model anti-seismic optimization method of road bridge structure - Google Patents
BIM model anti-seismic optimization method of road bridge structure Download PDFInfo
- Publication number
- CN116090078B CN116090078B CN202310361135.XA CN202310361135A CN116090078B CN 116090078 B CN116090078 B CN 116090078B CN 202310361135 A CN202310361135 A CN 202310361135A CN 116090078 B CN116090078 B CN 116090078B
- Authority
- CN
- China
- Prior art keywords
- unit
- bridge
- road
- seq
- domain
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Architecture (AREA)
- Bridges Or Land Bridges (AREA)
Abstract
The invention relates to the field of bridge earthquake resistance, and provides a BIM model earthquake resistance optimization method of a road and bridge structure. The method can effectively improve the shock resistance of the road and bridge structure, optimize the structural design scheme, improve the matching degree with the actual engineering requirement, and adopt the BIM technology to perform shock resistance optimization on the road and bridge structure, thereby realizing digital management and whole process visualization in the design and operation processes, improving the engineering design efficiency, reducing design errors and loopholes and providing reliable technical support for road and bridge engineering.
Description
Technical Field
The invention relates to the field of bridge earthquake resistance, in particular to a BIM model earthquake resistance optimization method of a road bridge structure.
Background
With the continuous development of computer technology and simulation technology, the combination of BIM technology and finite element analysis method for anti-seismic design has become an advanced and reliable design method. BIM technology is a design method based on a digital model, and by establishing a digital model of a structure, the collaborative design of each specialty such as building, structure and equipment is realized, and design conflicts and errors are reduced. In the BIM model, the information such as the geometric form, the material, the structure and the like of the structure are accurately expressed, and can be used for dynamic analysis and optimal design. The finite element analysis is a mechanical analysis method based on numerical calculation, and the mechanical response of the structure is obtained by dividing the structure into a plurality of finite small units and establishing equation solution in each unit. The finite element analysis can consider the complex form and the material mechanical property of the structure, and the dynamic response and the anti-seismic performance of the structure are evaluated more accurately and reliably.
Along with the continuous acceleration of the urban process, the road and bridge construction scale and the road and bridge construction quantity are continuously increased. However, the frequent occurrence of natural disasters such as earthquakes brings serious challenges to the safety and stability of road and bridge structures. Therefore, how to improve the anti-seismic performance of the road bridge structure becomes a problem to be solved in road bridge engineering.
Disclosure of Invention
The invention aims to provide a BIM model anti-seismic optimization method of a road bridge structure, which aims to solve one or more technical problems in the prior art and at least provides a beneficial selection or creation condition.
The invention relates to the field of bridge earthquake resistance, and provides a BIM model earthquake resistance optimization method of a road and bridge structure. The method can effectively improve the shock resistance of the road and bridge structure, optimize the structural design scheme, improve the matching degree with the actual engineering requirement, and adopt the BIM technology to perform shock resistance optimization on the road and bridge structure, thereby realizing digital management and whole process visualization in the design and operation processes, improving the engineering design efficiency, reducing design errors and loopholes and providing reliable technical support for road and bridge engineering.
To achieve the above object, according to an aspect of the present invention, there is provided a BIM model seismic optimization method of a road bridge structure, the method including the steps of:
s100, building a road bridge BIM model in BIM software according to the structural information of the road bridge;
s200, converting the road bridge BIM model into a road bridge finite element model through model conversion software;
s300, performing unit anti-seismic steady-state analysis on the road bridge finite element model to obtain a unit anti-seismic steady-state analysis result;
and S400, optimizing the road bridge BIM based on the unit anti-seismic steady-state analysis result.
Further, in step S100, the BIM software is any one of Revit, tekla, BIM5D, bentley, and the road and bridge is any one of a roadbed, a road surface, a bridge, a culvert, and a tunnel.
Further, in step S100, the structural information of the bridge at least includes span of the bridge, load requirement of the bridge, material specification of the bridge, and basic geometry and configuration of the bridge.
Further, in step S100, the method for building the bridge BIM model in the BIM software according to the structural information of the bridge specifically includes: in BIM software, an initial model of a road and a bridge is created according to basic geometric shapes and structures in the structural information of the road and the bridge, load requirements of the road and the bridge in the structural information of the road and the bridge and material specifications of the road and the bridge are added into the initial model of the road and the bridge, road and bridge surface materials and longitudinal and transverse reinforcing steel bars are added on the surface of the initial model of the road and the bridge, and the initial model of the road and the bridge is saved to be the BIM model of the road and the bridge; wherein the road and bridge surface material is reinforced concrete or masonry material.
Further, in step S100, a bridge BIM model is built in the BIM software according to the structural information of the bridge, and the method further includes:
s101, in BIM software, performing collision detection on a road bridge BIM model through a collision detection function, wherein the collision detection is used for detecting model components with collision in the road bridge BIM model, and modifying the collision components in the road bridge BIM model after the collision detection;
s102, circulating S101 until BIM software sends a prompt of no collision detection after collision detection is carried out on the road bridge BIM model.
Further, in step S200, the model conversion software is femtranfer or Hypermesh.
Further, in step S300, the method for performing unit anti-seismic steady-state analysis on the road bridge finite element model to obtain the unit anti-seismic steady-state analysis result specifically includes:
s301, finite element analysis is carried out on a road bridge finite element model, stress values received by each unit in the road bridge finite element model are solved (the road bridge finite element model converted through FEMTransfer comprises a plurality of discrete units which are gathered together to be capable of representing an actual continuous domain, so that the fact that the finite number of unknowns are used for approaching the infinite unknowns is achieved), the number of all units in the road bridge finite element model is recorded as N, and the N units are used as N structural units;
s302, representing the ith unit in the N structural units by sn (i), representing the stress value received by the ith unit in the N structural units by sts (i), wherein i is a sequence number, i=1, 2, … and N, creating a blank array sts, sequentially storing N values sts (1), sts (2), … and sts (N) into the array sts, recording sts (j) as the jth element in the array sts, j=1, 2, … and N, recording the element with the smallest element value in the array sts as sts (a), recording the a-th unit (sn (a)) in the N structural units as an inner unit, creating a blank sequence Seq, adding the sequence number a into the sequence Seq, and turning to S303;
s303, in the N structure units, forming a first unit domain by the inner units and all units connected with the inner units, using U1 (k 1) to represent the stress magnitude of the kth 1 unit in the first unit domain, wherein k1 is a sequence number, k1=1, 2, …, U1, U1 is the number of all units in the first unit domain, recording sv (k 1) =ABS (U1 (k 1) -sts (a)), ABS () represents taking absolute value to the number in () and traversing sequence number k1 in the formula sv (k 1) =ABS (U1 (k 1) -sts (a)) to obtain U1 values sv (1), sv (2), …, sv (U1), sv (2), …, sv (U1) to form a set U1{ } and recording m1 as the sequence number, m 1E [ 1] in the first unit domain, namely, the sequence number is set 304, and the sequence number is set outside the first unit domain;
wherein, the definition of the unit connected with the inner unit is: for any unit Un in the N structure units, when any side of the unit Un is overlapped with any side of the inner unit, the unit Un is called as a unit connected with the inner unit;
s304, if the number of elements in the current sequence Seq is not more than 2, taking the current external unit as an internal unit and turning to S303; if the number of elements in the current sequence Seq is greater than 2, then go to S305;
s305, representing the kth 2 element in the sequence Seq with Seq (k 2), k2 being the sequence number, k2=1, 2, …, N1 being the number of all elements in the current sequence Seq, and forming sn (Seq (1)), sn (Seq (2)), …, sn (Seq (N1-1)) in the N structural unit into an N1 structural unit;
if a unit sn (Seq (S)) connected to sn (Seq (N1)) exists in the N1 structural unit, the process proceeds to S306;
if no connection exists between sn (Seq (N1)) and any of the N1 structural units, taking the current sn (Seq (N1)) as an internal unit and proceeding to S303;
among these, the method for judging whether or not a unit sn (Seq (s)) connected to sn (Seq (N1)) exists in the N1 structural unit is as follows: setting an integer variable i1, wherein the initial value of the integer variable i1 is 1, the value range of the integer variable i1 is [1, N1-1], and starting traversing the variable i1 in the value range of the integer variable i1: when any one side of sn (Seq (N1)) overlaps with any one side of current sn (Seq (i 1)), the value of current variable i1 is recorded as s, sn (Seq (s)) is expressed as a unit connected with sn (Seq (N1)), s is a sequence number, s epsilon [1, N1-1];
s306, in the N structural unit, sn (Seq (S)), sn (Seq (s+1)), …, sn (Seq (N1)) are formed into a superimposed domain; when the number of all units in the superposition area exceeds half of the number of all units in the N structural units, calculating the earthquake resistance in the superposition area, and turning to S307; when the number of all units in the superimposed domain is less than half of the number of all units in the N structural units, then deleting the superimposed domain in the N structural units (i.e., removing sn (Seq (S)) in the N structural units, sn (Seq (s+1)), …, sn (Seq (N1)) the N1 units, taking the N structural units from which the superimposed domain was deleted as new N structural units and going to S302;
s307, taking the earthquake resistance value in the superposition domain as a unit earthquake resistance steady-state analysis result;
the method for calculating the anti-seismic degree in the superposition domain comprises the following steps: the anti-vibration degree in the superposition domain is recorded as basel= [ ave { overspixel }/max { overspixel } ] sum (Pixel); where ave { overspixel } represents an average value of stress values received by all cells in the overlay domain, max { overspixel } represents a stress value of a cell in the overlay domain that receives a maximum stress value, sum (Pixel) represents a sum of stress values received by all cells in the N2 structural unit, and the N2 structural unit is composed of all cells in the N1 structural unit after the current overlay domain is deleted (i.e., all cells in the overlay domain are deleted from the N structural unit, and all remaining cells in the N structural unit are composed of the N2 structural unit).
The beneficial effects of this step are: when the road bridge bears loads such as vehicles and pedestrians, internal stress and deformation can be generated, the stress and the deformation can cause deformation such as bending, torsion, elongation or shrinkage of components of the road bridge, so that the geometric shape and structural stability of the road bridge are affected, meanwhile, in the use environment of the road bridge, such as temperature change or structural strength change of the road bridge caused by internal force factors such as shearing force and axial force, the stress distribution of a model is obtained by utilizing a finite element analysis step, the stress intensity of each unit is calculated, a superposition domain in the road bridge structural model is screened out, the superposition domain is the area with the lowest stability in the road bridge structural model, the earthquake resistance in the superposition domain is used as a unit earthquake resistance steady-state analysis result, the unit earthquake resistance steady-state analysis result reflects the strength of different areas in the whole road bridge relative to the whole road bridge model, the area with lower structural strength in the road bridge model can be found out by utilizing the unit earthquake resistance steady-state analysis result, the position of the low structural strength can not exceed the structural strength range of the road bridge when the road bridge is subjected to load, the safety of the road bridge is greatly improved, meanwhile, the material with the appropriate stability is selected, the material is fully prolonged, and the service life is fully prolonged.
Optionally, in step S300, the step of performing finite element analysis on the road bridge finite element model specifically includes: loading the finite element model of the road bridge into finite element analysis software, setting load information and boundary conditions in the finite element analysis software, and outputting the stress received by each unit after solving the stress received by each unit in the finite element model of the road bridge.
Further, in step S400, the method for optimizing the bridge BIM model based on the unit anti-seismic steady-state analysis result specifically includes:
randomly selecting one unit from the N structural units and marking the unit as sn (x), forming a second unit domain by the sn (x) and all units connected with the sn (x), and carrying out reinforcement design on the position of the second unit domain in the road bridge BIM model when the earthquake resistance in the second unit domain is larger than the unit earthquake resistance steady-state analysis result (namely, the second unit domain is a region with lower structural strength and needs reinforcement treatment); wherein the degree of earthquake resistance in the second unit domain is equal to the sum of stress values experienced by all units in the second unit domain; the definition of the cell that interfaces with sn (x) is: for any one unit Um in the N structure unit, when any one side of the unit Um is overlapped with any one side of sn (x), the unit Um is called as a unit connected with the inner unit.
Because the low-strength position of the road bridge structure is generally in a regional block shape, the overall structural strength of the region is under the combined action of the stress magnitude suffered by each unit in the region, the unit anti-seismic steady-state analysis is often not accurate enough only from the whole stress, the analysis result is affected, the accuracy of the unit anti-seismic steady-state analysis is improved at the same time, and the method for calculating the anti-seismic degree in the superposition domain can be as follows:
s3011, using fie (j 1) to represent the stress value received by the jth 1 unit in the overlapped domain, where j1 is a sequence number, j1=1, 2, …, M is the number of all units in the overlapped domain, fie _a is an average value of the stress values received by all units in the overlapped domain, bell (j 1) = fie (j 1) -fie _a, traversing the sequence number j1 in the formula bell (j 1) = fie (j 1) -fie _a to obtain M values bell (1), bell (2), …, bell (M), sorting bell (1), bell (2), …, bell (M) in ascending order to obtain asc (1), asc (2), …, asc (M), and forming asc (1), asc (2), …, asc (M) into a set asc { and turning to S3012;
s3012, creating a blank array ord, initializing variables k3, wherein the initial value of k3 is 1, k3 epsilon [1, M ], traversing k3 from the initial value of k3, and turning to S3013;
s3013, using asc (n 1) to represent elements in the set asc { and the current bell (k 3) value equal to each other, adding the current n1 into the array ord, and turning to S3014;
s3014, if the value of the variable k3 is smaller than M, increasing the value of the variable k3 by 1, and turning to S3013; if the value of variable k3 is equal to M, go to S3015;
s3015, record the earthquake resistance BaseL in the superposition domain as:
in ln []Representation pair []The number in the method is obtained by natural logarithmic operation, fie (j 1) represents the stress value size received by the j1 st unit in the superposition domain, ord (j 1) represents the j1 st element in the array ord, sum (Pixel) represents the sum of the stress values received by all units in an N2 structural unit, and the N2 structural unit consists of all units in the N1 structural unit after deleting the current superposition domain; k (K) 0 Represents the sum of all elements with negative values in the set asc { }, K 0 I represents the pair K 0 Taking the absolute value.
The beneficial effects of this step are: in the analysis process of the unit earthquake-resistant steady state, the earthquake resistance in the superposition domain is influenced by the joint influence of the stress on each unit in the superposition domain, the stress on each unit can be balanced by the method of the step, the influence weights of different units in different positions in the superposition domain on the earthquake resistance are calculated, the coefficients in the earthquake resistance formula are modified and calculated, the screening accuracy of the low-structural-strength region is effectively improved, and the deflection resistance and the load bearing capacity of the road bridge structure are enhanced.
The beneficial effects of the invention are as follows: the method can effectively improve the shock resistance of the road and bridge structure, optimize the structural design scheme, improve the matching degree with the actual engineering requirement, and adopt the BIM technology to perform shock resistance optimization on the road and bridge structure, thereby realizing digital management and whole process visualization in the design and operation processes, improving the engineering design efficiency, reducing design errors and loopholes and providing reliable technical support for road and bridge engineering.
Drawings
The above and other features of the present invention will become more apparent from the detailed description of the embodiments thereof given in conjunction with the accompanying drawings, in which like reference characters designate like or similar elements, and it is apparent that the drawings in the following description are merely some examples of the present invention, and other drawings may be obtained from these drawings without inventive effort to those of ordinary skill in the art, in which:
FIG. 1 is a flow chart of a BIM model anti-seismic optimization method for road and bridge structures.
Detailed Description
The conception, specific structure, and technical effects produced by the present invention will be clearly and completely described below with reference to the embodiments and the drawings to fully understand the objects, aspects, and effects of the present invention. It should be noted that, in the case of no conflict, the embodiments and features in the embodiments may be combined with each other.
In the description of the present invention, a number means one or more, a number means two or more, and greater than, less than, exceeding, etc. are understood to not include the present number, and above, below, within, etc. are understood to include the present number. The description of the first and second is for the purpose of distinguishing between technical features only and should not be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated or implicitly indicating the precedence of the technical features indicated.
Fig. 1 is a flowchart of a method for optimizing earthquake resistance of a BIM model of a road bridge structure according to the present invention, and a method for optimizing earthquake resistance of a BIM model of a road bridge structure according to an embodiment of the present invention will be described with reference to fig. 1.
The invention provides a BIM model anti-seismic optimization method of a road and bridge structure, which comprises the following steps:
s100, building a road bridge BIM model in BIM software according to the structural information of the road bridge;
s200, converting the road bridge BIM model into a road bridge finite element model through model conversion software;
s300, performing unit anti-seismic steady-state analysis on the road bridge finite element model to obtain a unit anti-seismic steady-state analysis result;
and S400, optimizing the road bridge BIM based on the unit anti-seismic steady-state analysis result.
Further, in step S100, the BIM software is any one of Revit, tekla, BIM5D, bentley, and the road and bridge is any one of a roadbed, a road surface, a bridge, a culvert, and a tunnel.
Further, in step S100, the structural information of the bridge at least includes span of the bridge, load requirement of the bridge, material specification of the bridge, and basic geometry and configuration of the bridge.
Further, in step S100, the method for building the bridge BIM model in the BIM software according to the structural information of the bridge specifically includes: in BIM software, an initial model of a road and a bridge is created according to basic geometric shapes and structures in the structural information of the road and the bridge, load requirements of the road and the bridge in the structural information of the road and the bridge and material specifications of the road and the bridge are added into the initial model of the road and the bridge, road and bridge surface materials and longitudinal and transverse reinforcing steel bars are added on the surface of the initial model of the road and the bridge, and the initial model of the road and the bridge is saved to be the BIM model of the road and the bridge; wherein the road and bridge surface material is reinforced concrete or masonry material.
Further, in step S100, a bridge BIM model is built in the BIM software according to the structural information of the bridge, and the method further includes:
s101, in BIM software, performing collision detection on a road bridge BIM model through a collision detection function, wherein the collision detection is used for detecting model components with collision in the road bridge BIM model, and modifying the collision components in the road bridge BIM model after the collision detection;
s102, circulating S101 until BIM software sends a prompt of no collision detection after collision detection is carried out on the road bridge BIM model.
Further, in step S200, the model conversion software is femtranfer or Hypermesh.
Further, in step S300, the method for performing unit anti-seismic steady-state analysis on the road bridge finite element model to obtain the unit anti-seismic steady-state analysis result specifically includes:
s301, finite element analysis is carried out on a road bridge finite element model, stress values received by each unit in the road bridge finite element model are solved (the road bridge finite element model converted through FEMTransfer comprises a plurality of discrete units which are gathered together to be capable of representing an actual continuous domain, so that the fact that the finite number of unknowns are used for approaching the infinite unknowns is achieved), the number of all units in the road bridge finite element model is recorded as N, and the N units are used as N structural units;
s302, representing the ith unit in the N structural units by sn (i), representing the stress value received by the ith unit in the N structural units by sts (i), wherein i is a sequence number, i=1, 2, … and N, creating a blank array sts, sequentially storing N values sts (1), sts (2), … and sts (N) into the array sts, recording sts (j) as the jth element in the array sts, j=1, 2, … and N, recording the element with the smallest element value in the array sts as sts (a), recording the a-th unit (sn (a)) in the N structural units as an inner unit, creating a blank sequence Seq, adding the sequence number a into the sequence Seq, and turning to S303;
s303, in the N structure units, forming a first unit domain by the inner units and all units connected with the inner units, using U1 (k 1) to represent the stress magnitude of the kth 1 unit in the first unit domain, wherein k1 is a sequence number, k1=1, 2, …, U1, U1 is the number of all units in the first unit domain, recording sv (k 1) =ABS (U1 (k 1) -sts (a)), ABS () represents taking absolute value to the number in () and traversing sequence number k1 in the formula sv (k 1) =ABS (U1 (k 1) -sts (a)) to obtain U1 values sv (1), sv (2), …, sv (U1), sv (2), …, sv (U1) to form a set U1{ } and recording m1 as the sequence number, m 1E [ 1] in the first unit domain, namely, the sequence number is set 304, and the sequence number is set outside the first unit domain;
wherein, the definition of the unit connected with the inner unit is: for any unit Un in the N structure units, when any side of the unit Un is overlapped with any side of the inner unit, the unit Un is called as a unit connected with the inner unit;
s304, if the number of elements in the current sequence Seq is not more than 2, taking the current external unit as an internal unit and turning to S303; if the number of elements in the current sequence Seq is greater than 2, then go to S305;
s305, representing the kth 2 element in the sequence Seq with Seq (k 2), k2 being the sequence number, k2=1, 2, …, N1 being the number of all elements in the current sequence Seq, and forming sn (Seq (1)), sn (Seq (2)), …, sn (Seq (N1-1)) in the N structural unit into an N1 structural unit;
if a unit sn (Seq (S)) connected to sn (Seq (N1)) exists in the N1 structural unit, the process proceeds to S306;
if no connection exists between sn (Seq (N1)) and any of the N1 structural units, taking the current sn (Seq (N1)) as an internal unit and proceeding to S303;
among these, the method for judging whether or not a unit sn (Seq (s)) connected to sn (Seq (N1)) exists in the N1 structural unit is as follows: setting an integer variable i1, wherein the initial value of the integer variable i1 is 1, the value range of the integer variable i1 is [1, N1-1], and starting traversing the variable i1 in the value range of the integer variable i1: when any one side of sn (Seq (N1)) overlaps with any one side of current sn (Seq (i 1)), the value of current variable i1 is recorded as s, sn (Seq (s)) is expressed as a unit connected with sn (Seq (N1)), s is a sequence number, s epsilon [1, N1-1];
s306, in the N structural unit, sn (Seq (S)), sn (Seq (s+1)), …, sn (Seq (N1)) are formed into a superimposed domain; when the number of all units in the superposition area exceeds half of the number of all units in the N structural units, calculating the earthquake resistance in the superposition area, and turning to S307; when the number of all units in the superimposed domain is less than half of the number of all units in the N structural units, then deleting the superimposed domain in the N structural units (i.e., removing sn (Seq (S)) in the N structural units, sn (Seq (s+1)), …, sn (Seq (N1)) the N1 units, taking the N structural units from which the superimposed domain was deleted as new N structural units and going to S302;
s307, taking the earthquake resistance value in the superposition domain as a unit earthquake resistance steady-state analysis result;
the method for calculating the anti-seismic degree in the superposition domain comprises the following steps: the anti-vibration degree in the superposition domain is recorded as basel= [ ave { overspixel }/max { overspixel } ] sum (Pixel); where ave { overspixel } represents an average value of stress values received by all cells in the overlay domain, max { overspixel } represents a stress value of a cell in the overlay domain that receives a maximum stress value, sum (Pixel) represents a sum of stress values received by all cells in the N2 structural unit, and the N2 structural unit is composed of all cells in the N1 structural unit after the current overlay domain is deleted (i.e., all cells in the overlay domain are deleted from the N structural unit, and all remaining cells in the N structural unit are composed of the N2 structural unit).
Optionally, in step S300, the step of performing finite element analysis on the road bridge finite element model specifically includes: loading the finite element model of the road bridge into finite element analysis software, setting load information and boundary conditions in the finite element analysis software, and outputting the stress received by each unit after solving the stress received by each unit in the finite element model of the road bridge.
Further, in step S400, the method for optimizing the bridge BIM model based on the unit anti-seismic steady-state analysis result specifically includes:
randomly selecting one unit from the N structural units and marking the unit as sn (x), forming a second unit domain by the sn (x) and all units connected with the sn (x), and carrying out reinforcement design on the position of the second unit domain in the road bridge BIM model when the earthquake resistance in the second unit domain is larger than the unit earthquake resistance steady-state analysis result (namely, the second unit domain is a region with lower structural strength and needs reinforcement treatment); wherein the degree of earthquake resistance in the second unit domain is equal to the sum of stress values experienced by all units in the second unit domain; the definition of the cell that interfaces with sn (x) is: for any one unit Um in the N structure unit, when any one side of the unit Um is overlapped with any one side of sn (x), the unit Um is called as a unit connected with the inner unit.
Because the low-strength position of the road bridge structure is generally in a regional block shape, the overall structural strength of the region is under the combined action of the stress magnitude suffered by each unit in the region, the unit anti-seismic steady-state analysis is often not accurate enough only from the whole stress, the analysis result is affected, the accuracy of the unit anti-seismic steady-state analysis is improved at the same time, and the method for calculating the anti-seismic degree in the superposition domain can be as follows:
s3011, using fie (j 1) to represent the stress value received by the jth 1 unit in the overlapped domain, where j1 is a sequence number, j1=1, 2, …, M is the number of all units in the overlapped domain, fie _a is an average value of the stress values received by all units in the overlapped domain, bell (j 1) = fie (j 1) -fie _a, traversing the sequence number j1 in the formula bell (j 1) = fie (j 1) -fie _a to obtain M values bell (1), bell (2), …, bell (M), sorting bell (1), bell (2), …, bell (M) in ascending order to obtain asc (1), asc (2), …, asc (M), and forming asc (1), asc (2), …, asc (M) into a set asc { and turning to S3012;
s3012, creating a blank array ord, initializing variables k3, wherein the initial value of k3 is 1, k3 epsilon [1, M ], traversing k3 from the initial value of k3, and turning to S3013;
s3013, using asc (n 1) to represent elements in the set asc { and the current bell (k 3) value equal to each other, adding the current n1 into the array ord, and turning to S3014;
s3014, if the value of the variable k3 is smaller than M, increasing the value of the variable k3 by 1, and turning to S3013; if the value of variable k3 is equal to M, go to S3015;
s3015, record the earthquake resistance BaseL in the superposition domain as:
in ln []Representation pair []The number in the method is obtained by natural logarithmic operation, fie (j 1) represents the stress value size received by the j1 st unit in the superposition domain, ord (j 1) represents the j1 st element in the array ord, sum (Pixel) represents the sum of the stress values received by all units in an N2 structural unit, and the N2 structural unit consists of all units in the N1 structural unit after deleting the current superposition domain; k (K) 0 Represents the sum of all elements with negative values in the set asc { }, K 0 I represents the pair K 0 Taking the absolute value.
The invention provides a BIM model anti-seismic optimization method of a road and bridge structure, which comprises the steps of establishing a road and bridge BIM model in BIM software according to structural information of a road and bridge, converting the road and bridge BIM model into a road and bridge finite element model through model conversion software, carrying out unit anti-seismic steady-state analysis on the road and bridge finite element model to obtain a unit anti-seismic steady-state analysis result, and optimizing the road and bridge BIM model based on the unit anti-seismic steady-state analysis result. The method can effectively improve the shock resistance of the road and bridge structure, optimize the structural design scheme, improve the matching degree with the actual engineering requirement, and adopt the BIM technology to perform shock resistance optimization on the road and bridge structure, thereby realizing digital management and whole process visualization in the design and operation processes, improving the engineering design efficiency, reducing design errors and loopholes and providing reliable technical support for road and bridge engineering. Although the present invention has been described in considerable detail and with particularity with respect to several described embodiments, it is not intended to be limited to any such detail or embodiment or any particular embodiment so as to effectively cover the intended scope of the invention. Furthermore, the foregoing description of the invention has been presented in its embodiments contemplated by the inventors for the purpose of providing a useful description, and for the purposes of providing a non-essential modification of the invention that may not be presently contemplated, may represent an equivalent modification of the invention.
Claims (7)
1. The BIM model anti-seismic optimization method of the road and bridge structure is characterized by comprising the following steps of:
s100, building a road bridge BIM model in BIM software according to the structural information of the road bridge;
s200, converting the road bridge BIM model into a road bridge finite element model through model conversion software;
s300, performing unit anti-seismic steady-state analysis on the road bridge finite element model to obtain a unit anti-seismic steady-state analysis result;
s400, optimizing a road bridge BIM model based on a unit anti-seismic steady-state analysis result;
in step S300, a unit anti-seismic steady-state analysis is performed on the road bridge finite element model, and the method for obtaining the unit anti-seismic steady-state analysis result specifically includes:
s301, finite element analysis is carried out on the road bridge finite element model, stress values received by each unit in the road bridge finite element model are solved, the number of all units in the road bridge finite element model is recorded as N, and the N units are taken as N structural units;
s302, representing the ith unit in the N structural units by sn (i), representing the stress value received by the ith unit in the N structural units by sts (i), wherein i is a sequence number, i=1, 2, … and N, creating a blank array sts, sequentially storing N values sts (1), sts (2), … and sts (N) into the array sts, recording sts (j) as the jth element in the array sts, j=1, 2, … and N, recording the element with the smallest element value in the array sts as sts (a), recording the a th unit in the N structural units as an inner unit, creating a blank sequence Seq, adding the sequence number a into the sequence Seq, and turning to S303;
s303, in the N structure units, forming a first unit domain by the inner units and all units connected with the inner units, using U1 (k 1) to represent the stress magnitude of the kth 1 unit in the first unit domain, wherein k1 is a sequence number, k1=1, 2, …, U1, U1 is the number of all units in the first unit domain, recording sv (k 1) =ABS (U1 (k 1) -sts (a)), ABS () represents taking absolute value to the number in () and traversing sequence number k1 in the formula sv (k 1) =ABS (U1 (k 1) -sts (a)) to obtain U1 values sv (1), sv (2), …, sv (U1), sv (2), …, sv (U1) to form a set U1{ } and recording m1 as the sequence number, m 1E [ 1] in the first unit domain, namely, the sequence number is set 304, and the sequence number is set outside the first unit domain;
wherein, the definition of the unit connected with the inner unit is: for any unit Un in the N structure units, when any side of the unit Un is overlapped with any side of the inner unit, the unit Un is called as a unit connected with the inner unit;
s304, if the number of elements in the current sequence Seq is not more than 2, taking the current external unit as an internal unit and turning to S303; if the number of elements in the current sequence Seq is greater than 2, then go to S305;
s305, representing the kth 2 element in the sequence Seq with Seq (k 2), k2 being the sequence number, k2=1, 2, …, N1 being the number of all elements in the current sequence Seq, and forming sn (Seq (1)), sn (Seq (2)), …, sn (Seq (N1-1)) in the N structural unit into an N1 structural unit;
if a unit sn (Seq (S)) connected to sn (Seq (N1)) exists in the N1 structural unit, the process proceeds to S306;
if no connection exists between sn (Seq (N1)) and any of the N1 structural units, taking the current sn (Seq (N1)) as an internal unit and proceeding to S303;
among these, the method for judging whether or not a unit sn (Seq (s)) connected to sn (Seq (N1)) exists in the N1 structural unit is as follows: setting an integer variable i1, wherein the initial value of the integer variable i1 is 1, the value range of the integer variable i1 is [1, N1-1], and starting traversing the variable i1 in the value range of the integer variable i1: when any one side of sn (Seq (N1)) overlaps with any one side of current sn (Seq (i 1)), the value of current variable i1 is recorded as s, sn (Seq (s)) is expressed as a unit connected with sn (Seq (N1)), s is a sequence number, s epsilon [1, N1-1];
s306, in the N structural unit, sn (Seq (S)), sn (Seq (s+1)), …, sn (Seq (N1)) are formed into a superimposed domain; when the number of all units in the superposition area exceeds half of the number of all units in the N structural units, calculating the earthquake resistance in the superposition area, and turning to S307; when the number of all units in the superposition domain is less than half of the number of all units in the N structural units, deleting the superposition domain in the N structural units, taking the N structural units with the deleted superposition domain as new N structural units, and transferring to S302;
s307, taking the earthquake resistance value in the superposition domain as a unit earthquake resistance steady-state analysis result;
the method for calculating the anti-seismic degree in the superposition domain comprises the following steps: the anti-vibration degree in the superposition domain is recorded as basel= [ ave { overspixel }/max { overspixel } ] sum (Pixel); wherein ave { overspixel } represents an average value of stress values received by all cells in the superimposed domain, max { overspixel } represents a stress value of a cell having a largest stress value received in the superimposed domain, sum (Pixel) represents a sum of stress values received by all cells in an N2 structural unit, and the N2 structural unit is composed of all cells in the N1 structural unit from which the current superimposed domain is deleted.
2. The method for optimizing earthquake resistance of a BIM model of a road and bridge structure according to claim 1, wherein in step S100, the BIM software is any one of Revit, tekla, BIM5D, bentley, and the road and bridge is any one of a roadbed, a road surface, a bridge, a culvert and a tunnel.
3. The method for optimizing the earthquake resistance of a BIM model of a road and bridge structure according to claim 1, wherein in step S100, the structural information of the road and bridge at least includes a span of the road and bridge, a load requirement of the road and bridge, a material specification of the road and bridge, a basic geometric shape and a construction of the road and bridge.
4. The method for building a bridge BIM model in BIM software according to the structure information of the bridge in step S100 is specifically as follows: in BIM software, an initial model of a road and a bridge is created according to basic geometric shapes and structures in the structural information of the road and the bridge, load requirements of the road and the bridge in the structural information of the road and the bridge and material specifications of the road and the bridge are added into the initial model of the road and the bridge, road and bridge surface materials and longitudinal and transverse reinforcing steel bars are added on the surface of the initial model of the road and the bridge, and the initial model of the road and the bridge is saved to be the BIM model of the road and the bridge; wherein the road and bridge surface material is reinforced concrete or masonry material.
5. The method for seismic optimization of a BIM model of a road and bridge structure according to claim 4, wherein in step S100, a road and bridge BIM model is built in the BIM software according to the structural information of the road and bridge, and further comprising:
s101, in BIM software, performing collision detection on a road bridge BIM model through a collision detection function, wherein the collision detection is used for detecting model components with collision in the road bridge BIM model, and modifying the collision components in the road bridge BIM model after the collision detection;
s102, circulating S101 until BIM software sends a prompt of no collision detection after collision detection is carried out on the road bridge BIM model.
6. The method for optimizing earthquake resistance of a BIM model of a road bridge structure according to claim 1, wherein in step S200, the model conversion software is femtranfer or Hypermesh.
7. The method for optimizing the BIM model of the road and bridge structure according to claim 1, wherein in the step S400, the method for optimizing the BIM model of the road and bridge based on the unit anti-seismic steady-state analysis result is specifically as follows:
randomly selecting one unit from the N structural units and marking the unit as sn (x), forming a second unit domain by the sn (x) and all units connected with the sn (x), and carrying out reinforcement design on the position of the second unit domain in the road bridge BIM model when the earthquake resistance in the second unit domain is larger than the unit earthquake resistance steady-state analysis result; wherein the degree of earthquake resistance in the second unit domain is equal to the sum of stress values experienced by all units in the second unit domain; the definition of the cell that interfaces with sn (x) is: for any one unit Um in the N structure unit, when any one side of the unit Um is overlapped with any one side of sn (x), the unit Um is called as a unit connected with the inner unit.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310361135.XA CN116090078B (en) | 2023-04-07 | 2023-04-07 | BIM model anti-seismic optimization method of road bridge structure |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310361135.XA CN116090078B (en) | 2023-04-07 | 2023-04-07 | BIM model anti-seismic optimization method of road bridge structure |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116090078A CN116090078A (en) | 2023-05-09 |
| CN116090078B true CN116090078B (en) | 2023-06-09 |
Family
ID=86210423
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310361135.XA Active CN116090078B (en) | 2023-04-07 | 2023-04-07 | BIM model anti-seismic optimization method of road bridge structure |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116090078B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006195713A (en) * | 2005-01-13 | 2006-07-27 | Mitsui Eng & Shipbuild Co Ltd | Method for supporting design of aseismatic reinforcement for bridge |
| JP2010152857A (en) * | 2008-11-21 | 2010-07-08 | Ihi Corp | System, method and program for designing structure |
| CN111611634A (en) * | 2020-03-30 | 2020-09-01 | 南京震坤物联网科技有限公司 | Bridge health assessment system and health assessment method based on BIM-FEM |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108460231A (en) * | 2018-03-23 | 2018-08-28 | 中交公路长大桥建设国家工程研究中心有限公司 | A kind of bridge builds foster overall process intellectual monitoring assessment early warning decision system and method |
| EP3877600A4 (en) * | 2018-11-08 | 2022-08-03 | The Regents Of The University Of Michigan | Modular construction of load-bearing structures from reusable load-bearing elements |
| CN110263460A (en) * | 2019-06-26 | 2019-09-20 | 江苏工程职业技术学院 | A kind of assembled 4D construction simulation safety monitoring system based on BIM |
| CN110984366B (en) * | 2019-12-24 | 2021-05-11 | 北京华瑞国盾信息科技有限公司 | Construction method of assembly type building structure system of oil and gas pipeline engineering |
| US20210287138A1 (en) * | 2020-03-12 | 2021-09-16 | Autodesk, Inc. | Learning to simulate and design for structural engineering |
| CN114707377A (en) * | 2022-03-22 | 2022-07-05 | 大连民族大学 | Revit-OpenSEES-based reinforced concrete frame structure earthquake-resistant analysis method |
| CN115618479B (en) * | 2022-12-20 | 2023-03-28 | 广东卓智设计工程有限公司 | BIM-based curtain wall high-sealing installation method |
-
2023
- 2023-04-07 CN CN202310361135.XA patent/CN116090078B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006195713A (en) * | 2005-01-13 | 2006-07-27 | Mitsui Eng & Shipbuild Co Ltd | Method for supporting design of aseismatic reinforcement for bridge |
| JP2010152857A (en) * | 2008-11-21 | 2010-07-08 | Ihi Corp | System, method and program for designing structure |
| CN111611634A (en) * | 2020-03-30 | 2020-09-01 | 南京震坤物联网科技有限公司 | Bridge health assessment system and health assessment method based on BIM-FEM |
Non-Patent Citations (1)
| Title |
|---|
| 高速铁路系杆拱桥BIM技术研究;谢兴定;王振宇;;价值工程(第18期);第61-67页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116090078A (en) | 2023-05-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Oehlers et al. | Elementary behaviour of composite steel and concrete structural members | |
| Martínez et al. | Heuristic optimization of RC bridge piers with rectangular hollow sections | |
| JP3569899B2 (en) | Bridge damage estimation system and program | |
| Akgül et al. | Computational platform for predicting lifetime system reliability profiles for different structure types in a network | |
| CN110765534A (en) | Method for optimizing cable force of finished bridge of cable-stayed bridge | |
| Kurc et al. | A comparative study on different analysis approaches for estimating the axial loads on columns and structural walls at tall buildings | |
| Dang et al. | Bidirectional loading hybrid tests of square cross‐sections of steel bridge piers | |
| CN114329709A (en) | Method for rapidly diagnosing service performance and predicting development trend of shield tunnel | |
| CN116090078B (en) | BIM model anti-seismic optimization method of road bridge structure | |
| Islam et al. | An optimized design of network arch bridge using global optimization algorithm | |
| CN116108548B (en) | BIM-based road bridge structural strength analysis method and system | |
| Cheng et al. | Reliability analysis of a long span steel arch bridge against wind-induced stability failure during construction | |
| Tan et al. | A rapid evaluation method based on natural frequency for post-earthquake traffic capacity of small and medium span bridges | |
| Zhang et al. | Transfer relations: useful basis for computer-aided engineering of circular arch structures | |
| Xu et al. | Development of BP-based seismic behavior optimization of RC and steel frame structures | |
| Kazemi et al. | Effect of epsilon-based record selection on fragility curves of typical irregular steel frames with concrete shear walls in Mashhad city | |
| Sotelino et al. | Simplified load distribution factor for use in LRFD design | |
| Shittu et al. | A stochastic framework for the assessment of horizontally curved aluminium bridge decks on steel girders | |
| Suwal et al. | Probabilistic Evaluation of Modeling Parameters for Reinforced Concrete Moment Frame Building | |
| Ashraf et al. | Sway of semi-rigid steel frames: Part 1: Regular frames | |
| CN112629896A (en) | Beam structure damage identification method based on horizontal support reaction influence line | |
| CN117473620A (en) | Bridge reconstruction construction method integrating BIM technology and algorithm model | |
| Cîrdei et al. | Bridge dynamic displacement monitoring using geodetic measurement. | |
| Erduran et al. | Component damage functions for reinforced concrete frame structures | |
| CN113392450B (en) | Method, device and system for evaluating health state of railway pier |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |