Точное решение. Уравнение (5) представляет собой дифференциальное уравнение второго порядка, решаемое в квадратурах

Лекция 2

Уравнение (5) представляет собой дифференциальное уравнение второго порядка, решаемое в квадратурах, т. е. оно решается последовательным интегрированием. Проинтегрировав его дважды

(1.1)

мы получим решение с точностью до констант интегрирования С₁ и С₂.

Эти константы определяются из граничных условий (из условий закрепления торцов стержня). Граничные условия могут быть геометрическими (ГГУ), когда на торцах заданы перемещения, статическими (СГУ), когда на торцах заданы силы, или смешанными, когда на одном торце заданы ГГУ, а на другом – СГУ.

Рассмотрим примеры решения при различных граничных условиях.

Пусть, например, нагрузка изменяется по линейному закону

. (1.2)

Дифференциальное уравнение (5) примет вид (ЕА=const):

. (1.3)

Дважды интегрируем уравнение (7):

(1.4)

Рассмотрим геометрические ГУ (рис. 5 для варианта А).

Рисунок 5. Геометрические граничные условия (вариант А)

Подчиним решение (9) геометрическим граничным условиям

(1.5)

откуда

(1.6)

Внося полученные значения констант во второе соотношение (1.4), получаем выражения для перемещения

. (1.7)

Здесь введена безразмерная координата

. (1.8)

На основании соотношения (4) получаем выражение для продольного усилия

. (1.9)

Произведем ручной счет по формулам (1.7) и (1.9) в пяти равноотстоящих точках. Результаты представим в виде таблицы 1 и графиков (рис. 6 и 7).

Таблица 1

		0.25	0.5	0.75	1.0
		51/64 (0.797)	9/8 (1.125)	57/64 (0.891)
	4.0	37/64 (2.313)	1/4 (0.25)	-35/16 (-2.188)	-5

Рис 6. Изменение продольного перемещения по длине стержня (точное решение - вариант А)

Рис. 7. Изменение продольного усилия по длине стержня (точное решение - вариант А)

Запишем программу, численно реализующую выражения (1.7), (1.9), на алгоритмическом языке «Паскаль».

PROGRAM STTR;

uses crt;

const

qo=1.0; l=1.0; EA=1.0; m=4;

var i: integer;

qoln,qolu,dx,x,x3: real;

u,N: array[1..m+1] of real;

BEGIN

clrscr;

qoln:=qo*l/6;

qolu:=qoln*l/EA;

dx:=l/m;

writeln;

writeln(‘ Результат решения’);

writeln;

writeln(‘Координата Перемещение Усилие’);

for i:=1 to m+1 do begin

x:=dx*(i-1)/l;

x2:=sqr(x);

u[i]:={qolu*}x*(-x2-3*x+4);

N[i]:={qoln*}(-3*x2);

writeln;

writeln(‘ x=’,x:5:3,

‘ u=’,u[i]*6:7:4,’ N=’,N[i]*6:7:4);

end;

readln;

END.

Дата добавления: 2014-01-20