Summary of the invention
Purpose of the present invention is to provide a kind of ultra-large type cooling tower construction phase wind load obtaining value method.
Purpose of the present invention can reach by taking following technical scheme:
A kind of ultra-large type cooling tower construction phase wind load obtaining value method is characterized in that:
1) take Reliability Theory as theoretical foundation, utilize the period analytical approach of structural reliability set up respectively the structural design base period and therebetween the period relation of day part and construction time and during the fiduciary level relation of each time period, form day part relational expression and fiduciary level relational expression;
2) solution procedure 1) described relational expression, obtain construction time characteristi cvalue o fwindload and design reference period characteristi cvalue o fwindload ratio;
3) according to step 2) described ratio finds the solution the construction time wind load factor.
The present invention is take the basic theories of structural reliability design as the basis, has identical fiduciary level as prerequisite in cooling tower tube construction phase and design reference period, and the value typical problem of tower tube construction phase designed wind load is studied.According to the characteristics of the steady binomial random process model of wind load, utilize the period analytical approach of structural reliability to set up respectively structural design base period and the therebetween relation of day part and construction time and therebetween the fiduciary level relation of day part.Start with from the structure reactance of construction time and design reference period is proportional, requirement according to partial safety factor method for designing in the standard, obtain the same proportional conclusion of construction time characteristi cvalue o fwindload and design reference period characteristi cvalue o fwindload, and obtain accordingly the construction time wind load factor.
Purpose of the present invention can also reach by taking following technical scheme:
1) the structural design base period and therebetween the relation between the fiduciary level on the day part represent with following expression formula:
The form that the power function of structure all adopts structure reactance to compare with load effect, establish the power function of structure in arbitrary τ and be:
Zi=Ri-Si=R-Si,i=1,2,...,n, (1)
Wherein R is the structure reactance function, and Si is the load effect function;
The failure probability pfe of each period and reliability index β e are:
p
fe=Φ(-β
e)=P(Z
i≤0) (2)
Structural failure Probability p f and reliability index β are:
Perhaps be:
p
f=Φ(-β)≈1-Φ
n(β
e,ρ) (4)
In the expression formula (4), Φ n represents n dimension Standard Normal Distribution; β e is that element all is the n dimension reliability index vector of β e; ρ is n rank related coefficient square formations, and its nondiagonal element is ρ e.ρ e is Zi and the Zj (related coefficient of i ≠ j).
Further, if the equal Normal Distribution of R and S, the average μ R of R, standard deviation is σ R; The average of Se is μ S, and standard deviation is σ S, then
A perhaps β e through type mistake! Do not find Reference source.With the formula mistake! Do not find Reference source.Obtain.
Further, consider that R and Si or Sj are separate, (related coefficient of i ≠ j) is for Zi and Zj
Further, if R and S are Non-normal Variable, the cumulative distribution function of R is FR (R), and probability density function is fR (R); The cumulative distribution function of Se is FS (S), and probability density function is fS (S); Utilize the JC method to carry out the equivalent normalize of variable, establish with the corresponding equivalent normalize of R variable is R ' for this reason, and its average is μ
R ', standard deviation is σ
R 'With the corresponding equivalent normalize of S variable be S ', its average is μ
S ', standard deviation is σ
S 'An if power function mistake! Do not find Reference source.Designcheck point be (r
*, s
*), according to the equivalent normalize condition of JC method,
μ
R′=r
*-Φ
-1[FR(r
*)]σ
R′,
μ
S′=r
*-Φ
-1[F
S(s
*)]σ
S′,
In the formula,
Expression standard normal probability density function; So, have
2) the structure construction phase and during relation between the fiduciary level on the period represented by following expression formula:
If Tc inner structure reliability index is β c, failure probability is pf, c; τ c inner structure reliability index is β e, and c, failure probability are pfe, c, and corresponding structure reactance is Rc, and following expression is namely arranged:
Z
i,c=R
c-S
i=R
c-S
e=Z
e,
c, i=1,2,K,n
c (7)
p
fe,c=Φ(-β
e,c)=P(Z
i,c≤0)=P(Z
e,c≤0) (8)
β e wherein, c is that element all is β e, the n dimension reliability index vector of c; ρ c is n rank related coefficient square formations, and its nondiagonal element is ρ e, c.
Further, for construction time Tc and the wherein processing of the relation between the fiduciary level of period τ c, adopt structural design base period and the relational expression mistake between the fiduciary level on the day part therebetween! Do not find Reference source.~formula mistake! Do not find Reference source.If, Rc and S Normal Distribution, the average of Rc is
Standard deviation is
Then
If Rc and S are Non-normal Variables, the cumulative distribution function of Rc is
Probability density function is
Can utilize the JC method to carry out the equivalent normalize of variable, establish with the corresponding equivalent normalize of Rc variable is R ' for this reason
c, its average is
Standard deviation is
If power function Zi, the designcheck point of c=Rc-Se are (rc
*, s
*), the equivalent normalize condition according to the JC method obtains expression formula
The equivalent normalize condition of S is pressed expression formula (8), so, obtain expression formula
3) relation between construction time and the design reference period characteristic value of variable action is represented by following expression formula:
Construction time structure reactance Rc is the function of design reference period structure reactance R, and expression formula is
R
c=kR (15)
Wherein k is scale-up factor, is the construction time structure reactance factor.
Further, if the equal Normal Distribution of R and S, with formula (15) substitution formula (14) and formula (15),
If R and S are Non-normal Variables, will through corresponding parameter in the parameter substitution formula (16) after the equivalent normalize of JC method and the formula (17), get
Further, the design load of establishing structure reactance R and action effect S is respectively Rd, Sd, and standard value is respectively Rk, Sk, and partial safety factor is respectively γ R, γ S, requires γ 0Sd≤Rd during the ultimate limit states design, when being in ultimate limit state is
γ
0S
d=R
d,γ0γ
SS
k=R
k/γ
R (20)
Wherein
γ
R=R
k/R
d=R
k/r
*=(μ
R-α
Rσ
R)/r
* (21)
γ
S=S
d/S
k=s
*/S
k=s
*/(μ
S+α
Sσ
S) (22)
The desirable α R=of fraction during fraction coefficient 95% α S=1.645.
Further, the standard value of establishing required drag Rc of structure construction phase is Rck, and the standard value of action effect is Sck, presses equally the ultimate limit states design, namely requires γ 0Scd≤Rcd, when being in ultimate limit state is
γ
0S
cd=R
cd,γ
0γ
SS
ck=R
ck/γ
R (23)
By formula (15), formula (20) and formula (23),
S
ck=kS
k (24)
According to the structural design standard, the design load Sd of combination of load effect gets the least favorable value among the following combination and determines:
Combination by the variable load effect control
Combination by the permanent load effect control
In the formula, γ G is the dead load partial safety factor;
With
Be the 1st and i variable load partial safety factor; SGk is the characteristic value of load effect;
Be the 1st variable load standard value effect, this effect is greater than other any one variable load standard value effects;
Be i variable load standard value effect; ψ ci is the combined value coefficient of i variable load.
Further, the situation for the construction time has following expression
Consider first separately certain just like
Or
The characteristic value of variable action effect, and with Q1 as wind load W, make that other characteristic value of action effects are 0 in formula (25)~formula (28), recycling formula (24) can get
When structure was carried out theoretical analysis, action effect S and effect Q generally adopted linear relationship, i.e. S=CQ, and C is coefficient.If the front then has thinking that the coefficient of action effect of different times is constant at this moment
W
ck=kW
k (30)
4) construction time variable action factor k's finds the solution
Find the solution construction time variable action factor k according to following steps:
(1) gets initial value n=50, β=β c=2.7~4.2, γ S=1.4, γ R=1.1~1.4, μ S=1.Given a certain nc.
(2) initial value of supposition μ R, σ R and σ S, for example 5,0.5,0.5.
(3) separate the formula mistake! Do not find Reference source.With the formula mistake! Do not find Reference source., get β e, σ
R ', σ
S ', r
*And s
*Utilize the designcheck point method, wherein μ
R 'And σ
R 'Calculate μ with formula (7)
S 'And σ
S 'Calculate with formula (8).
(4) calculate ρ e, utilize formula (10).
(5) form system of equations, utilize the formula mistake! Do not find Reference source., formula (21) and formula (22).Solution μ R, σ R and σ S.
(6) calculate k, utilize formula (13), formula (18) and formula (19).
When not being the integral multiple in 1 year for construction time Tc, k calculates with following way:
(1) calculate first Tc integer year correspondence k, meter is made k1.
(2) regard remaining fate of Tc as new nc, T got into 365 days, and γ R is identical when calculating k1, calculates k, is denoted as k2.
(3)k=k1 k2。
About the formula mistake! Do not find Reference source.The calculating of middle n dimension normal distyribution function Φ n (β e, ρ), when n=1, Φ n (β e, ρ)=Φ (β e); When n=2, can utilize following formula to calculate: [3]
Wherein
When n>3, can utilize following formula: [3]
When R and S Normal Distribution, probability density function and the cumulative distribution function of normal random variable X are respectively
When the R obeys logarithm normal distribution, it is that Gumbel distributes that S obeys the distribution of extreme value I type, and probability density function and the cumulative distribution function of lognormal variable R are respectively
r>0 (36)
r>0 (37)
Parameter wherein
Probability density function and the cumulative distribution function of extreme value I type variable S are respectively
f
s(s)=αexp{-α(s-u)-exp[-α(s-u)]} (39)
F
S(s)=exp{-exp[-α(s-u)]} (40)
Wherein
The present invention has following outstanding beneficial effect:
1, the present invention relates to the numerical analysis method of complete ultra-large type cooling tower construction phase wind load value, comprise finding the solution of relation between fiduciary level relation, construction time and the design reference period characteristic value of variable action of cooling tower structure construction time and day part thereof and the wind load factor.The wind load value of construction time and cooling tower are built up final operating mode distinguish, can reduce construction investment, particularly for the ultra-large type cooling tower, determine that the value of construction time wind load can be reducingd the construction costs greatly.In addition, can according to method of the present invention, work up the designed wind load standard of reinforced concrete hyperbolic curve cooling tower tube construction phase at the relevant specification of construction work industry.
2, the present invention is take the basic theories of structural reliability design as the basis, in cooling tower tube construction phase and design reference period, has identical fiduciary level as prerequisite, value typical problem to tower tube construction phase designed wind load is studied, be conducive to determining of cooling tower construction phase designed wind load, stability and the cost of construction time cooling tower had decisive meaning.
Embodiment
The ultra-large type cooling tower construction phase wind load obtaining value method that present embodiment relates to comprises the steps:
1) sets up the steady binomial random process model of wind load;
2) according to the characteristics of the steady binomial random process model of wind load, utilize the period analytical approach of structural reliability to set up respectively structural design base period and the relation of day part therebetween, reach the fiduciary level relation of day part therebetween with the construction time, form day part relational expression and fiduciary level relational expression;
3) start with according to the structure reactance of construction time and design reference period is proportional, by solution procedure 2) in relational expression, obtain construction time characteristi cvalue o fwindload and design reference period characteristi cvalue o fwindload ratio;
4) according to step 3) obtain construction time characteristi cvalue o fwindload and the same proportional conclusion of design reference period characteristi cvalue o fwindload, obtain accordingly the construction time wind load factor.
In the present embodiment:
1) the structural design base period and therebetween the relation between the fiduciary level on the day part represent with following expression formula:
According to action effect, every change direct action time length structurally, design reference period T is divided into n period τ that equates, or think that action effect evenly changes n=T/ τ time in design reference period T, in each period τ, action effect S is stochastic variable, and upper Si of different period and Sj (probability distribution of i ≠ j) is identical, action effect amplitude stochastic variable on the different period τ is separate, and with whether to occur load on the period irrelevant;
The form that the power function of structure all adopts structure reactance to compare with load effect is Zi=Ri-Si=R-Si so can establish the power function of structure in arbitrary τ, i=1,2 ..., n, wherein R is the structure reactance function, and Si is the load effect function, and they are separate, when being linear relationship between load and the load effect, upper load effect Si of different periods and Sj, wherein i ≠ j, independent identically distributed, equivalent stochastic variable Se on the available period τ represents, Zi=R-Se, therefore
Z
i=R-S
i=R-S
e,i=1,2,K,n (1)
The failure probability pfe of each period and reliability index β e are
p
fe=Φ(-β
e)=P(Z
i≤0) (2)
Wherein Φ represents Standard Normal Distribution;
Like this, be exactly the failure event of n period series connection system in the failure event of T inner structure, i.e. structural failure Probability p f and reliability index β are
The computing method of cascade system failure probability have bound method and direct calculation method [3], and the latter's calculating formula is
p
f=Φ(-β)≈1-Φ
n(β
e,ρ) (4)
Wherein Φ n represents n dimension Standard Normal Distribution; β e is that element all is the n dimension reliability index vector of β e; ρ is n rank related coefficient square formations, and its nondiagonal element is ρ e.ρ e is Zi and the Zj (related coefficient of i ≠ j).
If the equal Normal Distribution of R and S, the average μ R of R, standard deviation is σ R; The average of Se is μ S, and standard deviation is σ S, then
β e also could unify the through type mistake! Do not find Reference source.With the formula mistake! Do not find Reference source.Obtain.Consider that again R and Si or Sj are separate, (related coefficient of i ≠ j) is for Zi and Zj
If R and S are Non-normal Variables, the cumulative distribution function of R is FR (R), and probability density function is fR (R); The cumulative distribution function of Se is FS (S), and probability density function is fS (S).Can utilize the JC method to carry out the equivalent normalize of variable.Establish with the corresponding equivalent normalize of R variable is R ' for this reason, and its average is μ
R ', standard deviation is σ
R 'With the corresponding equivalent normalize of S variable be S ', its average is μ
S ', standard deviation is σ
S 'An if power function mistake! Do not find Reference source.Designcheck point be (r
*, s
*), according to the equivalent normalize condition [3] of JC method,
μ
R′=r
*-Φ
-1[F
R(r
*)]σ
R′,
μ
S′=r
*-Φ
-1[F
S(s
*)]σ
S′,
In the formula,
Expression standard normal probability density function.So, have
Known that by formula (44) and formula (45) compare formula (42) and formula (43), formula (46) and formula (47) have more generality.
2) the structure construction phase and during relation between the fiduciary level on the period represented by following expression formula:
Suppose in the effect of structure construction phase and effect thereof identical in the characteristic of construction time and design reference period.The action effect of supposing again structure can be processed into steady binomial stochastic process, i.e. supposition: (1) is according to the every change direct action of action effect time length structurally, design reference period Tc is divided into nc period τ that equates, or thinks that the interior action effect of Tc evenly changes nc=Tc/ τ c time.(2) in each period τ c, action effect S is stochastic variable, and upper Si of different period and Sj (probability distribution of i ≠ j) is identical.(3) the action effect amplitude stochastic variable on the different period τ c is separate, and with whether to occur load on the period irrelevant.
Design reference period is the time parameter of selecting for the value of determining variable action etc.Shorter when the construction time, as less than 1 year, it is just improper for the variable action that has such as wind load (year maximum wind pressure) then Tc to be analogous to T.Therefore, the effect of structure construction phase and effect thereof being used as steady binomial stochastic process processes approximation is arranged.
If Tc inner structure reliability index is β c, failure probability is pf, c; τ c inner structure reliability index is β e, and c, failure probability are pfe, c, and corresponding structure reactance is Rc, and following expression is namely arranged:
Z
i,c=R
c-S
i=R
c-S
e=Z
e,c, i=1,2,K,n
c (48)
p
fe,
c=Φ(-β
e,c)=P(Z
i,c≤0)=P(Z
e,c≤0) (49)
β e wherein, c is that element all is β e, the n dimension reliability index vector of c; ρ c is n rank related coefficient square formations, and its nondiagonal element is ρ e, c.
It may be noted that for construction time Tc and the wherein processing of the relation between the fiduciary level of period τ c, be equivalent to the formula mistake! Do not find Reference source.~formula mistake! Do not find Reference source.If Rc and the equal Normal Distribution of S, the average of Rc is
Standard deviation is
Then
If Rc and S are Non-normal Variables, the cumulative distribution function of Rc is
Probability density function is
Can utilize the JC method to carry out the equivalent normalize of variable.Establish with the corresponding equivalent normalize of Rc variable is R ' for this reason
c, its average is
Standard deviation is
If power function Zi, the designcheck point of c=Rc-Se are (rc
*, s
*), according to the equivalent normalize condition [3] of JC method,
The equivalent normalize conditional (45) of S stands good.So, have
3) relation between construction time and the design reference period characteristic value of variable action is represented by following expression formula:
Construction time structure reactance Rc is the function of design reference period structure reactance R, and this function always can come match with polynomial expression, and when R=0 Rc=0 should be arranged, thus in the polynomial expression without constant term.Wherein simple and practical with proportionate relationship, namely
R
c=kR (56)
Wherein k is scale-up factor, can be called the construction time structure reactance factor.
If the equal Normal Distribution of R and S, with formula (15) substitution formula (51) and formula (52),
If R and S are Non-normal Variables, will through corresponding parameter in the parameter substitution formula (16) after the equivalent normalize of JC method and the formula (17), get
If the design load of structure reactance R and action effect S is respectively Rd, Sd, standard value is respectively Rk, Sk, and partial safety factor is respectively γ R, γ S, requires γ 0Sd≤Rd[4 during the ultimate limit states design], when just being in ultimate limit state be
γ
0S
d=R
d,γ
0γ
SS
k=R
k/γ
R (61)
Wherein
γ
R=R
k/R
d=R
k/r
*=(μ
R-α
Rσ
R)/r
* (62)
γ
S=S
d/S
k=s
*/S
k=s
*/(μ
S+α
Sσ
S) (63)
And the desirable α R=of the fraction α S=1.645 during fraction coefficient 95%.
If the standard value of required drag Rc of structure construction phase is Rck, the standard value of action effect is Sck, presses equally the ultimate limit states design, namely requires γ 0Scd≤Rcd, when just being in ultimate limit state is
γ
0S
cd=R
cd,γ
0γ
SS
ck=
Rck/γ
R (64)
By formula (15), formula (20) and formula (23),
S
ck=kS
k (65)
Structural design standard [4] regulation, the design load Sd of combination of load effect should get the least favorable value among following combination
Determine:
(1) by the combination of variable load effect control
(2) by the combination of permanent load effect control
In the formula, γ G is the dead load partial safety factor;
With
Be the 1st and i variable load partial safety factor; SGk is the characteristic value of load effect;
Be the 1st variable load standard value effect, this effect is greater than other any one variable load standard value effects;
Be i variable load standard value effect; ψ ci is the combined value coefficient of i variable load.
Situation for the construction time has
Consider first separately certain just like
Or
The characteristic value of variable action effect, and with Q1 as wind load W, make that other characteristic value of action effects are 0 in formula (25)~formula (28), recycling formula (24) can get
When structure was carried out theoretical analysis, action effect S and effect Q generally adopted linear relationship, i.e. S=CQ, and C is coefficient.If the front then has thinking that the coefficient of action effect of different times is constant at this moment
W
ck=kW
k (71)
Therefore, by formula (15) and formula (24) as can be known, according to the designing requirement of existing standard, structure should be identical at the ratio of the ratio of the drag standard value of construction time and design reference period and action effect value.Again according to formula (30), when effect and effect thereof were linear, also the ratio with variable action such as characteristi cvalue o fwindload was identical.This ratio k is hereinafter referred to as the construction time variable action factor.
4) construction time variable action factor k's finds the solution
According to above analysis, can find the solution k according to following steps:
(1) gets initial value n=50, β=β c=2.7~4.2, γ S=1.4, γ R=1.1~1.4, μ S=1.Given a certain nc.
(2) initial value of supposition μ R, σ R and σ S, for example 5,0.5,0.5.
(3) separate the formula mistake! Do not find Reference source.With the formula mistake! Do not find Reference source., get β e, σ
R ', σ
S ', r
*And s
*Utilize the designcheck point method, wherein μ
R 'And σ
R 'Calculate μ with formula (44)
S 'And σ
S 'Calculate with formula (45).
(4) calculate ρ e, utilize formula (47).
(5) form system of equations, utilize the formula mistake! Do not find Reference source., formula (21) and formula (22).Solution μ R, σ R and σ S.
(6) calculate k, utilize formula (50), formula (18) and formula (19).
When not being the integral multiple in 1 year for construction time Tc, k can calculate with following way:
(1) calculate first Tc integer year correspondence k, meter is made k1.
(2) regard remaining fate of Tc as new nc, T got into 365 days, and γ R is identical when calculating k1, calculates k, is denoted as k2.
(3)k=k1 k2。
About the formula mistake! Do not find Reference source.The calculating of middle n dimension normal distyribution function Φ n (β e, ρ), when n=1, Φ n (β e, ρ)=Φ (β e); When n=2, can utilize following formula to calculate: [3]
Wherein
When n>3, can utilize following formula: [3]
The simplest situation is the equal Normal Distribution of R and S.Probability density function and the cumulative distribution function of normal random variable X are respectively
More realistic is R obeys logarithm normal distribution [1], and S obeys extreme value I type distribution (Gumbel distribution) [4].Probability density function and the cumulative distribution function of lognormal variable R are respectively
r>0 (77)
r>0 (78)
Parameter wherein
Probability density function and the cumulative distribution function of extreme value I type variable S are respectively
f
s(s)=αexp{-α(s-u)-exp[-α(s-u)]} (80)
F
s(s)=exp{-exp[-α(s-u)]} (81)
Wherein
This paper is with the platform of MATLAB software as the above algorithm of realization.Form and solving equation group function:
5) result of calculation of construction time variable action factor k
According to " engineering structure reliability design unified standard ", design reference period T is 50 years, and 2 grades of the reliability index beta structure safe classes of structural-load-carrying capacity ultimate limit state design, destruction type are that ductile fracture gets 3.2.According to " industrial circulating water Cooling Design standard ", " loading code for design of building structures ", generally get the partial safety factor γ S=1.40 of load and effect thereof.Partial safety factor for resistance γ R is regulation not, and according to Related domestic documents to the analysis of ordinary reinforced concrete structure under the stresses such as uniaxial force, eccentric force, flexure, Shear, 1.10≤γ R≤1.40.1.30≤γ R≤1.39 when the member axial compression, 1.10≤γ R≤1.40 when eccentric compression, Shear, 1.10≤γ R≤1.15 when axial tension, eccentric tension, flexure.Structure is identical with the level of reliability of design reference period in the construction time.Accordingly, just can utilize this paper program to calculate variable action or wind load factor k under different partial safety factor for resistance γ R and the construction time Tc.
R, S obeyed the wind load factor k of the ultimate limit states design of different distributions when table 1 and table 2 had provided respectively β=β c=3.2.As expected, k increases with the increase of Tc, reduces with the increase of γ R.
Wind load factor k (β=β c=3.2) when table 1R and the equal Normal Distribution of S
Wind load factor k (β=β c=3.2) when table 2R obeys logarithm normal distribution, S obey the distribution of extreme value I type
According to the loading characteristic of cooling tower, get γ R=1.15, to β=β c=3.2, calculate the k value of Tc<1 year, and associative list 1, table 2 result, the calculated value by k provides recommendation, respectively shown in table 3, table 4.Can find out that k reduces gradually with the amplitude of variation of Tc.
Wind load factor k (β=β c=3.2, γ R=1.15) when table 3R and the equal Normal Distribution of S
Wind load factor k (β=β c=3.2, γ R=1.15) when table 4R obeys logarithm normal distribution, S obey the distribution of extreme value I type
In structural reliability problem, to only have when structural failure Probability p f 〉=0.001 (or reliability index β≤3.0902), the result of calculation of pf is just insensitive to the distribution form of variable.The β of the problems referred to above=3.2, so be necessary in the research of cooling tower wind load value, to consider the actual distribution form of R and S, and R obeys logarithm normal distribution, S are obeyed extreme value I type and are distributed more realisticly than the equal Normal Distribution of R, S, therefore have reason to think that the result of table 2, table 4 is more reasonable.
Use for convenient, utilize this paper program to calculate more k value, put the result of table 4 in order merger, as shown in table 5.
Table 5 cooling tower tube construction phase wind load factor k (β=β c=3.2, γ R=1.15)