книги / Цифровое моделирование вибраций в радиоконструкциях
..pdfудара. Наличие амортизаторов снижает частоту колеба ний платы и уменьшает ее деформации.
3.3.Расчет вибраций стойки
Вкачестве примера рассмотрим конструкцию стой
ки |
с четырьмя секциями |
(рис. 3.17). Материал стой |
ки |
алюминий: л= 5,42-1010 |
Н/м2, ц= 2,56• 1010 Н/м2, р= |
= |
2,65* 103 кг/м3. Размеры стойки 0,42X0,52X1,2 м. Чис- |
Рис. 3.17. К построению модели стойки.
ло шагов по длине каждого стержня равно десяти, поэтому величины шагов по координатам будут равны соответственно Iixg= 0fi5 м, hXh = 0,04 м и /ix& = 0,03 м. В качестве стержней используются уголки: вертикальные стержни — равнобокий уголок 25x25 мм с размерами эквивалентного прямоугольного сечения
Рис. 3.18. К вычислению масЬы и моментов инерции элемента пере крестия.
h'Vb = h'zb= 0,9• 10-2 |
м, hyb = hzb= 1,36-10- ? |
м, |
длинные горизонтальные |
стержни — неравнобокий уго |
|
лок 25X20 мм с размерами эквивалентного |
сечения |
hyg= h lg= 1,28-10-2 м, h'vg = 0 ,9 - \0~2 м, h'zg= 0,7 -10~2 м,
короткие горизонтальные стержни — неравнобокий уго-
Рис. 3.19. К вычислению моментов инерции частей элемента пере крестия.
лок 25X15 мм с размерами |
эквивалентного |
сечения |
|
hyk= hzH= 1,20-10-2 м> |
/;/уЛ= 0,8- 10 2 м, h'zk = 0.6-10-2 м. |
||
Основные размеры |
элемента |
перекрестия |
показаны |
на рис. 3.17. При вычислении массы элемента перекрес тия, его разбивают на части, изображенные на
рис. 3.18,а,б. Масса элемента перекрестия получается равной
m =2,65 -103 - 4 (30 •25+30 -21 +21 -21 +20 •7,7 + + 21-7,7 + 21-2,7+15-2,7)-10-9=2,37 10-2 кг.
При вычислении моментов инерции элемента перекрес тия его разбивают на те нее части и моменты инерции отдельных частей складывают. Вычисление моментов инерции частей производят по одной из двух формул:
/ = р Щ С 2(а, + а2) + (а31+ йп2) /3], |
(3.15) |
или |
|
/ = р(/г/3)[(а3, + а32) (6i + 62) + (bh + bh) (а, + о2)]. |
(3.16) |
Формула (3.15) получена для случая, когда |
ось, от |
носительно которой вычисляется момент инерции, не пе ресекает части элемента (рис. 3.19,a), a формула (3.16) — для случая, когда ось пересекает часть элемен та (рис. 3.19,6). В табл. 3.6 приведены размеры частей элемента и величина моментов инерции. Суммарные мо менты инерции элемента перекрестия имеют величину
Jg=3,14-10“ ® кг-м2, /6=4,57-10-* кг-м2.
Затем необходимо подготовить численные значения всех коэффициентов, входящих в расчетные соотноше ния. Для этого следует прежде всего выбрать шаг по времени т. Как указывалось выше, шаг по времени вы бирается из условия устойчивости разностного решения. Условие устойчивости будет выполняться, если сумма всех безразмерных коэффициентов, на которые умножа ются слагаемые в любой из формул (1.62), будет меньше четырех, деленных на максимальное значение М (1 + А ). Например, для первого уравнения (1.62) имеем
4А 0< - Æ(î + A)- или А 0 < - ц — Га)
Для второго уравнения получается неравенство
4А 1+2А 2<4/М (1+А ) и т. д. |
(3.17) |
Значения коэффициента массы М и коэффициента А за висят от шага по времени т, поэтому следует решить каждое из неравенств (3.17) относительно шага по вре мени и выбрать наименьшее значение шага. В практиче ских расчетах можно упростить эту задачу, если учесть, что значение коэффициента массы максимально в тех узлах модели, где нет дополнительной массы в виде при-
Т а б л и ц а 3.6
Вычисление моментов инерции перекрестия стержней
Номер |
Формула |
Л, мм |
|
л2, мм |
bi или |
Ьлили |
Момент |
|
части |
|
Ь, мм |
с, мм |
инерции части |
||||
|
|
|
|
|
|
|
/ кг-м* |
|
|
Момент инерции относительно оси короткого стержня |
|
||||||
1 |
(3.2) |
4 |
7,7 |
17,3 |
15 |
15 |
1,19.10-» |
|
2 |
(3.1) |
4 |
15 |
15 |
21 |
5,7 |
0,68 |
.10-» |
3 |
(3.1) |
4 |
5,7 |
15,3 |
21 |
5 |
0,40 |
-10-» |
4 |
(3.2) |
4 |
— 17,3 |
25 |
7 |
13 |
0,81 •10-о |
|
5 |
(3.1) |
4 |
-1 7 ,3 |
25 |
21 |
5 |
0,82•10~ » |
|
6 |
(3.1) |
4 |
5,3 |
15,3 |
2,7 |
5 |
0,05 |
-Ю -о |
7 |
(3.1) |
4 |
8 |
7 |
2,7 |
5,7 |
0,02 |
-Ю -о |
|
Момент инерции относительно оси длинного стержня |
|
||||||
1 |
(3.1) |
4 |
15 |
15 |
25 |
5,7 |
0,85*10-6 |
|
2 |
(3.2) |
4 |
3,7 |
17,3 |
15 |
15 |
1,05.10-6 |
|
3 |
(3.1) |
4 |
3,7 |
17,3 |
21 |
5 |
0,-50-10-6 |
|
4 |
(3.1) |
4 |
3 |
17 |
7,7 |
5,7 |
0,19-10-6 |
|
5 |
(3.1) |
4 |
3,7 |
17,3 |
7.7 |
5 |
0,18 |
-10-6 |
6 |
(3.1) |
4 |
— 17,3 |
20 |
21 |
5 |
0,22 |
-10-6 |
7 |
(3.2) |
4 |
— 17,3 |
20 |
7 |
8 |
0,15-10-6 |
|
|
Момент инерции относительно оси вертикального стер кня |
|||||||
1 |
(3.1) |
4 |
7,7 |
17,3 |
30 |
5,7 |
0,85-10-в |
|
2 |
(3.1) |
4 |
7,7 |
17,3 |
30 |
5,7 |
0,77 |
10-е |
3 |
(3.2) |
4 |
17,3 |
3,7 |
17,3 |
3.7 |
0,78-10-в |
|
4 |
(3.1) |
4 |
— 17,3 |
25 |
20 |
5,7 |
0,83-10-6 |
|
5 |
(3.2) |
4 |
-1 7 ,3 |
25 |
17,3 |
3,7 |
0,92-10-6 |
|
6 |
(3.2) |
4 |
— 17,3 |
20 |
17,3 |
3,7 |
0,26-10-6 |
|
7 |
(3.1) |
4 |
— 17,3 |
20 |
15 |
5,7 |
0,16-Ю -б |
крепленных к стойке блоков. Это максимальное значение равно 1 + Л^. Коэффициент А^ выбирается в зависи
мости от величины коэффициентов А и при учете потерь энергии в алюминии на низких частотах составляет при мерно одну десятую от максимального значения коэффи циентов А в уравнениях (1.62). Поэтому, если принять максимальный коэффициент А несколько меньшим еди-
104
ницы, например, равным 0,8, то условие устойчивости будет выполняться.
Величина коэффициента А6 будет максимальна для короткого стержня, у которого шаг h'Zk наименьший. Приняв значение этого коэффициента равным 0,8, опре делим величину шага по времени
т = 1 •10” 6 с.
Затем вычисляют коэффициенты в уравнениях дви жения. Значения коэффициентов приведены в програм ме задача «стойка» (Приложение 3).
При подготовке задачи необходимо задать начальные и граничные условия. В программе задача «стойка» вос производится удар стойки об абсолютно жесткую гори зонтальную плоскость. Для этого всем узлам моделисетки сообщается одинаковая начальная скорость, на правленная вниз. Узлам всех горизонтальных стержней и перекрестий в первый момент времени задается пере мещение, равное W0(X, Y, Z) =0,01. Узлам вертикальных стержней соответственно задается перемещение U0(X, Y, Z) = 0,01. Начальные значения всех остальных искомых функций в первый момент времени задаются равными нулю. Во второй момент времени все функции считаются равными нулю. Если например, начальная скорость стой ки составляла 50 м/с, то при принятом шаге по времени х = 1 •10~6 с единичное значение перемещений будет со ответствовать 0,2 мм.
Начиная с третьего шага по времени производится вычисление перемещений всех узлов, за исключением узлов нижних перекрестий, у которых перемещение W
вкаждом шаге приравнивается нулю. Воспроизведение
врасчете ударного воздействия на конструкцию позво ляет в первом же варианте расчета определить резо нансные частоты низших гармоник. Применяя при по следующих расчетах разложение перемещений интере сующих конструктора узлов в ряд Фурье, можно вычислить резонансные частоты и формы колебаний не
скольких гармоник.
В расчете можно воспроизвести и другие внешние воздействия, например испытания стойки на вибростен де. Для этого в узлах, соединенных с вибростендом, за дают перемещения по определенному, например гармо ническому, закону во времени.
В книге изложены основные идеи |
метода |
решения |
||||||
трехмерных |
нестационарных |
задач |
деформирования |
|||||
упругих тел — метода цифрового |
моделирования. Рас |
|||||||
смотрены |
расчеты вибраций |
достаточно |
сложных кон |
|||||
струкций, |
наиболее часто |
встречающихся |
в радиоэлек |
|||||
тронике. |
При |
постановке |
задач |
не |
делается |
никаких |
предварительных предположений о характере протека ния процессов деформирования во времени и пространст ве, т. е. не вводится никаких принципиальных допуще ний. Все примеры расчета доведены до числовых резуль татов. Приведены типовые программы для ЦВМ.
Дальнейшее развитие метода цифрового моделирова ния идет по пути совершенствования расчетных моделей с тем, чтобы увеличить точность отражения физики про цессов и расширить возможности программ для расчетов более сложных конструкций. В настоящее время ведут ся работы по решению нелинейных задач, задач расче та конструкций в криволинейных системах координат и задач, в которых учитывается влияние внешней среды (воздуха, воды) на процессы вибраций.
Несмотря на то, что метод цифрового моделирования является весьма универсальным, его применение ограни чено возможностями современных ЦВМ. Вычислитель ные машины среднего класса, такие как Минск-32, М-220, Урал-14М, позволяют производить расчеты про цессов достаточно быстро протекающих во времени, на пример процессов вибраций. С развитием вычислитель ной техники и совершенствованием метода расчета область его применения будет расширяться.
Метод цифрового моделирования в настоящее время является методом анализа, но не синтеза. Он дает воз можность априорно оценить работоспособность конструк ций при заданных внешних воздействиях, но не позволя ет ставить задачу оптимизации конструкции с точки зрения определенных требований, например минималь ного веса, разработки равнопрочных конструкций и т. п. Это связано с тем, что пока не существует надежных и экономичных алгоритмов такой оптимизации, а метод простого перебора приводит к непомерно большому объ ему вычислений.
ЗАДАЧА БЛОК
INTEGER X, V. Z, T
DIMENSION U0(7»e»I4)lui(7,e,l4),U2<7,e,i4)l
•VO(7*e»lA)>Vl<7#6#i4)iV2(7»e»i4),W0(7#e>14),
•Wl(T*e»14)#W2(7,e,;4),Do(7,s»l4),D2(7ie,i4),
«D4(7’ *»14> IH(7 » Ô» 14 )
присваивание ЗНАЧЕНИЙ КОЗФФИЦЙ^НТАМ
A= I , IOOOÎ B = O. O9I O I A U = O , OSOO
B13=0,0070*Cl3=Q,lB40iDl3sOtOllO
C15S0.0110*01330.О03О»А218О,вООО
В2230.02701В23=0.0270;С23*0.0110
|
023 = 0, ОНО »C23 = 0. 0Ц0 IO23S0. Olio |
|
|
A31=0,0300IB32*0.0070ÎB3350.0270 |
|
|
с 33= 0 , Olio ; 033 = 0 . 0 0 30 ; сз 3 = 0 л 6 ^ 0 |
|
|
0 3 5 = 0 , 0 1 1 0 1 0 1 2 = 0 . 0 2 7 0 |
|
|
ВВОД КОЭФФИЦИЕНТОВ МАССЫ |
|
1 |
FORMAT(7F7.4) |
|
|
J= 1 |
|
|
DO 21=1,112 |
|
|
REAO 1 , H(J) »HCJ♦1 ) , H( J+ 2 ) i H( J♦J) , H(J♦4 ) , |
|
2 |
• H(J♦5 ) , H( J + 6 ) |
|
J= J♦ 7 |
||
|
ВВОД начальных перемещении |
|
|
JS1 |
31=1*112- |
|
00 |
|
|
READ liU0<J)»U0(Jtl)»U0<J*2),U0U*3), |
|
|
• UO( J*4 ) ,U0 U*3 >*U0 ( J*6 > |
|
3 |
J=J47 |
|
|
J= 1 |
|
|
0 0 |
4 1 = 1 , 1 1 2 |
|
READ l , V O < J > » V O ( J M > » V O < J # 2 ) » V O ( J 4 3 ) , |
«*V0( Jt4) ,V0( J^5) ,V0( J*6>
4J= J4 7 J= 1
00 Ю01=1,Ц2
READ l , W O ( J ) » W O ( J ^ l > » W O ( J t 2 ) * W O U + 3 ) ,
*W0(J t4) ,W0 (J♦5 ï ,WO(J*6)
1 0 0 J = J ♦ 7
00 10lZ = l , 14 oo ioiv=i,ê
DO 101X= 1 , 7 v
иКХ^,2)=иС(Х,У*2>^1(Х,У,г)г^0СХ,У,1)
101wi(X, y,Z)=W0(X, V, Z> ЦИКЛЫПО X, V, Z, T
00 9 BT= 1 i 1000
00 962=1,14
0 0 96y = l , e<^---
300 96x = 1 , 7
ВЫЧИСЛЕНИЕ ДЕФОРМАЦИЙ
If(X=7> GO TO 6
É l l « A l l * ( U l < X M » V » Z > - U l O < # y # Z > ) > GO TO ?
6 |
Е Ц |
= 0 . |
|
|
|
|
|
|
|
|
7 I F ( X = 7 W = 1 V V = 8 ) GO TO 8 |
|
|
|
|||||||
|
Fllsbl2*(Vl<X«.bV*i»£>*Vl(X,V + l.Z>- |
|
||||||||
• V i ( X + 1 . y - 1 . Z ) - V i ( X i y - l . Z ) ) î GO TO 9 |
|
|||||||||
8 |
F1 x »0 • |
|
|
|
|
|
|
|
||
9 I F { X = 7 v Z = l v Z = i 4 ) GO TO 10 |
|
|
||||||||
|
6 1 1 5 B i 3 » ( w i ( X 4 . i , y , z ^ l ) * W l ( X , V * Z ^ X ) - |
|
||||||||
|
• W U M r y ^ Z - l l - W K X . V d - l ) ) ! |
GO TO Ц |
||||||||
10 |
G1 1 * 0 . |
|
|
|
|
|
|
|
||
11 I F ( X = 1 ) G О ТО i 2 |
|
|
|
|
|
|||||
12 |
Е 12 = А 1 1 # ( U l ( X , У , Z ) - U 1 ( х - 1 »V , Z> > » GO ТО 13 |
|||||||||
Е12"0 , |
|
|
|
|
|
|
|
|||
13 I F ( X r l v y = l v V s 8 ) GO ТО 14 |
|
|
||||||||
|
F l 2 s & l 2 * < V l ( X , y * i , z ) « . v i ( X . l , V * X » Z > - |
|||||||||
|
*V 1 (X»V-li Z) - VI(x-l>У-1; I )) |
\ GO ТО |
15 |
|||||||
14 |
F 12 = 0 , |
|
|
|
|
16 |
|
|
||
15 |
l F ( X s l v Z * l v Z s l4) GO |
ТО |
|
|
||||||
|
G l 2 = 8 l 3 » ( W l < X , y , Z ^ i ) 4 W l ( X - l lV , Z ^ l > - |
|||||||||
|
• W l ( X . y » Z - l > - W i ( X - l , V , Z - l > ) » GO TO 17 |
|||||||||
16 |
G 1 2 ’ 0 » |
|
|
|
|
|
|
|
||
17 |
I F ( V = 8 ) |
GO |
TO |
18 |
|
|
|
|
||
|
|
E l3 = C l 3 » i U l < X , V 4 - i , 2 ) - U l ( X , V , Z ) ) iGO |
TO 19 |
|||||||
18 |
E 13 = 0 , |
|
|
|
|
|
|
|
||
19 I F ( X s l v X s 7 v V 3 8 ) G O TO 20 |
|
|
||||||||
|
|
F l 3 s 0 l 3 # ( V l ( X ^ i , y * i , Z ) ^ V H X * l , V , Z > - |
||||||||
|
* V l ( X * l # y + l » Z ) - V l ( X - i , y t Z ) ) I GO TO 21 |
|||||||||
20 |
F13 3 0 « |
|
|
|
|
|
|
|
||
21 |
I F < y = l > |
GO |
TO |
22 |
|
|
|
|
||
22 |
E i ; 3 C l 3 M U l ( X lV , Z ) . u i ( X , y - l , Z ) ) J GO TO 23 |
|||||||||
E1430 , |
|
|
|
|
|
|
|
|||
23 |
I F { X s 1 V X = 7 v V = i )GO |
TO 24 |
|
|
||||||
|
|
F l 4 a 0 l 3 M V l ( X * i , y lz> + v i ( X + l , V * l , Z > - |
||||||||
|
* V l ( X < - l , V , Z > - V i l X - l » V - i , Z > ) I GO TO 25 |
|||||||||
24 |
F 1 4 5 0 » |
|
|
|
|
|
|
|
||
25 |
I F (X a 7 Î GO |
TO |
26 |
|
|
|
|
|||
|
|
E 2 5 e C 2 5 M V l < X * i lv , z > - V l ( X , V , Z > ) Î GO TO 27 |
||||||||
26 |
E 25* 0 • |
|
|
|
|
|
|
|
||
27 ! F ( X a 7 v V 3 i v V a e f G O TO 28 |
|
|
||||||||
|
|
|
F253025^(U1(X^I I V^1I Z)^U1{.X,V4];,Z>- |
|||||||
|
|
* u i ( x n » y - i » z > - u i < x , v - i » z > j ; oo T O 29 |
||||||||
28 |
|
F2 53 0 , |
|
|
|
|
|
|
|
|
|
29 |
|
I F {X = 1 )GO |
TO |
30 |
|
|
|
|
|
|
|
|
E 2 6 e P 2 5 * ( V l ( X fy , Z ) - v l ( X - l , V , Z ) ) I GO TO 31 |
|||||||
|
30 |
|
£26 = 0 • |
|
|
^ |
|
|
|
|
|
31 |
|
I F ( X = lvV = l W = 8)G0 T O |
)2, |
|
|
||||
|
|
|
F2 6 =0 2 5 * ( Ul ( X, V4 i , z ) ^Ui ( X - l , y - l , Z> - |
|||||||
|
|
•U2(X»V-l»Z)-lli(X-liV-liZ> ) * 00 TO 33 |
||||||||
|
32 |
F 1 6 ï 0 , |
|
|
|
|
|
|
|
|
|
33 I F ( V =8 ) G О TO 34 |
|
|
|
|
|||||
|
|
|
E 2 i = A 2 1 * ( V I ( X , y * l , z ) - v И X , V , Z ) ).l GO TO 35 |
|||||||
|
34 |
E 2 1 = 0 . |
|
|
|
|
|
|
|
|
|
35 I F ( V » 8 v Z = l v Z = l 4 ) G 0 TO 36 |
|
|
|||||||
|
|
|
F 2 i = B 2 2 M W l < X , v * b Z * l M W l ( X , y , Z * l > - |
|||||||
|
|
• W n X » V M * Z - l > - W l < X ; y , Z - l > ) • GO TO 37 |
||||||||
|
36 F 2 l = 0 , |
|
|
|
|
|
|
|
||
|
37 |
I F ( X = lvX = 7vV38)GO |
TO |
38 |
|
|
||||
|
|
|
G2i = &23MUl<X4i#Vn ,z>*UUXfbV,Z>- |
|||||||
|
|
• U l ( X - l i V M i Z ) - U l ( X - i , V . Z > ) i GO TO 39 |
||||||||
|
38 |
G2 1 - 0 , |
|
|
|
|
|
|
|
39 1F( V5 1 ) G0 ТО 40
”, 2, ! Aftn *tVl<X' V-Z>-VUX.V.l.Z))if'0 т0 41 40 ЬZ 2•u*
41 iF(V = lvZ= lv2*14)G0 Tо 45
F22»®22*(Wl(X,v,Z*p^wvy ,,j*i)*
i2‘ r22U!’ Z' 1>' Wl<X,V- l *^î''> 00 T° °
43! F(X=lvXr7vV= i ) GO TO 44
G2 2 sB2 î*«ui(Xti,y,z)+ut(Xlv. l , Z)-
u , « ] / 0! , v , n ' UlU- l ' v |
|
BO 10 |
45 |
- |
|
45 |
tF(Z»14)G0 TO 4 6 |
|
|
|
|
46 |
E23«0î î *<V1(X' V’ Z*1>”VllX,V*Z)1' |
G° |
T |
||
47 |
IF(V*lvV=8 vZ=1 4 )G0 то 48 |
|
|
|
|
|
r 2 3 = D 2 3 M „ l ( X . V T 1 . z n > t W M X | V t l . Z >: o |
||||
4I ‘ F23!0 l , 2 ‘ ‘ l ' W l U , 'l ' 1 ' n ' 1 |
G° T° |
|
|||
47 |
IflIMIGO TO 50 |
|
|
|
|
50 |
£24-023*tntX' y’ Z,*n<X,V' Z' 11 ’ ‘ |
G° |
T0 |
||
31 |
1F(VS1VZ=1VZÏ 14)G0 TO 52 |
v + l#z-l)- |
|
||
|
F24s023»(WUX,y.+ i , z)*wl(x |
|
|||
52 |
•Wl(X»y-l,Z>-WKX»y-i,z-i))| |
GO TO 33 |
|||
F24a0 , |
|
|
|
|
|
53 |
1F(Vs 8 >GO TO 3 <, |
v |
Z) , ; |
GO TO 33 |
|
|
E35ïC35.twHX,y*llZ). Wl(X |
54E35*0 •
55IF<VsôvZ=ivZ=l4)G0 Ю 56
F338D35.<V1<X,V+1.Z+1>*V1(X,V.Z*1>-
• V U X » V * 1 , Z- 1 )-VI (X , v , z- 1 ) ) ÎGO то 57
56 F35*0 ♦
57 I F(V= 1 )G0 TO 58
5 8 |
E34.C35*tWl(X,v,Z).wi(X.V-l,Z)): 00 T0 59 |
E36«0 , |
|
39 |
IF(VslvZsivZ=14)G0 TO 60 |
|
f36«035«(Vl(X,y,Ztl>tVl(X,y-l»Z*X>- |
60 |
•V1(X(y,Z4i).Vl(X.y.j,z-1 ))i GO TO «4 |
F36* 0 . |
61IF(Z3 H)GO TO 62
E3;«A31*(Ml(XfvvZ41).Hl(X»y,Z) >• G0 T° 63
62E3l«0,
63l F<X3 1 v Xs 7 v Z= 1 4 ) GO TO 64
F3I «B32«(UUX*1#V.Z^I)*UI(X*1*VI E>*
*Ul(X*l#V,Z^l)-Ui(X«i»ViZ)) » GO то 63
64 |
F3 1 «0 » |
|
65 |
I F(V5 1vVs8 vZ = 14)GO T0*66 |
|
|
G3i»B33#(vi(X,yti,Zfl)^VUX#V^l»z>][ |
|
66 |
*V1(X.M,Z»1)-V1(X,V-1.Z))5 |
GO TO 67 |
G31 »0 # |
TQ df |
|
67 |
I F(Za i ) G0 TO 68 |
E32«A3i*(Wl(Xry,Z)-WltX»VrZ-l)î‘ 50
68E32■0 ,
69I F( XaivXs?vZ=l)GO TO 70 F32*B32*(Ul<X^l,y,z>*Ul(X«bV»*-l;, ; i
70 |
•Ul(X*l,y,Z)-Ui(X-l»V»Z-lï)» |
GO то 71 |
F32*0• |
|
|
71 |
I F(y = ivy= 8 vZ = l ) GO TO 72 |
|
|
6 3 2 * В э з М У 1 ( Х | У + 1 , (2 > + У П Х , У * 1 » * : * , " з |
|
|
•У1 (х»у.1 ,г)-У1 (х»у,ь2 т1 )1 ; |
oo T0 |
72GJ2-0 ,
73I F(X37 )G0 TO 74
74 |
H r ; co’ 3Mvn()Ubv' z) - wi()<' v' z),îGo |
T° 75 |
73 |
I F(X3 7vz*ivZ»i 4 )GO TO 76 |
|
76 |
. U l ( 5 D* î 3 V , 7 i r ) ( 4 l ' V , Z * l M U l ( ) ( , V , I 4 ) e |
|
F33-01’ V,Z*l),’Ul(X' V,Z’ l) 1! °° T° |
77 |
77I F(Xa1 )GО TO 7 ô
E3A,C33MwKXlv,7>.Wl(X.bV,Z))î GO TO 79
78E34“0 ,
79 |
I F( X= lvZ= lvZ= 14)G0 TO 80 |
8 0 |
•U1 (X»V»Z-i)-Ui(X-i,v,z-l))î GO TO ex |
F34* 0 , |
81IF(I=14)G0 ТО 82 El3sC15«(Ul(X,v,2*j)»ui(X,V,Z)1i GО ТО 83
82EI 5s0 ,
83 I F(Xs1vX* 7vZ * 14>GO ТО 84
F13= Dl3.(Wi(X*i,V.z + l)+WUX*l,V,Z)-
•W1(X-l, V, Z*l)-W1 (X-i, у,2 )) >GO ТО 85
84 F15*0,
85IMZs1> GO ТО 86
El6ïC15»(UUX , V.z )-Ul(X ,V ,z-i) ) 1GO те 87
86E l 6 * 0 ,
87 |
1F(X=1VX=7VZ=14)GO ТО 88 |
|
Fi6 ïDl6 »(Wl(X*l,y,z)+Wl<X6 l,y#Z-i)- |
88 |
• W K X - l * V i Z >-H*i l x - l . V#Z- 1 >) » GO ТО 89 |
F l 6*0 , |
|
89 |
вычисление напряжений |
ail*EiltFl l +Giiiai2sEl2*Fi2+Gl2 |
ai3*El3*Fiïiül*sEl4*Fi4îqi5sElS*Fl5
a i 6 s E i 6 * F i 6 î a , 2 1 s E 2 l * F 2 l * 6 j l a22*E22*F22 + G22!a23*E23*F;23
a26sE2**F2*:a23sE23*F23:a2<sE26+F26 a3iaE31*Fïi + G3i;ü32 = E32*F.J2+G32
0 3 3 = E 3 3 * F 3 3 î a 3 4 * E 3 4 * F 3 4 t a j 3 a 8 3 5 * F 3 5
Q 3 6 * E 3 6 * F 3 6
ВЫЧИСЛЕНИЕ-УСКОРЕНИЙ 0i«au-ai2*ai3-au*ai5-ai6
03*021-02 2*823*02.4*025-026
05*031-032*833-аз4*035-036
ВЫЧИСЛЕНИЕ ПЕРЕМЕЩЕНИЙ и -ЗАКРЕПЛЕНИЕ пояса
IF({X*1AZ*1>V<X*1AZ*14>V(X«?AZ»14>V
.(Х*7лг*1))GО ТО 90
ü2(X»y»ï>*H(X,V,Z>»A*(Dl-6»00tX.-y»Z)>* •2.*Ul<X,y,Z>-U0<X,y,Z>*G0 То 91
90U 2 ( X > V . Z » * 0 .
91IF(X*1VX*7VZ*IVZ*14»60 ТО 9j
V2(X.y.Z>*H(X,V,Z)*4*<O3-6»02CX,VU) >♦
. 2 .*VHX,y,Z)-V0 (X,V,Z) 100 То 9J
92U2(X»V»Z)*0«
93IF((X?1AZ*1)V(XS1AZ*14>V(X«7AZ«14>V
*(X*7AZ*1))GО ТО 94
W2(X.y.Z)*H(X,V,Z)*A*t05-B*02CX,V»Z)>*
.2 .*N1(X,V,z>-W0tx.v.Z); GO то '95
94W2(X•V>Z>*0 • ПЕРЕЗАПИСЬ УСКОРЕНИЙ
93 DO<X*y.Z)*01iD2(X,y.Z>*03f04IX,v,Z»*PS
06 CONTINUE