Научная статья на тему 'RELIABILITY ESTIMATION OF INDUSTRIAL BUILDING STRUCTURES'

RELIABILITY ESTIMATION OF INDUSTRIAL BUILDING STRUCTURES Текст научной статьи по специальности «Строительство и архитектура»

CC BY
123
25
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Magazine of Civil Engineering
Scopus
ВАК
RSCI
ESCI
Ключевые слова
СООРУЖЕНИЕ / ЗДАНИЕ / КОНСТРУКЦИЯ / НАДЕЖНОСТЬ / ОТКАЗ / ВЕРОЯТНОСТНЫЙ МЕТОД / НАГРУЗКА / СТАТИЧЕСКИ НЕОПРЕДЕЛИМАЯ СИСТЕМА / CONSTRUCTION / BUILDING / STRUCTURE / RELIABILITY / FAILURE / PROBABILISTIC METHOD / LOAD / REDUNDANT SYSTEM

Аннотация научной статьи по строительству и архитектуре, автор научной работы — Pichugin S.F.

The article is devoted to solving the problem of the reliability of industrial building structures: crane and roof beams, columns, redundant building systems. For this aim, the probabilistic method of structural reliability estimation is developed after the criterion of bearing strength, in which the main component of structure reliability is faultlessness. The method takes account of the random loads and material strength, loads joint action, the specific character of work and failure elements, nodes and the whole structure as well. For the ground of method, the large amount of statistic results on crane load is examined for the bridge cranes of the different types. A large amount of wind and snow meteorological data is collected for the territory of Ukraine. Stationary probabilistic model of crane load and quasi-stationary model of the snow and mean wind load are substantiated. The most widely spread probabilistic presentations of random loads are observed. They are as follows: stationary random process and its absolute maxima, random sequence of independent and correlated loads, discrete presentation and extreme model. It was deduced that the redundant structure reliability estimation is a very complicated problem as depends upon the system complexity. The method of states, a probabilistic method of ultimate equilibrium and logic and probabilistic method are developed for solving this problem. On the base of the determined method, the numerical reliability computations of a wide range of industrial building structures are realized. It is shown that the structures have quite different levels of reliability. In particular, the light roof structures are not reliable enough being under the great influence of snow load. At the same time, the Design Code allows over-estimation of reliability for industrial columns. The estimation of industrial redundant structures with a different degree of redundancy is obtained on the base of developed approach. It gave the possibility to evaluate the high safety level of redundant structures in comparison with separate members and statically determined structures. With regard to mentioned results, it is recommended to correct some load factors, a combination factor and a factor for model uncertainties of the Design Codes of structures and loads.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

ОЦЕНКА НАДЕЖНОСТИ КОНСТРУКЦИЙ ПРОМЫШЛЕННЫХ ЗДАНИЙ

В статье рассмотрены вопросы оценки надежности конструкций промышленных зданий: подкрановых и стропильных балок, колонн, статически неопределимых конструктивных систем. Для этого разработан вероятностный метод оценки надежности конструкций по критерию несущей способности, основой которого является безотказность конструкции. Метод учитывает случайные нагрузки и прочность материалов, совместное действие нагрузок, специфический характер работы и отказов элементов, узлов и сооружений в целом. Для обоснования метода был сформирован обширный объем статистических данных по нагрузкам мостовых кранов различных типов. Большое количество метеорологических данных по ветру и снегу было собрано для территории Украины. Обоснованы стационарная вероятностная модель для крановой нагрузки и квазистационарная модель для снеговой и ветровой нагрузок. Рассмотрены наиболее распространенные вероятностные представления случайных нагрузок: стационарный случайный процесс и его абсолютные максимумы, случайная последовательность независимых и коррелированных нагрузок, дискретное представление и экстремальная модель. Было подтверждено, что оценка надежности статически неопределимых систем – достаточно сложная проблема, поскольку она зависит от сложности таких систем. Для решения данной проблемы разработаны метод состояний, вероятностный метод предельного равновесия, логико-вероятностный метод. На основе разработанного метода выполнены численные расчеты надежности широкого круга строительных конструкций. Показано, что различные конструкции имеют разный уровень надежности. В частности, недостаточно надежными могут быть легкие стропильные конструкции при действии значительных снеговых нагрузок. В то же время нормы проектирования закладывают в колонны промышленных зданий излишний запас надежности. На основе разработанного подхода выполнены расчеты надежности статически неопределимых систем разной степени сложности, которые дали возможность оценить повышенный уровень надежности таких систем по сравнению с отдельными элементами и статически неопределимыми системами. В результате были обоснованы коэффициенты норм проектирования, такие как коэффициенты надежности по нагрузке, коэффициент сочетания нагрузок, коэффициент условий работы.

Текст научной работы на тему «RELIABILITY ESTIMATION OF INDUSTRIAL BUILDING STRUCTURES»

S.F. Pichugin,

Poltava National Technical Yuri Kondratiuk University, Poltava, Ukraine

Key words: construction; building; structure; reliability; failure; probabilistic method; load; redundant system

doi: 10.18720/MCE.83.3

Reliability estimation of industrial building structures Оценка надежности конструкций промышленных зданий

Д-р техн. наук, заведующий кафедрой С.Ф. Пичугин,

Полтавский национальный технический университет имени Юрия Кондратюка, г. Полтава, Украина

Ключевые слова: сооружение; здание; конструкция; надежность; отказ; вероятностный метод; нагрузка; статически неопределимая система

Abstract. The article is devoted to solving the problem of the reliability of industrial building structures: crane and roof beams, columns, redundant building systems. For this aim, the probabilistic method of structural reliability estimation is developed after the criterion of bearing strength, in which the main component of structure reliability is faultlessness. The method takes account of the random loads and material strength, loads joint action, the specific character of work and failure elements, nodes and the whole structure as well. For the ground of method, the large amount of statistic results on crane load is examined for the bridge cranes of the different types. A large amount of wind and snow meteorological data is collected for the territory of Ukraine. Stationary probabilistic model of crane load and quasi-stationary model of the snow and mean wind load are substantiated. The most widely spread probabilistic presentations of random loads are observed. They are as follows: stationary random process and its absolute maxima, random sequence of independent and correlated loads, discrete presentation and extreme model. It was deduced that the redundant structure reliability estimation is a very complicated problem as depends upon the system complexity. The method of states, a probabilistic method of ultimate equilibrium and logic and probabilistic method are developed for solving this problem. On the base of the determined method, the numerical reliability computations of a wide range of industrial building structures are realized. It is shown that the structures have quite different levels of reliability. In particular, the light roof structures are not reliable enough being under the great influence of snow load. At the same time, the Design Code allows over-estimation of reliability for industrial columns. The estimation of industrial redundant structures with a different degree of redundancy is obtained on the base of developed approach. It gave the possibility to evaluate the high safety level of redundant structures in comparison with separate members and statically determined structures. With regard to mentioned results, it is recommended to correct some load factors, a combination factor and a factor for model uncertainties of the Design Codes of structures and loads.

Аннотация. В статье рассмотрены вопросы оценки надежности конструкций промышленных зданий: подкрановых и стропильных балок, колонн, статически неопределимых конструктивных систем. Для этого разработан вероятностный метод оценки надежности конструкций по критерию несущей способности, основой которого является безотказность конструкции. Метод учитывает случайные нагрузки и прочность материалов, совместное действие нагрузок, специфический характер работы и отказов элементов, узлов и сооружений в целом. Для обоснования метода был сформирован обширный объем статистических данных по нагрузкам мостовых кранов различных типов. Большое количество метеорологических данных по ветру и снегу было собрано для территории Украины. Обоснованы стационарная вероятностная модель для крановой нагрузки и квазистационарная модель для снеговой и ветровой нагрузок. Рассмотрены наиболее распространенные вероятностные представления случайных нагрузок: стационарный случайный процесс и его абсолютные максимумы, случайная последовательность независимых и коррелированных нагрузок, дискретное представление и экстремальная модель. Было подтверждено, что оценка надежности статически неопределимых систем - достаточно сложная проблема, поскольку она зависит от сложности таких систем. Для решения данной проблемы разработаны метод состояний, вероятностный метод предельного равновесия, логико-вероятностный метод. На основе разработанного метода выполнены численные расчеты надежности широкого круга строительных конструкций. Показано, что различные конструкции

имеют разный уровень надежности. В частности, недостаточно надежными могут быть легкие стропильные конструкции при действии значительных снеговых нагрузок. В то же время нормы проектирования закладывают в колонны промышленных зданий излишний запас надежности. На основе разработанного подхода выполнены расчеты надежности статически неопределимых систем разной степени сложности, которые дали возможность оценить повышенный уровень надежности таких систем по сравнению с отдельными элементами и статически неопределимыми системами. В результате были обоснованы коэффициенты норм проектирования, такие как коэффициенты надежности по нагрузке, коэффициент сочетания нагрузок, коэффициент условий работы.

1. Introduction

The building structure reliability is the problem of a high priority today. It can be treated as a scientific trend, which does not need both material and financial expenses.

The eminent scientists created the classic works in building reliability [1-3], numerous foreign and native scientists were involved in the study of this problem [4-12]. The general probabilistic models of wind, snow and crane loads can be considered as a statistic basis of reliability estimation of industrial building structures. Loads of bridge cranes were considered in [13], the description of the wind load was summarized in [14]. Active researches of the wind load proceed the last years [15-17]. The analysis of snow load was resulted in [18], the influence of global temperature rise on the snow load is analyzed in work [19]. In different countries, systematic researches of climatic influences are conducted with the purpose of Codes correcting [20]. The perspective model is a correlated casual sequence of loads (based on the method of the generalized covariance, developed by A.P. Kudzis [21]). The linearizing of the functions of random parameters worked out in work [22]. However, these approaches lack both systematic analysis and model comparison, which can cause different results in structure reliability design.

The probabilistic computation of structures has not been developed with regard to joint loads application, their real distributions, and frequency characteristics. At the same time, active studies of probabilistic models for concrete elements proceed the last years [23], in particular, the initial reliability of reinforce-concrete elements was defined [24]. However, the estimation of reliability is yet absent for many buildings and structures. The problem of reliability continues to be probed [25-30] and expects further investigations and solutions. The exploitation experience demonstrated the existence of quite different reliability of various structures. The reliability of existing structures has to be examined together with the reliability of elements (sections) as well as nodes [31-33]. The assessment of a redundant structure safety is rather a difficult task. In spite of many studies that have been done the problem of reliability analysis of these structures has not been solved yet.

As follows from the preceding, the evolution of reliability computation theory and design Codes of building structures are still of interest because of the complexity of a problem on the one hand and on the other of the ignoring the random loads.

The article is devoted to solving the problem related to the assessment of the reliability of industrial building structures: crane and roof beams, columns, redundant building systems. This task deals with the developing of the general method of reliability estimation of industrial building structures after the criterion of bearing strength, in which the main component of structure reliability is faultlessness. At such interpretation the estimation of structure reliability includes the decision of the following questions: probabilistic description of loads; mechanical descriptions of materials, joints and other random parameters; reliability estimation of elements of building structures; calculation of reliability of the structural systems taking into account the possible character of their failure; quantitative estimation of reliability of buildings and structures of the different setting.

2. Methods

It was put and decided the complicated and difficult task of formation of the generalized probabilistic models of crane and atmospheric loads. The decision was based on sufficient statistical material, took into account the specific features of the different loads, type of their distributions and frequency structure.

The stochastic model of the crane loads was built as a result of experimentally statistical researches. The vertical and horizontal loads were investigated for bridge cranes with rigid and flexible hanged cargo of different carrying capacity and regimes which were exploited in the conditions of metallurgical and machine-building shops from 10 to 30 years of service. The generalized numerical and frequency descriptions of the crane loads were grounded and the possibility of application for them of normal low was substantiated (Figure 1, d). The special attention in model experiments was spared lateral forces of bridge cranes.

The systematic information about the wind velocity measurements done with ten minutes average at 70 Ukrainian meteorological stations were used as an initial data. The mean value of wind load is of quasi-stationary origin with the constant frequent parameters and normalized distribution. For the description of ordinate density of wind mean value Veibull's law was used, for wind dynamic value the normal distribution was based. Having integrated some initial data all necessary mean parameters of probabilistic wind models were determined for territory of Ukraine.

The results of regular snow measurements for 15-40 years at 62 Ukrainian meteorological stations have been taken as the reference statistic material for the snow probabilistic model. It is substantiated that the ground snow load is of a quasi-stationary origin. Its mathematical expectation and standard have a seasonal trend. At the same time snow frequent characteristics and normalised ordinate distribution remain constant during the season. For its description it was first applied polinomo-exponential distributing (Figure 1, b). In future in PoltNTU there were executed the researches of deposits of snow on coverages with the overfalls of heights, the estimation of influence of roofs heating descriptions on the value of the snow load, the calculation of the snow load on the cold roofs of buildings with the positive internal temperatures of air.

composition p=O.S composition p-1.0

Figure 1. Distributions of random arguments: a - durability of steel; b - snow load; c - wind load; d - crane load; e - composition of polinomo-exponential and normal distributions; f - difference of polinomo-exponential and normal distributions; g - composition of Veibull's and normal distributions; h - difference of Veibull's and normal distributions

The systematic analysis of random loads was realized for the six most commonly used probabilistic models. The main one is presented in the form of a stationary (crane load) or quasi-stationary (wind and snow loads) random processes (Figure 2, a, pos.1); their parameters are an effective frequency m and

the coefficient of trend Kr which accounts the atmospheric load season change.

The absolute maxima of random process are one of simpler model; they are determined by the tail part of the distribution of outliers and are higher then characteristical maximum level y0 . The letter is a

solution of the following equation N+ (yo; 0 <t< t) = 1, where N+ (•) - the number of outliers of random process.

The model in the form of a random sequence of independent random loads with the intensity X is widely spread (X - the number of loads in per-unit time tx). For discrete presentation of loads (Figure 2,a,

pos. 2) the frequent parameter of which is the mean duration of overloading A connected with the intensity by the ratio A = tx / X .

The analysis of the problem has demonstrated that the load values sampling can be classified as the exponential type. That's why their maximum values can be presented correctly by the extreme double exponential Gumbel distributions of a normalized type:

y = an (Y- un) ,

(1)

where un - characteristic extremum; an - extreme intensity; y = - ln[- ln F(t)] - Gumbel distribution argument.

The normalised load level was taken into account y = (x - x)/x, where x - the load ordinate, x -mathematical expectation, x - standard deviation.

The choice of models of loads depends on a specific of solving probabilistic problems: the more complicated ones are solved in the manner of random processes which, however, are difficult for description and take much more computation time. More accessible and simple models, mentioned earlier, are based on the random values and corresponding frequent characteristics and provide not less exact solution if they are proper grounded. In particular, all the examined models are close to its sense of the following evaluation:

Q(t )=

= -Q = Xt[l - F (y)] = thFtii = [r (1 -p) + p][l - F (y)] ,

(2)

where Q(t) - probability of exceeding y - level during t; F(t) and f (y) - accordingly integral and differential functions of load distribution; r - amount of calculation cuts of casual sequence of loads; p - generalized coefficient of correlation of sequence.

Load probabilistic comparison can be well performed at the extreme scale which is illustrated in Figure 2, b. On the axis of ordinate of a scale the standardized load is laid off, on the axis of abscissa the

Gumbel's distribution argument y is laid of which is connected with the load return period T. Gumbel's distribution (1) is described on the scale in the form of straight lines. The models of a random process, a random sequence of independent loads and discrete presentation are introduced as different curves. The main advantage of this scale is its visual effect of the tail parts of the load distributions which have rather small distinctions in the usual form of presentation. It enables to present the visual comparison and correspondence of parameters of different load models.

b) a)

a - realization of random process: 1 - model of continuous process; 2 - discrete model;

b - general form of load presentation

A model of random processes, without regard to its relative complication, is for today most studied, the model of absolute maximums of random processes also enabled to get interesting scientific results. It should be noted that other perspective models, the especially random sequences of independent and correlated overloads, and also discrete presentation of loads, are less studied and recommended for the use in subsequent researches of loads.

The worked out random parameters give possibility to develop the reliability estimation of building structures.

3. Results and Discussion

3.1. Reliability estimation of building structures in the technique of random values

The failure of an element takes place when a stochastic stress under the joint applied loads S (t) (t-time) exceeds the resistance of an element R(t). The failure of the element is defined by the equation

Y(t) = R(t)- S(t) < 0, (3)

where Y (t) - margin function of load carrying capacity.

The probabilistic technique of random values without the account of time factor is applied in this section. The technique is grounded at the action of loads with little change in time (dead and some technological loads) or those which have non-permanent character. Then the margin function will be written down as

Y = R - S < 0 . (4)

Mathematical expectation and standard mean deviation of margin function are determined as for a linear function:

Y = R - S ; Y =V R2 + S2 . (5)

The safety characteristic is of great importance and is derived from equation

p = Y / Y = (R - S)/(R2 + S2)V2, (6)

where Y,R,S - the corresponding mathematical expectations, Y,R,S - the corresponding standard mean deviations.

This characteristic determines probability of structural failure

Q(Y < 0) = Fy (0) = FY (Y -pY), (7)

where Fy (.) - an integral function of margin function distributing.

In the case of normal distribution f (Y) the safety characteristic is very convenient:

Q(Y < 0) = 0.5-0(p) , P(Y > 0) = 0.5 + 0(p) . (8)

Here o(p) - the well-known function of Laplace contained in widespread statistical tables.

The element resistance and margin function of the load carrying capacity (3) are generalized by the function of a few arguments:

Y = f (X i, X 2,..., Xn ).

In general case this function can be nonlinear, and for simplification of calculations, its linearizing is often executed with replacement of initial expression the linear function:

Y = f (Xl5X2,...,)+£ ^L(Xl5X2,.., J„fc -X) = Y + £Di(x, - Xt).

i=ldX i

(9)

i=1

In the calculations of reliability the linearizing is executed near the mathematical expectation of function, here the angular coefficient of entered line is determined by the partial differentiation of initial

function on the proper argument. Error in determination of Y as a result of linearizing for (8) it is possible to estimate the following expression:

^ 1 n d2Y ^ = ^

AY =

i=i

2 ^ dX i

■(xl, X2,...XW)Xi .

Here Xi - dispersion of i - argument.

At insignificance of this error a formula for determination of mathematical expectation (9) is simplified and taken to the substitution of mathematical expectations of arguments in an initial function:

Y ( Xi, X 2,..., Xn ) = f (Xi, X2,..., Xn ) . (10)

Standard deviation of the margin function is determined, in the case of independence of arguments, from equation:

Y = #:

I

i=i

dY

~dXi

( X ^ X 2,...., Xn )

X 2 =

I d2 X 2

i=i

(11)

Error of this expression as a result of linearizing is estimated as follows

n2

AY = A

)=i

=21

2 i=1

d2Y i- -

dX 2

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(x 1, X 2,..., Xn )

XX 2 +I

i > j

d2 Y - - -(X 1, X2,...., Xn )

dXj dX j

2

Xi X

ij

(12)

3.2. Practical reliability estimation of structure elements

In quality an example we will consider the reliability calculation of reinforced concrete beams with carbon-plastic external strengthening. For the receipt of reliability estimations we will use the worked out reception with the substitution of probabilistic parameters in the deterministic decisions of durability of reinforced concrete beams. Taken into account thus, that most random arguments of the margin function of load carrying capacity of reinforced concrete beams can be reasonably described by a normal law, in particular, durability of concrete, armature, carbon-plastic, and also row of loads (dead, technological, crane etc.).

Random value of maximum bend moment which is perceived by a beam with a double armature:

2

Mult = f (<b<> < ) = ¿¿.A (ho - 05x) + <rscAs (0.5X - a'), (13)

where: <jb - random value of concrete resistance to the compression for the limit states of the first group; ¿~s - random value of the armature durability to tension; As - section area of the stretched armature; ¿7sc -random value of concrete resistance to tension; A.- section area of the compressed armature; a' - distance from resultant effort in the compressed armature to the compressed verge of element; ho - calculation height of section; x - height of the compressed zone of concrete, equal

7 = (<sAs -<scAs Xbb ■

where: b - width of section.

We put expression for x in a formula (13):

0 5 2

Mult =&sAsh0 -<?scA'sa'-^r(àsAs -<?scAl) . (14)

<bb

We have the mathematical expectation of limit moment putting the expected values of casual arguments in this common expression.

We will define coefficients for the calculation of standard of limit moment:

Dsc = dMulL = A- [- <bba' + (<sAs - <scAS)] (15)

d<sc <bb

A =drL = < ¿b-(<A-¿Ac )] ; (16)

D dMult 0 5 (< A < A)2

Db "-= —27 (<sAs - <cAs ) . (17)

d<b <bb

The numerical values of coefficients we get, putting the expected values of random arguments in these expressions.

The standard of limit bend moment is determined as

Mult = ^(Pb«b )2 +Ds<. )2 +(DSc<JSc )2 . (18)

For the reliability estimation of beams we determine a safety characteristic, having in this case next

kind:

P = (Mult - Mcal )/Mult, (19)

ult cal}/ ult

where: Mcai - calculation value of external bend moment in a beam.

Farther we will consider a beam strengthened by carbon-plastic (FAP). Random value of the limit bend moment, perceived by an increased beam:

Mult = f (<b <s > ¿sc <fu ) = <fuAf (h - 05X ) + <sAs ( h - 05X ) + <sA ( 05X - a'), (20)

where: <f - random value of durability on tension of FAP; Af - area of section of carbon-plastic armature; h - height of section; x - height of the compressed zone of concrete, equal

7 = <fuAf + <sAs - <scA's)<bb . We put expression for x in a formula (20):

Mult =<fuAfh+<AK{ôfA+<A-<A) . (21)

Mathematical expectation of limit moment we get, putting the expected values of random arguments in this expression.

We will define coefficients for the calculation of standard of limit moment:

(25)

(22)

(23)

(24)

The numerical values of coefficients we get, putting the expected values of random arguments in these expressions.

The standard of limit bend moment is determined as

For the reliability estimation of beams, we determine a safety characteristic utilizing formulae (8), (19) and the function of Laplace.

Application of the got results is illustrated on a numerical example with a reinforce-concrete beam. A beam has section sizes b = 300 mm, h = 700 mm, external calculation bend moment is Mcal = 500 kNm. After the above-mentioned formulae, numerical descriptions of maximum moment of

beam are Mult = 762.7 kNm, Mult = 29.92 kNm. The safety characteristic in relation to the calculation value of external bend moment in a beam makes / = (762.7 - 500) / 29.92 = 8.78, that answers the low probability of beam failure Q (// = 4.0 • 10-19. Farther on the terms of example, a beam is loaded with the increase of calculation moment to M'cal = 750 kNm. The safety characteristic here is diminished substantially / = (762.7 - 750) / 29.92 = 0.42, the probability of failure Q (/)= 0.337 becomes too high, that testifies to the accident rate of beam and necessity of its strengthening.

The strengthening of the beam is one layer of carbon-plastic fibre with width 250 mm, thickness 1.4 mm, area section of strengthening Af = 350 mm2 =3.5 cm2. Calculated durability of carbon-plastic fibre is

Rfu = 1071 MPa. Numerical descriptions of maximum moment of beam are increased -

Mult = 991.0 kNm, Mult = 35.5 kNm. The safety characteristic in relation to the calculation value of

heightened external bend moment in a beam makes / = (991.0 - 750) / 35 / 5 = 6 / 79, the probability of

failure of the increased beam Q (/)= 13 • 10-11 is substantially below than failure probability of beam

without strengthening. The got estimation of failure probability testifies to sufficient reliability of beam and efficiency of its strengthening of carbon-plastic fibre.

The analysis of results of reliability calculation of the increased beam confirmed the validity of application for this purpose the described higher linearized probabilistic model, as an insignificant error is here assumed: - 0.12 % for the mathematical expectation of maximum moment and 0.43 % for its standard.

For elements which perceive a multicomponent dead load a stochastic model as a sum of random values allows to enter the reducing coefficient of combination y = 0.90...0.95. It allowed also estimating the reliability of linear part of main pipeline and faultlessness of bearings and non-load-bearing structures from steel of the thin-walled types.

The compositions of sums and differences of loads and strength should be done for several loads and for design of margin of carrying capacity (4). Convolution formulae are used for this purpose.

Mult =^(Dfu°fu )2 + (Db°b )2 + (Ds°s )2 + (Dsc°sc )2

(26)

Y

Y = X1 + X 2, f(Y)=Jf1(X1)f2(Y-X1)dX1 ;

0

Y = Xi — X 2, f(Y )= J fl(Xl )f2 X - Y )dX\ . (27)

—X

They cannot be determined in a closed form for the mentioned distributions and therefore numerical integration is adopted. The convolutions are transformed (presented in an unified comparable form) for the design simplification and their parameters are combined. From 13 received analytical expressions of convolutions, only the formula for difference of normal and arbitrary f2 (z) density distributions is suggested

т~\ z2

f (т) = ~^ f exp (-0.5E2 ) f2 (z) dz , (28)

V2n z.

where: y = (y - Y )/Y- standardised variable; E = D -y + z / p ; D =yj 1 + p2 ; p = x n / x a ; vn = xn / xn; xn; xn - parameters of the normal distribution, xa - standard deviation of arbitrary

distribution, zi = yV 1 + p2 - p/vn - lower limit of integration, z2 - upper limit of integration which is determined with the account of necessary calculation accuracy.

The examples of joint distributions based on the proposed formulae are presented on the Figure 1, d-g. They can be non-symmetrical and differ considerably from the normal distribution in solutions of structural reliability problems. Without the account of this factor the estimation of structure reliability can have noticeable inaccuracy.

3.3. Reliability estimation of structure elements with the account of time factor

The loads are presented in the form of stationary and quasistationary random processes. Hence the Y (t) -function is a random process as well. Therefore the element failure is a projection of random process

of the load carrying capacity margin (3) into the negative region. In this case the probability of failure Q(t) is estimated as follows

Q(t )=0qfy (P)t/(P№42n), (29)

where: pq and f (ft) - the effective frequency and density of the distribution of the random process Y(t)

ordinate; - the coefficient of structural complexity of the random process Y(t) which takes into account

the spectrum of frequencies of the real loads; t - the work time of element (different from his term of service).

The solving of the problem (29) demands the presentation of the frequency characteristics of summary random processes of loads and Y (t). Frequency description of load connection, presented in the form of quasistationary random processes is determined as

(pi • Kri )2 +(®2 • K ■ Kr 2 )2 ]2 , (30)

1

®12 =

Vl + K2

where: PP2P12 - effective frequencies of separate random processes and their sum; K = X2/X1 - a relation of standards of loads which are composed; Ktri, Ktr2 - trend coefficients that take into account the slow seasonal changes of the atmospheric loads.

The frequency analysis of the composition of the random loads are constantly investigated by scientists [34, 35].

It is necessary to underline that the investigated random loads are substantially differed on frequency, it is evidently shown their realization (Figure 3, a,b) and also substantially different spectral densities (Figure 3, d). Effective frequencies of random processes of various loads are differed considerably: for crane loads p = 1700 - 5160 1/day, for wind ones = 5.4 - 6.6 1/day, for snow

loads < = 0.073 - 0/141 1/day. Therefore the joint action of several loads becomes multifrequent (Figure 3, c).

This peculiarity can be defined in the formula (29) by the coefficient of structural complexity of the random process which equals the ratio of mean frequency on maximum <am to the effective frequency

on zero a>q :

<m _ KIV (0)•k(0) 1/2 = A • (1 + k2 • ©4)(1 + k2)f2

<q - K (0)] (1 + k2 ©2)

(31)

where: K (0), K (0), KIV (0) - random process correlation function and its derivatives with zero argument; - structure complexity coefficient of separate loads; k = x2 / x1 ; © = <2/< - standard deviations and effective frequency proportion of summed up random loads.

In accordance with equation (31) the frequency composition analysis can be realised if the load random process is differentiated sufficiently many times. For the atmospheric and crane loads the coefficient of structure complexity differs negligibly from Pt = 3 . As a result of substantial difference of frequencies of loads, this coefficient of their composition can be considerable.

If the effective frequency of one of the loads prevails and < is large, then formula (31) is simplified as follows:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

P<=Pl (1 + k 2)1/2/ k,

(32)

Crane beams as well as crane trestle structures can be treated like elements under only crane load. We will examine general stress of these elements without local effects and fatigue. In this case the probability of structure failure Q(t) depends on the ratio of a load stress of one bridge crane XM1 to the load stress of two bridge cranes XM 2, then on the ratio of a load carrying capacity of cranes to their mass and it also depends on crane work condition. The analysis showed that elements have a deficient reliability in the case of the one crane loads dominance. In the rest cases the elements reliability under crane loads is sufficient. In some cases the design crane load on structures can be decreased: if the approach of a loaded crane trolley to a crane track is limited or if the structure service term T is also limited.

The reliability of structures (beams, trusses) under snow load was estimated by a developed method. The calculation demonstrated the lack of reliability of these structures. That justifies the idea of understating of existing snow loads in the norms of SNiP, which operated on territory of Ukraine to 2006 year. Besides this fact validates the causes of steel truss failures. It applies to steel structures with lightweight roofs in the southern districts of Ukraine when there is much snow in winter. The increasing of snow design load for 1.5-2 times for our region can solve the reliability problem of steel structures under snow load. It was executed in national Code in Ukraine DBN V.1.2:2-2006 "Loads and influences". The reliability of steel roofs under snow load is constantly probed by the specialists of different countries [36].

«1

«■■MM

m

X/iD

X,ft)

01 liiav

Figure 3. To the estimation of frequency of load composition a - high-frequency component; b - low-frequency component; c - total randoml process; d - comparison of the load spectral densities

Structure elements under wind load designed in accordance with the existing code (glass elements, wind protection screens etc.) are of sufficient reliability. So if the service term of these elements will be

limited it is possible to introduce the temporary coefficient yT < 1 for the design wind load. It is

recommended to introduce the increasing coefficients of combination y = 0.7 - 0.9into the structure

design under snow, wind and crane loads.

The special attention was spared the evaluation of reliability of steel beam-column structures. Time factor, existing loads, random steel strength were taken into account during computation process of these elements. Existing steel columns of industrial buildings in a broad range of parameters were examined. The general conclusion is as follows: the reliability of steel columns of industrial buildings is quite sufficient. In quality tasks for subsequent researches it is possible to name the search of direct probabilistic decisions of tasks of compressed elements and subsequent development of probabilistic version of method of finite elements.

3.4. Reliability analysis of structure nodes and redundant systems

It was founded out that reliability of typical structure nodes could be compared with the schemes of successive connections of correlated elements. Thus, the node reliability depends upon the number of engaged independent load carrying elements.

Some beams and simple frames, as well as multi-storey and multi-span structures of industrial and residential buildings present the redundant structures. Redundant structure failures occur after some member failures in the form of transmission to different workable states. Thus, the redundant structure failure estimation is a very complicated problem as depends upon the system complexity. Illustration of it is a multibranch graph of the states of twice statically indefinite system (Figure 4, a). The complete estimation of the failure probability of elastic-plastic redundant systems with random durability and loading was got as a result of application of probabilistic method of limit equilibrium. A method allows, walking around the transient states of system, directly to examine the reliability of the last (true) mechanism of failure of the system (Figure 4, b). Corresponding formulae are derived; the algorithms and computing programs are developed. The estimation of a wide range of industrial redundant structures with different degree of redundancy was obtained on the base of this approach. It gave possibility to evaluate the safety level of redundant structures in comparison with separate members and statically determined structures. This level can be taken into account introducing the additional coefficient of work condition yc = 1.1 - 1.4.

a - graph of the states of twice statically indefinite system; b - true mechanism of frame failure

Among questions, connected with the estimation of reliability of redundant systems, which remained unsolved, it follows to name formalization and translation on personal computer of method of the states and logic-probabilistic method, and also classification of redundant system elements, after the criterion of reliability. This direction becomes especially actual in connection with development of scenarios of failures, questions of vitality and progressive failure of buildings [37-39].

4. Conclusions

1. The reliability analysis of wide range of structures designed in accordance with existing Design Codes is examined in the article. Considered structures of industrial buildings are as follows: crane and roof beams, columns, redundant building systems.

2. This analysis exposes that the structures have quite different levels of reliability. It is shown that the light roof structures are not reliable enough being under the great influence of snow load. At the same time, the Design Code allows over-estimation of reliability for steel columns and redundant structures.

3. For the achievement of the put aim, the practical method of reliability analysis of structural elements is worked out. This method takes into account real load distributions and frequency characteristics.

4. The general probabilistic models of wind, snow and crane loads were designed and taken as a statistic basis of reliability estimation. These models are as follows: stationary and quasi-stationary random processes, its absolute maxima, random sequences of independent and correlated loads, discrete presentation and extreme model.

5. The simplest method of structure reliability estimation is realized in the technique of random values. Reliability estimation of reinforced concrete beams with carbon-plastic external strengthening is presented as an example.

6. Substantial results are got in the probabilistic techniques of random processes. It is important to notice that correlation functions and effective frequencies of random processes of various loads differ considerably. Therefore the joint action of several loads becomes multi-frequent. This peculiarity can be defined by the structure complexity coefficient (3m which equals the ratio of mean frequency on maximum to the effective frequency on zero.

7. There is confirmed that redundant structure failure estimation is a very complicated problem as depends upon the system complexity. The method of states, a probabilistic method of ultimate equilibrium and logic, as well as a probabilistic method is developed for solving this problem.

8. With regard to mentioned above, it is recommended to correct some load factors, a combination factor and a factor for model uncertainties of the Design Codes of steel structures and loads.

References

1. Streletskiy N.S. Osnovyi statisticheskogo ucheta koeffitsienta zapasa prochnosti sooruzheniy [Bases of the statistical account of constructions safety factor]. Moscow: Stroyizdat. 1947. 92 p. (rus)

2. Rzhanitsyin, A.R. Teoriya rascheta stroitelnyih konstruktsiy na nadezhnost [The theory of calculation of building structures reliability]. Moscow: Stroyizdat. 1978. 239 p. (rus)

3. Bolotin, V.V. Statistical methods in structural mechanics. San Francisco: Holden Day. 1969. 276 p.

4. Freudenthal, A.M. The Safety of Structures. Proc. ASCE. 1947. Vol. 112. No. 1. Pp. 125-180.

5. Benjamin, J.R. Cornell, C.A. Probability, Statistics and Decision for Civil Engineers. New York: McGraw-Hill. 1970. 684 p.

6. Thoft-Christensen, P., Baker, M.J. Structural reliability theory and its applications. Berlin: Springer. 1982. 267 p.

7. Augusti, G., Baratta, A., Casciati, F. Probabilistic methods in structural engineering. London: Chapman & Hall.1984. 583 p.

8. Ditlevsen, O., Madsen, H. Structural reliability methods. Chichester: John Wiley. 1996. 372 p.

9. Elishakoff, I. Probabilistic Theory of Structures. New York: Dover Publications, 1999. 502 p

10. Raizer, V. Reliability of Structures. Analysis and Application. USA: BACKBONE Publishing Company. 2009. 146 p.

11. Perelmuter, A.V. Izbrannyie problemyi nadezhnosti i bezopasnosti stroitelnyih konstruktsiy [The selected problems of building structures reliability and safety]. Moscow: Izdatelstvo ASV, 2007. 256 p. (rus)

Литература

1. Стрелецкий Н.С. Основы статистического учета коэффициента запаса прочности сооружений. М.: Стройиздат, 1947. 92 с.

2. Ржаницын А.Р. Теория расчета строительных конструкций на надежность. М.: Стройиздат, 1978. 239 с.

3. Bolotin V.V. Statistical methods in structural mechanics. San Francisco: Holden Day.1969. 276 p.

4. Freudenthal A.M. The Safety of Structures // Proc. ASCE, 1947. Vol. 112. № 1. Pp. 125-180.

5. Benjamin J.R. Cornell C.A. Probability, Statistics and Decision for Civil Engineers. New York: McGraw-Hill. 1970. 684 р.

6. Thoft-Christensen P., Baker M.J. Structural reliability theory and its applications. Berlin: Springer. 1982. 267 p.

7. Augusti G., Baratta A., Casciati F. Probabilistic methods in structural engineering. London: Chapman & Hall.1984. 583 р.

8. Ditlevsen O., Madsen H. Structural reliability methods. Chichester: John Wiley. 1996. 372 р.

9. Elishakoff I. Probabilistic Theory of Structures. New York: Dover Publications. 1999. 502 р.

10. Raizer V. Reliability of Structures. Analysis and Application. USA: BACKBONE Publishing Company. 2009. 146 p.

11. Перельмутер А.В. Избранные проблемы надежности и безопасности строительных конструкций. 3-е изд., перераб. и доп. М.: Изд-во АСВ, 2007. 256 с.

12. Пичугин С.Ф. Надежность стальных конструкций производственных зданий. М: Изд-во АСВ, 2011. 456 с.

12. Pichugin, S.F. Nadezhnost stalnyih konstruktsiy proizvodstvennyih zdaniy [Reliability of industrial building steel structures]. Moscow: Izdatelstvo ASV, 2011. 456 p. (rus)

13. Pichugin, S. Analysis of bridge crane loads on industrial buildings. XLIV Konferencja Naukowa KILiW PAN i KN PZITB. Tom VII. Poznan-Krynica, 1998. Pp. 171-178.

14. Pichugin, S. Probabilistic analysis on wind Load and reliability of structures. Proc. of the 2 EACWE. Vol. 2. Genova. Italy. 1997. Pp. 1883-1890.

15. Kinash, R.I.; Huk, J.S. Geographical altitude coefficient utilisation for wind gust loading determination in the mountainous terrain of the Transcarpathian region. Technical Transactions, Civil Engineering. 2015. No. 2. B/2015. Politechnika Krakowska: Krakow. Poland. Pp. 83-89.

16. Kwon, D.K., Kareem, A., Stansel, R., Ellington, B.R. Wind load factors for dynamically sensitive structures with uncertainties. Engineering Structures. 2016. Vol. 103. Pp. 53-62.

17. Baravalle, M., Kohler, J. On the probabilistic representation of the wind climate for calibration of structural design standards. Structural Safety. 2018. Vol. 70. Pp. 115-127.

18. Pichugin, S. Probabilistic description of ground snow loads for Ukraine. Snow Engineering. Recent Advanced and Developments. Rotterdam: A.A.Balkema. 2000. Pp. 251-256.

19. Ledovskoy, I.V. Trend snegovyih nagruzok na grunt pod vliyaniem globalnogo potepleniya klimata [Trend of the snow loading on ground under influence of the global warming of climate] Stroitelnaya mehanika i raschet sooruzheniy. 2012. No. 4. Pp. 71-76.

20. Fikke, S.M., Nygaard, B.E.K. Assessing environmental actions from modern meteorology. Technical Transactions. Civil Engineering. 2015. No. 2. B/2015. Politechnika Kracowska: Krakow. Poland. Pp. 405-415.

21. Kudzys, A.P. Otsenka nadezhnosti zhelezobetonnyih konstruktsiy [Reliability estimation of reinforced concrete structures]. Vilnius: "Mokslas". 1985. 156 p.

22. Cornell, C.A. Bounds on the reliability of structural systems. Journal of the Structural Division. ASCE. 1967. Vol. 93. No. ST. Pp. 171-200.

23. Liang, J., Ding, Z., Li, J. A probabilistic analyzed method for concrete fatigue life. Probabilistic Engineering Mechanics. 2017. Vol. 49. Pp. 1-32.

24. Moiseenko, R.P. Nachalnaya nadezhnost zhelezobetonnyih elementov [Initial reliability of reinforce-concrete elements] Stroitelnaya mehanika i raschet sooruzheniy. 2016. No. 2(265). Pp. 66-70.

25. Gromatskiy, V.A. O nekotoryih metodah veroyatnostnogo rascheta nadezhnosti konstruktsiy [About some methods of probabilistic calculation of reliability of constructions] Stroitelnaya mehanika i raschet sooruzheniy. 2016. No. 2(265). Pp. 61-65.

26. Krakovskiy, M.B. Povyishenie nadezhnosti zhelezobetonnyih konstruktsiy pri raschetah po neskolkim programmnyim kompleksam [Increase of reliability of reinforce-concrete constructions at calculations on a few programmatic complexes]. Stroitelnaya mehanika i raschet sooruzheniy. 2016. No. 5(268). Pp. 45-49.

27. Rayzer, V.D. Analiz nadezhnosti konstruktsiy pri iznose nesuschih konstruktsiy [Analysis of structure reliability at the wear of load carrying structures]. Stroitelnaya mehanika i raschet sooruzheniy. 2013. No. 6(251). Pp. 16-20.

28. Perelmuter, A.V., Pichugin, S.F. Ob otsenke uyazvimosti stroitelnyih konstruktsiy [About an estimation of building structures vulnerability]. Magazine of Civil Engineering. 2014. No. 5(49). Pp. 5-14. doi: 10.5862/MCE.49.1

13. Pichugin S. Analysis of bridge crane loads on industrial buildings // XLIV Konferencja Naukowa KILiW PAN i KN PZITB. Tom VII. Poznan-Krynica. 1998. Pp. 171-178

14. Pichugin S. Probabilistic analysis on wind load and reliability of structures // Proc. of the 2 EACWE. Vol. 2. Genova. Italy. 1997. Pp. 1883-1890.

15. Kinash R.I., Huk J.S. Geographical altitude coefficient utilisation for wind gust loading determination in the mountainous terrain of the Transcarpathian region // Technical Transactions, Civil Engineering. 2015. № 2. B/2015. Politechnika Krakowska: Krakow. Poland. Pp. 83-89.

16. Kwon D.K., Kareem A., Stansel R., Ellington B.R. Wind load factors for dynamically sensitive structures with uncertainties // Engineering Structures. 2016. Vol. 103. Pp. 53-62.

17. Baravalle M., Köhler J. On the probabilistic representation of the wind climate for calibration of structural design standards // Structural Safety. 2018. Vol. 70. Pp. 115-127.

18. Pichugin S. Probabilistic description of ground snow loads for Ukraine // Snow Engineering. Recent Advanced and Developments. Rotterdam: A.A.Balkema., 2000. Pp. 251-256.

19. Ледовской И.В. Тренд снеговых нагрузок на грунт под влиянием глобального потепления климата // Стр. мех и расчет сооруж. 2012. № 4. С. 71-76.

20. Fikke S.M., Nygaard B.E.K. Assessing environmental actions from modern meteorology // Technical Transactions. Civil Engineering. 2015. № 2. B/2015. Politechnika Kracowska: Krakow. Poland. Pp. 405-415.

21. Кудзис А.П. Оценка надежности железобетонных конструкций. Вильнюс: Мокслас, 1985. 156 с.

22. Cornell C.A. Bounds on the reliability of structural systems // Journal of the Structural Division, ASCE. 1967. Vol. 93, № ST. Pp. 171-200.

23. Liang J., Ding Z., Li J. A probabilistic analyzed method for concrete fatigue life // Probabilistic Engineering Mechanics. 2017. Vol. 49. Pp. 1-32.

24. Моисеенко Р.П. Начальная надежность железобетонных элементов // Стр. мех и расчет сооруж. 2016. № 2(265). С. 66-70.

25. Громацкий В.А. О некоторых методах вероятностного расчета надежности конструкций // Стр. мех и расчет сооруж. 2016. № 2(265). С. 61-65.

26. Краковский М.Б. Повышение надежности железобетонных конструкций при расчетах по нескольким программным комплексам // Стр. мех и расчет сооруж. 2016. № 5(268). С. 45-49.

27. Райзер В.Д. Анализ надежности конструкций при износе несущих конструкций // Стр. мех и расчет сооруж. 2013. № 6(251). С. 16-20.

28. Перельмутер А.В., Пичугин С.Ф. Об оценке уязвимости строительных конструкций // Инженерно-строительный журнал. 2014. № 5(49). С. 5-14.

29. Климович В.Н., Парьев К.В.Надежность конструкции защиты от размывов дна вблизи платформы «Приречная» // Инженерно-строительный журнал. 2017. № 7(75). С. 52-65.

30. Индейкин А.В., Чижов С.В., Шестакова Е.Б., Антонюк

A.А., Кулагин Н.И., Смирнов В.Н., Карпов В.В., Голицынский Д.М. Приближенные методы оценки надежности балочных пролетных строений железнодорожных мостов // Инженерно-строительный журнал. 2017. № 7(75). С. 150-160.

31. Исхаков Ш.Ш., Ковалев Ф.Е., Васкевич В.М., Рыжиков

B.Ю. Оценка надёжности эксплуатации зданий и сооружений по методикам возникновения риска их неработоспособных состояний // Инженерно-строительный журнал. 2012. № 7(33). С. 76-88.

29. Klimovich, V.I. Pariev, K.V. Reliability of scour protection design near the platform «Prirazlomnaja». Magazine of Civil Engineering. 2017. No. 7(75). Pp. 52-65. doi: 10.18720/M CE. 75.5

30. Indeykin, A.V., Chizhov, S.V., Shestakova, E.B., Antonyuk, A. A., Kulagin, N.J., Smirnov, N.V., Karpov, V.V., Golitsynsky, D.M. Approximated method of estimation of the reliability of framed railway structures of railway bridges. Magazine of Civil Engineering. 2017. No. 7(75). Pp. 150160. doi: 10.18720/MCE.75.15.

31. Ishakov, Sh.Sh., Kovalev, F.E., Vaskevich, V.M., Ryizhikov, V.Yu. Otsenka nadezhnosti ekspluatatsii zdaniy i sooruzheniy po metodikam vozniknoveniya riska ih nerabotosposobnyih sostoyaniy [Estimatiing reliability of buildings and structures according to the methods of the risk of unserviceability]. Magazine of Civil Engineering. 2012. No. 7(33). Pp. 76-88.

32. Pryadko, Yu.N., Muschanov, V.F., Bartolo, H., Vatin, N.I., Rudneva, I.N. Improved numerical methods in reliability analysis of suspension roof joints. Magazine of Civil Engineering. 2016. No. 5(65). Pp. 27-41. doi: 10.5862/MCE.65.3

33. Kowal, Z. Influence of the bearing capacity and reliability of connections and nodes on the bearing capacity and reliability of complex bar structures. Proc. of the XIII Int. Conf. on Metal Structures (ICMS2016). Rotterdam: Balkema. Pp. 202-203.

34. Mori, Y., Kato, T., Murai, K. Probabilistic models of combinations of stochastic loads for limit state design. Structural Safety. 2003. Vol. 25. No. 1. Pp. 69-97.

35. Dimitrov, N. Probabilistic model for multi-axial dynamic load for wind turbines. Engineering Structures. 2016. Vol. 117. Pp. 239-249.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

36. Derek, L., Kozak, A.B.L. Reliability of steel roof structures under snow load. Structural Safety. 2015. Vol. 54. Pp. 46-56.

37. Thai, H.T., Uy, B., Kang, W.H., Hikc, S. System reliability evaluation of steel frames with semi-rigid connections. Journal of Constructional Steel Research. 2016. Vol. 121. Pp. 29-39.

38. Tao, Y.R., Cao, L., Cheng, G.Q., Huang, Zh.H. Safety analysis of structures with probability and evidence theory. International Journal of Steel Structures. 2016. Vol. 16. No. 2. Pp. 289-298.

39. Dehtyar, A.S. O vozmozhnosti progressiruyuschego razrusheniya [About possibility of progressive failure]. Stroitelnaya mehanika i raschet sooruzheniy. 2017. No. 6 (275). Pp. 2-6.

32. Прядко Ю.Н., Мущанов В.Ф., Бартоло Х., Ватин Н.И., Руднева И.Н. Усовершенствование числовых методов и расчета узлов висячих конструкций // Инженерно-строительный журнал. 2016. № 5(65). С. 27-41.

33. Kowal Z. Influence of the bearing capacity and reliability of connections and nodes on the bearing capacity and reliability of complex bar structures // Proc. of the XIII Int. Conf. on Metal Structures (ICMS2016). Rotterdam: Balkema. Pp. 202-203.

34. Mori Y., Kato T., Murai K. Probabilistic models of combinations of stochastic loads for limit state design // Structural Safety. 2003. Vol. 25. № 1. Pp. 69-97.

35. Dimitrov N. Probabilistic model for multi-axial dynamic load for wind turbines // Engineering Structures. 2016. Vol. 117. Pp. 239-249.

36. Derek L., Kozak A.B.L. Reliability of steel roof structures under snow load // Structural Safety. 2015. Vol. 54. Pp. 46-56.

37. Thai H.T., Uy B., Kang W.H., Hikc S. System reliability evaluation of steel frames with semi-rigid connections // Journal of Constructional Steel Research. 2016. Vol. 121. Pp. 29-39.

38. Tao Y.R., Cao L., Cheng G.Q., Huang Zh. H. Safety analysis of structures with probability and evidence theory // International Journal of Steel Structures. 2016. Vol. 16. № 2. Pp. 289-298.

39. Дехтярь А.С. О возможности прогрессирующего разрушения // Стр. мех и расчет сооруж. 2017. № 6 (275). С. 2-6.

Sergii Pichugin, Сергей Федорович Пичугин,

+3(8050)591-77-28; pichugin.sf@gmail.com +3(8050)591-77-28;

эл. почта: pichugin.sf@gmail.com

© Pichugin, S.F., 2018

i Надоели баннеры? Вы всегда можете отключить рекламу.