ПРОГРАММА РАСЧЕТА ГЕОМЕТРИЧЕСКИХ ХАРАКТЕРИСТИК ПОПЕРЕЧНЫХ СЕЧЕНИЙ С УЧЕТОМ ДВИЖЕНИЯ ФРОНТА ПОВРЕЖДЕНИЯ
Предполагается, что при числе точек излома контура поперечного сечения, равном четырем /m=5/, вычисления по приводимым в данных методических указаниях форму .там должны производиться с помощью калькулятора. Однако при использовании калькулятора не исключается возможность появления не только погрешностей, но и грубых ошибок в вычислениях, которые сделают результаты качественно неверными. Кроме того, уже при т=5 весьма большое время тратится на ручной счет, а при больших значениях т, использование ПЭВМ просто необходимо.
В связи с этим величины геометрических характеристик плоских сечений вычисляются с помощью программы, составленной для ПЭВМ ГВМ РК 386 на языке Quick BASIC 4.5. При использовании данной программы необходимо помнить, что точки излома контура нумеруются против движения часовой стрелки.
В программе приняты следующие обозначения:
М- число, на единицу превышающее число точек излома контура, U(I), V(I) (I=1,2,…,M) - декартовы координаты точек излома контура, соответствующие X и Y /СМ/,
G(I), (I=1,2,... ,М) -угол, составляемый вектором X с вектором скорости распространения повреждения, 0≤G(I)≤360,
H - время в годах, прошедшее с момента начала процесса разрушения ,
W(I), (I=1,2,...,М) - скорости распространения повреждений /СМ/год/,
F1 - величина площади-сечения /СM2 /,
Х1 - координата центра тяжести сечения /СМ /,
Y1- координата центра тяжести сечения /СМ /,
J1(1)=JX1- величина осевого момента инерции /СМ4/,
J1(2)=JY1- величина осевого момента инерции /СМ4/,
J1(3)=JX1Y1- величина центробежного момента инерции /СМ4/,
E1=JU - величина первого из экстремальных моментов инерции/СМ4/,
E2=JV - величина второго из экстремальных моментов инерции /СМ4/,
E3=αº - угол наклона главных центральных осей к осям X, Y /градусы/,
SW- указатель. При задании величины SW< 0 величинам W(I), G(I), (I=1,2,…,M) присваиваются нулевые значения, то есть не учитывается воздействие агрессивной среды. Если же SW>0, то величины W(I), G(I), (I=1,2,…,M) вводятся с клавиатуры,
WS- указатель. При задании величины WS<0 программа прекращает определение характеристик сечения и выходит на ввод координат точек излома контура нового поперечного сечения. Если же WS>0, то программа запрашивает величину времени с начата процесса разрушения уже введенного сечения.
Для работы с программой необходимо подготовить и ввести с клавиатуры следующие исходные данные:
I. М - число, на единицу превышающее число точек излома контура,
2,3. U(I), V(I) (I=1,2.....М) - декартовы координаты точек излома контура /СМ /,
4. SW величину указателя, описания которого приведены в первых пяти строках данной страницы методических указаний.
5,6. Только при SW>0 вводятся W(I), G(I), (I =1,2,... ,М) - скорости /СМ/год / и углы /градусы/ распространения повреждений.
7. H- время с начала процесса разрушения /годы/.
8. Только после печати результатов для конкретного поперечного сечения для определенного времени H вводится величина указателя WS определяющего дальнейший ход работы программы. Описания указателя WS приведены в строках 6 - 10 данной страницы методических указаний.
На приводимой ниже схеме алгоритма цифрами обозначены:
2 - ввод величины М .
5- ввод величин U(I), V(I) в цикле по I,
6 - ввод величины указателя SW,
8 - ввод величин W(I), G(I) в цикле по I,
10 - печать величин I, U(I), V(I), W(I), G(I) в цикле по I,
11 - ввод величины Н,
14 - вычисление величин F1, X1, Y1 в цикле по I,
I7 - вычисление величин J1(1), J1(2), J1(3) в цикле по I,
19 - вычисление величин E1, E2, E3,
21 - ввод величины WS.