Программа расчета геометрических характеристик поперечных сечений с учетом движения фронта повреждения

Предполагается, что при числе точек излома контура попереч­ного сечения, равном четырем /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,...,М) - скорости распространения поврежде­ний /СМ/год/,

F₁ - величина площади-сечения /СM² /,

Х₁ - координата центра тяжести сечения /СМ /,

Y₁- координата центра тяжести сечения /СМ /,

J₁(1)=J_X1- величина осевого момента инерции /СМ⁴/,

J₁(2)=J_Y1- величина осевого момента инерции /СМ⁴/,

J₁(3)=J_X1Y1- величина центробежного момента инерции /СМ⁴/,

E₁=J_U - величина первого из экстремальных моментов инер­ции/СМ⁴/,

E₂=J_V - величина второго из экстремальных моментов инер­ции /СМ⁴/,

E₃=αº - угол наклона главных центральных осей к осям 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 - вычисление величин F₁, X₁, Y₁ в цикле по I,

I7 - вычисление величин J₁(1), J₁(2), J₁(3) в цикле по I,

19 - вычисление величин E₁, E₂, E₃,

21 - ввод величины WS.