Оценка возможности транспортировки специальных изделий воздушным транспортом. Транспортировка системы “носитель-трал” военно-транспортной авиацией, страница 6

Первое из этих уравнений определяет форму, а второе характеризует распределение скоростей точек троса в начальный момент времени t=0.

Граничные условия:

U(0, t)= 0;

U(l, t)= 0 (8).

показывают, что перемещение точек закрепления троса равны нулю в любой момент времени.

2.5.1.  Решение волнового уравнения методом Фурье.

 - исходное уравнение собственных колебаний троса.

Рассмотрим свободные колебания.

Будем искать решение уравнения (9) в виде:

U(t, x)= T(X)*X(x) (10).

Следовательно, в уравнение (9):

T(t)*X(x)= a2*T(t)*X(x).

Дальше делим переменные:

Утверждение: поскольку левая часть может зависеть только от t, а правая от x, то ни левая, ни правая части не зависят ни от t, ни от x и являются постоянными.

Почему: пусть левая часть зависит от t. Зафиксируем правую часть (x). Тогда в правой части при фиксировании (x), если менять t, то правая часть (константа) равняется нескольким числам слева, а это не может быть. Аналогичные утверждения можно провести с x.

Из равенства (11):

где С= const.

Последнее соотношение показывает, что функции X(x) и T(t) удовлетворяют следующим дифференциальным уравнениям:

T-c*a2*T= 0;

X-c*X= 0.

Рассмотрим X-c*X= 0 (12).

Из граничных условий (8) следует:

X(0)*T(t)= 0,

X(l)*T(t)= 0.

Так как T(t)0 хотя бы в одной точке (нас интересует тривиальное решение), то

X(0)= 0, X(l)= 0 (13).

Найдем решение уравнения (13).

1.Пусть С= 0, следовательно X= 0,

μ2= 0 – характеристическое уравнение, μ1,2= 0.

X= C1*eμ1*x+C2*x*eμ2*x= C1+C2*x.

Определим С1 и С2:

X(0)=C1+C2*0= 0

X(l)= C1*C2*l

cледовательно,

C1=0, C2= 0, т. е. X(x)= 0.

2. Пусть c= λ2, следовательно, X- λ2*X= 0.

μ22= 0

μ1,2= λ

X(x)= C1*eλ*x+C2*e-λ*x

C1+C2= X(0)= 0

C1*eλ*e+C2*e-λ*e= 0

C1= C2= 0

X(x)= 0.

3. C=-λ2

X+ λ2*X= 0

μ2+ λ2= 0

μ1,2=i*λ.

Т. е. Характеристическое уравнение имеет комплексно сопряженные корни μ1 и μ2.

Указанным корням соответствуют действительные решения дифференциального уравнения eμ1  и eμ2.

Из этих решений линейно-независимым является только eμ1*x. В противном случае определитель Вронского равен нулю, т. е. получим нулевое решение.

X(x)= eμ1*x *cosλ*x+C2*sinλ*x

X(0)= C1*1+C2*0= 0

X(l)= C1*cosλ*l+C2*sinλ*l= 0

C1= 0

C2*sinλ*l= 0

Нас интересует случай, когда C2 0 (нетривиальное решение).

Получаем:

sinλ*l= 0, λ*l= π*k

λ= λ*k=  π*k/l ; k= 0,

Каждому  в λk соответствует:

Xk(x)= Ak*sinλkx= Ak*sin(k*π*x/l).

Число λk= k*π/l называются собственными числами, задачи (12), (13), а функции Xk(x)= Ak*sinλkx – собственными функциями.

Найдем T(t):

T(t)-c*a2*T= 0 при c= -λ2*k

Tk2k*a2*T= 0

μ2+ λ2k*a2= 0

μ1,2=i*λk*a

Tk(t)= Bk*cosλk*a*t+Dk*sinλk*a*t= Bk*cos(k*π*a*t/l)+Dk* sin(k*π*a*t/l)

Заметим, что k= 1, 2, 3. При k0 функции Xk(x) будут отличаться от соответственных функций для k>0 только значком, и можно включить в константу.

Т. е. каждому k соответствует функция:

Uk(t, x)= T(t)*X(x)

Uk(t, x)= [ak*cos(k*π*a*t/l)+bk* sin(k*π*a*t/l)]*sin(k*π*x/l) (14).

ak= Bk*Ak

bk= Dk*Ak

Функция (14) называется собственными функциями задачи о свободных колебаниях троса.

Определим  ak и bk.

Составим ряд:

U(t, x)=

 


Будем предполагать, что этот ряд сходится, его можно почленно дифференцировать.

Подберем числа ak и bk так, чтобы выполнялись начальные условия (7).

Заметим, что функция U(t, x), определяемая равенством (15) удовлетворяет граничным условиям (8).

Мы получили решение (14) и используем (8) при определении значений констант C1 и C2.

Кроме того, эта функция является решением дифференциального уравнения (9).

Из равенства (15) при t= 0 получаем:

U(t= 0)=

При t= 0 Ut= F(x)

Для того, чтобы выполнялись начальные условия (7) необходимо:

U(t= 0)=

U(t= 0)=

Для того, чтобы выполнялись (16) и (17) разложим f(x) и F(x) в ряд Фурье на промежутке [0, l] по sin.

Отступление: ряд Фурье по [-l, l].

Ряд Фурье для нечетных функций (по sin) имеет вид:

Равенство (16) будет разложением функции f(x) по sin, если:

Равенство (17) будет разложением функции F(x) по sin, если:

 

Следовательно, получаем:

Т. е. функция U(t, x), удовлетворяющая уравнению (9), начальным условиям (7) и условию (8), определяется следующим образом: