Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Фролов, С. А. Кибернетика и инженерная графика

.pdf
Скачиваний:
15
Добавлен:
20.10.2023
Размер:
7.75 Mб
Скачать

точки Ац\).

Затем на образе 6 берем точку

А\

<2> и выпол­

няем для нее все отмеченные

выше

построения.

После

точки Л! (2)

переходим к Ац3),

A \ { i ) ,

. . .,

А\

( п )

и т. д.,

пока не будут исчерпаны все точки образа 6'. В резуль­

тате этих построений

получим

преобразованный вид кри-

->

1 — I

1

1—I

I I I г

Рис. 66

вой 2 (на рис. 66 кривая показана штриховой линией).

Находим точку

К2 пересечения

преобразованной

кривой

с образом / (операторы XI,

I ,

I ,

V). На

образе 2

опреде­

ляем точку К2

(операторы

X,

I ,

V), по /С3 находим гори­

зонтальную проекцию

 

К. г

(операторы

X, I ,

V).

Схема

счета решения задачи

имеет следующий вид:

 

 

 

Ш

Л

;

 

W,V*Xa

 

 

 

tioVnVI.Ji8VIai

VI6n,

X 1 6 / 1

7 / 1

8 y i 9 X 2 0 / 2

1 V 2 2 X 2 3

/ 2 4 V 2 5 . (49)

100

Взод исходных данных

С л е ж е н и е

Формирование переменных команд, зависящих от адресов первых точек образов

Обращение к операторам X, 1,Y (определение точки В? на образе 1;

Обращение

к

оператору

I

 

 

(определение

коэффициента

прямой

B 2 - S 2 )

Выборка на

рабочие

адреса

6

 

первой (очередной) точки

образа

 

Проба:на

конец

 

 

 

 

 

точек

образа

6

к

нет

 

I

 

 

Обращение

оператору

 

 

Определение коэффициента прямой,

проходящей

через точку

S,

и

т о ч к и образа 6

- Ацц)

Обращение к операторам X . I . Y

 

(определение

точки Сц^браза

4}

 

Обращение

к

операторам

 

Y

 

(опредепение

точки А г (

образа

2

)

Обращение

к о п е р а т о р у Ш а

 

С 2 )

(определение

фронтальной

проекции

Обращение

к

оператору

I

C7-S2)

(определение

коэффициентов

прямой

Сг имеет

координаты: х С г = х с , ; у с , =

у в 2

Обращение

к

опеоатооу

YJa

 

 

(определение

точки

А ' г ш )

 

 

Обращение

н

оператору

YIj

 

 

(определение

точки

А

ад)

 

 

Определение

точки

Д г , „ х д г ( „ - <д'2

и Ул г ( 1 ) = уд г

Формирование

массива

Г

Агтл

 

Рис. 67

1

Обрашенне

к операторам

H.1.I.Y

(определение точки Кг!

Обращение

к операторам

Х,1,¥

((определение точки Кг)

Обращение

коператорам

X . I . Y

(определениетточки К,)

Стоп

101

На рис. 67 приведена блок-схема программы для реше­ ния задачи по нахождению точки встречи произвольной

кривой с поверхностью

эллиптического

параболоида.

 

 

Приведенная

 

программа

 

 

может быть использована для

 

 

решения

задач по

определе­

 

 

нию линии

пересечения двух

 

 

поверхностей

с

криволиней­

 

 

ными

образующими,

если се­

 

 

чения

параллельными

пло­

 

 

скостями хотя

бы

одной по­

 

 

верхности

будут

 

подобны и

 

 

подобно

расположены.

При­

 

 

мером таких

задач

служит

 

 

нахождение

линии

пересече­

 

 

ния

кольцевого

 

трубопро­

 

 

вода

с

патрубком,

который

в,

 

представляет

произвольную

 

 

поверхность канала (рис. 68).

 

 

Для

решения

этой

задачи

 

 

достаточно

после

 

выполне­

 

 

ния массива

( Х 2 3 ,

/ 2 4 ,

V2&)

 

 

предусмотреть

операцию про­

 

Рис. 68

ведения

прямой,

параллель­

 

ной

заданной,

проходящей

через

следующую (из

записанных

в

б-массиве

обра­

за /)

точку кривой

/ х (оператор

/ / ) .

 

 

 

 

 

§ 11. РЕШЕНИЕ ЗАДАЧ С ИСПОЛЬЗОВАНИЕМ СВОЙСТВ ЦЕНТРА ТЯЖЕСТИ

Использование некоторых понятий о центре тяжести позволяет получить изящные графические решения мно­ гих задач. Определение положения центра тяжести лю­ бого геометрического образа может быть сведено к задаче по нахождению объединения двух материальных точек. Под объединением двух материальных точек подразуме­ вается такая новая материальная точка, носителем которой является центр тяжести данных материальных точек, а масса равна сумме масс этих материальных точек. Дей­ ствительно, пусть требуется определить центр тяжести пластинки криволинейного контура (рис. 69). Эту пла­ стинку можно с некоторой степенью точности рассматри­ вать как «многоугольник», составленный из прямоуголь-

102

ников. Положение центра тяжести пластинки может быть определено путем последовательного объединения мате­ риальных точек, носителями которых являются центры тяжести прямоугольников /, // , / / / — точки Z1 , Z1 1 , Z 1 1 1 . Если сосредоточить массы прямоугольников, которые чис­ ленно равны их площади, в центре тяжести, то центр тя­

жести фигуры, состоящей из двух

прямоугольников / +

+

//.. может

быть

найден

объединением

материальных

точек:

 

 

 

 

 

 

 

 

 

 

 

Z1 (S!);Z" (Sll)

= Z±

(S1 +

SU).

 

 

 

 

 

 

Точка

Z x

лежит

на

от­

 

 

 

 

 

резке Z'Z1 1 , ее положение

 

 

 

 

 

определяется

условием

 

 

 

 

 

 

 

 

 

Ztzl

 

=

 

 

 

 

 

 

 

 

 

 

 

zxzn

 

s1

'

 

 

 

 

 

 

 

 

Для

определения

центра

 

 

 

 

 

тяжести

фигуры,

образован­

 

 

 

 

 

ной

прямоугольниками

/ -+-

 

 

 

 

 

+

/ / -4- / / / , достаточно найти

 

 

 

 

 

объединение

материальных

 

Рис. 69

 

 

точек

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Zj (S1

+ S»);

Z I U (Sl ») = Z2

(S1 + S! I

+ S1 1 1 ).

 

 

Так же находим точку Z 3

(S1 +

S" + Sm +

SI V ), ко­

торая

будет

 

объединением материальных

точек

Z 2

(S1 +

+

S n + S I n )

и Z I V

(SI V ). Аналогично можно найти

точки

Z4 ,

Z 5 . . .

и, наконец,

точку Z0 ,

которая

будет центром

тяжести

пластинки.

 

 

 

 

 

 

 

 

Описанный в § 3 метод считывания чертежа

дает для

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

тура

(рис. 70) будет представлена

в виде значений коор­

динат

точек /, 2,3,4,5,6,7...,

т. е. в процессе считы­

вания происходит расчленение площади, ограниченной замкнутой кривой, на следующие один за другим прямо­ угольники, у которых высота равна шагу квантования по уровню 0,4 мм, а длина — величине разности абсцисс точек считывания.

103

Рассматривая эти прямоугольники как «отрезки», легко найти положение их центров тяжести:

 

 

 

 

^л, п+1 — '

 

 

 

 

 

 

 

 

Вычисление значения

Z„ „ + 1

легко

осуществляется на

ЭЦВМ. Зная положение

центра тяжести каждого из этих

 

 

 

 

 

 

 

«отрезков»,

можно опре­

 

 

 

 

 

 

 

делить

объединение

ма­

 

 

 

 

 

 

 

териальных

точек,

т. е.

 

 

 

 

 

 

 

центр

 

тяжести

фигуры,

 

 

 

 

 

 

 

образованной

«отрезка­

 

 

 

 

 

 

 

ми» 12 и 3—4.

 

Для

 

 

 

 

 

 

 

этого

достаточно

в точ­

 

 

 

 

 

 

 

ках Z1

2 и

Z 3

4

поме­

 

 

 

 

 

 

 

стить

 

массы «отрезков»

 

 

 

 

 

 

 

( m l i 2

= х% — хх

и m3 i 4 =

 

 

 

 

 

 

 

=

х4—х3)

и на отрезке

 

 

 

 

 

 

 

ZB

2 Z 3 i 4

найти точку Z,

 

 

 

Рис.

70

 

удовлетворяющую усло-

 

 

 

 

 

ZZX,

 

 

 

Зная

 

 

 

 

 

 

 

вию ZZ,

 

 

 

положение

точки

Z и ее массу Ь2

+

m 3 i 4 ) , определяем

положение

точки Z x

как объединение материальных точек

Z (mx 2

+

tfz3,4) и 2 5 i в

(m6 i

в ) . Следуя таким путем,

можно

найти

Z 0

— положение центра

тяжести'^пластинки.

 

 

Рассмотренный

способ позволяет

определить

положе­

ние центра тяжести

и в том случае,

если пластинка

имеет

отверстие (рис. 71). Также находится центр тяжести от­ резка /—2 — точка Zx 2 , затем отрезка 3—4 (определяе­ мого двумя точками, записанными в следующих порядко­

вых номерах ячеек Нф)

— точка

Z 3 4 .

Определяется

объ­

единение

материальных

точек:

 

 

 

 

Zi,2

[ЩлУ, Z3A

(Щ,ь)

= z

(ть

2 + Щл)-

 

Зная

Z (mh

2 + Щ,4 )

и Z5 , e

(m5 j „), находим Z x

2 +

+Щ, 4 + т5, б) И Т. Д.

Решая задачу по определению центра тяжести пла­ стинки, мы одновременно находим и площади любой пло­

ской

фигуры, для чего полученную разность абсцисс

точек,

расположенных в одной строке и следующих друг

за

другом, умножаем на 2, а не на 0,5, как это делается

 

* Мы рассматриваем случай, когда линейная плотность постоянна

по

всей

длине «отрезка».

104

при определении центра тяжести. То же делаем со сле­ дующими разностями и полученные результаты сумми­ руем. Зная положение центров тяжести пластинок, на

Рис. 71

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

торых находят) берем центры

 

 

 

 

тяжести

пластинок

(носи­

 

 

 

 

тель) и

их площадь

(массы).

 

 

 

 

На рис. 72 показан геометри­

 

 

 

 

ческий

образ

произвольной

 

 

 

 

формы.

Рассекая

его парал­

 

 

 

 

лельными

.

плоскостями

а 1 ,

 

 

 

 

а " , а

ш ,

.

.,

а " ,

получим

 

 

 

 

п +

1-сечение. Если

расстоя­

 

 

 

 

ние

между

секущими

пло­

 

 

 

 

скостями

будет

достаточно

 

 

 

 

малым,

то

части

со1,

со11,

 

 

 

 

со1 1 1 ,. ..

с

определенной

сте­

 

 

 

 

пенью точности можно отож­

 

 

 

 

дествлять

 

с

пластинками,

 

 

 

 

положение

 

центра

тяжести

 

 

 

 

которых

 

Z\,

Z",

Zo \

. • •

 

 

 

 

легко определяется. Сосредо­

 

 

 

 

точив массы

отсеченных

ча­

 

 

 

 

стей

со1

и

со"

в

их

центрах тяжести Z\

 

и

Z0

находим

объединение

Z\

(со1);

Z " (со")

=

Z x

(со' + со11).

Для определения центра тяжести геометрического об­

раза,

состоящего

из

пластинок

(частей)

со1

+ со11 + со1 1 1 ,

достаточно

 

найти

 

 

 

 

 

 

 

 

 

2 2

( с ' +

о," +

со ш ) =

Zi (со1 + со11);

Zl11

(со1 1 1 )-

Аналогично находим Z 3

и т. д. до Z0 .

 

 

 

105

Суммируя величины со(-, можно подсчитать объем гео­ метрического образа.

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

Проба на конец точек контура

г

Замещение в рабочих ячейках координат точки Z,>2 координатами точки Z, (далее Z 2 , Z 3 . . . j

 

 

г

 

 

Вычисление (xN + 4-xN + 3 )

*

 

H( XN + 2 -Xn + ,)+(XN-Xn-i)

 

4

 

 

 

Вычислениех, **

ч г

 

 

z l 2

 

 

Т

>

. f

 

Вычисление х ? м ; X Z s 6 , . . Выборкау Z j 4 ; y Z 5 6 ; . .

Вычисление расстояний

между точками Z,2 и Z 3 4 ; Z , и Z 5 4 ; Z 2 h Z 7 8

Определение координат точек Z, (Z?;Z3;...)

делящих отрезки в нужном отношении

*) N=0;2;4;б — если N=0,то вторая скобка считается нулем

**) выполняется только при первом цикле

Рис. 73

необходимо лишь определить положение центра тяжести каждого из «отрезков» по отношению к центру тяжести всей пластины (сечения). На рис. 73 приведена блок-схема

106

программы для решения задачи по определению центра тяжести пластины. Программа может быть легко при­ способлена для решения задачи по подсчету величины площади пластинки или объема любого геометрического тела.

Использование свойств центра тяжести и, в частности, теорем Гюльдена позволяет решать на ЭЦВМ задачи по нахождению площади поверхности и объема тел враще­ ния, если известны характер линии или форма пластинки, образующих при своем вращении эти поверхности или тела. Алгоритм для решения таких задач выражается простыми формулами, которые вытекают непосредственно из формулировок теорем.

Первая теорема Гюльдена. Если поверхность образо­ вана вращением некоторой линии вокруг оси, причем линия лежит в одной плоскости с осью и целиком по одну сторону от оси, то площадь поверхности равна произве­ дению длины линии на длину окружности, описанной центром тяжести линии. Поэтому, если S — площадь по­ верхности, / — длина линии и R — расстояние центра тя­ жести линии от оси, то S = 2nRl.

Вторая теорема Гюльдена. Если тело образовано вра­ щением некоторой плоской фигуры (пластинки) вокруг некоторой оси, причем пластинка лежит в одной плоскости с осью и целиком по одну сторону от оси, то объем тела равен произведению площади пластинки на длину окруж­ ности, описанной центром тяжести пластинки, следова­ тельно, V = 2nRS. Выше было показано, как машина самостоятельно может определить положение центра тя­ жести и величину площади пластинки, пользуясь только чертежом. Нахождение длины линии сводится к много­ кратному выполнению оператора VII (определить рас­ стояние между двумя точками) и суммированию получен­ ных результатов. Определение величины R (при задан­

ных — положении центра тяжести

и

а — оси вращения)

реализуется с помощью операторов

/

/ / , V, VII.

К РАСЧЕТУ П Л А В У Ч Е С Т И , ОСТОЙЧИВОСТИ

ИНЕПОТОПЛЯЕМОСТИ К О Р А Б Л Я

ЭЦВМ целесообразно использовать для решения задач, в которых приходится выполнять большое количество однотипных вычислений. К таким задачам относится ста­ тический расчет корабля.

107

Основные мореходные качества, исследуемые в теории

корабля (в разделе

статика): плавучесть,

остойчивость

и непотопляемость. Существующие

методы

расчета этих

характеристик

требуют

вычисления

следующих величин:

1)

площадь

ватерлинии

 

 

 

 

 

 

 

 

S = 2

J

ydx;

 

(50)

 

 

 

 

 

 

 

т

 

 

2)

площадь

шпангоута

 

т

 

 

 

 

 

 

 

 

 

 

 

 

 

со =

2 J

г/dz;

 

(51)

 

 

 

 

 

о

 

 

 

3)

статический момент площади ватерлинии

 

 

Мк

= 2

J

xydx;

 

(52)

4)

статические моменты площади

шпангоутов

 

 

 

 

 

1

г

 

 

 

 

 

 

b =

 

 

 

(53)

 

 

 

\\t?bdz;

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

С

\yzdz;

 

(54)"

 

 

 

 

 

о

 

 

 

5)

моменты

инерции

площади ватерлинии

 

 

J x

=

T J

 

y 3 d x '

 

( 5 5 )

 

 

 

 

 

 

L _

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

+

2

 

 

 

 

 

Jy = 2

\

x2ydx;

 

(56)

2

* Мы сохраняем обозначения, принятые в теории корабля, в ча­ стности, в курсе «Статика корабля». В. В. Семенов-Тянь-Шаньский, Судпромгиз, 1940.

108

6) два выражения для водоизмещения корабля:

 

V = + j4

о» dx\

(57)

L

 

 

2

 

 

Z

 

 

V=\

Sdz-

(58)

о

7) статический момент погруженного объема относи­ тельно основной плоскости

Z

 

Муг = j Szdz;

(59)

о

 

8) два выражения статического момента погруженного

объема относительно плоскости

миделя:

 

+

4

 

 

МУ2=

\

coxdx-

(60)

_

L _

 

 

 

2

 

 

 

Z

 

 

Myz=\Sxdz.

 

(61)

о

 

 

Все приведенные величины выражаются определен­ ными интегралами, вычисление которых связано с труд­ ностями, так как корпус корабля образует, как правило, трансцендентную поверхность, которая не может быть за­ дана аналитически. Сведения о ней можно получить лишь по теоретическому чертежу. Поэтому при вычислении интегралов приходится пользоваться приближенными спо­ собами, используя для расчета исходные данные, снятые непосредственными измерениями с чертежа. Насколько трудоемким является процесс статического расчета ко­ рабля, можно судить по приведенным выше величинам, которые приходится определять много раз для всех сече­ ний по шпангоутам, всех уровнях ватерлиний и для раз­ личных равнообъемных наклонениях корабля. Указан­ ные интегралы необходимы для того, чтобы определить водоизмещение корабля (V), координаты центра тя­ жести (ЦТ) и центра величины (ЦВ) при различном уровне и положении ватерлинии.

109