Контрольная работа по сопромату Методика решения задач Дополнительные задачи на сдвиг Сложное сопротивление Действие динамических нагрузок

Контрольная работа по сопромату Методика решения задач

ПРИМЕНЕНИЕ ЭВМ ДЛЯ РЕШЕНИЯ ЗАДАЧ ПО СОПРОТИВЛЕНИЮ МАТЕРИАЛОВ

  Введенные во всех высших и средних технических учебных заведениях новые учебные планы и программы создают необходимые объективные условия для широкого использования ЭВМ. Рациональность использования ЭВМ особо ощутима при расчете статически неопределимых систем. Однако и при расчете некоторых статически определимых систем могут быть использованы ЭВМ. Это в первую очередь относится к таким задачам, решение которых состоит из большого числа аналогичных последовательных операций.

9.1. Вычисление моментов инерции плоских составных сечений

 Геометрические характеристики плоских сечений рассматривались в главе 2. В разделе 2.3 предлагается порядок расчета для сложных составных сечений. Эту методику легко реализовать на ЭВМ. Ниже приведена программа на алгоритмическом языке Фортран-IV. Ввод числовых данных осуществляется самым простым способом – способом «присвоения».

 В качестве образца взят числовой пример расчета, рассмотренный в п.2.3 (рис. 2.3.1).

  Ввод числовых данных для конкретного числового примера осуществляется в программе от метки (5) до метки (8) включительно.

 ЭВМ выдает на печать координаты центра тяжести, осевые и центробежный моменты инерции относительно центральных осей хс, ус, а также главные моменты инерции и тангенс двойного угла наклона главных осей.

 PROGRAM AXE Построчные пояснения Мебель от производителя / производство / кухни из массива дерева дуба, ясеня, бука.. Компания inverter24 поможет установить кондиционер сплит систему и выполнит заправку фреоном. Монтаж кондиционеров производится мастерами, имеющих большой опыт работы. Установка осуществляется в удобное для вас время.

C  Геометрические характеристики плоских сечений 

DIMENSION X(5), Y(5), A(5), FIX(5), FIY(5), FIXY(5), AB(5), BA(5)

 5 N=4

X(1)=25. х1 = 25 см

Y(1)=24.8  у1 = 24,8 см

A(1)=50.*1.6 А1 = 50·1,6 (см2)

FIX(1)=50.*1.6**3/12.  Ix1 = 50·1,63/12 (см4)

FIY(1)=1.6*50.**3/12. Iy1 = 1,6·503/12 (см4)

FIXY(1)=0.  Ix1y1 = 0

X(2)=43.42 x2 = 43,42 см

Y(2)=12. y2 = 12 см

A(2)=30.6  A2 = 30,6 см2

FIX(2)=2900. Ix2 = 2900 см4

FIY(2)=208. Iy2 = 208 см4

FIXY(2)=0. Ix2y2 = 0 см4

X(3)=36.11

Y(3)=4.89

A(3)=42.19

FIX(3)=1316.62

FIY(3)=1316.62

FIXY(3)=776.5

X(4)=5.32

Y(4)=21.64

A(4)=30.04

FIX(4)=238.75

FIY(4)=784.22

  8 FIXY(4)=249.2

C Определение координат центра тяжести

SY=0.

SX=0.

AA=0

FIXC=0.

FIYC=0.

FIXYC=0.

DO 10 I=1, N, 1

SY=SY+A(I)*X(I) Sy = ΣAixi (см. формулы (2.1.5))

SX=SX+A(I)*Y(I)  Sч = ΣAiyi (см. формулы (2.1.5))

AA=AA+A(I) A = ΣAi

 10 CONTINUE

XC=SY/AA xc = Sy /A(см. формулы (2.1.7))

YC=SХ/AA yc = Sx /A(см. формулы (2.1.7))

WRITE (7,15) XC, YC

 15 FORMAT (5X, ‘Координаты центра тяжести’,// 7X, 3HXC=, F5.2, 3X, 3HYC=, F5.2)

C Вычисление моментов инерции относительно центральных осей

DO 20 I=1, N, 1

AB(I)=Y(I)–Y(C)  ai = yi – yc

BA(I)=X(I)–XC bi = xi – xc

FIXC=FIXC+ FIX(I)+AB(I)**2*A(I)  Ixc = Σ(Ixi + ai2Ai)

FIYC=FIYC+FIY(I)+BA(I)**2*A(I) Iyc = Σ(Iyi + bi2Ai)

FIXYC=FIXYC+FIXY(I)+AB(I)*BA(I)*A(I) Ixcyc = Σ(Ixiyi + aibiAi)

  20 CONTINUE

С Вычисление главных моментов инерции

FIMAX=(FIXC+FIYC)/2.+0.5*SQRT((FIXC–FIYC)**2+4.*FIXYC**2)  (см.(2.2.11))

FIMIN=(FIXC+FIYC)/2.–0.5*SQRT((FIXC–FIYC)**2+4.*FIXYC**2)  (см.(2.2.11))

TG=2.*FIXYC/(FIYC–FIXC) (см.(2.2.12))

WRITE (7,25) FIXC, FIYC, FIXYC, FIMAX, FIMIN, TG

 25 FORMAT (5X,’Осевые и центробежный моменты инерции относительно центральных

 * осей’//7X, 4HIXC=, F9.2, 3X, 4HIYC=, F9.2, 5HIXYC=, F9.2,// 5X, ‘Главные

 * моменты инерции’//7X, 5HIMAX=, F9.2, 3Х, 5HIMIN=, F9.2//5X, ‘Тангенс двойного

 * угла наклона главных осей’//7X, 3HTG=, F10.5)

STOP

END

Результаты расчета, выдаваемые на печать:

 Координаты центра тяжести

 Хс=27.41 Ус=17.54

  Осевые и центробежный моменты инерции относительно центральных осей

 Iхс=16884.53  Iус=45135.47 Iхус=-10452.02

 Главные моменты инерции

 Imax=48581.96  Imin=13438.04

 Тангенс двойного угла наклона главных осей

 TG= –.73994

  Задача 9.1.1. Найти координаты центра тяжести и вычислить главные моменты инерции для составного сечения, показанного на рис. 2.1.11.

 Решение. Для получения результатов используем приведенную программу для ЭВМ на алгоритмическом языке Фортран-IV. Ввод числовых данных для примеров в этой программе осуществляется от метки (5) до метки (8) включительно. Введем нумерацию прокатных профилей и обозначим швеллер № 20 – 3, двутавр № 20 – 2 и горизонтальную полосу – 1. Используя принятые в программе обозначения, осуществим ввод числовых данных (в программе запись осуществляем в один столбец):

 5 N=3 FIX(2)=1840.

X(1)=10  FIY(2)=115.

Y(1)=0.5 FIXY(2)=0.

A(1)=20.*1. X(3)=10.

FIX(1)=20.*1.**3/12.  Y(3)=19.45

FIY(1)=1.*20.**3/12. A(3)=23.4

FIXY(1)=0. FIX(3)=113.

X(2)=10.  FIY(3)=1520.

Y(2)=11. 8 FIXY(3)=0.

A(2)=26.8

 Координаты центров тяжести швеллера, двутавра и горизонтальной полосы записаны относительно случайных осей х (рис. 2.1.11) и у/. Ось у/ проходит на расстоянии 10 см от оси у влево.

 Ответ: результаты, выдаваемые ЭВМ на печать: хс = 10 см;

  ус = 10,83см; Iхс = 5828,34 см4; Iус = 2301,67 см4; Iхсус = 0;

 Imax = 5828,34 см4; Imin = Iу = 2301,67 см4; tg2α = 0.

 Задача 9.1.2. Вычислить главные моменты инерции для составного поперечного сечения, представленного на рис. 2.1.12. Найти положение главных осей инерции.

 Ответ: результаты, выдаваемые ЭВМ: хс = 11,7 см; ус = 10,83 см;

 Iхс = 3710,75 см4; Iус = 2065 см4; Iхсус = –382 см4; Imax = 3795 см4;

 Imin = 1981 см4; tg2α = 0,4639.

 Задачи 9.1.3–9.1.10. Применяя приведенную программу на языке Фортран-IV, проверить расчетом на ЭВМ ответы, данные для примеров 2.3.3–2.3.10.

  Задача 9.1.11. Вычислить главные центральные моменты инерции поперечного сечения, показанного на рис. 2.1.13. Найти положение главных осей инерции.

 Ответ: результаты, выдаваемые ЭВМ: х = 7,74 см; ус = 6,76 см;

 Iхс = 2241,75 см4; Iус = 545,46 см4; Iхсус = –480,95 см4;

 Imax = 2368,63 см4; Imin = 418,58 см4; tg2α = 0,56706.

9.2. Построение эпюр прогибов упругой оси балки

  В разделе 4.4 приводится дифференциальное уравнение изгиба упругой оси балки (4.4.1), интегрируя которое можно найти прогиб произвольного поперечного сечения балки. Удовлетворив граничным условиям, находят произвольные постоянные, в результате чего уравнение упругой оси балки можно записать в виде уi = уi(х), где i – число участков, на которые разбивается балка.

 Построим эпюру прогибов балки, которая рассматривается в задаче 4.4.1, рис. 4.4.2 (глава 4). Балка содержит три участка ,  и . Для каждого участка балки получено уравнение изогнутой оси балки (4.4.8).

 Реализуем процесс построения эпюры прогибов упругой оси балки (рис. 4.4.2) и нахождения максимального прогиба в виде программы для ЭВМ на алгоритмическом языке Фортран-IV.

 Для числового примера примем балку из двутавра № 20 (Iz = 1840 см4; Е = 2,1·106 кг/см2) с а = 1 м, l = 2 м, q = 200 кг/м. Общая длина однопролетной балки L = 4 м.

 Ввод числовых данных осуществляется способом «присвоения» с метки (5) до метки (8) включительно. От метки (9) до метки (16) описывается первый участок (), от метки (20) до метки (26) – второй (), а от метки (30) до метки (36) – третий ().

 В программе используются следующие идентификаторы:

текст

а

l

L

x

y

q

E

Iz

C

D

программа

А

В

FL

X

Y

Q

FE

FI

C

D

 Машина будет выдавать на печать значения прогибов с шагом Δх = =20 см. Изменяя значения параметров AN, AL (см. программу) можно увеличить или уменьшить шаг счета Δх.

  PROGRAM BEAM

 5 A=100.

B=200.

FL=400.

Q=2.

FE=2.1E6

FI=1840.

G=FE*FI

C=Q*((FL–A)**4/(24.*FL)–B*FL*FL/12.-A**4/(24.*FL))/G

  8 D=0.

 9 AN=5.

 X=–A/AN

 12 X=X+A/AN

 IF (X–A) 15, 15, 20

 15 Y=Q*B*X**3/(12.*G)+C*X+D

WRITE (7, 50) X,Y

 16 GO TO 12

 20 AL=10.

 X=A–B/AL

 22 X=X+B/AL

 IF (X–A–B) 25, 25, 30

 25 Y=Q*(B*X**3/12. – (X–A)**4/24.)/G+C*X+D

 WRITE (7,50) X,Y

 26 GO TO 22

 30 AN=5

 X=A+B–A/AN

 32 X=X+A/AN

  IF (X-FL) 35, 35, 60

 35 Y=Q*(B*X**3/12.– (X–A)**4/24.+(X–A–B)**4/24.)/G+C*X+D

WRITE (7, 50) X,Y

 36 G0 TO 32

 50 FORMAT (2HX=, F5.1, 3X, 2HY=, E12.4)

  60 STOP

 END

Результаты, выдаваемые ЭВМ на печать:

Х= .0 У= .0000Е+00

Х= 20.0 У= -.1891Е-01

Х= 40.0 У=  -.3741Е-01

Х= 60.0 У= -.5507Е-01

Х= 80.0 У= -.7150Е-01

Х=100.0  У= -.8627Е-01

Х=100.0 У= -.8627Е-01

Х=120.0 У= -.9897Е-01

Х=140.0 У=  -.1092Е+00

Х=160.0 У= -.1168Е+00

Х=180.0 У= -.1214Е+00

Х=200.0  У= -.1229Е+00 (ymax)

Х=220.0 У= -.1214Е+00

Х=240.0 У= -.1168Е+00

Х=260.0 У=  -.1092Е+00

Х=280.0 У= -.9897Е-01

Х=300.0 У= -.8627Е-01

Х=300.0 У= -.8627Е-01

Х=320.0 У= -.7150Е-01

Х=340.0 У=  -.5507Е-01

Х=360.0 У= -.3741Е-01

Х=380.0 У= -.1891Е-01

Х=400.0  У= -.1483Е-07

 Задача 9.2.1. Построить эпюру прогибов консольной балки, нагруженной сосредоточенным моментом m = 300 кг·м (рис. 4.4.6). Балка представляет собой двутавр № 10 (Iz = 198 см4; Е = 2,1·106 кг/см2) с l = 1 м.

 Решение. Используем алгоритм, примененный для составления программы для ЭВМ, рассмотренной в качестве образца (PROGRAM BEAM). Для нашего примера эта программа будет иметь вид:

 PROGRAM BEAM1 Построчные пояснения:

 5 A=100. l = 100 см

  FM=30000. m = 30000 кг·см

 G=2.1E6*198. G = EIz

 C=0. С – произвольная постоянная

 8 D=0. D – произвольная постоянная

 9 U=5. . разбиение 1-го участка на 5 участков

 X=–A/U Δx = l/u – шаг вычислений

  12 X=X+A/U xi = xi-1 + Δx

 IF (X–A) 15, 15, 20

 15 Y=(–FM*X*X/2.+C*X+D)/G  y1=(-mx2/2+Cx+D)/(EIz) – прогиб на 1-м участке

 WRITE (7, 50) X,Y

  16 G0 TO 12

 20 V=5. разбиение 2-го участка на 5 участков

 X=A–A/V

  22 X=A+A/V xi = xi-1+ Δx

 IF (X–2*A) 25, 25, 60

 25 Y=(–FM*X*X/2.+FM*(X–A)**2/2.+C*X+D)/G  y2 = [–mx2/2+m(x–l)2/2+Cx+D]/(EIz)

 WRITE (7, 50) X,Y

 26 GO TO 22

 50 FORMAT(2HX=, F5.1, 2HY=, E12.4)

 60 STOP

 END

Результаты, выдаваемые ЭВМ на печать:

X= .0 Y= .0000E+00 X=100.0 Y= -.3608E+00

X=  20.0 Y= -.1443E-01 X=120.0 Y= -.5051E+00 

X= 40.0 Y= -.5772E-01 X=140.0  Y= -.6494E+00

X= 60.0 Y= -.1299E+00 X=160.0 Y= -.7937E+00

X= 80.0  Y= -.2309E+00 X=180.0 Y= -.9380E+00

X=100.0 Y= -.3608E+00 X=200.0 Y= -.1082E+01

  В вышеприведенной программе применяются размерности: см, кг·см, кг/см2, см4.

  Задача 9.2.2. Составить программу для ЭВМ и построить эпюру прогибов однопролетной балки, показанной на рис. 4.4.7. Принять, что m = =300 кг·м, Е = 2,1·106 кг/см2, l = 1 м Балка – из двутавра № 18.

 У к а з а н и е. Уравнение упругой оси балки взять из задачи 4.4.6.

 Задача 9.2.3. Составить программу для ЭВМ и построить эпюру прогибов консольной балки, изображенной на рис. 4.4.8. Принять q = 1 кН/м, а= 1 м, b = с = 2 м. Балка изготовлена из двутавра № 18. Уравнения изогнутой оси балки для каждого участка взять из ответа к примеру 4.4.7.

  Задача 9.2.4. Составить программу для ЭВМ и построить эпюру прогибов однопролетной балки, показанной на рис. 4.1.17. Пусть F = 1 кН, l = =1м, балка изготовлена из двутавра № 20. Уравнения изогнутой оси балки для двух участков взять из ответа к примеру 4.4.8.

 Задача 9.2.5. Построить эпюру прогибов консольной балки, нагруженной сосредоточенными силами F = 1 кН. Пусть l = 1 м (рис. 4.2.4). Балка круглого поперечного сечения с d = 20 см, Е = 0,1·105 МПа (сосна). Полученные результаты сравнить с ответом в задаче 4.4.10.

 Задача 9.2.6. Построить эпюру прогибов стальной однопролетной балки из двутавра № 18, показанной на рис. 4.4.10. При составлении программы для ЭВМ использовать уравнения изогнутой оси балки, приведенные в ответе к задаче 4.4.12. Пусть l = 4 м, F = 1 кН. Имеется ли симметрия эпюры прогибов относительно оси, проходящей вертикально через сосредоточенную силу 2F?

 Задача 9.2.7. Построить эпюры прогибов и углов поворота сечений  стальной балки из двутавра № 20, показанной на рис. 4.4.11, где а = 1 м, b= 0,8 м; F = 1 кН. Результаты сравнить с ответом к задаче 4.4.13.

 Задача 9.2.8. Имеется стальная балка из двутавра № 22, нагруженная сосредоточенной силой F = 30 кН (рис. 4.5.1). Удовлетворяет ли сортамент балки условию жесткости (4.5.1), если [1/n0] = 1/250? Материал балки – сталь С255.

 У к а з а н и е. Предварительно необходимо построить эпюру прогибов и определить ymax.

 Ответ: балка удовлетворяет условию жесткости (4.5.1).

 

 

Аналитический расчет кривых брусьев малой кривизны


Расчет кривых брусьев малой кривизны рассматривался в разделе 5.4. Предложенная в примере 5.4.1 методика построения эпюр изгибающих моментов, поперечных и нормальных сил легко реализуется в виде программы для ЭВМ.

 Например, составим программу на алгоритмическом языке ПЛ-1 для расчета круговой трехшарнирной арки, изображенной на рис. 9.3.1.

ARCA: PROCEDURE OPTIONS (MAIN);

 /*КРУГОВАЯ АРКА РАДИУСОМ R*/

GET LIST (F,FL,A,B,C,D,Q1,Q2,Q3,Q4,P1,P2,P3,H);

PUT SKIP EDIT (‘КРУГОВАЯ АРКА РАДИУСОМ R’)(X(10),A);

PUT SKIP DATA (F,FL,A,B,C,D,Q1,Q2,Q3,Q4,P1,P2,P3,H);

DECLARE  F,FL,A,B,C,D,Q1,Q2,Q3,Q4,P1,P2,P3,H,VA,VB,HH,R,X,Y,TGFI,SINFI,COSFI,

  FMO, FQO,FM,FQ,FN;

VA=(Q1*A*(2.*FL–A)+P1*(FL+2.*B)+Q2*B*(FL+B)+P2*FL+Q3*C*(FL–C)+2.*P3*D+

 Q4*D*D)/(2.*FL);

VB=(Q1*A*A+2.*P1*A+Q2*B*(FL–B)+P2*FL+Q3*C*(FL+C)+2.*P3*(FL–D)+

  Q4*D*(2.*FL–D))/(2.*FL);

HH=(VB*FL/2. –P3*C–Q4*D*(C+D/2.) –Q3*C*C/2.)/F;

R=FL*FL/(8.*F)+F/2.;

PUT SKIP DATA (VA,VB,HH,R);

DO X=0. TO A BY H;

FMO=VA*X–Q1*X*X/2.;

FQO=VA–Q1*X;

  CALL TR; CALL REZ; END;

DO X=A TO A+B BY H;

FMO=VA*X-Q1*A*(X–A/2.)-P1*(X–A) –Q2*(X–A)**2/2.;

FQO=VA–Q1*A–P1–Q2*(X–A);

 CALL TR; CALL REZ; END;

DO X=FL/2. TO FL/2.+C BY H;

FMO=VA*X–Q1*A*(X–A/2.) –P1*(X–A) –Q2*B*(X–AB/2.)-P2*(X-FL/2.)-Q3*(X-FL/2.)**2/2.;

FQO=VA-Q1*A-P1-Q2*B-Q3*(X-FL/2.)-P2;

  CALL TR; CALL REZ; END;

DO X=FL–D TO FL BY H;

FMO=VB*(FL–X) –Q4*(FL–X)**2/2.;

FQO=–VB+Q4*(FL–X);

  CALL TR; CALL REZ; END;

TR: PROCEDURE;

Y=F–R+SQRT(R*R– (X-FL/2.)**2);

TGFI=(FL/2. –X)/SQRT(R*R– (X–FL/2.)**2);

COSFI=SQRT(1./(1.+TGFI**2));

SINFI=TGFI*COSFI;

END TR;

REZ: PROCEDURE;

FQ=FQO*COSFI–HH*SINFI;

FM=FMO–HH*Y;

FN=-FQO*SINFI–HH*COSFI;

PUT SKIP DATA (X,FM,FQ,FN);

END REZ;

END ARCA;

 В программе применены следующие идентификаторы:

Текст

l

f

q

φ

H

Mz

Qy/

N

Δx

Программа

FL

F

Q

FI

FMO

FQO

HH

FM

FQ

FN

H

 Алгоритм построения эпюр изгибающих моментов, поперечных и нормальных сил можно реализовать также на алгоритмическом языке Бейсик. Например, для параболической трехшарнирной арки, изображенной на рис. 9.3.1, программа на этом языке для персональной ЭВМ примет вид:

 5 OPEN “ARCA05.DAT” FOR OUTPUT AS FILE#1

 10 PRINT ‘Расчет параболической трехшарнирной арки’

  20 PRINT ‘Ввести F,FL,A,B,C,D,Q1,Q2,Q3,Q4,P1,P2,P3’

 30 INPUT F,FL,A,B,C,D,Q1,Q2,Q3,Q4,P1,P2,P3

  40 VA=(Q1*A*(2.*FL–A)+2*P1*(FL–А)+Q2*B*(FL+B)+P2*FL+Q3*C*(FL–C)+2.*P3*D+

  Q4*D*D)/(2.*FL)

 50 VB=Q1*A+P1+Q2*B+P2+Q3*C+P3+Q4*D-VA

 60 H=(VB*FL/2. –P3*C-Q4*D*(C+D/2.) –Q3*C*C/2.)/F

 70 PRINT#1, “ Расчет параболической трехшарнирной арки”

 80 PRINT#1, “VA=”,VA,”VB=”,VB,”H=”,H

 100 PRINT#1, “Таблица значений изгибающих моментов, поперечных”

 110 PRINT#1, “ и нормальных сил”

 120 X=-1.0

 125 X=X+1.

 130 TG=4.*F*(FL–2.*X)/(FL*FL)

  140 COS=SQRT(1./(1.+TG*TG))

 150 SIN=TG*COS

 160 Y=4.*F*X*(FL–X)/(FL*FL)

  170 IF X>A GO TO 205

 180 MO=VA*X–Q1*X*X/2.

 190 Q0=VA–Q1*X

  200 GO TO 320

 205 AA=0.5*FL

 210 IF X>AA GO TO 245

 220 МО=VA*X-Q1*A*(X–A/2.) –P1*(X–A)-Q2*(X–A)**2/2.

 230 QO=VA–Q1*A–P1–Q2*(X–A)

  240 GO TO 320

 245 AA=A+B+C

 250 IF X>AA GO TO 290

 260 MO=VA*X–Q1*A*(X–A/2.) –P1*(X–A) –Q2*B*(X–A–B/2.) –P2*(X–FL/2.) –

 Q3*(X–FL/2.)**2/2.

  270 QO=VA–Q1*A–P1–Q2*B–Q3*(X–FL/2.) –P2

 280 GO TO 320

 290 IF X>FL GO TO 370

 300 MO=VB*(FL-X)-Q4*(FL-X)**2/2.

 310 QO= -VB+Q4*(FL-X)

  320 M=MO-H*Y

 330 Q=QO*COS–H*SIN

 340 N=–QO*SIN–H*COS

 350 PRINT#1, “X=”,X,”Y=”,Y,”M=”,M,”Q=”,Q,”N=”,N

 360 GO TO 125

 370 STOP

 380 END

 Здесь применены следующие идентификаторы:

Текст

l

f

q1

cosφ

sinφ

tgφ

Mz

Qy/

Программа

FL

F

Q1

COS

SIN

TG

MO

QO

M

Q

 Затем необходимо дополнительно вычислить усилия Q и N в местах приложения сосредоточенных сил, причем определять Q и N следует в сечениях справа от сил. ЭВМ выдает на печать значения Q, N в сечениях слева от сосредоточенной силы.

 Порядок построения эпюр изгибающих моментов, поперечных и нормальных сил для эллиптической арки, изображенной на рис. 9.3.2, реализуем на алгоритмическом языке Фортран-IV.

  PROGRAM ARCA01

C Расчет эллиптической трехшарнирной арки 

 TYPE ¤, ‘введите стрелу подъема F=’

 ACCEPT ¤, F

 TYPE ¤, ‘введите пролет арки FL=’

 ACCEPT ¤, FL

 TYPE ¤, ‘введите длину первого участка А=’

 ACCEPT ¤, А ………… и т.д. для B,C,D,Q1,Q2,Q3,Q4,P1,P2,P3

 WRITE (6,102)

102 FORMAT (5X, ‘Расчет эллиптической трехшарнирной арки’/)

  WRITE (6,103)

103 FORMAT (30X, ‘Опорные реакции’/)

 WRITE (6,106)

106 FORMAT (5X, ‘X’,5X,’Y’,10X,’M’,13X,’Q’,15X,’N’/)

 VA=(Q1*A*(2.*FL–A)+P1*(FL+2.*B)+Q2*B*(FL+B)+P2*FL+Q3*C*(FL–C)+2.*P3*D+

 Q4*D*D)/(2.*FL)

 VB=Q1*A+P1+Q2*B+P2+Q3*C+P3+Q4*D–VA

 H=(VB*FL/2. –P3*C–Q4*D*(C+D/2.) –Q3*C*C/2.)/F

 WRITE (6,104) VA, VB, H

104 FORMAT (5X,’VA=’,F8.4,5X,’VB=’,F8.4,’H=’,F8.4/)

 X=0.

 2 X=X+1.

  IF (X.EQ.FL) GO TO 10

 Y=F/FL*SQRT(FL*FL–4.*(X–0.5*FL)**2)

 TGFI=(F/FL)**2*4.*(FL/2. –X)/Y

 COSFI=SQRT(1./(1.+TGFI**2))

 SINFI=TGFI*COSFI

 IF (X-A) 5,5,6

 5 FMO=VA*X–Q1*X*X/2.

 FQO=VA–Q1*X

 G0 TO 12

  6 IF (X–A–B) 7,7,8

 7 FMO=VA*X-Q1*A*(X–A/2.)-P1*(X–A)-Q2*(X–A)**2/2.

  FQO=VA–Q1*A–P1–Q2*(X–A)

 G0 TO 12

 8 IF (X-A-B-C) 9,9,10

  9 FMO=VA*X–Q1*A*(X–A/2.) –P1*(X–A) –Q2*B*(X–A–B/2.) –P2*(X–FL/2.) –

 Q3*(X–FL/2.)**2/2.

  FQO=VA–Q1*A–P1–Q2*B–Q3*(X–FL/2.) –P2

 GO TO 12

 10 IF (X–FL) 11,14,14

  FMO=VB*(FL–X) –Q4*(FL–X)**2/2.

 FQO= –VB+Q4*(FL–X)

 12 FM=FMO–H*Y

  FQ=FQO*COSFI–H*SINFI

 FN=–FQO*SINFI–H*COSFI

 PRINT ¤, X,Y,FM,FQ,FN

  G0 TO 2

 14 STOP

 END

 В сечении арки х = 0 м, т.е. на опоре А имеем у = 0, tgφ =, φ = π/2, cosφ = 0, sinφ = 1. Следовательно, по формулам (5.4.3) находим М(х=0) = 0, Q(х = 0) = –Н, N(х = 0) = –VА. Аналогично в сечении х = l, т.е. на опоре В имеем у = 0, φ = –π/2, cosφ = 0, sinφ = –1, и по формулам (5.4.3) находим М(х = l) = 0, Q(х = l) = H, N(х = l) = –VВ.

 Затем необходимо вычислить усилия Q и N в местах приложения сосредоточенных сил, причем определять Q и N следует в сечениях справа от сил. ЭВМ выдает на печать значения Q, N в сечениях слева от сосредоточенной силы.

 Идентификаторы для программы на языке Фортран аналогичны идентификаторам для программы на языке ПЛ-1.

 Задача 9.3.1. Построить эпюры изгибающих моментов Мz, поперечных Qу/ и нормальных N сил для трехшарнирной круговой арки, показанной рис. 5.4.1, a.

 У к а з а н и е. Для расчета можно использовать любую из трех предложенных программ. Программы на языке ПЛ-1 применять без каких-либо изменений. В программах на языках Бейсик и Фортран необходимо заменить уравнение оси арки на уравнение окружности (5.4.4), а значение tgφ дать по формуле (5.4.5).

 Ответ: эпюры Мz, Qу/, N приведены на рис. 5.4.1, г.

 Задача 9.3.2. Построить эпюры изгибающих моментов Мz, поперечных Qу/ и нормальных N сил для трехшарнирной параболической арки, по-

казанной рис. 5.4.3. Ось параболической арки очерчена по кривой

y = 4fx(l – x)/l2,  а tgφ = dy/dx = 4f(l – 2x)/l2.

 У к а з а н и е. Для расчета можно использовать без каких-либо изменений программу для ЭВМ на языке Бейсик. При применении предложенных программ на языках ПЛ-1 или Фортран необходимо заменить в них уравнение оси арки на уравнение параболы, данное в условии задачи и, кроме того, поставить соответствующее значение tgφ.

 Ответ: эпюры Мz, Qу/, N приведены на рис. 5.4.3.

 Задача 9.3.3. Построить эпюры изгибающих моментов Мz, поперечных Qу/ и нормальных N сил для трехшарнирной эллиптической арки, показанной рис. 5.4.4. Ось эллиптической арки очерчена по кривой

  а 

 У к а з а н и е. Для расчета можно использовать без каких-либо изменений предложенную программу для ЭВМ на языке Фортран.

 Ответ: эпюры Мz, Qу/, N приведены на рис. 5.4.4, которые построены на основании таблицы значений изгибающих моментов, поперечных и нормальных сил, выданной ЭВМ на печать.

Таблица значений изгибающих моментов, поперечных и нормальных сил для эллиптической арки, показанной на рис. 5.4.4

x

y

M

Q

N

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

2.397916

3.316625

3.968627

4.472136

4.873397

5.196152

5.454356

5.656854

5.809475

5.916080

5.979130

6.000000

5.979130

5.916080

5.809475

5.656854

5.454356

5.196152

4.873397

4.472136

3.968627

3.316625

2.397916

-25.00936

-24.17419

-21.13823

-18.20665

-16.28125

-15.82497

-15.60995

-14.30865

-12.03476

-8.863564

-4.843040

-0.000000

-4.343040

-7.863564

-10.53476

-12.30865

-13.10995

-12.82497

-21.28125

-28.20665

-33.13823

-35.17419

-32.00936

-1.388797

2.037357

2.778986

2.308956

1.174283

-0.3642907

0.7613568

1.775328

2.710070

3.589154

4.430793

5.250000

-3.931229

-3.090930

-2.214185

-1.282963

-0.2739916

0.8446751

-7.291819

-5.504004

-3.213944

-4.0983200E-02

5.003451

-28.07817

-25.95499

-23.99145

-22.39378

-21.14080

-20.19108

-20.18001

-20.11618

-20.01170

-19.87286

-19.70230

-19.50000

-19.68141

-19.83075

-19.94768

-20.02914

-20.06832

-20.05241

-23.33756

-23.82265

-24.23805

-24.45017

-23.93278

 Задача 9.3.4. Построить эпюры изгибающих моментов Мz, поперечных Qу/ и нормальных N сил для трехшарнирной параболической арки, показанной на рис. 9.3.3. Уравнение параболической оси арки, значение tgφ и указания к расчету приведены в задаче 9.3.2.

 Ответ: VA = VВ = 120 кН; Н = 120 кН; Mс = Qс = 0; Nс = –Н.

 Задача 9.3.5. Построить эпюры изгибающих моментов Мz, поперечных Qу/ и нормальных N сил для трехшарнирной круговой арки, показанной на рис. 9.3.4. Уравнение круговой оси арки задано в виде (5.4.4), значение tgφ вычислить как dy/dx. Указания к расчету приведены в задаче 9.3.1.

 Ответ: V = 5кН; Н = 7,5 кН; Мс = 0, R = 6,5 м; Nс = –7,5 кН; Qс=5 кН.

 Задача 9.3.6. Построить эпюры изгибающих моментов, поперечных и нормальных сил для трехшарнирной параболической арки, показанной на рис. 9.3.5. Уравнение параболической оси арки, значение tgφ и указания к расчету приведены в задаче 9.3.2.


Ответ: VA = 90 кН; VB = 30 кН; Н = 60 кН; Qс = –30 кН; Nс = –Н.

Непосредственное интегрирование

Отыскание неопределенного интеграла с помощью таблицы, правил и тождественных преобразований называют непосредственным интегрированием.

Пример 1.

При интегрировании использованы правила 2 и 3, а также табличные формулы 2,4,6,11.

Пример 2.

При интегрировании использованы тождественные преобразования подынтегральной функции, правила 2 и 3, табличная формула 2.


Лабораторный практикум по сопромату