Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Мастаченко, В. Н. Надежность моделирования строительных конструкций монография

.pdf
Скачиваний:
19
Добавлен:
19.10.2023
Размер:
3.51 Mб
Скачать

На основе формул (58), (59) с учетом результатов, получаемых по формулам (129) — (131), можноопределить коэффициенты корреляции индикаторов. I\, h, h-

w2 (Е)

(132)

Sj3 = г 2з — 0 .

s (Е)s U2)

 

Для иллюстрации методики расчета ряд

значений

коэффициентов вариации принимается на основе прак­ тических соображений (в рассматриваемом опыте они не определялись). Коэффициент вариации испытатель­ ной нагрузки принимается равным w (Р) =0,03 исходя из условия, что при использовании обычных маномет­ ров, но в^лабораторных условиях, возможна погреш­ ность порядка ± (9 —10) %.

Коэффициент вариации модуля упругости стали марки СтЗ принят равным w (£) = 0,015. Коэффициент вариации геометрических размеров принимается рав­ ным: w(L) =0,02.

Коэффициент вариации модуля деформаций каната

определяется

на

основе проведенных

испытаний. Под­

счеты дают

при

п=12: Ек= 1,7■ 107

Па; s(£ K)= 7,6X

ХЮ5 Па; w(EK) =0,0445.

в той же последо­

Выполняя

вычисления для w'(E)

вательности, как и в предыдущем примере при а*=0,95

и п = 12, можно найти 1 ^ 0 ,5 5 ; /^^0,08;

=0,60;

w' (Ея) =0,071.

k

Коэффициент вариации начального натяжения опре­ деляется возможной погрешностью при натяжении за­ тяжки. Эта погрешность может быть порядка ±10%. Учитывая, что в данном случае усилия натяжения кон­ тролировались по манометру и по относительным удли­ нениям затяжки, можно считать, что предельная погреш­ ность не превышала ± (4—5) %; это дает основание при­ нять небольшой коэффициент вариации ш(ао) =0,015. Предполагается, что законы распределения всех рас­ смотренных величин нормальны.

Численные значения Sj, / = 1, 2, 3 при учете получен­ ных и принятых значений коэффициентов вариации оп­ ределяются по формулам (129) —(131):

s (Е) = (0,032 + 0,015а+

4-0,022)1/2 = 0,0525;

s (/а) = (0,0152 + 0,0712)1/2 =

0,0725; s( / 3) = 0,015.

70

Коэффициент корреляции определяется по формуле

(132):

Т12

0,0152

0,059.

0,0525-0,0725

Коэффициент г12

мал, поэтому

дальнейший подсчет

ведется как для некоррелированных индикаторов. При значении г>0,10-1-0,20 учет корреляции индикаторов обязателен. В. этом случае расчеты следует выполнять по формуле (67).

Заданному значению Н = 0,95 отвечает погрешность подобия А, определяемая из уравнения (49)

Подбором находим А «0,15.

Погрешность в оценке напряжений для оригинала оп­ ределяется по формулам (103). Для напряжений полу­ чаем возможные относительные погрешности в диапазо­ не от +13 до —17,5%; для перемёщений (см. формулу

(125) — от + 7 до —8,5%.

Результаты оценки точности, приведенные в § 13 и 14, хотя и базируются на ряде условных статистических характеристик свойств материалов, геометрии и нагру­ зок, дают возможность сказать, что моделирование в уп­ ругой области при надлежащем контроле параметров эксперимента может проводиться с использованием од­ ной модели при условии, что погрешность результатов порядка 10—20% является приемлемой.

Анализируя числовые данные, использованные в при­ мерах, нетрудно заключить, что в области неупругих де­ формаций, когда коэффициенты вариации механических характеристик резко возрастают, погрешность конечных результатов будет существенно выше. Из-за недостатка места примеры по оценке точности результатов модели­ рования в неупругой области работы конструкций не рассматриваются, но дается иллюстрация планирования испытаний.

§15. ПРИМЕР ОПРЕДЕЛЕНИЯ ЧИСЛА МОДЕЛЕЙ

ИОБЪЕМА ВЫБОРОК ПРИ ИССЛЕДОВАНИИ

МОДЕЛИ СТРУКТУРЫ

Проведем планирование эксперимента для случая моделирования структуры, рассмотренной в § 13.

71

Пусть требуемый уровень вероятности подобия задан равным Н = 0,95, а предельная относительная погреш­ ность результатов оценки напряжений (для оригинала) равна ± 15%- Предварительная информация о свойст­ вах материала и прочих характеристиках по § 13.

Планирование проводится в соответствии с указани­

ями §11.

На основе формул (103) находим при а/*.,= 1 (см. структуру индикаторов в § 13)

откуда Д=0,13.

Предположим, что коэффициент вариации известен лишь по результатам испытания семи образцов и равен 0,067. При этом значении найдем s(I) =0,078.

По формуле (109) при Л=0,13 и s(/) =0,078 можно найти, что при р—2 Н = 0,981 >0,95. Принимается р = 2.

Далее по указаниям § 13 находим: А!н =0,0327; Н '= =0,95; 6= 2,34; 6 '= 1,96; Л.б =0,162; X- =0,194.

По формуле (115) при |а*| = 1 получим Х ^ ^ =0,195, а по формуле (116) при р = 2 найдем Х^р(е| ) =0,277.

Пользуясь графиком (см. рис. 9), определяем, что

необходимо испытать не менее 38 образцов.

поэтому

Реально было испытано лишь 7 образцов,

испытания необхвдимо продолжить, уточнить

w(E) и

снова проверить правильность планирования

экспери­

мента.

 

§ 16. ПРИМЕР ПЛАНИРОВАНИЯ ИСПЫТАНИЯ БЕТОННОЙ КОНСТРУКЦИИ НА ДЕЙСТВИЕ СОБСТВЕННОГО ВЕСА

В ЦЕНТРОБЕЖНОМ ПОЛЕ

Некоторая бетонная конструкция должна быть иссле­ дована на действие собственного веса с помощью моде­ лирования в центробежном поле. Требуемые параметры моделирования: [Н*] =0,90; [А] =0,08; а* = 0,90. Нуж­ но определить число моделей и объемы выборок для образцов материала моделей .

Система индикаторов, обеспечивающих подобие мо­ дели и оригинала, имеет структуру:

л

Л Л Л

л л

 

gr PrLr

 

Д=

Л

= 1; /2—гг —1.

(133)

<у,

72

Здесь

л

л

масштабы напряжений

<ту, ег— соответственно

л

л

о и деформаций е;

л

масштабы ускорений

gn рг и Ьг— соответственно

свободного падения, плотности мате­ риалов и геометрических размеров.

В первом предположении принимается, что все вели­ чины, входящие в индикаторы (133), некоррелированы

между собой. Генеральные средние для всех величин А Л Л

оригинала известны и равны gH, рн, Гн. Предполагается, что q-*-оо, поэтому коэффициенты вариации для указан­ ных средних равны нулю.

На основе предварительных исследований ускорений центрифуги, материала модели и точности изготовления моделей было установлено, что щ(£м)=0,03; ю(рм) — =0,02; w(LM) =0,01; w(oM) = w (E м) =0,08; ш(ем) =0,08.

Определены также средние выборочные для модели: gm

рм, Гм, 0 М, ем. Объем каждой выборки п =\0 . Все числа условны.

Планирование эксперимента. Прежде всего на основе уравнения (109) необходимо установить число моделей, а для этого требуется вначале определить средние квад­ ратические отклонения величин индикаторов 1\ и /2. По формуле (77) при р = 1 и q^-oo имеем:

s (/р = [ю2 (gM) + ш2 (рм) + ш2 (LM) +

w2 (ам)]1/2 =

=

[0,032 +

0,022 +

0,012 +

0,082]1/2 = 0,0872;

 

 

s (/2) =

w(sM) =

0,08.

 

Примем р — 6,

тогда найдем

 

 

Н =

22Ф V 6-0,08

/ 1^6-0,08

= 0,975-0,986 > 0,90.

 

0,0872

V

0,08

 

Значения Н = 0,90, Н! =

0,975 и Н2=0,986 есть неко­

торые оценки надежности, которые будут действительно обеспечены, если коэффициенты вариации не случайны. Если необходимо обеспечить Н ^0,90 при а* = 0,90, то при р = 6 абсолютная погрешность может быть допуще­ на до 0,96—0,90=0,06, при этом относительная погреш­ ность будет равна: Лн =0,06:0,96=0,0625. Распреде­

лим эту погрешность между Hi и Н2 следующим обра­ зом: 7,Hi =0,0325, ЯНг =0,0290. Условие (ПО) при этом

удовлетворяется,

6-72

73

Далее рассмотрим планирование по каждому Hj, / = = 1,2 отдельно:

/ = 1; Hi =

0,975;

б ,

= Ф - 1

/ 0 ,9 7 5 \

( — 1 = 2 , 2 4 ;

Hj =

(1 — 0,0325) 0,975 = 0,944;

б] ■ ф- 1 0,944

1,91.

 

 

,

2

 

По формулам

§ 11

находим

Аст=0,147, А.-=0,172,

Хш(а)=0,17. Здесь

вместо

w(X)

введено обозначение

до (а), поскольку именно величина о обладает наиболь­

шей изменчивостью. В соответствии

с формулой (118)

имеем: Ха =0,416—1,416 А^. Зададим

А^ =0,0363 и най­

дем из выражения (119) отношение

 

а

ta.

о .о з б з

 

 

=

0.08 = 0,453.

ir„

По табл. 12.2 из работы [41] при а*=0,90 и —гг —

Vn

=0,453 находим п = 15. Если А^ =0,0363, то Аа =0,369.

Представим выражение (120) в виде условия

1 %1—а*

(1+А0)2> — •

Подставив сюда полученное значение А0, по табл. 5 из

работы [44]

устанавливаем

при а* =

0,90, что записан­

ное условие

удовлетворяется

при

[=13, т. е. «=14;

 

 

 

0 986\

 

 

 

 

(- у -

=2,44;

н 2 = О — 0,029) 0,986 =

0,956;

б2= 2,02,

Далее все действия повторяются

в '

соответствии с

предыдущим. По формуле

(118)

Аа =0,515—1,515 А^

Если принять такое значение А^, чтобы получить п = 15

(как в предыдущем случае), то можно найти, что для обеспечения соответствующего значения Ас необходимо

иметь п=10. Разумеется, результат п = 15 удовлетворит всем требованиям, но целесообразно выяснить, нет ли возможности уменьшить число образцов. Для этого пе­ рераспределим погрешности: Ан =0,0345, Ан =0,0270,

74

Повторяя все выкладки, можно найти: / = 1, n^=V2, па = 13; /= 2 , = 12, и0 = 12. Следовательно, число об­

разцов можно сократить, по крайней мере, до 13.

В* начале примера было указано, что_ выборки для

предварительного определения значений 1 и s были при­ няты по п=10. Как показал расчет, это число образцов для определения механических характеристик бетона моделей недостаточно, поэтому объем выборок следует довести до 13 образцов и затем уточнить расчет. Общее

число образцов при испытаниях будет равно:

N — pn=

= 6-13=78.

 

числа

Интересно проследить, как влияет изменение

моделей на требуемое число образцов. Если

принять

р — 5,

то при тех же исходных данных можно

было бы

найти

yv=5-39=195. Дальнейшее уменьшение

числа

моделей при любом сколько угодно большом п не обес­ печивает удовлетворения условия Н^0,90. Поскольку число образцов и число моделей связаны между собой с точки зрения оценки Н при заданном а*, то в преды­ дущем случае, видимо, можно оставить прежний объем выборок, но увеличить число моделей до 7, если такое увеличение не представит больших трудностей, чем уве­ личение числа образцов.

§ 17. ПРИМЕР ПЛАНИРОВАНИЯ ИССЛЕДОВАНИЯ НА МОДЕЛЯХ ЖЕЛЕЗОБЕТОННЫХ КОНСТРУКЦИЙ

Железобетонные балки, предназначаемые, для серий­ ного изготовления, предполагается исследовать на моде­ лях. Оригинал запроектирован из бетона марки 300 с ра­

бочей арматурой класса A-I. Нагрузки

на балки — со­

средоточенные, собственный вес балок

по сравнению

с внешними нагрузками мал. В соответствии с табл. 29*

СНиП П-В.1-62* имеем:

Дпрн = R npH =2100 Па; Дрн =

= Яр =210 Па. Среднее

значение предела

текучести

# а= 30 000 Па.

 

несущей

Основная

задача моделирования — оценка

способности,

геометрический масштаб модели

Lr— 5/1.

Путем подбора состава бетона на образцах 2X2X8 см3

л

л

 

необходимо обеспечить # пр.м=2100 Па,

Яр.м=210 Па

при одинаковых предельных деформациях

с бетоном

оригинала. Арматура моделей — класса

А-1.

Для оцен-

6*

75

ки X и s (А) взяты выборки; для бетона моделей /г=50,

для арматуры моделей п — 40. Коэффициент

вариации

испытательной нагрузки равен:

w (Р) =0,02

на основе

весьма большого п.

для бетона

моделей

Коэффициенты вариаций:

w (#пр.м)=йи(еб.м)=0.05; длй

арматуры

моделей

ш(/?а.м) =йУ(еа.м) ==0,03. Коэффициент вариации для геометрических размеров моделей, определенный на до­ статочно большой выборке, равен w (L) =0,206 (все цифры условные).

Структура индикаторов подобия для данных условий:

л

л л

, л о

 

h — г К

'

’, =

Pr R7rL7 2 -

1'.

 

 

А

Л Л ,

_ 1

(134)

 

 

/ 3 =

6 а г 8 бг

 

 

 

 

 

 

Л

Л Л

 

 

сосредоточенных

Обозначения:

р г = р нр м- 1- - масштаб

л

Л

Л

л

л

 

нагрузок; Р§гRupm/Rupm—^рн/^рм — масштаб сопро-

лл л

тивлений бетона; Rar=RanlRaM — масштаб сопротивле-

лл

ния арматуры; г&г, г^г— масштабы предельных дефор­ маций арматуры бетона.

Анализируя структуру индикаторов подобия, прини­ маем, что коррелированы лишь индикаторы / 1 и /2 (при наличии опытных данных следует учитывать корреля­ цию с третьим индикатором).

Требуется определить число моделей для испытаний и проверить достаточность числа образцов в выборках для определения механических свойств бетона и армату­

ры моделей при условии, что необходимая

надежность

моделирования должна быть не менее

Н= 0,94 при

погрешности подобия Д=0,05.

 

 

 

Планирование эксперимента. Прежде всего определя­

ются D(Ij)

и г12 по формулам (77) и (58):

 

 

D (А) =

(Р) +

ш2 (Рам) + 4ш2 (L) = 30 -10-*;

sv = 0,0548

D (/,) =

(Рам) +

(Рбм) =

34-10-4;

s, =

0,0584

D (/,) = ш2 (8ам) +

(8бм) =

34-10-4;

s3 =

0,0584

ri2 = (— 1) (— 1)

(7?ам)/sis2 =

0,282.

Расчеты проводятся по формулам

(61) —(67). С уче­

том полученных значений

Sj, / = 1, 2, 3 и г12 по форму-

76

Лам (65) и (66)

можно найти п =

8,05 / 1+7,55 /2 — 15,60;

t2= —10,75 /i+10,10 /г+0,65. По

этим уравнениям, за­

давая /= 1 ± А ,

можно получить

значения

ila и 12адля

а — 1,

2, 3,

4. Эти значения представлены в табл. 1.

 

 

 

 

 

 

 

 

Т а б л и ц а 1

V

 

 

а

V

 

а

 

 

1

2

3

1

2

3

 

 

4

4

h

1,05

1,05

0,95

0,95 h

+0,80

0,00

—0,80

0,00

и

1,05

0,95

0,95

1,05

12

0,00

—1,00

0,00

+

 

О

Можно было бы ограничиться подсчетами ila

и

i2a

при а=1,2. Подсчет при а = 3,4

целесообразно

прово­

дить лишь для контроля: iu a = —+ a+2.

 

 

По формулам

(61— (63) с учетом данных табл. 1 по­

лучены результаты, сведенные в табл. 2.

 

2*1

 

 

 

 

Т а б л и ц а

а

 

Р = 1

 

р = 9

р= 1б

 

йа

" а

а и

 

 

 

 

 

 

 

1

0,625

1,25

—0,80

1,875

2,500

 

2

0,625

0,80

—1,25

1,875

2,500

 

 

П р и м е ч а н и е .

Значения

Ьа и аа

не зависят от р.

 

 

По данным табл. 2 по формуле (67) с использовани­ ем таблиц [38] получены оценки надежности: р = 1, Н = = 0,1185; р = 9, Н = 0,885; р = 16, Н = 0,974> [0,90]. Ме­ тодика подсчета проводится для последнего случая. Число индикаторов т — 3, из них два коррелированы, по­ этому

Н = {1 — 2 [Г (2,5;

1,25)

Т (2,50; — 0,8) + Г (2,50; 0,80) —

Т (2,50; — 1,25)]} 2 Ф (

°5 ) = 0,974-0,9994 = 0,974.

1

\

0,0584 /

Функции Т(Ae; b j определяются в соответствии с ука­ заниями к таблицам, приведенным в [38].

77

Для испытаний принимаем р=16. Теперь необходимо проверить достаточность выборок, взятых для определе­

ния свойств бетона и арматуры.

а* ==0,90

по

Число бетонных образцов и = 50. При

табл.

12.2 из [41]

имеем ta, / V n = 0,237.

По формулам

§ 11

получаем =

0,012. Используя табл. 5, из [44]

при

/= 5 0 —1=49 и <х* = 0,90 определяем Хо = 0,152. Погреш­

ность в определении да(^пр) равна:

KW(r) — 0,194. На­

ибольшее возможное при данных

.условиях значение

w'(Rnv) определяется по формуле (121): w'(Rnp) =0,0597. Аналогично выполняются подсчеты для определения

w'(Ra) =0,0362. После уточнения'

значений

да(/?Пр) и

w(Ra)

произведен

перерасчет для

Sf, новые

значения

равны:

sj=0,0585;

s2= s 3=0,0698;

ri2=0,32.

Значения

j'i и t2 для а = 1 , 2, 3, 4 соответственно равны: + 0,70; 0,05; —0,70; -0,05; t2: —0,10; —0,95; +0,10; +0,95. Да­ лее подсчитаны: fti = 0,617; bi= 1,175; ai = —0,576; h2= =0,512; &2=0,955; a2= —1,570. С учетом этих значений

при р — 16

получаем

]//7 /ii= 2,468;

]/рЛ 2=

2,048 и

Н =

{1 — 2 (2,468; 1,175) — Т (2,468; — 0,576) +

 

+ Т (2,048; 0,955) -

Т (2,048; -

1,570)]} 2 Ф ( ^ ^ )

=

 

= 0,946-0,996 =

[0,94].

 

 

Следовательно, выборки для определения свойств бе­ тона и арматуры моделей достаточны.

Для уменьшения вычислений при планировании ис­ пытаний железобетонных конструкций, условия подобия которых могут быть выражены индикаторами, приведен­ ными в этом параграфе, можно использовать графики (см. Труды МИИТ, 1971, вып. 369), позволяющие при за­ данном числе образцов материалов и известных выбо­ рочных коэффициентах вариации прочности бетона оп­ ределять необходимое число моделей исходя из задан­ ных значений Н и А.

ЗАКЛЮЧЕНИЕ

Результаты, представленные в этой работе, позволя­ ют высказать суждение о том, что теория моделирования, построенная на идее сохранения подобия лишь по пер­ вым начальным моментам с оценкой вероятности подо­ бия в каждом отдельном случае, дает практическую воз­ можность рассмотреть ряд задач, для которых ранее не было удовлетворительного решения. Среди них можно назвать задачи по оценке точности и надежности резуль­ татов исследования на моделях и задачи по планирова­ нию эксперимента, строящегося на принципах подобия. Последняя группа задач решается на основе обеспечения заданных значений точности и надежности результатов. Как следствие из этих решений появляется возможность ставить и решать-задачи по оптимальному планированию эксперимента на моделях исходя из экономических пос­ ледствий решений, принимаемых на основе результатов исследований, содержащих погрешности. Все эти задачи рассмотрены применительно к случаю моделирования математического ожидания поведения натурного объек­ та и в определенной степени проиллюстрированы на при­ мерах.

Приведенный в работе формульный аппарат дает воз­ можность решать и ряд других задач, например прово­ дить сопоставление осредненных результатов испытаний ряда моделей и оригиналов при любом их числе с уче­ том случайных свойств; давать оценку точности резуль­ татов испытаний натурных конструкций; оценивать адек­ ватность расчетных моделей конструкций их реальной работе и устанавливать систематическую погрешность расчетных моделей [24, 25]. Хотя эта группа задач и не иллюстрирована на примерах из-за ограниченности ме­ ста, их решение не должно представлять больших труд­ ностей.

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

79

Соседние файлы в папке книги из ГПНТБ