книги / Цифровое моделирование вибраций в радиоконструкциях

DO 99 Z= i #14

DO 99 У= 1,8

DO 99 Х=Д,7

СПЕРЕЗАПИСЬ ПЕРЕМЕЩЕНИЙ

U0 U*V»Z)sUl(Xiy,Z)ÎV0 (X»ViZ)=VUX,V,2 )

W0(X.V»Z)=Wl(X,ViZ)ïUl(X,V»Z)=U2(XrV»2)

Vi(X#V,Z>=V2<X»V»Z)

99Wl(X#V»Z)=W2<X,V,Z>

СПЕЧАТЬ ПЕРЕМЕЩЕНИЙ ЦЕНТРА БЛОКА

X= 4 JV=4 î Z= 7

U4s(î.23e<U2<X,y,Z)«.U2<X#V*l,Z)*

• U2 <X»У»Z41ï4U2(X»Vf 1.Z*1) )

V4r:0.25»<V2<x»V»z)4,V2(X#V*liZ>t

•V2(X»yiZ*l>*V2tX,y4l,Z*l>)

WArfl.25*(W2tX»ViZ)4W2<X»V4l,Z)i

•W2 (X*y»ZM)*W2 <X.y*l»ZM>)

97ГORMAT(3F8#A) WRITE 97.U4.V4.W4

98CONTINUE ST0 P»END

Приложение 2

ЗАДАЧА ПЛАТА IMTEGER X, У»T, H

DIMENSION U1<1A,12),U2(14I.12)»U3(14,12),

V1 ( 1 A»J2 ) » #V2(l<>»l2),v3(i4,i2)»W0(l4li2),Wl(l4,i2)> W2 ( i4 , 12 >*

•D0(l4,i2) , D1 <14 #1 2 ) >D2(14 y1 2) * H( 1 4,12 ) , A( 14,1 2 )

СПРИСВОЕНИ ЗНАЧЕНИЙ КОЭФФИЦИЕНТАМ

AU«0 . 0091 JB11 = 0 . 0095 ;A13e0 .0027îB13no, ooo? А15в0.вооо;В15=0'0135;А21з0.0091'В21гО'037|

А23вО,0027;В23=0.0007»А23яО,ёООО»В25=0,0027

Û3I *0.0027;B3I =0,0^05;A33«0,00271833=0,0403 8 = о»looo;

с ВВОД ПРИЗНАКОВ 1 ГORMA T(14I2 >

J= 1

DO 21=1.12

*H( J♦5 )»H( J ♦6 >»H( J♦7 ) * H(J♦• ) .M(j *9 ) » H ( J ♦10 ) •

•H( Jtll >»HU*12 Î *H<J*13)

2J=J*14

СВВОД КОЭФФИЦИЕНТОВ A

3FERMAT!7F7.4)

JM

DO 41=1.24

READ 3,A<J)IAU*1>,A<J*2>IA<JO>»A<J+ 4>»

*A(Jf5)» A( J + 6)

и J = Jf7

С4 ЗАДАНИЕ НАЧАЛЬНЫХ ЗНАЧЕНИИ ПЕРЕМЕЩЕНИЙ DO 51=1*168

5 WlСI>= 0♦01

С ЦИКЛЫПО КООРДИНАТАМ ИВРЕМЕНИ 00 59Т=1.1000

*' 13) 60 T 0 3 4

D 2 3 ? О, » 60 T 0 3 5

34 D 2 3 - & 2 3 * ( U 2 ( X + 1 I V+>I ) + U 2 ( X | V + 1 ) »

•U2 (X♦1 »v-l)-U2 (X»V- 1 ) )

35 i r ( H ( X » y ) s 0 l v H ( X , V ) s 0 4 v H ( X i V ) = 0 3 ) G O : T 0 36

42 I F - ( H U i V ) s 0 2 v H ( X » V ) a o 3 v H « X » y ) s 0 7 ) G O T 04 3

C 3 1 5 A 3 1 * ( V U < X 4 1 , V > - W 1 < X , V ) ) ; Q 0 T 0 4 4

43 C31s0 ,

44 I F ( H ( X I V) = 0 1 V H ( X » ‘V ) = 0 2 V H ( X I V ) * 0 3 V H ( X # V ) S

* 0 4v H ( . X # y ) s 0 3 v H ( X » V ) = 07)GO То 4 5

D 3 i s B 3 1 * ( U 2 ( X ^ i l V ) - u 2 ( X * l , y ) ) J G 0 T 0 46 45 0 3 1 = 0 ,

46 I F ( H ( X » У ) = 0 1 v H ( X » V ) = 0 4 v H ( X » У ) = 05 ) GO T 04 7 C 3 2 5 A 3 1 « ( W I ( X , V ) - W I ( X - 1 , V ) ) ; G 0 T O 4 6

47 C3 2 = 0 ,

48 I F ( H ( X »У ) = 0 3 v H ( X . У ) = 0 4 v H ( X »V ) = 08 ) GO T Q49 C33 = A j 3 M V U < X , y « . l ) - W l ( X » V } > ; G O T 0 5 0

49C3 3 = 0 ,

50I F ( H ( X » V ) s o i v H ( X i y ) s 0 2 v H ( X i - V > s 0 3 v H ( X # y ) *

* 0 4 v H l X » y ) = 0 6 v H ( X i V ) a 0 8 ) G O T o 3 l

0 3 3 = B 3 3 « ( V 2 ( X , y 4 i ) . v 2 ( X i y - l ) ) i G 0 T052

31 0 3 3 = 0 .

52 l F ( H ( X » y > = 0 1 v H ( X i y > 3 0 2 v H ( X » y ) = 0 6 ) G O T 0 5 3 C 3 4 5 A 3 3 * <W 1 f X »■V )■•W 1 ( X »У• 1 ) ) ;G O T 0 54

53C34= 0,

СВЫЧИСЛЕНИЕ НАПРЯЖЕНИЙ

34 F l i 5 C i l 4 D l l } F i 2 s C l 2 ^ D l 2 » F l 3 = 0 1 3 4 0 1 3 F 14 = C l 4 4 D l 4 » F i 5 =C 1 5 4 0 1 5 î F g l s C 2 1 4 0 2 1

F 2 2 = C 2 2 * D 2 2 » F 2 3 = C 2 3 4 D 2 3

F245C244D24<F23=023^023?F31=0314031

F325C32^F335C334D33»F34aC34

СВЫЧИСЛЕНИЕ УСКОРЕНИЙ

E 4 = F 1 1 - F 1 2 4 F 1 3 - F 1 4 - F 1 5

E 3 s p 2 i - F 2 2 4 F 2 3 - F 2 4 - F 2 3

E 3 * F 3 i - F 3 2 4 F 3 3 - F 3 4

СВЫЧИСЛЕНИЕ ПЕРЕМЕЩЕНИИ

СЗАЩЕМЛЕННЫЙ КРАЙ

IF(V5ivV312)GO TО6 о

U3(X#y)£A<X»y)4(E4-B*ElUfy))42.4U2<X#y>

* U 1 ( X » V )

УЗ(Х»У)5А(Х,У)4(Е5-В4Е2(Х|У))42|*У2(Х#У>

•VllX.y )

^ 2 ( Х # У ) я А ( Х # У ) 4 ( Е з . В 4 Е 0 ( Х , У ) ) 4 2 . 4 ^ 1 ( Х » У )

•W0(X»V) i GО T061

60 U 3 ( X » y ) 3 0 . » V ( X , y ) = o , ) W 2 t X , y > 5 o #

СПЕРЕЗАПИСЬ УСКОРЕНИЙ

41EltX*V>*E4)E2(X»V)cE5iÊ0<X*y>*E3

p * « A 0 l M B u x M , V . Z ) . 2 ; . B H X , y , Z I - » l C l . W . * , z

#A 4 l . ( W l ( X * 1 , v , Z ) - U l ( l , y - * , Z ) ) - A 4 1 , e i ( X . y lZ)

p 5 î A 0 1 » < C l ( X * l , y , Z ) . 2 > # c l ( X , v , Z ) * - » l ( b y . 4 , z , ) .

U 2 ( X » y * Z ) = H ( X , V , Z ) . ( D O - A » | 0 ( X , V i Z ) ) ^ 2 . » U l ( X , V , 2 )

-1 (Cl<2>v'Zl*CJ«l>y>Zl)

»c 3 1 <

h;Y sВА1 *(W1(2 »V. Z»-W1 ( 1 »V, Z>) ♦ #CÂ1#(&1(2*V,Z)^BK1,V,Z))JG0 T0 2 7

QXi-Gll*<Ul(2.ViZ)-U1t1 »V, Z)ï 26 uVi=B21*(Vl(2»ViZUVUl»V|Z)î*

eCU*lCl<2»V,Z)*Cl<l,V,Z)> a*i::B2lMWi<2.V,Z>-Wl(l.V,Z>>*

# C2 1 * I B 1 ( 2 . V , Z ) ^ B 1 U » V , Z ) ) 1 F ( V = 3 >GO T O2 9

2 I F ( V = 1> GO T O2 8

ÛX2-BU*<-V1<1».V»Z)-U1<10, V-l» Z) ) QV2-B2l*(Ul(l»V»Z)-Vl(10»V-l»Z))+ eCU*(Cl(bV.Z)*Cl(lO.V-bZ))

û Z ' » * B 2 1 * ( W l ( l i V » Z ) - W l ( i O » V - l » Z ) ) +

#C2 2 *<-^1 <1 »V»Z)^B1 (1 0 »V-1 I Z))#gO ТОЗО

* A flX2?B31»(-Vl(l»V»Z)-Ul(10|V*3#Z))

Z ûV2-B/>1#(ul ( l »v»Z,- vl tl0' V*3,Z))* #c3i»<ci<bv,z>*ci<t.o*vo»z> >

ûZ2'B4i*(wl (i ' v»Z)-^lll0,V+3,Z))* •+C41#<-A1( l#y»Z>*Bl(lO»V*3»Z) ) »GО ТОЗО

59 û X 2 = B 3 l * ( - V l ( l » V » Z ) - U l ( 1 0 i ^ - l » Z ) )

ûy2-p^ l MIJl (l*V,z).Vl(10»V-l»Z))^

•C3l*(Cl<1»V,Z)4C1(10»V-1,Z)> ûZ2«B41*<W1<1,V»Z)-wi(1 0 »V-liZ))*

# C 4 1 » ( - A 1 ( l , V i Z ) ^ B l ( l O » V - l , Z ) )

ВЫЧИСЛЕНИЕ ПЕРЕМЕЩЕНИЙ перекрестия

3 0 U2U* ViZ)=aXi-QV2^Z3^QZ4^2^UlCl#X,Z> -UOU»V,Z)

V 2 ( X » V » Z ) = Q V l ^ Q X 2 ^ a V 3 - Q V A* 2 4» V l ( X » V | Z )

- V 0 U » V »2)

W2( 1»V>Z)=QZ1-QZ2-GIX3*,QX4*2 »* W1 ( 1 #У» z *

-и о(î #v#z )

ВЫЧИСЛЕНИЕ МОМЕНТОВ

IF(VslW = 3) GO ТО л RXis-Ell*C-Al(2 #V,2 >-AiCl»VlZ))

Ryi=-E21*(Bl(2»V»Z)-Bl(l»V»Z)) RZl5E31#(Cl(2»V,Z)-ClU,VrZ))JG0 T032

31RXi5-Dll*(Al<2iV*Z)-AKliV»Z>> RVlsD2l#(B1^2»V,Z)-Bl(bV|Z)î

RZl503le<CX(2,V,Z)-Cl(l*V,Z))

RX350, ÎRV3 = 0. ÎRZ3= 0 »»G0 ТО 35

33 RX3--E71*(CKi,V,Z)-AlU0»VtA,Z-l))

R V 3 s E e i * < - B l ( l # V » Z ) - B l ( l O # V # A » 2 - D )

RZ33E9l*<Al( 1 ,V»Z)-Cl(10 »V^4»Z»1)>» GO ТО 33

34RX33-D7lMCi<i,V,z>-Al(l0rV*A'Z-i>)

RV3sD81#(-Bl(irV»Z)-Bl(lO»V^A#Z-i))

RZ3= D91*(Al(bv,Z) . ci (10,V + 4»Z-l)ï

RX4 5 O. ÎRV4= 0. ; R24SQ , ; GO TQ3B

36RX45-E7l«(Ai(2*y*4fZ>-Cl<l»V»Z1) RV4SE61MBI (2,V+ 4,Z>*6I <1I V,Z))

RZ45E9l»<Cl(2,V44,Z>-Al(l,VfZ) );GO TO -.fi

37RX4s-D7i<Ai<2,y*4,z>-Cl<bV,Z) ) RW4«081*(Bl(2,y^4,z)^bl<l»V,Z))

RZ4sD91MCl(2*V*W>-AiU#V,Z>)

RX2î-HAi*<ei<l#V#Z)-Al<lO»V-l»Z))

RV2iE5H<-Al(l»V»Z)-Bl( 10, V-l»Z) ) RZ25E6J*{Cl(l#V>Z).CltlO,V*X»ZniGO T04l

39 RX2S*D41MB1<1*V#Z)-A1(10,V4-J»Z) )

RV2SD314<-A1(1*V»Z)*B1(10,V*JI Z))

RZ2-D6H<Cl(bViZ>-CltlO#V*3»Z))JGO TQ4X

4n RX2S^04Î MBl ( l * V» Z ) - A l U 0 , V-1*Z> )

° RV2-03iA(*Al(i#V#Z)-Bi<lO,V-i#Z))

RZ2-06i*<Cl(l*ViZ)«CltlO#V-l#Z)) вычисление УГЛОВ ПОВОРОТА ПЕРЕКРЕСТИЯ

L 4 ! Г(y = 2vV*4)00 Т042

1 A2(i#y#Z)a-RXi^RV2-RZ3^RZA>F01* #ûZ2-PUeCV3-ril«ûy4^2.^Am,V,Z)-A0<l»V»Z)

B2(l»V»Z)aRV14RX2*RV3-RV4-F2l«aZl^ #F3l*ÛZ3^F3X#eZ442«#BX-Cl»V,Z)-B0<l»V»Z)

C2Cl»V»Z)*RZl-RZ2^RX3-RX4.F4i*QVl-F31* #ûV2*Z.#cl (l»V' z>eC0U»V#Z)iGO T0 4 3

A2(1»V»Z>=-RXI *RV2-RZ3*RZA*G00*ÛZ2-

42#-С10жвУЗ-010жОУ4^2,*А1(1»у#2>-АО(1»У»2)

B2(l»V»Z)eRyi^RX2^RV3-RV4.G2l»eZX'*‘ •G3l#e23^G3X#GZ442.*0X(l»V»Z)*BO(X»V»Z)

C2(X»V»Z)=RZX-RZ24RX3-RX4.G4X*QV1- •G3X*Gy2^2.#CX(X#y»Z^eO°(Xry#Z)

,, CONTINUE

ПЕРЕЗАПИСЬ ПЕРЕМЕЩЕНИИ и углов

00. 442з X»з

DO 4 4 у я х » 8

DO 44 Х*х » ХО uO(XiViZ)sux(x,v.z);vo(x,y,z)=vi<x»v»z)

WO(X»V,Z)*KX(X, V,Z ); AO(X, V, Z)SAX(X»V»Z)

BO(Xiv»z)=Bx(x,y,zj;co(x,y,z)*cx(x»y»2)

ül(X,y.Z)*U2(X,V,Z)ÎVX<X#V,Z)*V2<xiV#Z)

WXîx»y»Z)*W2(X,y,Z))AX(X,V,z)=A2(x»y»Z)

ВX( x » V#Z) *B2 ( x., У, Z) lL CX(X»V»Z)*C2(X,y,Z)

ПЕЧАТЬ ПЕРЕМЕЩЕНИЙ

иЬ ГORMAT(6FX0 . 4>

00 46z * X» 3

00 46V* 1 » 6

и6 VERITE 43 ,U1 ( 6 , У,1 ) , Vx ( 6 , y , Z)•,WX ( 6 #V.Z) »AX ( * 1

•Bx(6iV#Z)>Ci<6»V»Z)

u7 C O N T I N U E

STOP : Er-':>

1.Бабушка И., Витасск Э., Прагер М. Численные процессы реше­ ния дифференциальных уравнений Пер. с чешек. М., «Мир», 1969.

2.Вазов В., Форсайт Дж, Разностные методы решения дифферен­ циальных «уравнений 'в частных «производных. Пер. с англ. М., ИЛ, 1963.

3.Годунов С. К., Рябенький В. С. Разностные схемы. М., «Наука»,

1973.

4.Годфри Д. Е. Р. Теория упругости и пластичности. Киев, «Буднвельник», 1969.

5.Карпушин В. Б. Вибрации и удары в радиоаппаратуре. М., «Сов. радио», 1971.

6.Маквецов Е. Н. Математические основы цифрового моделирова­

ния вибраций в радноконструкциях — В кн.: Вопросы проекти­ рования специальных радиоэлектронных устройств. Под ред.

Е.Н. Маквецова. Изд. Пеизси. политехи, ин-та, 1972, вып. 2.

7.Рихтмайер Р. Д. Разностные методы решения краевых задач. Пер. с англ. М , ИЛ, 1960.

8.Самарский А. А. Введение в теорию разностных схем. М., «Пау­ ка», 1971.

9.Самарский А. А., Гулин А. В. Устойчивость разностных схем. М., «Наука», 1973.

10.Слепян Л. И. Нестационарные упругие волны. Л., «Судострое­ ние». 1972.

11.Угодчиков А. Г., Длугач М. И., Степанов А. Е. Решение крае­

12.Филиппов А. П. Колебания деформируемых систем. М., «Маши­ ностроение», 1970.