
книги из ГПНТБ / Мастаченко, В. Н. Надежность моделирования строительных конструкций монография
.pdfНа основе формул (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