
книги из ГПНТБ / Мастаченко, В. Н. Надежность моделирования строительных конструкций монография
.pdfПолагая в этом выражении |
) = ^ х Ау |
с учетом |
равенства (115) найдем, что |
|
|
|
|
(118) |
Погрешность в оценке коэффициента вариации w(X) зависит от погрешности оценки а и ц, причем представля ет интерес погрешность оценки для а сверху, а для р,— снизу. Вследствие этого далее используются односторон ние оценки.
Известно, что если на основе выборки объема п (из нормально распределенной совокупности) определены
средние выборочные X и среднее квадратическое откло нение s случайной величины X, то оценка для генераль
ного среднего снизу будет ц ^ Х 'J/' п . где *а* — кван
тиль (-распределения, а * —доверительная вероятность. |
||
Исходя из этого выражения, можно найти, что |
||
|
а>(Х) |
|
|
(119) |
|
Аналогично из условия |
У п |
|
|
||
а < s |
f |
|
Х?-а« |
||
|
где f = n —1 — число степеней свободы, |
%?_«* — квантиль |
%2-распределения, можно установить, что |
|
1 . |
( 120) |
Значения погрешностей, определяемые выражениями (119) и (120), относятся, строго говоря, лишь к тем ве личинам Х{, которые обладают наибольшей дисперсией, поэтому естественно предположить, что для оценки ста
тистических характеристик других |
величин |
можно |
||
брать несколько меньшее число образцов. |
|
|||
К сожалению, условие (118) не позволяет однознач |
||||
но выбрать |
и Ка исходя |
из |
Известной величины |
|
JУр^юр\). |
Следовательно, |
здесь |
также |
приходится |
распределять |
погрешность между |
и %а, причем при |
одной и той же доверительной вероятности и погрешно
сти для оценки а требуется значительно больше образ цов. Оптимальное решение будет таким, при котором число образцов будет одним и тем же как для оценки генерального среднего р, так и генерального среднего квадратического отклонения а. Чтобы обеспечить это ус-
Рис. 7. Номограмма для решения уравнения %= \'-\-'к"+'к'Х"
ловие, следует |
в уравнении |
(118) принимать |
значи |
тельно меньше |
(примерно 0,1 |
от 1/р ^ м х )), |
чем Ха . |
Приходится решать эту задачу методом подбора; при не котором навыке этот подбор занимает не очень много
времени.
Наибольшее значение коэффициента вариации изу чаемого свойства Xi может быть приближенно опреде лено по формуле
» ' ( * ) = [ ! + * • ( * , ) ] « ( * / ) . |
( 1 2 1 ) |
61
где ‘kw(Xi) определяется по формуле (117), a ifl(A'j)
выборочное значение.
С целью уменьшения несколько утомительной рабо ты по подбору Ян. построена номограмма (рис. 7) для
решения этой задачи.
Рис. 8. Графики для определения объема выбор ки гс <35 в зависимости от относительной погреш ности %Ю(х) при известном выборочном значении w(X) и заданном значении а*
62
Во многих практических случаях число испытывае мых образцов можно сразу определять по графикам, представленным на рис. 8 и 9. Для использования графи ка необходимо найти значение относительной погрешно сти Kw(x) по формуле (116) и, приняв а*, найти число об разцов п.
Рис. У. Графики для определения объема выборки 2 0 |
< я < 8 0 |
в зависимости от относительной погрешности %w(x) |
при из |
вестном выборочном значении w(X) и заданном значении а * |
Изложенный путь решения задачи может оказаться удобным лишь для случая некоррелированных индикато ров и при отсутствии информации о коэффициентах ва риации свойств материала и геометрии модели, а также испытательных нагрузок.
§ 12. ПЛАНИРОВАНИЕ ИСПЫТАНИИ ПРИ КОРРЕЛИРОВАННЫХ ИНДИКАТОРАХ ПОДОБИЯ
В общем случае, когда индикаторы коррелированы и когда предварительная статистическая информация о свойствах модели должна быть получена на основе экс перимента, можно рекомендовать следующий путь.
Прежде всего, исходя из вероятности подобия Н, нуж но выбрать доверительную вероятность а*« Н. Затем, используя таблицы [41] в зависимости от предполагае мого типа распределения исследуемого свойства Х{, нуж-
63
но найти число испытаний п для каждой из величин мо дели, входящей в индикаторы подобия, используя хотя бы априорные сведения о коэффициенте вариации w(X). В частности, при допущении о нормальности распреде ления изучаемого свойства для оценки объема предвари тельных испытаний образцов материалов, замеров гео метрических размеров и пр. может быть использована формула
_ (а* |
( 122) |
|
w ( X ) |
||
■ |
||
Относительная погрешность |
%х в оценке среднего |
значения не должна быть больше предельной погрешно сти конечного результата и не более 0,05. Определив от ношение в левой части равенства (122), по табл. 12.2 из книги [41] при заданной доверительной вероятности а* устанавливается объем предварительных испытаний п.
Если найдено и проведено п испытаний и при этом выборочное значение коэффициента вариации окажет ся большим, чем принятое при использовании формулы (122), то необходимо найти новое число испытаний и продолжить испытания. Этот процесс повторяется до тех пор, пока выборочные значения коэффициентов ва риации и фактически испытанное число образцов не совпадут с определяемыми по формуле (122).
Если же на основе первого подсчета числа п прове дены испытания и w(X) оказались меньше принятого, то испытания заканчиваются.
Далее операции выполняются в такой последователь ности:
1) в соответствии с формулами (121), (116), (118) и
(119) определяются при доверительной вероятности а * « Н коэффициенты вариации w'(X) с учетом принято го для предварительной оценки числа испытаний для наиболее изменчивой величины X (или двух-трех наибо лее изменчивых величин). Однако в любом случае эти коэффициенты не должны превышать коэффициента вариации соответствующей величины определенного для генеральной совокупности (если он известен);
2) с использованием коэффициентов вариации w'(X) для всех величин модели, входящих в индикаторы подобия, определяются при р = 1 средние квадратиче ские отклонения индикаторов и коэффициенты корреля ции [см. формулы (77) (при q-^-oо), (91), (58)];
64
3) по |
формулам (61) — (66) определяются величи |
ны h a , ba |
и аа\ |
4) принимая последовательно р = 1, 2, 3, . . по фор муле (90) с учетом формул (79) или (87) определяют р, при котором Н ^[Н *] при принятом значении А.
В том случае, когда индикаторы некоррелированы, подсчет ha, ba и аа не нужен и в формуле (90) следует
принять k = 0 или использовать непосредственно фор мулу (31).
По известным выборкам для всех величин модели, входящих в условия подобия, устанавливаются необ ходимые статистические характеристики, в том числе и средние арифметические значения. Для соответствую щих величин оригинала математические ожидания или их оценки принимаются на основе литературных дан ных. Знание математических ожиданий для величин оригинала и средних арифметических значений для со ответствующих величин модели позволяет найти оцен ки для множителя преобразования величин оригинала в величины модели. Далее необходимо обеспечить ра венство единице индикаторов подобия для рассматри
ваемой задачи. Это |
обеспечение производится путем |
надлежащего выбора |
I — т множителей преобразова |
ния, здесь I— число различных множителей, входящих во всю систему индикаторов подобия, т —число связей, т. е. число индикаторов подобия. По известным значени ям множителей преобразования уточняются численные значения величин нагрузок или других воздействий, прикладываемых к модели для обеспечения подобия в состояниях модели и оригинала. Результаты испытания моделей оцениваются в соответствии с формулами (103) для оценки точности результатов.
Примеры оценки точности результатов для некото рых случаев реально проведенных испытаний и приме ры планирования испытаний даны в гл. 4.
Глава 4
ПРИМЕРЫ ОЦЕНКИ ТОЧНОСТИ РЕЗУЛЬТАТОВ
ИПЛАНИРОВАНИЯ ИСПЫТАНИЙ НА МОДЕЛЯХ
§13. ОЦЕНКА ТОЧНОСТИ РЕЗУЛЬТАТОВ, ПОЛУЧЕННЫХ НА УПРУГОЙ МОДЕЛИ СТРУКТУРЫ, ИЗГОТОВЛЕННОЙ ИЗ ОДНОГО МАТЕРИАЛА
Модель |
структуры [45] |
изготовлена |
в |
масштабе |
1/70. Нагружение модели |
осуществлялось |
сосредото |
||
ченными грузами по £ = 2 0 Н*. Материал |
модели — по |
|||
листирол. |
Модуль упругости полистирола |
|
для модели |
определялся на семи образцах^ по результатам испыта
ний которых |
можно найти £ = 346 000 |
Па; s(E) = |
= 23 300 Па; |
w(E) = 0,067. Других статистических дан |
|
ных в работе [45] не имеется. |
определялись |
|
Деформации в элементах «структуры» |
в предположении линейного напряженного состояния в упругой области. Условие подобия в этом случае выра жается одним индикатором подобия [32].
л |
л л . л „ |
= 1 |
(123) |
I x = |
Pr E~l L f |
||
Л Л- А Л |
А Л Л |
А |
Л |
при а,.= £ г, ur= L r, где Pr, Er, Lr, crr |
и ur— соответст |
венно множители преобразования для сосредоточенных нагрузок, модуля упругости, геометрических размеров, напряжений и перемещений.
лл л л
Учитывая, что ат= Е г и ur— Lr, индикатор подобия
А л
(123) можно записать, используя значения аг и ur. Интервальная оценка для напряжений и перемеще
ний в оригинале, подобном данной модели, будет иметь вид:
1 + Д 1 |
< 1— |
0“зм 0 |
(124) |
|
|
||
Т |
“Г ^ < “н < ( ] Д Г 1 ) 2 < Ч - |
(125) |
|
1 А |
В индикатор (123) входят значения нагрузок, моду ля упругости, геометрических размеров. Ввиду отсутст-
* Здесь и далее сделан приближенный перевод единиц измерения системы МКГСС в единицы системы СИ: 1 кгс=9,81 Н » 1 0 Н = =0,01 кН; 1 кгс/см2=9,81 Н/см2=9,81 П а » 10 Па,
66
Вия Данных полагаем, что при условиях эксперимента w(P) = 0, w(L) = 0,02.
Модуль упругости определен на семи образцах, его
значение случайно. |
В соответствии с формулой |
(121) |
|||||
найдем w'(E). Пусть а*=0,95. По табл. [44] |
при |
/ = |
|||||
= 7—1 = 6 |
находим Xi_a* =%о.о5 = |
1.64. По табл. 12.2 в |
|||||
книге [41] |
найдем при п — 7 |
|
|
|
|
|
|
|
1д* _ ^0,95 = 0,89. |
|
|
|
|||
|
v~a |
|
|
|
|
|
|
По формуле (120) |
получаем Л.а <'0,92, |
а по формуле |
|||||
(119) ^«£0,06= |
|
|
|
|
|
|
|
По выражению |
(117) |
возможная |
погрешность в |
||||
оценке w(E) равна: |
(£) = |
1,04. |
Далее |
по |
формуле |
(121) находим w'(E) = 0,136.
Среднее квадратическое отклонение индикаторов по
формуле (77) при q-*-оо и р = 1 |
|
|
|
\_ |
|
s = s (/) = [ш2 (Р) + |
w2 (£) + 4ai2 (L)] 2 |
= |
= [0 -f- 0 ,1362 + |
4-0,022] 2 = 0,142. |
|
Предположим, что плотности распределения Е и L |
||
нормальны, тогда Ех (7 )= 0 и по формуле |
(49) |
Н = 2Ф (б) = 2Ф
или
^(■ fb
Здесь Ф-1 (....) —обратная функция Лапласа.
Приняв Н= 0,95, по |
таблицам функции Лапласа |
|
легко найти Д = 0,278. |
|
(124) можно устано |
В соответствии с выражением |
||
вить интервалы для стн: |
|
|
л |
л |
л |
0,78схг (JM< стн < 1,39 ar aM,
л
а по формуле (125) найдем интервал для ии:
0,885 ul™Lr < и и < 1,18 < ам£ г.
Приняв найденные значения ам и им в качестве ма тематических ожиданий, по формулам (103) можно
67
найти с риском ошибки 5%, что погрешность в оценке напряжений не превышает +22% и —39%, а в оценке
перемещений +11,5% и —48%.
Более высокий уровень вероятности Н здесь не был бы оправдан, учитывая, что моделирование в данном случае проводилось с целью проверки метода расчета.
§ 14. ОЦЕНКА ТОЧНОСТИ РЕЗУЛЬТАТОВ ИССЛЕДОВАНИЯ УПРУГОЙ РАБОТЫ
МОДЕЛИ ПРЕДВАРИТЕЛЬНО-НАПРЯЖЕННОЙ ФЕРМЫ, ИЗГОТОВЛЕННОЙ ИЗ ДВУХ МАТЕРИАЛОВ
Схема фермы показана на рис. 10 [46]. Все элемен ты модели изготовлены из стали марки СтЗкп, а затяж-
Рис. 10. Схема опытной фермы
/ — стержневая конструкция; 2 — затяжка
ка — одна прядь каната типа ТК 7X37, диаметр пря ди — 17 мм (модель Ф-2)„ Соединения сварные (элект роды типа Э-42).
Задачей исследования модели была проверка допу стимости расчета стальной предварительно-напряжен ной фермы с затяжкой из стального каната как упругой статически неопределимой системы по методу сил с од ним лишним неизвестным — усилием в затяжке. Модель испытывалась как в упругой, так и пластической об ластях работы. В данном примере будет рассмотрена оценка погрешности результатов лишь в области упру гих деформаций без учета влияния работы соединений на поведение фермы. Оценка точности результатов по зволяет обоснованно найти зону возможных значений
для математического ожидания определяемых в экспе рименте величин.
При испытании модели фермы использовались два монтажных толкающих домкрата : грузоподъемностью по 25 т. Для синхронной работы рабочие цилиндры обо их домкратов были объединены в общую маслосистему
68
трубками, присоединенными к общему манометру. Дав ление домкратов на ферму передавалось в виде сосре доточенных сил в узлах фермы. Натяжение затяжки контролировалось по манометру и по относительным удлинениям. В процессе испытаний регистрировались прогибы и относительные деформации.
При испытании определялись механические свойства стали марки СтЗкп: предел текучести, временное сопро тивление, относительное удлинение. Данных по моду лям упругости нет. Приведены данные о модуле дефор маций каната-затяжки в области рабочих напряжений. Для определения его статистических характеристик ис пользованы результаты испытаний по образцам АЗ, А4, Б1—Б5, В5, Вб, В8—В10. Образец В7 не принимается во внимание, поскольку для этого образца значение мо дуля деформации выходит за пределы случайного от клонения. Общее число образцов для каната-затяж ки — 12.
Индикаторы подобия для данного случая записыва-
ются в виде:
|
|
|
А Л 4А л |
|
|
(126) |
|
|
i |
^ |
p r E 7 l |
= >; |
|
|
|
|
|
|
Л Л |
, |
|
|
(127) |
|
|
h |
= E r E7 |
r = V' |
|
|
|
|
|
|
|
|
|||
|
|
|
А |
|
|
|
(128) |
|
|
|
1 з = в о г ^ 1 , |
|
|
||
где Е, Ек— модули деформаций для стали |
СтЗ |
и ка |
|||||
ната; |
|
|
|
(в затяжке и |
дру |
||
а0 — начальное напряжение |
|||||||
гих элементах фермы). |
непосредственно не |
||||||
Силы тяжести |
в |
эксперименте |
|||||
моделировались. |
Подобие по |
коэффициенту |
Пуассона |
при исследовании линейного напряженного состояния
может быть опущено. |
(77) при р = 1 и q-*~oo |
||||
В соответствии |
с формулой |
||||
для |
приведенных |
индикаторов |
можно записать |
значе |
|
ния |
средних |
квадратических |
отклонений индикаторов |
||
S], / = 1, 2, 3,..., определяемых на основе выборочных ха |
|||||
рактеристик. |
Эти |
значения будут определяться |
выра |
||
жениями: |
|
|
_i_ |
|
|
|
|
|
[ д а 2(Р) + д а 2( £ ) |
|
|
|
|
S (/j) = |
+ 4ш2 (L)] 2 ; |
(129) |
|
|
|
s(I2) = [w*(E) + w*(Ek)]2 ; |
(130) |
||
|
|
|
s(I3) = w(a0). |
(131) |
69