CN109558680B - Bridge multi-target equivalent static wind load calculation method based on POD technology - Google Patents

Bridge multi-target equivalent static wind load calculation method based on POD technology Download PDF

Info

Publication number
CN109558680B
CN109558680B CN201811451886.6A CN201811451886A CN109558680B CN 109558680 B CN109558680 B CN 109558680B CN 201811451886 A CN201811451886 A CN 201811451886A CN 109558680 B CN109558680 B CN 109558680B
Authority
CN
China
Prior art keywords
wind load
equivalent static
static wind
lrc
matrix
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
Application number
CN201811451886.6A
Other languages
Chinese (zh)
Other versions
CN109558680A (en
Inventor
董锐
李狄钦
翁祥颖
韦建刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fuzhou University
Original Assignee
Fuzhou University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fuzhou University filed Critical Fuzhou University
Priority to CN201811451886.6A priority Critical patent/CN109558680B/en
Publication of CN109558680A publication Critical patent/CN109558680A/en
Application granted granted Critical
Publication of CN109558680B publication Critical patent/CN109558680B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

The invention relates to a bridge multi-target equivalent static wind load calculation method based on a POD technology. Firstly, calculating an equivalent static wind load matrix FLRCIntrinsic mode matrix ΦLRCEquivalent static wind load base vector { phi ] based on eigenmode matrixLRC}N×kAnd a combined coefficient column vector C; then, according to the equivalent static wind load basic vector { phi ] obtained by calculationLRC}N×kAnd the column vector C of the combination coefficient is summed to obtain a multi-target equivalent static wind load matrix
Figure RE-DDA0001963508130000011
Finally, for the multi-target equivalent static wind load matrix
Figure RE-DDA0001963508130000012
The calculation accuracy and the distribution rationality of (2) are evaluated. If the computational error of the response of the key part of the bridge structure and the rationality of the equivalent static wind load distribution simultaneously meet the conditions, the solved equivalent static wind load matrix is considered to be
Figure RE-DDA0001963508130000013
And the engineering design requirement is met. The invention comprehensively considers the characteristics of the LRC method, the Universal method and the PSWL method, and realizes the multi-target equivalence of equivalent static wind load of the bridge structure.

Description

Bridge multi-target equivalent static wind load calculation method based on POD technology
Technical Field
The invention belongs to the field of civil engineering, and particularly relates to a bridge multi-target equivalent static wind load calculation method based on a POD (platform-oriented programming) technology.
Background
Due to the complex calculation of the buffeting response of the bridge structure, the design specification of the current mainstream wind load adopts equivalent static wind load to replace the buffeting response. Wind engineering researchers have proposed a plurality of equivalent static wind load calculation methods aiming at the characteristics of different types of structures, wherein the original and representative important methods mainly include a gust wind load factor method (GLF), a load-response correlation method (LRC), an inertia force wind load method (IWL), an IWL + LRC combination method, a uniform equivalent static wind load method (Universal) and a basic wind load method (PSWL).
The GLF method is to take the product of the gust load factor G and the average wind load as the equivalent static wind load of the structural wind vibration response. The GLF method has clear concept and simple and convenient operation, but can only ensure that the wind vibration response extreme values at the equivalent positions of the structure are equal, and the responses at other positions have errors of different degrees. The LRC method is an accurate method for calculating the equivalent static wind load of the structural background, but actually, the extreme values of wind vibration response at equivalent positions can only be guaranteed to be equal, and the response at non-equivalent positions has errors of different degrees. The IWL method is a method for calculating the structural resonance equivalent static wind load, is mainly applied to a high-rise building structure with a control function of a first-order vibration mode, and is difficult to consider the influence of the coupling of a high-order vibration mode and the vibration mode. The IWL + LRC combination method integrates the advantages of the IWL method and the advantages of the LRC method, but the method is difficult to consider the influence of high-order mode shape and mode shape coupling, and is actually a single-target equivalent method, and the response at a non-equivalent position has errors of different degrees. The Universal method is a multi-target equivalent static wind load calculation method, but the selection of the equivalent static wind load basic vector is usually completed by the judgment of an engineer. The PSWL method is a method for constructing an equivalent static wind load envelope value, the method needs to use the premise that all response envelope values are not overestimated, although the obtained equivalent static wind load distribution can guarantee the calculation accuracy of the structure, the rationality of the equivalent static wind load distribution is not evaluated, and the local response distortion of the structure can be caused in some cases. At present, an IWL method considering a first-order vibration mode of a structure is adopted in 'building structure load specification' (GB50009-2012) of China, and a static gust load method is adopted in 'highway bridge wind resistance design specification' (JTG/TD 60-01-2004). The static gust load method adopts a force equivalence principle, and compared with the main method based on structural response equivalence, the buffeting response calculation error is larger.
The invention provides a bridge structure multi-target equivalent static wind load vector method based on an intrinsic orthogonal decomposition (POD) technology on the basis of the existing research, which comprehensively considers the respective characteristics of an LRC method, a Universal method and a PSWL method, realizes the multi-target equivalent of the bridge structure equivalent static wind load, and simultaneously considers the distribution rationality.
Disclosure of Invention
The invention aims to provide a POD (platform POD) technology-based bridge multi-target equivalent static wind load calculation method, which comprehensively considers the characteristics of an LRC (line replaceable channel) method, a Universal method and a PSWL (particle swarm optimization) method, realizes multi-target equivalent of bridge structure equivalent static wind loads, and considers the distribution rationality.
In order to achieve the purpose, the technical scheme of the invention is as follows: a bridge multi-target equivalent static wind load calculation method based on POD technology includes the steps of firstly, calculating an equivalent static wind load matrix FLRCEigenmode matrix phiLRCEquivalent static wind load base vector { phi ] based on eigenmode matrixLRC}N×kAnd a combined coefficient column vector C; then, according to the equivalent static wind load basic vector { phi ] obtained by calculationLRC}N×kAnd the column vector C of the combination coefficient is summed to obtain a multi-target equivalent static wind load matrix
Figure BDA0001886622100000021
Finally, for the multi-target equivalent static wind load matrix
Figure BDA0001886622100000022
The calculation accuracy and the distribution rationality of (2) are evaluated. If the computational error of the response of the key part of the bridge structure and the rationality of the equivalent static wind load distribution simultaneously meet the conditions, the solved equivalent static wind load matrix is considered to be
Figure BDA0001886622100000023
And the engineering design requirements are met. In an embodiment of the present invention, the method specifically includes the following steps:
step S1, calculating an equivalent static wind load matrix FLRCEigenmode matrix phiLRCBased on intrinsic modeEquivalent static wind load base vector { phi of matrixLRC}N×kAnd a combined coefficient column vector C; wherein,
equivalent static wind load matrix FLRC：
According to the relation between load and response, the extreme value is responded by pulsating wind
Figure BDA0001886622100000024
For an equivalent target, obtaining an equivalent static wind load column vector by an LRC method,
Figure BDA0001886622100000025
in the formula,
Figure BDA0001886622100000026
is a response at the ith node of the structure
Figure BDA0001886622100000027
Calculating an equivalent static wind load column vector as an equivalent target according to an LRC method;
Figure BDA0001886622100000028
RMS value column vector of fluctuating wind load;
Figure BDA0001886622100000029
for fluctuating wind loads and responses
Figure BDA00018866221000000210
A load-response correlation coefficient sequence vector between; the symbol "", indicates that the operation between the matrices is multiplication of the corresponding elements;
forming an N x m-dimensional equivalent static wind load matrix F by using the obtained equivalent static wind load column vectorsLRC，
Figure BDA00018866221000000211
Eigenmode matrix phiLRC：
Using intrinsic orthogonal decomposition (POD) technique, called singular value decomposition technique for discrete data structure, equivalent static wind load matrix FLRCThe decomposition is carried out, and the decomposition is carried out,
FLRC＝UΣVT (2)
in the formula, U is an AND matrix FLRCIn an N × N orthogonal matrix with the same number of rows, Σ ═ diag (λ)1…λN) Is a non-negative diagonal matrix of dimension Nxm, and1＞λ2＞…＞λNnot less than 0, V is AND matrix FLRCThe number of columns of the orthogonal matrix is the same; t represents the conjugate transpose of the matrix, and transpose operation is performed on the real matrix;
let phiLRC＝U，Q＝ΣVTThen F isLRCCan be converted into a standard form of POD decomposition to obtain an eigenmode matrix phiLRC；
FLRC＝ΦLRCQ (3)
In the formula, Q is a corresponding main coordinate matrix;
equivalent static wind load base vector { phi ] based on eigenmode matrixLRC}N×k: in the eigenmode matrix phiLRCThe front k-order column vector is used as a basic vector for constructing the multi-target equivalent static wind load;
combined coefficient sequence vector C:
c is a combination coefficient column vector of the equivalent static wind load base vector, and is a k multiplied by 1 order column vector; defining a function f (C) with C as an unknown number,
Figure BDA0001886622100000031
wherein | × | non-conducting phosphor2Expressing a two-norm, and min (#) represents the minimum value;
step S2, obtaining an equivalent static wind load basic vector { phi ] according to calculationLRC}N×kThe bridge buffeting response extreme value is taken as an equivalent target to obtain
Figure BDA0001886622100000032
In the formula,
Figure BDA0001886622100000033
the vector is a buffeting response extreme value column vector (not containing an average wind load response), wherein the sigma r is an RMS value column vector corresponding to the buffeting response r; i is an influence function matrix of the response r;
Figure BDA0001886622100000034
the equivalent static wind load base vector is formed by the first k-order column vectors of the eigenmode matrix; c is a combination coefficient column vector of the equivalent static wind load base vector, and is a k multiplied by 1 order column vector;
then, solving an optimal numerical solution of the combination coefficient column vector C according to a least square method criterion and the formula (4); the multi-target equivalent static wind load of the structure can be obtained after the combination coefficient column vector C is obtained,
Figure BDA0001886622100000035
in the formula,
Figure BDA0001886622100000036
the multi-target equivalent static wind load column vector is obtained;
step S3, obtaining a multi-target equivalent static wind load matrix according to the obtained multiple target equivalent static wind load matrix
Figure BDA0001886622100000037
Calculating the equivalent static wind load precision, and taking the response of the key part of the bridge as a judgment target; the computational error of the response at the ith critical site of the structure is defined as,
Figure BDA0001886622100000041
in the formula, epsiloniCalculating an error for the buffeting response at the ith key location of the structure;
Figure BDA0001886622100000042
calculating a buffeting response extreme value at the ith key point of the structure according to a random buffeting theory, and simply referring to the accurate value;
Figure BDA0001886622100000043
the buffeting response extreme value at the ith key position of the structure is calculated according to the formula (6), and is simply called a calculated value; epsiloncriDetermining according to engineering experience for calculating an error control value; and if the calculation errors of the responses at the key parts of the structure all satisfy the formula (7), the calculation accuracy of the equivalent static wind load is considered to satisfy the requirement.
In an embodiment of the present invention, in step S3, on the premise that the calculation accuracy is ensured, the distribution rationality of the equivalent static wind load needs to be checked:
defining the mean value of equivalent static wind load
Figure BDA0001886622100000044
And standard deviation of
Figure BDA0001886622100000045
Are respectively as
Figure BDA0001886622100000046
Figure BDA0001886622100000047
If the equivalent static wind load distribution at two adjacent points of the structure meets the relation (10), the load at the node is considered to have no sudden change, and the rationality of the equivalent static wind load distribution meets the requirement;
Figure BDA0001886622100000048
in the formula, alpha is an empirical coefficient and is determined according to engineering experience;
and if one of the equivalent static wind load and the equivalent static wind load is not satisfied, the equivalent static wind load needs to be calculated again until the two are satisfied.
In one embodiment of the present invention, the epsiloncriCan be 10 percent, and alpha can be 10 to 20 percent.
In one embodiment of the present invention,. epsilon.criAnd alpha should not be too small at the same time.
Compared with the prior art, the invention has the following beneficial effects:
(1) the basis vector method integrates the respective characteristics of the LRC method, the Universal method and the PSWL method; the method has the advantages that the LRC method can reflect main information of pulsating wind load distribution and structural response distribution, the defect that the Universal method selects the basic vector of the equivalent static wind load through personal judgment of engineers is overcome, and the hypothesis that the rationality of the equivalent static wind load distribution and the calculated value are not overestimated are not considered in the PSWL method is increased;
(2) the computing precision and the rationality of the equivalent static wind load distribution are comprehensively considered, and the wind load in the bridge design is more reasonably guided to be computed.
Drawings
FIG. 1 is a flow chart of the present invention.
Fig. 2 is an overall layout view (unit: m) of the cable-stayed bridge.
FIG. 3 is a section (unit: cm) of a standard girder.
Fig. 4 is a finite element model of a cable-stayed bridge.
Fig. 5 is a wind load acting on the main beam structure.
FIG. 6 is a graph showing buffeting displacement response RMS of a main beam of the cable-stayed bridge.
FIG. 7 is an equivalent static wind load eigenmode matrix ΦLRCSchematic (first 5 stages).
FIG. 8 is a multi-target equivalent static wind load distribution.
FIG. 9 shows a comparison of buffeting shift response extremes.
FIG. 10 is a multi-target equivalent static wind load distribution change rate.
In the figure: FIG. 6 is a plot of buffeting response RMS value (a) and main beam position (b); in fig. 7, (a), (b) and (c) are schematic diagrams of eigenmode matrixes in three directions of vertical, horizontal and torsional rotation, respectively; in fig. 8, (a) is a vertical equivalent static wind load distribution diagram, (b) is a horizontal equivalent static wind load distribution diagram, and (c) is a torsional equivalent static wind load distribution diagram; in fig. 9, (a), (b), and (c) are buffeting displacement response extrema in three directions of vertical, horizontal, and torsional, respectively; in fig. 10, (a), (b) and (c) are respectively the multi-target equivalent static wind load change rate in the vertical direction, the horizontal direction and the torsion direction.
Detailed Description
The technical scheme of the invention is specifically explained below with reference to the accompanying drawings.
As shown in FIG. 1, the method for calculating the multi-target equivalent static wind load of the bridge comprises 4 parameters which are respectively an equivalent static wind load matrix FLRCEigenmode matrix phiLRCEquivalent static wind load base vector { phi ] based on eigenmode matrixLRC}N×kAnd a combined coefficient vector C. Wherein:
equivalent static wind load matrix FLRC: according to the relation between load and response, the extreme value is responded by pulsating wind
Figure BDA0001886622100000051
For an equivalent target, obtaining an equivalent static wind load column vector by an LRC method,
Figure BDA0001886622100000052
in the formula,
Figure BDA0001886622100000053
is a response at the ith node of the structure
Figure BDA0001886622100000054
Calculating an equivalent static wind load column vector as an equivalent target according to an LRC method;
Figure BDA0001886622100000055
RMS value column vector for fluctuating wind load;
Figure BDA0001886622100000056
for fluctuating wind loads and responses
Figure BDA0001886622100000057
A load-response correlation coefficient sequence vector between; the symbol "", indicates that the operation between the matrices is multiplication of the corresponding elements;
forming an N x m-dimensional equivalent static wind load matrix F by using the obtained equivalent static wind load column vectorsLRC，
Figure BDA0001886622100000061
Eigenmode matrix phiLRC: the load matrix F is processed using an intrinsic orthogonal decomposition (POD) technique, referred to as a Singular Value Decomposition (SVD) technique for discrete data structuresLRCThe decomposition is carried out, and the decomposition is carried out,
FLRC＝UΣVT (2)
in the formula, U is an AND matrix FLRCIn an N × N orthogonal matrix with the same number of rows, Σ ═ diag (λ)1…λN) Is a non-negative diagonal matrix of dimension Nxm, and1＞λ2＞…＞λNnot less than 0, V is AND matrix FLRCM × m-dimensional orthogonal matrices having the same number of columns; the symbol "T" represents the conjugate transpose of the matrix, for real matrices the transpose operation;
let phiLRC＝U，Q＝ΣVTThen F isLRCCan be converted into a standard form of POD decomposition to obtain an eigenmode matrix phiLRC；
FLRC＝ΦLRCQ (3)
In the formula, Q is a corresponding main coordinate matrix;
equivalent static wind load base vector { phi ] based on eigenmode matrixLRC}N×k: in the eigenmode matrix phiLRCThe front k-order column vector is used as a basic vector for constructing the multi-target equivalent static wind load.
Combining coefficient vector C: c is a combination coefficient column vector of the equivalent static wind load base vector and is a k multiplied by 1 order column vector. Defining a function f (C) with C as an unknown number,
Figure BDA0001886622100000062
wherein | | circum-calry2The two norms are shown, and min (x) is the minimum.
In the invention, an equivalent static wind load basic vector { phiLRC}N×kThen, the bridge buffeting response extreme value is taken as an equivalent target to obtain
Figure BDA0001886622100000063
In the formula,
Figure BDA0001886622100000064
the vector is a buffeting response extreme value column vector (not containing an average wind load response), wherein the sigma r is an RMS value column vector corresponding to the buffeting response r; i is an influence function matrix of the response r;
Figure BDA0001886622100000065
the equivalent static wind load base vector is formed by the first k-order column vectors of the eigenmode matrix; c is a combination coefficient column vector of the equivalent static wind load base vector, and is a k multiplied by 1 order column vector;
and then, solving the optimal numerical solution of the base vector combination coefficient C according to the formula (4) according to the least square method criterion. The multi-target equivalent static wind load of the structure can be obtained after the combination coefficient C is obtained,
Figure BDA0001886622100000071
in the formula,
Figure BDA0001886622100000072
for multi-target equivalent static windA load column vector.
In the invention, the equivalent static wind load is obtained
Figure BDA0001886622100000073
And then calculating the equivalent static wind load precision, and taking the response of the key part of the bridge as a judgment target. The error in the computation of the response at the ith critical site of the structure is defined as,
Figure BDA0001886622100000074
in the formula, epsiloniCalculating an error for the buffeting response at the ith key position of the structure;
Figure BDA0001886622100000075
calculating a buffeting response extreme value at the ith key point of the structure according to a random buffeting theory, and simply referring to the accurate value;
Figure BDA0001886622100000076
the buffeting response extreme value at the ith key position of the structure is calculated according to the formula (6), and is simply called a calculated value; epsiloncriFor calculating the error control value, the error control value is determined according to engineering experience and can be generally 10%; and if the calculation errors of the responses at the key parts of the structure all satisfy the formula (7), the calculation accuracy of the equivalent static wind load is considered to satisfy the requirement.
In the invention, the distribution rationality of the equivalent static wind load needs to be checked on the premise of ensuring the calculation precision of the buffeting response of the structure. The equivalent static wind load is obtained equivalently according to the overall response of the structure, and is not local response. If the distribution of the equivalent static wind load has violent sudden change, the phenomenon of local damage or instability of the component is likely to occur in the static load combination calculation. In order to avoid the above false results, the rationality of the equivalent static wind load distribution must be ensured. Defining the mean value of equivalent static wind load
Figure BDA0001886622100000077
And standard deviation of
Figure BDA0001886622100000078
Respectively, are as follows,
Figure BDA0001886622100000079
Figure BDA00018866221000000710
if the equivalent static wind load distribution at two adjacent points of the structure meets the relation (10), the load at the node is considered to have no sudden change, and the rationality of the equivalent static wind load distribution meets the requirement.
Figure BDA00018866221000000711
In the formula, alpha is an empirical coefficient, and is determined according to engineering experience, and can be generally 10-20%.
And if one of the calculation accuracy and the rationality of the equivalent static wind load distribution is not satisfied, the equivalent static wind load needs to be calculated again until the calculation accuracy and the rationality of the equivalent static wind load distribution are satisfied. It is to be noted thatcriAnd the value of alpha should not be too small at the same time.
The present invention is described in detail below with reference to the drawings, but the present invention is not limited thereto.
Referring to fig. 1-10, and tables 1-4, tables 1-4 are as follows:
TABLE 1 buffeting calculation of key parameters
Figure BDA0001886622100000081
TABLE 2 Multi-target equivalent static wind load calculation error (unit:%)
Figure BDA0001886622100000082
TABLE 3 Multi-target equivalent static wind load root mean square
Figure BDA0001886622100000083
TABLE 4 maximum value of multiple-target equivalent static wind load change rate
Figure BDA0001886622100000084
Example 1:
the main channel bridge of the east-sea bridge is connected with Shanghai and Yangshan deep-water harbors, and is a single-tower single-cable-plane cable-stayed bridge with a span combination of 73+132+420+132+73 as 830m, and the whole arrangement is shown in figure 2. The main beam is a steel-concrete composite box beam, the beam width is 33m, the central beam height is 4m, and the section of the standard main beam is shown in figure 3; the height of the inverted Y-shaped concrete bridge tower is 148 m. The finite element program used in this example was ANSYS, and a finite element model of a cable-stayed bridge is shown in fig. 4.
And obtaining the main vibration mode and frequency of the front 20 orders of the main girder of the east-sea bridge through modal analysis, and using the main vibration mode and frequency as the basis of buffeting calculation. In practical analysis, the wind speed acting on the main beam mainly takes the average wind speed U, the horizontal pulsating wind speed U (t) in the same direction as the average wind and the vertical pulsating wind speed w (t) into consideration. The wind load acting on the main beam mainly considers the resistance D, the lift L and the lift moment M, and the deformation in the corresponding direction is p, h and alpha respectively. The positive direction of each parameter is shown in fig. 5.
And (3) performing coupled buffeting frequency domain analysis on the main beam of the cable-stayed bridge by adopting a calculation model considering self-exciting force and buffeting force simultaneously on the basis of a finite element model of the cable-stayed bridge. The self-excitation force adopts a self-excitation force expression based on 18 flutter derivatives proposed by Scanlan; the buffeting force is expressed in terms of buffeting force by Davenport considering the pneumatic admittance, which is a simplified expression of Liepmann of the Sears function. The pneumatic force parameters in the self-excitation force and buffeting force expressions are obtained through a main beam segment model wind tunnel test; derivative of flutter
Figure BDA0001886622100000091
And
Figure BDA0001886622100000092
the flutter derivative is obtained through a segment model vibration measurement wind tunnel test
Figure BDA0001886622100000093
And P1 *～P6 *The static force three-component force coefficient and the change rate thereof are obtained by a section model force measurement test, and the data at a wind attack angle of 0 degree are shown in a table 1. Horizontal pulsation wind spectrum S in buffeting calculationuu(n) selecting a Kaimal spectrum; vertical pulsating wind spectrum Sww(n) selecting Lumley-Panofsky to correct the wind spectrum; cross-spectra of horizontal and vertical pulsating wind are according to the literature (Simiu E, Scanlan R H. wind Effects on Structures [ M ] M].JohnWiley&Sons, inc.1996); the spatial correlation takes the form suggested by the Highway bridge wind resistance design Specification (JTG/T D60-01-2004). Other major parameters in the buffeting calculation are shown in table 1, using a 0 ° wind attack angle as an example. Considering the front 20 orders of vibration type of the main beam, the buffeting displacement response RMS value of the main beam is obtained according to the CQC combination mode, as shown in FIG. 6. According to fig. 6, the main span midspan, main span quartet, side span-L2 and side span-R2 midspan of the cable-stayed bridge are respectively taken as key points, as shown in fig. 2. Node1, Node3, Node5 refer to edge-L2, main span, and edge-R2 mid-span locations, respectively, and Node2 and Node4 refer to main-span quarter-point locations.
Determining an equivalent static wind load matrix FLRC. Taking the peak factor g as 3.5 according to
Figure BDA0001886622100000094
Obtaining an equivalent static wind load matrix F corresponding to the main beam buffeting displacement response extreme valueLRC。
Determining the eigenmode matrix phiLRC. Determining an equivalent static wind load matrix FLRCThen obtaining an equivalent static wind load eigenmode matrix phi through SVDLRCThe first 5 th order eigenmode distribution is shown in fig. 7.
Determining equivalent static wind load basis vector [ phi ] based on eigenmode matrixLRC}N×k. Taking the eigenmode matrix phiLRCThe front k-order column vectors form an equivalent static wind load base vector.
According to the method, the multi-target equivalent static wind load when the base vector k is 1, 3, 6, 10 and 20 can be obtained respectively
Figure BDA0001886622100000101
The distribution of (c) is shown in fig. 8. And (3) respectively acting the multi-target equivalent static wind load in the graph 8 on the main girder of the cable-stayed bridge to obtain a buffeting response extreme value calculation value, as shown in the graph 9. The calculation error of the dither response is small when k is 10 and 20, and for clarity, only the extreme values of the accurate dither response and the calculated values when k is 1, 3 and 6 are shown in fig. 9.
According to the method, the multi-target equivalent static wind load accuracy is calculated, and when k is 1, 3, 6, 10 and 20, the buffeting response calculation error at the key position is shown in the table 2. When k is 6, the response error at all key positions is not more than 7%, wherein the calculation error at the main span key point is not more than 5%, and the requirement of engineering design can be met.
According to the method provided by the invention, the rationality of the distribution of the equivalent static wind load needs to be checked on the premise of meeting the calculation precision. As can be seen from fig. 8, as the number of basis vectors participating in calculation increases, the discreteness of the multi-target equivalent static wind load distribution gradually increases.
Figure BDA0001886622100000102
The root mean square of (d) is shown in table 3. To further evaluate the rationality of the equivalent static wind load distribution, the rate of change of the wind load distribution at different positions was calculated according to equation (10), as shown in fig. 10. The maximum values of the change rates of the equivalent static wind load distributions in each case are shown in table 4. As can be seen from FIG. 10, the change rate of the equivalent static wind load distribution at the side span is obviously greater than that of the main span. When k is not more than 6, the distribution change rate of the equivalent static wind load in the vertical direction, the horizontal direction and the torsion direction is not more than 9%, and the requirement of engineering design can be met.
The above are preferred embodiments of the present invention, and all changes made according to the technical solutions of the present invention that produce functional effects do not exceed the scope of the technical solutions of the present invention belong to the protection scope of the present invention.

Claims (2)

1. A bridge multi-target equivalent static wind load calculation method based on POD technology is characterized by comprising the following steps: firstly, calculating an equivalent static wind load matrix FLRCEigenmode matrix phiLRCEquivalent static wind load base vector { phi ] based on eigenmode matrixLRC}N×kAnd a combined coefficient column vector C; then, according to the equivalent static wind load basic vector { phi ] obtained by calculationLRC}N×kAnd combining the column vectors C of the coefficients to obtain a multi-target equivalent static wind load matrix
Figure FDA0003579781150000011
Finally, for the multi-target equivalent static wind load matrix
Figure FDA0003579781150000012
Evaluating the calculation precision and the distribution rationality; if the computational error of the response of the key part of the bridge structure and the rationality of the equivalent static wind load distribution simultaneously meet the conditions, the solved equivalent static wind load matrix is considered to be
Figure FDA0003579781150000013
The engineering design requirements are met;
the method comprises the following concrete implementation steps:
step S1, calculating an equivalent static wind load matrix FLRCEigenmode matrix phiLRCEquivalent static wind load base vector { phi ] based on eigenmode matrixLRC}N×kAnd a combined coefficient column vector C; wherein,
equivalent static wind load matrix FLRC：
According to the relation between load and response, the extreme value is responded by pulsating wind
Figure FDA0003579781150000014
For an equivalent target, obtaining an equivalent static wind load column vector by an LRC method,
Figure FDA0003579781150000015
in the formula,
Figure FDA0003579781150000016
is a response at the ith node of the structure
Figure FDA0003579781150000017
Calculating an equivalent static wind load column vector as an equivalent target according to an LRC method;
Figure FDA0003579781150000018
RMS value column vector for fluctuating wind load;
Figure FDA0003579781150000019
for fluctuating wind loads and responses
Figure FDA00035797811500000110
Load-response related coefficient sequence vectors between; the symbol "", indicates that the operation between the matrices is multiplication of the corresponding elements;
forming an N x m-dimensional equivalent static wind load matrix F by using the obtained equivalent static wind load column vectorsLRC，
Figure FDA00035797811500000111
Eigenmode matrix phiLRC：
Adopting POD technology, for discrete data structure called singular value decomposition technology, equivalent static wind load matrix FLRCThe decomposition is carried out, and the decomposition is carried out,
FLRC＝UΣVT (2)
in which U is an AND matrix FLRCIn the same row number, i.e., diag (λ)1…λN) Is a non-negative diagonal matrix of dimension Nxm, and1＞λ2＞…＞λNnot less than 0, V is AND matrix FLRCM × m-dimensional orthogonal matrices having the same number of columns; t represents the conjugate transpose of the matrix, and transpose operation is performed on the real matrix;
let phi toLRC＝U，Q＝ΣVTThen F isLRCConverting into standard form of POD decomposition to obtain eigenmode matrix phiLRC；
FLRC＝ΦLRCQ (3)
In the formula, Q is a corresponding main coordinate matrix;
equivalent static wind load base vector { phi ] based on eigenmode matrixLRC}N×k: in the eigenmode matrix phiLRCThe front k-order column vector is used as a basic vector for constructing the multi-target equivalent static wind load;
combined coefficient sequence vector C:
c is a combination coefficient column vector of the equivalent static wind load base vector, and is a k multiplied by 1 order column vector; defining a function f (C) with C as an unknown number,
Figure FDA0003579781150000021
wherein | × | non-conducting phosphor2Expressing a two-norm, and min (#) represents the minimum value;
step S2, obtaining an equivalent static wind load basic vector { phi ] according to calculationLRC}N×kThe bridge buffeting response extreme value is taken as an equivalent target to obtain
Figure FDA0003579781150000022
In the formula,
Figure FDA0003579781150000023
is an extreme column vector of buffeting response, containing no mean wind load response, where σrA column vector of RMS values corresponding to dither response r; i is an influence function matrix of the response r;
Figure FDA0003579781150000024
the equivalent static wind load base vector is formed by the first k-order column vectors of the eigenmode matrix; c is a combination coefficient column vector of the equivalent static wind load base vector, and is a k multiplied by 1 order column vector;
then, solving an optimal numerical solution of the combination coefficient column vector C according to a least square method criterion and the formula (4); the multi-target equivalent static wind load of the structure is obtained after the combination series vector C is obtained,
Figure FDA0003579781150000025
in the formula,
Figure FDA0003579781150000026
the multi-target equivalent static wind load column vector is obtained;
step S3, obtaining a multi-target equivalent static wind load matrix according to the obtained multiple target equivalent static wind load matrix
Figure FDA0003579781150000027
Calculating the equivalent static wind load precision, and taking the response of the key part of the bridge as a judgment target; the computational error of the response at the ith critical site of the structure is defined as,
Figure FDA0003579781150000028
in the formula, epsiloniCalculating an error for the buffeting response at the ith key position of the structure;
Figure FDA0003579781150000029
to pressCalculating a buffeting response extreme value at the ith key point of the structure according to a random buffeting theory;
Figure FDA0003579781150000031
the buffeting response extreme value at the ith key position of the structure is calculated according to the formula (6); epsiloncriDetermining an error control value according to engineering experience for calculating the error control value; if the calculation errors of the responses at the key parts of the structure all satisfy the formula (7), the calculation accuracy of the equivalent static wind load is considered to satisfy the requirement;
in step S3, the distribution rationality of the equivalent static wind load must be checked on the premise of ensuring the calculation accuracy:
defining the mean value of equivalent static wind load
Figure FDA0003579781150000032
And standard deviation of
Figure FDA0003579781150000033
Are respectively as
Figure FDA0003579781150000034
Figure FDA0003579781150000035
If the equivalent static wind load distribution at two adjacent points of the structure meets the relation (10), the load at the node is considered to have no sudden change, and the rationality of the equivalent static wind load distribution meets the requirement;
Figure FDA0003579781150000036
in the formula, alpha is an empirical coefficient and is determined according to engineering experience;
and if one of the equivalent static wind load and the equivalent static wind load is not satisfied, the equivalent static wind load needs to be calculated again until the two are satisfied.
2. The POD technology-based bridge multi-target equivalent static wind load calculation method according to claim 1, wherein the epsiloncriTaking 10 percent and taking 10 to 20 percent of alpha.
CN201811451886.6A 2018-11-30 2018-11-30 Bridge multi-target equivalent static wind load calculation method based on POD technology Active CN109558680B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811451886.6A CN109558680B (en) 2018-11-30 2018-11-30 Bridge multi-target equivalent static wind load calculation method based on POD technology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811451886.6A CN109558680B (en) 2018-11-30 2018-11-30 Bridge multi-target equivalent static wind load calculation method based on POD technology

Publications (2)

Publication Number Publication Date
CN109558680A CN109558680A (en) 2019-04-02
CN109558680B true CN109558680B (en) 2022-06-14

Family

ID=65868149

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811451886.6A Active CN109558680B (en) 2018-11-30 2018-11-30 Bridge multi-target equivalent static wind load calculation method based on POD technology

Country Status (1)

Country Link
CN (1) CN109558680B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110083798A (en) * 2019-05-16 2019-08-02 中国十九冶集团有限公司 Combined analysis method for three-direction wind load components of high-rise building
CN114459726A (en) * 2021-12-29 2022-05-10 北京市建筑设计研究院有限公司 Building wind pressure measuring method utilizing wind tunnel test and random field decomposition technology
CN116244793B (en) * 2022-12-30 2023-09-22 重庆科技学院 Method for calculating design wind load of bent torsion column spiral Liang Jingguan tower based on GLF method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1595163A (en) * 2004-06-25 2005-03-16 武汉理工大学 Intelligent health observation system for roof network frame construction with large-scale complex body type
KR20100026133A (en) * 2008-08-29 2010-03-10 한양대학교 산학협력단 Method for topology optimization design using equivalent static loads
CN103020471A (en) * 2012-12-27 2013-04-03 黑龙江大学 Block Ritz vector generation method for fluctuating wind-induced response calculation of long-span roof structure

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1595163A (en) * 2004-06-25 2005-03-16 武汉理工大学 Intelligent health observation system for roof network frame construction with large-scale complex body type
KR20100026133A (en) * 2008-08-29 2010-03-10 한양대학교 산학협력단 Method for topology optimization design using equivalent static loads
CN103020471A (en) * 2012-12-27 2013-04-03 黑龙江大学 Block Ritz vector generation method for fluctuating wind-induced response calculation of long-span roof structure

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于修正LRC的响应分组法求解等效静力风荷载;周晅毅等;《万方数据期刊库》;20100612;第1-9页 *

Also Published As

Publication number Publication date
CN109558680A (en) 2019-04-02

Similar Documents

Publication Publication Date Title
CN109558680B (en) Bridge multi-target equivalent static wind load calculation method based on POD technology
CN105843073B (en) A kind of wing structure aeroelastic stability analysis method not knowing depression of order based on aerodynamic force
CN107577870B (en) Power distribution network voltage power sensitivity robust estimation method based on synchronous phasor measurement
CN104484502A (en) Finite element model correction method based on positive substructure
CN111353910A (en) Distributed intelligent power grid economic dispatching method based on finite time consistency under directed topology
CN111709085A (en) Topological optimization design method for constrained damping sheet structure
CN107168281A (en) Multi-agent system method for diagnosing faults based on finite time observer
KR102578263B1 (en) Multi-agent system consistency distributed control method for agriculture based on sampling data
CN110504676A (en) A kind of State Estimation for Distribution Network based on APSO-BP
CN111988786B (en) Sensor network covering method and system based on high-dimensional multi-target decomposition algorithm
CN114925462B (en) Thin-wall part machining deformation prediction method based on cutting force and rigidity correlation evolution
CN115455793A (en) High-rise structure complex component stress analysis method based on multi-scale model correction
CN105391056A (en) Power system distributed economic dispatching method taking unbalanced communication network into consideration
CN111539138A (en) Method for solving time domain response sensitivity of structural dynamics peak value based on step function
CN109657301B (en) Structural topology optimization method containing pathological load based on double-aggregation function
CN102682175A (en) Method for analyzing reliability of construction error of grid structure based on buckling mode combination
CN107069708B (en) Extreme learning machine-based transmission network line active safety correction method
CN113328467B (en) Probability voltage stability evaluation method, system, terminal device and medium
CN111008686A (en) Grid structure optimization method for realizing reliable splitting based on node correlation
CN104965506A (en) Adjustable parameter-based distributed flight control system real-time fault diagnosis method
CN113298376B (en) Given time consistency control method for economic dispatching of smart power grid with valve point effect
CN114609909A (en) Design method of random multi-intelligence system control protocol under switching topology
CN106599541B (en) A kind of structure and parameter on-line identification method of dynamic power load model
CN114221334B (en) Rapid state estimation method based on graph neural network
CN113110321B (en) Distributed estimation method of networked industrial control system based on event trigger

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