CN116561875A - Bridge network vulnerability analysis method considering bridge seismic response correlation - Google Patents
Bridge network vulnerability analysis method considering bridge seismic response correlation Download PDFInfo
- Publication number
- CN116561875A CN116561875A CN202310827938.XA CN202310827938A CN116561875A CN 116561875 A CN116561875 A CN 116561875A CN 202310827938 A CN202310827938 A CN 202310827938A CN 116561875 A CN116561875 A CN 116561875A
- Authority
- CN
- China
- Prior art keywords
- bridge
- vulnerability
- road section
- correlation
- network
- 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.)
- Granted
Links
- 230000004044 response Effects 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000012038 vulnerability analysis Methods 0.000 title claims abstract description 19
- 238000006073 displacement reaction Methods 0.000 claims abstract description 34
- 238000004364 calculation method Methods 0.000 claims description 43
- 238000000611 regression analysis Methods 0.000 claims description 5
- 102220093746 rs876661027 Human genes 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 description 27
- 238000010586 diagram Methods 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 239000000203 mixture Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
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)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Structural Engineering (AREA)
- Civil Engineering (AREA)
- Architecture (AREA)
- Computer Networks & Wireless Communication (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to the technical field of traffic network vulnerability analysis, in particular to a bridge network vulnerability analysis method considering bridge seismic response correlation, which comprises the following steps: calculating the maximum drift ratio of bridge piers of each bridge; calculating the correlation between any two bridges in the same road section through a correlation formula; establishing the association between the correlation of the bridge and the maximum drift ratio of the corresponding bridge piers, and selecting a preset number of bridge seismic displacement responses through the association; comparing the maximum drift ratio of the bridge piers of each bridge with the limit value of the maximum drift ratio of the bridge piers, and calculating the vulnerability of each road section according to the comparison result; judging the magnitude relation between the vulnerability of each road section in the bridge network and the set value, and calculating the vulnerability of the bridge network according to the judging result; the method can accurately calculate the seismic vulnerability of the bridge network with correlation among the bridges.
Description
Technical Field
The invention relates to the technical field of traffic network vulnerability analysis, in particular to a bridge network vulnerability analysis method considering bridge seismic response correlation.
Background
The bridge is used as a ring in the traffic network and can play a role in connecting roads. However, bridges are more easily damaged under the action of an earthquake than roads, and the post-earthquake state of the bridge directly affects the traffic capacity of the traffic network. Therefore, it is necessary to evaluate the post-earthquake state of the bridge and further determine the actual traffic capacity of the traffic network.
When evaluating the post-earthquake state of a bridge, the vulnerability of the bridge is generally used as a measurement index. The traditional vulnerability calculation method needs to carry out fine modeling on the bridges, and for bridges with a small number, the accuracy of the method is high. However, in a practical traffic network, the number of bridges is huge, the constructed bridge network is complex, and there are a large number of similar bridges in the complex network, which results in that a large cost is spent and accuracy is lowered by adopting the method.
The reason that a large number of similar bridges exist is that the design and construction methods are similar, the reference specifications are the same, so that the bridges have quite similar structural characteristics in aspects of types, bridge spans, seismic reduction and isolation measures and the like, and correlation is established between the bridges. The method is affected by regions, similar traffic flows exist on each road section or bridge on the bridge network, and the same extreme events are experienced under similar environmental conditions, so that the inter-bridge compatibility is more obvious when an earthquake happens. The traditional vulnerability calculation method does not consider the correlation between bridges, which is one of the important reasons for low calculation accuracy.
Disclosure of Invention
In order to avoid and overcome the technical problems in the prior art, the invention provides a bridge network vulnerability analysis method considering the correlation of bridge seismic response. The method can accurately calculate the seismic vulnerability of the bridge network with correlation among the bridges.
In order to achieve the above purpose, the present invention provides the following technical solutions:
a bridge network vulnerability analysis method considering bridge seismic response correlation comprises the following steps:
s1, calculating the maximum drift ratio of bridge piers of each bridge;
s2, calculating the correlation between any two bridges in the same road section through a correlation formula;
s3, establishing association between the correlation of the bridge and the corresponding maximum drift ratio of the bridge piers, and selecting a preset number of maximum drift ratios of the bridge piers through the association; comparing the maximum drift ratio of the bridge piers of each bridge with the limit value of the maximum drift ratio of the bridge piers, and calculating the vulnerability of each road section according to the comparison result;
and S4, judging the magnitude relation between the vulnerability of each road section in the bridge network and the set value, and calculating the vulnerability of the bridge network according to the judging result.
As still further aspects of the invention: the specific steps of step S1 are as follows:
s11, firstly, acquiring a calculation formula of the bridge seismic displacement response, wherein the calculation formula is as follows:
wherein, μrepresenting the median value of the seismic displacement response of the bridge;Trepresenting the fundamental period of the bridge;PGVrepresenting peak ground speed;a,b,cthe parameters are parameters of a bridge seismic displacement response calculation formula;
s12, the basic period of each bridge is shortenedT、Peak ground speedPGVCorresponding median value of bridge seismic displacement responseμIntroducing into a pier maximum drift ratio calculation formula, and solving parameters by regression analysisa,b,cThe specific value of (3);
s13, median value of earthquake displacement response through bridgeμCalculating the maximum drift ratio of bridge piersM,MThe calculation formula of (2) is as follows:
wherein, Mrepresenting the maximum drift ratio of the bridge pier;Hthe height of the bridge pier is set;
s14, acquiring the basic period of each bridge in the step S12T、Bridge pier heightHPeak ground speed of the bridge locationPGVAnd leading the values into a bridge pier maximum drift ratio calculation formula to calculate so as to obtain the bridge pier maximum drift ratio of each bridge.
As still further aspects of the invention: the specific steps of the step S2 are as follows:
s21, acquiring the basic period of each bridge in the same road sectionTRespectively comparing the basic periods of any two bridges in the road section, and generating a basic period ratio;
s22, calculating covariance between seismic displacement responses of any two bridges in the road section according to a covariance calculation principle;
s23, substituting the calculated basic cycle ratio and the corresponding covariance into a correlation calculation formula, fitting the calculation formula, and solving the specific numerical value of each parameter in the calculation formula through regression analysis;
s24, substituting the calculated basic period ratios in the road section into a correlation calculation formula determined by parameters to deduce the corresponding correlation between the bridges.
As still further aspects of the invention: the specific steps of step S3 are as follows:
s31, selecting a road section, and calculating the ground speed of each bridge in the road section at each peak valuePGVMaximum drift ratio of bridge pierM;
S32, selectingN L Group same peak ground speedPGVMaximum drift ratio of bridge pier of each lower bridgeM;
S33, selecting the maximum drift ratio of each pierMLimit value of maximum drift ratio of bridge pierM max Comparing the sizes;
s34, counting maximum drift ratio of bridge pierMIs larger than the limit value of the maximum drift ratio of the bridge pierM max The number of groups of (a)n L Will beN L Andn L introducing a road section vulnerability probability function representing vulnerability of the road section to obtain the ground speed of the road section at the peak valuePGVThe vulnerable probability of the lower road section;
s35, calculating the ground speed of each road section in the bridge network at each peak value by using the road section vulnerable probability functionPGVThe probability of vulnerability of the lower road section.
As still further aspects of the invention: the specific steps of step S4 are as follows:
s41, intersecting and connecting all road sections through the end points of the road sections, and forming corresponding paths, wherein all paths are combined to form a bridge network;
s42, randomly selecting each road section from the set rangeN W Group random numberP W As an observation value, the link vulnerability probability of the link is comparedP L And observed valueP W Is of a size of (2);
s43, if the vulnerable probability of the road section exists in the pathP L Greater than the corresponding observed valueP W The path does not exist, otherwise it exists;
s44, statisticsN W Group number of group random numbers such that no path exists between the start point and the end point of bridge networkn W And willN W Andn W and (5) leading the probability of the vulnerability of the bridge network into a calculation formula of the probability of the vulnerability of the bridge network to calculate the probability of the vulnerability of the bridge network.
As still further aspects of the invention: the basic cycle ratio is calculated as follows:
wherein, R k,i,j represent the firstkIn road segmentsiBridge seat and the firstjThe basic cycle ratio of the seat bridge;T k,i represent the firstkIn road segmentsiThe basic period of the bridge seat;T k,j represent the firstkIn road segmentsjBasic cycle of the bridge.
As still further aspects of the invention: the correlation calculation formula is expressed as follows:
wherein, Cov[i,j]represent the firstiBridge seat and the firstjCovariance between base bridge fundamental periods, the covariance characterizing the firstiBridge seat and the firstjCorrelation between seat bridges;eandfare parameters of the correlation calculation formula.
As still further aspects of the invention: the calculation formula of the road section vulnerable probability function is as follows:
wherein, P L the link vulnerability probability of the link is represented.
As still further aspects of the invention: the calculation formula of the probability of vulnerability of the bridge network is as follows:
wherein, P W representing the probability of vulnerability of the bridge network.
As still further aspects of the invention: the setting range in step S42 is [0,1].
Compared with the prior art, the invention has the beneficial effects that:
1. the method and the device can accurately calculate the vulnerability of the bridge network with correlation among the bridges.
2. The method and the device can calculate the earthquake displacement response of the bridge pier more quickly and obtain the correlation between the bridge displacement responses more quickly.
Drawings
FIG. 1 is a schematic diagram of the operation flow structure of the present invention.
FIG. 2 is a scatter plot of median values of seismic displacement responses of 20 bridges in accordance with the present invention.
FIG. 3 is a spatial curved surface diagram of median composition of seismic displacement responses of 20 bridges in the invention.
FIG. 4 is a correlation data distribution diagram of 20 bridges in the present invention.
Fig. 5 is a network diagram of a bridge formed by 20 bridges according to the present invention.
Fig. 6 is a network diagram of a bridge formed by 20 bridges after earthquake in the invention.
Fig. 7 is a graph of vulnerability of the bridge network from point a to point E in the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Because bridges are more easily damaged than roads in an earthquake, the invention assumes that the traffic capacity of the roads remains intact when studying the earthquake vulnerability of road sections and traffic networks, i.e. does not consider the influence of the damage probability of the roads. The vulnerability of a road segment may be obtained from bridges on the road segment. Likewise, the vulnerability of the bridge network may be obtained from a combination of vulnerabilities of road segments. Generally expressed as:
(1)
(2)
wherein, P road section Representing the vulnerable probability of the road section;P mbridge-shaped structure Represent the firstmThe probability of the easy damage of the seat bridge,Mrepresenting the total number of bridges in the road section;P network system Representing the probability of vulnerability of the bridge network;P nroad segment- Represent the firstnThe probability of vulnerability of the road segment,Nrepresenting the total number of road segments in the bridge network.
In order to achieve the above object, the present invention thus employs the operation steps as shown in fig. 1:
s1, calculating the maximum drift ratio of bridge piers of each bridge.
Median seismic displacement response of bridge based on assumptions of Probabilistic Seismic Demand Model (PSDM)μAnd earthquake intensity indexIMThe logarithmic linear relation exists, and the logarithmic linear relation is a bridge seismic displacement response calculation formula, and the calculation formula is as follows:
or->(3)
Wherein, a,cthe parameters are parameters of a bridge seismic displacement response calculation formula.
By peak ground speedPGVRepresenting the intensity of earthquake motionIM. The bridge seismic displacement response calculation formula becomes the following form:
(4)
and carrying out nonlinear time-course analysis on bridges with 20 different basic periods to obtain the median value of the seismic displacement response. Research finds that the median value of earthquake displacement response and the fundamental period of bridgeTSince there is a quantitative relationship, the formula (4) can be rewritten as:
(5)
wherein, band calculating parameters of a formula for the bridge seismic displacement response.
The data obtained by nonlinear time-course analysis of the bridges with 20 different basic periods are shown in fig. 2, and the data in fig. 2 are subjected to statistical regression analysis according to a formula (5), and the result is shown in fig. 3. In combination with fig. 2 and 3, givenPGVUnder the condition of (1), obtaining the quantitative relationship between the median value of the seismic displacement response of the bridge and the fundamental period of the bridge, and determining the parametersaTaking the weight of the mixture to be 0.36,btaking the sample of 0.82 of the total weight of the mixture,ctaking 1.21, equation (5) can be rewritten as:
(6)
maximum drift ratio of bridge pierMPDMedian value of seismic displacement response to bridgeμThe relation of (2) is:MPD=μ/Husing charactersMAlternate representationMPD. The bridge pier height is also added to the formula (6)HEquation (6) translates to equation (7), equation (7) being as follows:
(7)
according to formula (7), at the basic period of the known bridgeTAnd pier heightHOn the premise of the method, the median value of the earthquake displacement response of the bridge can be rapidly obtained. The parameters of the bridge sections for 20 different fundamental periods are shown in table 1.
Table 1 bridge parameters
S2, calculating the correlation between any two bridges in the same road section through a correlation formula.
The covariance reflects the cooperative relation between the two variables, and when the covariance is positive, and the larger the covariance is, the larger the degree of the same direction of the two variables is; the two variables change inversely when the covariance is negative. The invention uses covariance to describe the correlation of bridge seismic displacement response. In PSDM, it is generally assumed that the seismic displacement response and seismic capacity of the bridge are log-normal distribution-compliant. Thus, when considering the correlation of the bridge seismic displacement response, the natural logarithm of the bridge seismic displacement response is taken first, and then the covariance is calculated. The 20 models were analyzed for 400 sets of covariance and the data obtained are shown in fig. 4. The following quantization model is given to determine the relationship between covariance and fundamental period ratio.
(8)
Wherein, R k,i,j represent the firstkIn road segmentsiBridge seat and the firstjThe basic cycle ratio of the seat bridge;T k,i represent the firstkIn road segmentsiThe basic period of the bridge seat;T k,j represent the firstkIn road segmentsjBasic cycle of the bridge.
(9)
Wherein, Cov[i,j]represent the firstiBridge seat and the firstjCovariance between base bridge fundamental periods, the covariance characterizing the firstiBridge seat and the firstjCorrelation between seat bridges;eandfare parameters of the correlation calculation formula.
Substituting the data into equation (9) can be done as follows from the data scatter and fit curve represented in fig. 4: parameters (parameters)eTaking 0.58, parametersfTaking-0.50, then equation (9) is updated to equation (10), equation (10) being expressed as follows:
(10)
through the formula (10), the correlation between bridge seismic displacement responses can be obtained quickly when the basic period of the bridge is known.
S3, establishing association between the correlation of the bridge and the corresponding maximum drift ratio of the bridge piers, and selecting the maximum drift ratio of the bridge piers with a preset number through the association; and comparing the maximum drift ratio of the bridge piers with the limit value of the maximum drift ratio of the bridge piers, and calculating the vulnerability of each road section according to the comparison result.
When calculating the vulnerability of the road sections containing the bridge, the correlation between the bridge seismic displacement responses needs to be considered.
The invention acquires a bridge structure response sample considering correlation by means of a multi-element normal sampling function in a Python platform.
The Monte Carlo method is used for simulating the vulnerability of the road section by considering the completely uncorrelated unconnected probability calculation principle with reference to a series system.
The detailed steps for calculating the vulnerability of the road section are as follows:
1. for a certain road section, selecting one road sectionPGV。
2. Obtaining the maximum drift ratio of the bridge pier in the step S1MAnd that obtained in step S2Cov[i,j]As a parameter of the multiple normal sampling function, by selectingN L Group same peak ground speedVThe lower pier is the mostLarge drift ratioM. In practice, a bridge is at a peak ground speedPGVMaximum drift ratio of bridge pierMIs a value, which can be recorded as a numberM 1 . However, if the earthquake happens again under the same condition, the bridge pier maximum drift ratio of the bridgeMIs marked asM 2 。M 1 And (3) withM 2 Are generally not equal, the same happensN L Group same peak ground speedPGVDuring the earthquake, the bridge will correspondingly generateN L Personal (S)M. Among the seismic vulnerability analysis systems we consider thisN L Personal (S)MObeys normal distribution, and because the invention considers a plurality of bridges, each bridgeN L Personal (S)MFor example, 20 bridges have 20×N L Personal (S)MIs subject to a multivariate normal distribution of mean and covariance.
Sequentially sorting according to the measured time to obtain 20 bridges, wherein each bridge is respectively measured into 10000 groups, namelyN L Each group comprises the maximum drift ratio of the corresponding bridge piers generated by each bridgeM。
3. Judging whether the selected data sets have the maximum drift ratio of the bridge pierMIs larger than the limit value of the maximum drift ratio of the bridge pierM max If so, the bridge damage is indicated, and the variable is countedcount 1 The number of the two-way valve is added with one,count 1 is 0; if not, count variablecount 1 Is unchanged. Up toN L The group data were all analyzed at this timecount 1 =n L RepresentingN L Among the group data aren L Maximum drift ratio of bridge piersMIs larger than the limit value of the maximum drift ratio of the bridge pierM max Is the case in (a).count 1 /N L The vulnerable probability of the road section is expressed as follows by adopting a formula (11):
(11)
wherein, P L the link vulnerability probability of the link is represented.
4. Modification ofPGVRepeating the steps to obtain the vulnerable probability of the road section under each peak ground speed, wherein the corresponding discrete data point is #PGV,P L ) A smooth line is used to connect the discrete data points together to form a vulnerability curve.
And S4, judging the magnitude relation between the vulnerability of each road section in the bridge network and the set value, and calculating the vulnerability of the bridge network according to the judging result.
In the connectivity analysis of a bridge network, it is generally necessary to use adjacency matrices to represent connection relationships between nodes in the bridge network by means of graph theory. The road sections are connected by the end points of the road sections in a crossing way, and form corresponding paths, and each path is combined to form a bridge network. The nodes are the end points of each road section and are also junction points.
In the theory of diagrams, the representation of a diagram is generallyG= (Y,E),YIs a set of points that are to be selected,Y= {y 1 ,y 2 , …,y p ,…,y q ,…,y m },Eis an edge set of the two-dimensional object,E= {e 1 ,e 2 , …,e n for adjacent matrixA m m× =[a pq ]The representation is made of a combination of a first and a second color,a pq is referred to as formula (12):
(12)
the detailed steps for simulating the vulnerability of the bridge network are as follows:
1. for a certain road section, selecting one road sectionPGV。
2. The vulnerable probability of each road section obtained in the step S3 is calculatedP L As a known value, at [0,1]Internally acquiring for each road segmentRandom valueP W As an observation. Comparison ofP L And (3) withP W Is of a size of (a) and (b). If it isP W Less thanP L Representing the road segment damage, the road segment corresponding in the adjacency matrixa pq Marked as "0", otherwiseP W Not less thanP L Indicating that the road section works normally, the road section corresponding to the adjacent matrixa pq And is designated as "1".
3. After the states of all road sections are determined, the adjacency matrix is determined, and whether paths exist from the starting point to the destination point can be judged by analyzing the adjacency matrix according to the graph theory principle.
Randomly selecting each road section from the set rangeN W Group random numberP W As an observation value, the link vulnerability probability of the link is comparedP L And observed valueP W Is of a size of (a) and (b).
If there is vulnerable probability of road section in the pathP L Greater than the corresponding observed valueP W Then the path does not exist, the variable is countedcount 2 The number of the two-way valve is added with one,count 2 is 0; if not, count variablecount 2 Is unchanged.
StatisticsN W Group number of group random numbers such that no path exists between the start point and the end point of bridge networkn W At this timecount 2 =n W ,count 2 /N W The probability of easy damage of the bridge network is obtained,N W =10000. expressed using equation (13), equation (13) is expressed as follows:
(13)
wherein, P W representing the probability of vulnerability of the bridge network.
4. Modification ofPGVRepeating the steps to obtain the probability of easy damage of the bridge network under each peak ground speed, wherein the corresponding discrete data points are @PGV,P W ) A smooth line is used to connect the discrete data points together to form a vulnerability curve.
Examples:
the example object is a bridge network composed of 5 nodes and 7 road sections, the bridge network is shown in FIG. 5, 20 bridges in step S1 are adopted, and the bridge pier height of each bridgeHAnd basic periodTAs shown in table 2.
Table 2 bridge parameter table
The calculation flow is as follows:
1. according to the basic period of the bridge in Table 2TObtaining the maximum drift ratio of the bridge piers of each bridge by the bridge pier height H and the formula (7)MPDThe results are shown in Table 3:
TABLE 3 expression of maximum drift ratio of pier
According to research, bridge pier maximum drift ratio limit valueM max The values of (2) are shown in Table 4:
TABLE 4 maximum drift ratio limit for piers
2. According to the basic period ratio of the bridgeT r As can be obtained from equation (10), the covariances are combined to form a covariance matrix, the covariance matrix generated by 20 bridges is 20×20, and the covariance matrix is a symmetric matrix, i.e. each column or each row respectively represents data of one bridge. Listing first column data of covariance matrix, e.g. tableShown at 5.
3. For a givenPGVFor example, initiallyPGV=0.03 m/s. Sampling the samples considering the correlation of the seismic displacement response by using a plurality of normal sampling functions in a Python platform, and designatingN L =100000, the sample obtained is oneN L X 20 data matrix, thiColumn data represents the firstiBridge seatN L And response samples.
4. The damaged state, such as the medium damaged state, of a road segment, such as road segment AB, is analyzed. As shown in fig. 5, the section AB in the embodiment includes the bridge 1 and the bridge 2.
When (when)N L When=1, the bridge pier maximum drift ratio of two bridgesMPDIn a data matrix, usingDr[1,1]、Dr[1,2]And (3) representing. The maximum drift ratio of bridge piers of two bridgesMPDLimit value of maximum drift ratio of bridge pier in slight damage stateMPD max Comparison is performed as long as there isMPD>MPD max The road section is considered to be in a slight damage state, and the variable is countedcount 1 Adding one; if not, continuing to judge the second group of data.
When (when)N L When=2, the data matrix is formedDr[2,1]、Dr[2,2]The corresponding data are respectively compared with the limit value of the maximum drift ratio of the bridge pierMPD max Comparison is performed as long as there isMPD>MPD max The road section is considered to be in a slight damage state, and the variable is countedcount 1 Adding one; if not, continuing to judge other groups of data.
When (when)N L When=3, 4 …, the above procedure is repeated untilN L =100000。N L After all of the group samples have been analyzed,count 1 /N L the road section vulnerability probability in the slight damage state of the road section AB is the road section vulnerability probability. Changing limit value of maximum drift ratio of bridge pierMPD max The link vulnerability probability of the link AB in the other three damaged states can be obtained.And the road section vulnerability probability calculation process of other road sections is the same as that of the road section AB.
TABLE 5 covariance matrix portion data
5. Modification ofPGVFor example, the nextPGV=0.06 m/s, which can be obtained by repeating the above stepsPGVRoad segment vulnerability probability of all road segments below. In the present embodimentPGVThe values are from 0.03m/s to 1.50m/s, 50 groups are total, and partial data of the medium probability of occurrence of all road sections are shown in the following table 6.
Table 6 moderate probability of vulnerability
Each is to bePGVThe vulnerable probability of 7 road sections forms a matrixU 1×7 For example:
U 1×7 (PGV=0.30m/s)=[ 0.20278, 0.23943, 0.34123, 0.41875, 0.39404, 0.46120, 0.46119 ];
U 1×7 (PGV=0.60m/s)=[ 0.70286, 0.73475, 0.84276, 0.89045, 0.86740, 0.90796, 0.89630];
U 1×7 (PGV=0.90m/s)=[ 0.90655, 0.92385, 0.96950, 0.97995, 0.97310, 0.98578, 0.98148]。
5. from the graph theory knowledge, the adjacency matrix of the bridge network in the embodiment is written directly, which is a 5×5 matrix:
(14)
for matrix A 5×5 The values of (2) are explained as follows: rows 1 through 5 represent A, C, D, E, F five nodes, and columns 1 through 5 represent these five nodes, respectively; a value of 0 at row 1 and column 3 locations indicates a node in the bridge network of an embodimentNo road segment is directly connected between the point A and the node C, and the value at the position of the 2 nd row and the 3 rd column isa BC Indicating that there is a link between node B and node C when no earthquake has occurred, the state parameter is used here because the damaged state of the link BC after the earthquake is unknowna BC Instead, the value of the state parameter refers to formula (12).
The matrix is determined as followsaIs the value of (1):
to be used forPGVFor example, apply np.random.rand function in Python to generate random numbers between 0 and 1 for 7 road segments, and form a random number matrix:
P W 1×7 (PGV=0.30m/s)=[ 0.03153, 0.45179, 0.90813, 0.03868, 0.78933, 0.39648, 0.26016 ]. Comparing the known matrix using a random number matrix:
U 1×7 (PGV=0.30m/s)=[ 0.20278,0.23943, 0.34123, 0.41875, 0.39404, 0.46120, 0.46119 ]。
judging according to the description in the formula (11) and the step S4aIs the value of (1):
0.03153<0.20278,a AB =0;
0.45179>0.23943,a AD =1;
the remaining results were as follows:
a BC =1;a BE =0;a CD =1;a CE =0;a DE =0. The adjacency matrix in equation (14) is rewritten as:
(15)
after the adjacent matrix is updated, the corresponding bridge network is also updated synchronously, and a bridge network diagram corresponding to the adjacent matrix expressed by the formula (15) is shown in fig. 6.
The network function library in Python can be used to obtain whether there is an undamaged path between two points, for example, in this example, A is used as the starting pointE is the endpoint. Calculating that no path exists between A and E according to the above adjacency matrix, counting variablescount 2 And adding one.
Regenerating random number matrixU 1×7 Repeating the above operation untilN W The analysis of the fractions was complete. 50 groups were analyzedPGVThe process is carried outPGV,P W ) And the scattered points are connected smoothly, so that the vulnerability curve of the bridge network can be obtained. Other lesion status analysis steps are as above. The vulnerability curves of the four types of damage states of the final example are shown in fig. 7 below.
The foregoing is only a preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art, who is within the scope of the present invention, should make equivalent substitutions or modifications according to the technical scheme of the present invention and the inventive concept thereof, and should be covered by the scope of the present invention.
Claims (10)
1. A bridge network vulnerability analysis method considering bridge seismic response correlation is characterized by comprising the following steps:
s1, calculating the maximum drift ratio of bridge piers of each bridge;
s2, calculating the correlation between any two bridges in the same road section through a correlation formula;
s3, establishing association between the correlation of the bridge and the corresponding maximum drift ratio of the bridge piers, and selecting a preset number of maximum drift ratios of the bridge piers through the association; comparing the maximum drift ratio of the bridge piers of each bridge with the limit value of the maximum drift ratio of the bridge piers, and calculating the vulnerability of each road section according to the comparison result;
and S4, judging the magnitude relation between the vulnerability of each road section in the bridge network and the set value, and calculating the vulnerability of the bridge network according to the judging result.
2. The bridge network vulnerability analysis method considering bridge seismic response correlation according to claim 1, wherein the specific steps of step S1 are as follows:
s11, firstly, acquiring a calculation formula of the bridge seismic displacement response, wherein the calculation formula is as follows:wherein, μrepresenting the median value of the seismic displacement response of the bridge;Trepresenting the fundamental period of the bridge;PGVrepresenting peak ground speed;a,b,cthe parameters are parameters of a bridge seismic displacement response calculation formula;
s12, the basic period of each bridge is shortenedT、Peak ground speedPGVCorresponding median value of bridge seismic displacement responseμIntroducing into a pier maximum drift ratio calculation formula, and solving parameters by regression analysisa,b,cThe specific value of (3);
s13, median value of earthquake displacement response through bridgeμCalculating the maximum drift ratio of bridge piersM,MThe calculation formula of (2) is as follows:
wherein, Mrepresenting the maximum drift ratio of the bridge pier;Hthe height of the bridge pier is set;
s14, acquiring the basic period of each bridge in the step S12T、Bridge pier heightHPeak ground speed of the bridge locationPGVAnd leading the values into a bridge pier maximum drift ratio calculation formula to calculate so as to obtain the bridge pier maximum drift ratio of each bridge.
3. The bridge network vulnerability analysis method considering bridge seismic response correlation according to claim 2, wherein the specific steps of step S2 are as follows:
s21, acquiring the basic period of each bridge in the same road sectionTRespectively comparing the basic periods of any two bridges in the road section, and generating a basic period ratio;
s22, calculating covariance between seismic displacement responses of any two bridges in the road section according to a covariance calculation principle;
s23, substituting the calculated basic cycle ratio and the corresponding covariance into a correlation calculation formula, fitting the calculation formula, and solving the specific numerical value of each parameter in the calculation formula through regression analysis;
s24, substituting the calculated basic period ratios in the road section into a correlation calculation formula determined by parameters to deduce the corresponding correlation between the bridges.
4. The bridge network vulnerability analysis method considering bridge seismic response correlation according to claim 3, wherein the specific steps of step S3 are as follows:
s31, selecting a road section, and calculating the ground speed of each bridge in the road section at each peak valuePGVMaximum drift ratio of bridge pierM;
S32, selectingN L Group same peak ground speedPGVMaximum drift ratio of bridge pier of each lower bridgeM;
S33, selecting the maximum drift ratio of each pierMLimit value of maximum drift ratio of bridge pierM max Comparing the sizes;
s34, counting maximum drift ratio of bridge pierMIs larger than the limit value of the maximum drift ratio of the bridge pierM max The number of groups of (a)n L Will beN L Andn L introducing a road section vulnerability probability function representing vulnerability of the road section to obtain the ground speed of the road section at the peak valuePGVThe vulnerable probability of the lower road section;
s35, calculating the ground speed of each road section in the bridge network at each peak value by using the road section vulnerable probability functionPGVThe probability of vulnerability of the lower road section.
5. The bridge network vulnerability analysis method considering bridge seismic response correlation as claimed in claim 4, wherein the specific steps of step S4 are as follows:
s41, intersecting and connecting all road sections through the end points of the road sections, and forming corresponding paths, wherein all paths are combined to form a bridge network;
s42, randomly selecting each road section from the set rangeN W Group random numberP W As an observation value, the link vulnerability probability of the link is comparedP L And observed valueP W Is of a size of (2);
s43, if the vulnerable probability of the road section exists in the pathP L Greater than the corresponding observed valueP W The path does not exist, otherwise it exists;
s44, statisticsN W Group number of group random numbers such that no path exists between the start point and the end point of bridge networkn W And willN W Andn W and (5) leading the probability of the vulnerability of the bridge network into a calculation formula of the probability of the vulnerability of the bridge network to calculate the probability of the vulnerability of the bridge network.
6. The bridge network vulnerability analysis method considering bridge seismic response correlation according to any one of claims 3-5, wherein the basic cycle ratio is calculated as follows:
wherein, R k,i,j represent the firstkIn road segmentsiBridge seat and the firstjThe basic cycle ratio of the seat bridge;T k,i represent the firstkIn road segmentsiThe basic period of the bridge seat;T k,j represent the firstkIn road segmentsjBasic cycle of the bridge.
7. The bridge network vulnerability analysis method considering bridge seismic response correlation as claimed in claim 6, wherein the correlation calculation formula is as follows:
wherein, Cov[i, j]represent the firstiBridge seat and the firstjCovariance between base bridge fundamental periods, the covariance characterizing the firstiBridge seat and the firstjCorrelation between seat bridges;eandfare parameters of the correlation calculation formula.
8. The bridge network vulnerability analysis method considering bridge seismic response correlation according to claim 4 or 5, wherein the calculation formula of the road section vulnerability probability function is as follows:
wherein, P L the link vulnerability probability of the link is represented.
9. The bridge network vulnerability analysis method considering bridge seismic response correlation as claimed in claim 5, wherein the bridge network vulnerability probability calculation formula is as follows:wherein, P W representing the probability of vulnerability of the bridge network.
10. The bridge network vulnerability analysis method considering bridge seismic response correlation as claimed in claim 5, wherein the setting range in step S42 is [0,1].
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310827938.XA CN116561875B (en) | 2023-07-07 | 2023-07-07 | Bridge network vulnerability analysis method considering bridge seismic response correlation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310827938.XA CN116561875B (en) | 2023-07-07 | 2023-07-07 | Bridge network vulnerability analysis method considering bridge seismic response correlation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116561875A true CN116561875A (en) | 2023-08-08 |
CN116561875B CN116561875B (en) | 2023-09-15 |
Family
ID=87488304
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310827938.XA Active CN116561875B (en) | 2023-07-07 | 2023-07-07 | Bridge network vulnerability analysis method considering bridge seismic response correlation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116561875B (en) |
Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010144487A (en) * | 2008-12-22 | 2010-07-01 | Railway Technical Res Inst | Method and system for evaluating level of damage to rc member |
JP2015078554A (en) * | 2013-10-18 | 2015-04-23 | 公益財団法人鉄道総合技術研究所 | Bridge pier soundness evaluation method and bridge pier soundness evaluation device |
US20160217226A1 (en) * | 2015-01-28 | 2016-07-28 | Brian Westcott | Methods and systems for infrastructure performance: monitoring, control, operations, analysis and adaptive learning |
JP2017020172A (en) * | 2015-07-07 | 2017-01-26 | 公益財団法人鉄道総合技術研究所 | Bridge dynamic response evaluation method |
CN106897510A (en) * | 2017-02-16 | 2017-06-27 | 西南交通大学 | A kind of bridge structure 3-D seismics vulnerability analysis method |
JP2018031676A (en) * | 2016-08-24 | 2018-03-01 | 公益財団法人鉄道総合技術研究所 | Method for inspecting structural performance of railroad bridge |
CN109214119A (en) * | 2018-10-16 | 2019-01-15 | 中铁二院工程集团有限责任公司 | Bridge Earthquake Resistance Design method based on response surface model |
CN109918819A (en) * | 2019-03-15 | 2019-06-21 | 哈尔滨工业大学 | A kind of extensive bridge network estimation method based on Bayesian network |
WO2019149271A1 (en) * | 2018-02-05 | 2019-08-08 | 横琴共轭科技有限公司 | Assembled concrete pier system having mixed configuration of frp ribs and regular steel ribs |
JP2019212143A (en) * | 2018-06-07 | 2019-12-12 | 株式会社日立製作所 | Damage prediction method, damage prediction system, and program |
CN110807562A (en) * | 2020-01-07 | 2020-02-18 | 杭州鲁尔物联科技有限公司 | Regional bridge risk prediction method and system |
CN112364426A (en) * | 2020-11-19 | 2021-02-12 | 华东交通大学 | Railway pier damage assessment method and system based on driving safety and dynamic response, terminal device and readable storage medium |
WO2021179350A1 (en) * | 2020-03-13 | 2021-09-16 | Dalian University Of Technology | Method of damage detection for decks of girder bridges using an actively excited vehicle |
CN115130175A (en) * | 2022-06-15 | 2022-09-30 | 大连理工大学 | Bridge structure nonlinear random vibration and earthquake vulnerability analysis method |
US20230106956A1 (en) * | 2021-10-01 | 2023-04-06 | Kinetica Dynamics Inc. | System for building-specific multi-peril risk assessment and mitigation |
DE202023101791U1 (en) * | 2023-04-10 | 2023-04-18 | Banda Sai Chaitanya | A system for RC bridge pier damage assessment using a simplified damage index |
CN116127818A (en) * | 2023-03-20 | 2023-05-16 | 安徽建筑大学 | Bridge seismic vulnerability assessment method based on influence line correction model |
-
2023
- 2023-07-07 CN CN202310827938.XA patent/CN116561875B/en active Active
Patent Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010144487A (en) * | 2008-12-22 | 2010-07-01 | Railway Technical Res Inst | Method and system for evaluating level of damage to rc member |
JP2015078554A (en) * | 2013-10-18 | 2015-04-23 | 公益財団法人鉄道総合技術研究所 | Bridge pier soundness evaluation method and bridge pier soundness evaluation device |
US20160217226A1 (en) * | 2015-01-28 | 2016-07-28 | Brian Westcott | Methods and systems for infrastructure performance: monitoring, control, operations, analysis and adaptive learning |
JP2017020172A (en) * | 2015-07-07 | 2017-01-26 | 公益財団法人鉄道総合技術研究所 | Bridge dynamic response evaluation method |
JP2018031676A (en) * | 2016-08-24 | 2018-03-01 | 公益財団法人鉄道総合技術研究所 | Method for inspecting structural performance of railroad bridge |
CN106897510A (en) * | 2017-02-16 | 2017-06-27 | 西南交通大学 | A kind of bridge structure 3-D seismics vulnerability analysis method |
WO2019149271A1 (en) * | 2018-02-05 | 2019-08-08 | 横琴共轭科技有限公司 | Assembled concrete pier system having mixed configuration of frp ribs and regular steel ribs |
JP2019212143A (en) * | 2018-06-07 | 2019-12-12 | 株式会社日立製作所 | Damage prediction method, damage prediction system, and program |
CN109214119A (en) * | 2018-10-16 | 2019-01-15 | 中铁二院工程集团有限责任公司 | Bridge Earthquake Resistance Design method based on response surface model |
CN109918819A (en) * | 2019-03-15 | 2019-06-21 | 哈尔滨工业大学 | A kind of extensive bridge network estimation method based on Bayesian network |
CN110807562A (en) * | 2020-01-07 | 2020-02-18 | 杭州鲁尔物联科技有限公司 | Regional bridge risk prediction method and system |
WO2021179350A1 (en) * | 2020-03-13 | 2021-09-16 | Dalian University Of Technology | Method of damage detection for decks of girder bridges using an actively excited vehicle |
CN112364426A (en) * | 2020-11-19 | 2021-02-12 | 华东交通大学 | Railway pier damage assessment method and system based on driving safety and dynamic response, terminal device and readable storage medium |
US20230106956A1 (en) * | 2021-10-01 | 2023-04-06 | Kinetica Dynamics Inc. | System for building-specific multi-peril risk assessment and mitigation |
CN115130175A (en) * | 2022-06-15 | 2022-09-30 | 大连理工大学 | Bridge structure nonlinear random vibration and earthquake vulnerability analysis method |
CN116127818A (en) * | 2023-03-20 | 2023-05-16 | 安徽建筑大学 | Bridge seismic vulnerability assessment method based on influence line correction model |
DE202023101791U1 (en) * | 2023-04-10 | 2023-04-18 | Banda Sai Chaitanya | A system for RC bridge pier damage assessment using a simplified damage index |
Non-Patent Citations (6)
Title |
---|
PANG Y等: "Effects of dynamic fluid-structure interaction on seismic response of multi-span deep water bridges using fragility function method", ADVANCES IN STRUCTURAL ENGINEERING * |
UMA S R等: "Development of probabilistic framework for performance-based seismic assessment of structures considering residual deformations", JOURNAL OF EARTHQUAKE ENGINEERING * |
兰日清;李小军;丰彪;王玉石;: "近场强震作用下城市桥梁结构灾害评价方法研究", 自然灾害学报, no. 04 * |
王其昂;吴子燕;贾兆平;: "桥梁***地震多维易损性分析", 工程力学, no. 10 * |
王帅;宋帅;: "桥梁构件易损性曲面的地震动强度参数向量优选", 地震工程与工程振动, no. 04 * |
陈亮;李建中;盛光祖;: "考虑强地面运动持时的概率分布对桥梁结构概率地震需求的影响", 公路交通科技, no. 02 * |
Also Published As
Publication number | Publication date |
---|---|
CN116561875B (en) | 2023-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Di Febbraro et al. | Expert-based and correlative models to map habitat quality: Which gives better support to conservation planning? | |
Boots | Developing local measures of spatial association for categorical data | |
Pontius Jr et al. | GEOMOD modeling | |
Opdam | Assessing the conservation potential of habitat networks | |
Raimbault et al. | Space matters: Extending sensitivity analysis to initial spatial conditions in geosimulation models | |
Makido et al. | Assessing alternatives for modeling the spatial distribution of multiple land-cover classes at sub-pixel scales | |
CN117035465B (en) | Method and device for evaluating landslide susceptibility | |
CN115222196B (en) | Vulnerability assessment and driving force detection method for mountain road network | |
CN115345069A (en) | Lake water volume estimation method based on maximum water depth record and machine learning | |
CN116663960A (en) | Ecological network construction method based on accessibility and availability of habitat under climate change | |
CN112347670A (en) | Rockfill material creep parameter prediction method based on neural network response surface | |
Hagen-Zanker et al. | Map comparison methods for comprehensive assessment of geosimulation models | |
CN116561875B (en) | Bridge network vulnerability analysis method considering bridge seismic response correlation | |
Weaver et al. | Modelling land cover transitions: a solution to the problem of spatial dependence in data | |
Lin et al. | Identifying the spatial mixture distribution of bird diversity across urban and suburban areas in the metropolis: A case study in Taipei Basin of Taiwan | |
Wagner et al. | Realistic simulation of the effects of abundance distribution and spatial heterogeneity on non-parametric estimators of species richness | |
CN116401962A (en) | Method for pushing optimal characteristic scheme of water quality model | |
CN115618611A (en) | Method for establishing spatial multi-scale gradient lifting tree model and electronic equipment | |
CN113689048A (en) | Method, system and computer-readable storage medium for predicting refined spatial distribution of future population | |
Leitao et al. | Modeling overland flow from local inflows in almost no-time, using Self Organizing Maps | |
Pérez et al. | Model testing and assessment: perspectives from a swarm intelligence, agent-based model of forest insect infestations | |
Chen et al. | Disentangling multi‐species aggregate versus overlapping distributions | |
Jardón et al. | Spatial Markov chains implemented in GIS | |
CN116777300B (en) | Urban road traffic accident risk assessment method based on gradient lifting decision tree | |
Lott et al. | Developing and testing TernCOLONY 1.0: an individual-based model of Least Tern reproduction |
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 |