atapin
.pdfЗатем поочередно приложим единичные силы, заменяющие X2 и X3, и определим еще две функции изгибающих моментов (рис. 10.26 и 10.27).
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
2a |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
M2 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
M3 |
1 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3a |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
Рис. 10.26. Эпюра М2 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 10.27. Эпюра М3 |
||||||
qa2/2 |
|
|
2qa |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Нам потребуется также функция изги- |
|||||||||||||
3qa2/2 |
|
|
|
|
MP |
|
|
|
бающего момента, возникающего в основной |
||||||||||||||||||
|
|
|
|
|
|
системе |
от действия внешней |
|
нагрузки |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
(рис. 10.28). Значения функций моментов по |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
участкам |
удобно |
свести в |
|
таблицу |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
(табл. 10.1). |
|
|
|
|
|
|
|
|||||||
Рис. 10.28. Эпюра МP |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Т а б л и ц а 10.1 |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Силовой |
|
|
1-й силовой |
|
2-й силовой |
|
|
|
3-й силовой |
|||||||||||||||||
|
|
|
|
|
участок |
|
участок |
|
|
|
участок |
||||||||||||||||
|
фактор |
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
0 x a |
|
0 x 2a |
|
|
|
0 x 3a |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
M1 |
|
|
|
|
|
x |
|
0 |
|
|
|
|
|
|
|
–a |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
M2 |
|
|
|
|
0 |
|
|
|
0 |
|
|
|
|
|
|
|
x |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
M3 |
|
|
|
|
0 |
|
|
|
x |
|
|
|
|
2a |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
MP |
|
|
|
|
|
q |
x2 |
|
|
|
q |
x2 |
|
|
|
|
|
|
|
3 |
qa2 |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
224 |
|
|
|
|
|
|
|
|
|
|
|
|
a |
1 |
|
|
|
|
q |
1000 |
|
|
|
|
|
|
|
|
|
M11 (x) |
x |
|
|
|
M12 (x) |
0 |
|
M13 (x) |
a |
|
|
|
|
|
|||
|
a |
1 |
|
|
|
|
q |
1000 |
|
|
|
|
|
|
|
|
|
M21 (x) |
0 |
|
|
|
M22 (x) |
0 |
|
M23 (x) |
x |
|
|
|
|
|
|||
M11 (x) |
x |
|
|
|
M12 (x) |
0 |
|
M13 (x) |
a |
|
|
|
|
|
|||
M31 (x) |
0 |
|
|
|
M32 (x) |
x |
|
M33 (x) |
2 a |
|
|
|
|
|
|||
M21 (x) |
0 |
|
|
|
M22 (x) |
0 |
|
M23 (x) |
x |
|
|
|
|
|
|||
M31 (x) |
0 |
|
|
x2 |
M32 (x) |
x |
|
x2 |
|
2 a |
3 |
|
|
|
|||
|
|
|
|
M33 (x) |
|
|
|
|
|
||||||||
MP1 (x) |
|
q |
|
MP2 (x) |
|
q |
|
MP3 (x) |
|
q a |
|
|
|||||
|
2 |
|
2 |
2 |
|
|
|||||||||||
MP1a(x) |
|
q |
|
x2 |
MP2 (x) |
|
q |
x2 |
MP3 (x) |
3 |
|
|
|
||||
|
2 |
2 a |
2 |
|
3 qa a |
|
|
||||||||||
11 |
|
|
|
|
|
|
|
|
|
2 |
M13 (x) M13 (x) dx |
|
|||||
|
M11 (x) M11 (x) dx |
|
M12 (x) M12 (x) dx |
|
|
|
|||||||||||
|
0a |
|
|
|
|
|
|
02 a |
|
|
|
|
03 a |
|
|
|
|
11 |
|
M11 (x) M11 (x) dx |
|
M12 (x) M12 (x) dx |
|
|
M13 (x) M13 (x) dx |
|
|||||||||
|
0 |
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
|
|
|
|
a |
|
|
|
|
|
|
2 a |
|
|
|
|
|
3 a |
|
|
|
12 |
a M11 (x) M21 (x) dx |
2 a M12 (x) M22 (x) dx |
|
3 a M13 (x) M23 (x) dx |
|
||||||||||||
12 |
0 M11 (x) M21 (x) dx |
0 |
M12 (x) M22 (x) dx |
0 |
M13 (x) M23 (x) dx |
|
|||||||||||
|
0 |
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
|
|
|
|
a |
|
|
|
|
|
|
2 a |
|
|
|
|
|
3 a |
|
|
|
13 |
a M11 (x) M31 (x) dx |
2 a M12 (x) M32 (x) dx |
|
3 a M13 (x) M33 (x) dx |
|
||||||||||||
13 |
0 M11 (x) M31 (x) dx |
0 |
M12 (x) M32 (x) dx |
0 |
M13 (x) M33 (x) dx |
|
|||||||||||
|
0 |
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
|
|
|
|
a |
|
|
|
|
|
|
2 a |
|
|
|
|
|
3 a |
|
|
|
22 |
a M21 (x) M21 (x) dx |
2 a M22 (x) M22 (x) dx |
|
3 a M23 (x) M23 (x) dx |
|
||||||||||||
22 |
0 M21 (x) M21 (x) dx |
0 |
M22 (x) M22 (x) dx |
0 |
M23 (x) M23 (x) dx |
|
|||||||||||
|
0 |
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
|
|
|
|
a |
|
|
|
|
|
|
2 a |
|
|
|
|
|
3 a |
|
|
|
|
a |
|
|
|
|
|
|
2 a |
|
|
|
|
|
3 a |
|
|
|
23 |
|
M21 (x) M31 (x) dx |
|
|
) |
M32 |
( ) |
|
|
(x) |
(x) dx |
|
|||||
|
|
M22(x) |
(x) dx |
|
|
M23 (x) M33 (x) dx |
|
||||||||||
|
0 |
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
|
|
|
|
a |
|
|
|
|
|
|
a |
|
|
|
|
|
3 a |
|
|
|
|
|
|
|
|
|
|
2 a |
|
|
|
|
|
3 a |
|
|
|
|
33 |
|
M31 (x) M31 (x) dx |
|
M32 (x) M32 (x) dx |
|
|
M33 (x) M33 (x) dx |
|
|||||||||
33 |
0 M31 (x) M31 (x) dx |
0 |
M32 (x) M32 (x) dx |
0 |
M33 (x) M33 (x) dx |
|
|||||||||||
|
0 |
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
|
|
|
21 |
12 |
|
|
31 |
13 |
32 |
23 |
|
|
|
|
|
|
||||
21 |
12 |
|
|
31 |
13 |
32 |
23 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
|
|
|
|
2 a |
|
|
|
|
3 a |
|
|
|||
1P |
a M11 (x) MP1 (x) dx |
2 a M12 (x) MP2 (x) dx |
|
3 a M13 (x) MP3 (x) dx |
|||||||||||||
1P |
0 M11 (x) MP1 (x) dx |
0 |
M12 (x) MP2 (x) dx |
0 |
M13 (x) MP3 (x) dx |
||||||||||||
|
0 |
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
|
|
|
|
a |
|
|
|
|
|
2 a |
|
|
|
|
3 a |
|
|
|||
2P |
a M21 (x) MP1 (x) dx |
2 a M22 (x) MP2 (x) dx |
|
3 a M23 (x) MP3 (x) dx |
|||||||||||||
2P |
0 |
M21 (x) MP1 (x) dx |
0 |
M22 (x) MP2 (x) dx |
0 |
M23 (x) MP3 (x) dx |
|||||||||||
|
0 |
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
|
|
|
|
a |
|
|
|
|
|
2 a |
|
|
|
|
3 a |
|
|
|||
3P |
|
M31 (x) MP1 (x) dx |
|
M32 (x) MP2 (x) dx |
|
|
M33 (x) MP3 (x) dx |
||||||||||
|
0a |
|
|
|
|
|
02 a |
|
|
|
|
03 a |
|
|
|||
|
|
|
|
|
|
|
|
|
|||||||||
3P |
|
M31 (x) MP1 (x) dx |
|
M32 (x) MP2 (x) dx |
|
|
M33 (x) MP3 (x) dx |
||||||||||
|
|
|
|
|
Рис. 10.29. Фрагмент программы |
|
|
||||||||||
|
0 |
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
|
|
225
Переходим теперь к программированию задачи в MathCad. Примем, что размер a = 1 м, а интенсивность распределенной нагрузки q = = 1000 н. Так как изгибная жесткость всех стержней одинакова EJ = const, то она не будет участвовать в нахождении коэффициентов канонических уравнений. На рис. 10.29 приведен фрагмент программы реализующей вычисление коэффициентов.
Для вычисления коэффициентов используются единообразные выражения, которые многократно копируются. Остается только внимательно исправить значения индексов. На следующем фрагменте (рис. 10.30) приведено решение системы канонических уравнений и найденные значения неизвестных реакций.
X1 |
1 |
|
X2 |
|
1 |
X3 |
|
1 |
|
|
|
|
|
|
|
|
|
|
Given |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
11 |
X1 |
12 |
X2 |
13 |
X3 |
1P |
|
|
0 |
|
|
|
|
|
||||
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
||||||||||||||
|
21 |
X1 |
|
22 |
X2 |
|
23 |
X3 |
2P |
|
|
0 |
|
|
||||
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
||||||||||||||
|
|
31 |
X1 |
|
32 |
X2 |
|
33 |
X3 |
|
|
3P |
|
|
0 |
|
||
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
X1 |
|
|
|
|
|
|
|
|
|
|
|
182.143 |
|
|
|
|
|
|
X2 |
|
Find(X1 X2 |
X3) |
42.857 |
|||||||||
|
|
|
|
|
X3 |
|
|
|
|
|
|
|
|
|
|
|
798.214 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 10.30. Фрагмент решения системы уравнений
Теперь нужно вычислить изгибающие моменты, действующие в эквивалентной системе, и построить эпюры. На рис. 10.31 приведен соответствующий фрагмент. Графики функций моментов копируются затем в текстовый редактор Word – меню “Правка”, “специальная вставка”, вставить как “Рисунок (метафайл Windows)”. К выделенному рисунку применяется команда “Разгруппировать”, при этом он преобразуется в векторные примитивы редактора Word. Необходимые элементы (кривые изображающие функции) переносятся затем на расчетную схему. Значения функций в характерных точках сняты с графиков с помощью инструмента “Trace” (рис. 10.31).
226
318,6 Н · м
85,7 Н · м
317,9 Н · м |
403,6 Н · м |
Рис. 10.31. Окончательная эпюра изгибающих моментов
Последний фрагмент на рис. 10.32 относится к проверке правиль-
ности нахождения неизвестных реакций – X1, X2, X3. Здесь с помощью интегралов Мора вычисляются перемещения в эквивалентной системе по направлению отброшенных “лишних связей”. Полученные значения с полным основанием считаем равными нулю, тем более, что эти значения должны быть разделены на изгибную жесткость EJ.
227
M1 (x) |
MP1 (x) |
M11 (x) X1 |
M21 (x) X2 |
|
M31 (x) X3 |
||||||
M2 (x) |
MP2 (x) |
M12 (x) X1 |
M22 (x) X2 |
|
M32 (x) X3 |
||||||
M3 (x) |
MP3 (x) |
M13 (x) X1 |
M23 (x) X2 |
|
M33 (x) X3 |
||||||
x |
0 0.01 |
a |
|
|
|
|
|
|
|
|
|
|
|
|
|
M1(x) |
100 0 |
0.2 |
0.4 |
0.6 |
0.8 |
1 |
|
|
|
|
|
200 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
300 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
x |
0 0.01 |
2 a |
|
|
|
|
|
|
|
|
200 |
|
|
|
|
|
M2(x) |
0 |
0.5 |
1 |
1.5 |
2 |
200 |
|||||
400 |
|
|
|
|
|
|
|
|
x |
|
|
x 0 0.01 3 a
|
|
M3(x) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
1 |
2 |
|
3 |
|
|
|
||
|
|
|
|
50 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
2 a |
|
|
|
3 a |
|
|
11 |
||||
1P |
M11 (x) M1 (x) dx |
|
|
M12 (x) M2 (x) dx |
|
M13 (x) M3 (x) dx |
5.551 |
10 |
|||||
|
0 |
|
|
0 |
|
|
|
|
0 |
|
|
|
|
|
a |
2 a |
|
|
|
3 a |
|
|
12 |
||||
2P |
M21 (x) M1 (x) dx |
|
|
M22 (x) M2 (x) dx |
|
M23 (x) M3 (x) dx |
1.179 |
10 |
|||||
|
0 |
|
|
0 |
|
|
|
|
0 |
|
|
|
|
|
a |
2 a |
|
|
|
3 a |
|
|
10 |
||||
3P |
M31 (x) M1 (x) dx |
|
|
M32 (x) M2 (x) dx |
|
M33 (x) M3 (x) dx |
4.823 |
10 |
|||||
|
0 |
|
|
0 |
|
|
|
|
0 |
|
|
|
|
228
11
РЕШЕНИЕ ЗАДАЧ С ПРИВЕДЕНИЕМ ИСХОДНЫХ ДАННЫХ
К БЕЗРАЗМЕРНОМУ ВИДУ
Как правило, формулы сопротивления материалов предназначены для так называемого ручного счета и результаты решения задачи могут быть получены без применения вычислительной техники.
Однако в некоторых случаях (оптимизация конструкции, ее параметрическое исследование), когда требуется большое количество однотипных вычислений, целесообразно производить расчеты с помощью ЭВМ. При этом приходится разрабатывать специализированные программы, которые потом могут использоваться многократно для расчета подобных конструкций. При написании программ можно использовать как универсальные языки программирования, такие как Паскаль, Фортран, так и универсальные среды – MathCAD, EXCEL и др.
При разработке программ иногда приходится модифицировать расчетные формулы сопротивления материалов, приводить их к безразмерному виду и к виду, удобному для программирования.
В настоящем разделе будет рассмотрена задача параметрического исследования двутавровой стальной балки с применением среды
EXCEL.
11.1. Расчетные формулы
При нахождении опорных реакций для балки удобно использовать уравнения двух моментов, но в отличие от “ручного” счета эти уравнения удобно записывать не относительно опор, а для некоторых двух фиксированных сечений балки, например, ее начала и конца.
Рассмотрим случай нагружения балки, имеющей N силовых участков (рис 11.1). На границах силовых участков могут быть приложены
сосредоточенные усилия и моменты P , M , а также равномерно рас- i i
пределенная нагрузка qi , которая должна заканчиваться на правом конце балки.
229
Ошибка! |
|
Pi |
|
|
y |
RA |
|
qi |
RB |
|
|
Mi |
|
x |
|
|
|
|
xA |
xi |
xB |
xN |
Рис. 11.1. Общий случай нагружения балки
Выражения для моментов от внешней нагрузки для левого M 0P и
правого M NP концов балки будут иметь вид
M P |
N |
|
|
|
|
|
|
|
|
|
xN |
xi |
|
|
|
|
|
M |
i |
P x x |
|
q x |
N |
x |
|
x ; |
|
||||||||
|
|
|
|
|
|
||||||||||||
0 |
|
i |
i |
0 |
i |
i |
2 |
|
|
0 |
|
||||||
i |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
(1) |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
M P |
|
|
|
|
|
|
|
|
|
|
xN |
xi |
|
|
|
||
M |
i |
P x x |
N |
q x |
N |
x |
|
x |
N |
. |
|||||||
|
|
|
|
||||||||||||||
N |
|
i |
i |
|
i |
|
i |
2 |
|
|
|
|
|||||
i |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Записав условия равенства нулю моментов для левого и правого концов балки
R |
A |
x |
A |
|
R x |
B |
M P |
0, |
|
|
|
|
|
|
|
|
B |
|
0 |
|
|
|
|
(2) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
M P |
|
R |
A |
x |
N |
x |
A |
R x |
N |
x |
B |
0, |
|||
|
|
|
|
B |
|
N |
|
получим выражения для определения опорных реакций. Решать систему (2) удобнее всего численно. Так, по методу Крамера будем иметь:
R |
A |
DA |
; |
R |
DB |
; |
D x |
N |
x |
A |
x |
B |
; |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|||||||||||||||||
|
D |
|
B |
D |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
D |
A |
M P |
x |
N |
x |
B |
M P x |
B |
; D |
|
M P |
x |
N |
x |
A |
M P x |
A |
. |
||||||
|
|
0 |
|
|
N |
|
|
B |
|
|
0 |
|
|
N |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
230 |
|
|
|
|
|
|
|
|
|
|
|
При нахождении перерезывающих сил и изгибающих моментов удобно использовать метод перемещений. При этом, определив прогиб балки v из универсального уравнения упругой линии (см. задачу 7, II), силовые факторы можно вычислить из дифференциальных зависимостей
dv |
, M EI |
d |
EI |
d 2v |
, Q |
dM |
. |
(3) |
|
dx |
dx |
dx2 |
dx |
||||||
|
|
|
|
|
Универсальное уравнение упругой линии для балки, показанной на рис 11.1, имеет вид
|
|
v j |
|
v0 |
|
0 x j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
1 |
|
N |
|
|
|
|
|
x j |
|
|
xi |
2 |
|
|
|
|
x j xi |
3 |
|
|
x j |
xi |
4 |
|
, |
(4) |
|||||
|
|
|
|
|
|
|
tij |
Mi |
|
|
|
|
|
|
|
Pi |
|
|
|
|
|
|
qi |
|
|
|
|
|
|
|
|||||
|
|
|
|
EI i |
|
|
2! |
|
|
|
|
3! |
|
|
4! |
|
|
|
|
||||||||||||||||
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
где |
tij |
1 |
|
при x j |
|
xi |
, а v j − прогиб в сечении x j . |
|
|
|
|
(5) |
|||||||||||||||||||||||
0 |
при |
x j |
|
xi |
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
Дифференцируя выражение (4), согласно (3) будем иметь: |
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
j |
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
x j |
1 |
|
|
|
|
x j |
2 |
|
|
x j |
xi |
3 |
|
|
||||||
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
xi |
|
|
xi |
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
tij Mi |
|
|
|
|
|
|
|
|
Pi |
|
|
|
|
|
qi |
|
|
|
|
; |
|
|||
|
|
|
|
|
|
|
EJ i |
0 |
|
|
|
1! |
|
|
|
2! |
|
|
|
3! |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
|
x j xi |
1 |
|
|
x j |
|
xi |
2 |
|
|
|
|
|
|
(11.6) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
M j |
|
|
tij |
|
Mi |
|
Pi |
|
|
|
|
|
|
|
qi |
|
|
|
|
; |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
1! |
|
|
|
|
|
2! |
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
i |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x j |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
|
xi |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
Q j |
|
tij |
Pi |
|
qi |
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
1! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
i |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
231