При поиске функции зеленого, как определить отклонение u(x)u(x)u(x)

В PDE ниже, ты описывает прогиб жесткой балки с плотностью длины р л , модуль упругости Е , момент инерции я и плотность внешней силы ф ( Икс , т ) .

Вопрос состоит в том, чтобы определить функцию Грина стационарного прогиба балки, когда ф ( Икс , т ) "=" ф ( Икс ) , то есть когда плотность внешней силы не зависит от времени.

ПДЭ:

р л ты т т Е я ты Икс Икс Икс Икс "=" ф ( Икс , т )
Граничные условия:
ты ( 0 , т ) "=" ты Икс ( 0 , т ) "=" ты ( л , т ) "=" ты Икс ( л , т ) "=" 0

К моей проблеме: мне помог учитель, который сказал мне, что я должен начать с определения

ты ( Икс ) "=" 0 л г ( Икс , с ) ф ( с ) Е я г с

С этого момента у меня есть идея, как действовать дальше. Но я понятия не имею, как он знал, что поставить как ты ( Икс ) , и было бы здорово, если бы я мог получить объяснение этого. Например, как он узнал, что ф ( с ) следует разделить на Е я ? Спасибо!

Ответы (1)

Когда луч неподвижен, ваш PDE становится

Е я ты Икс Икс Икс Икс "=" ф ( Икс ) ,
или
ты Икс Икс Икс Икс "=" ф ( Икс ) Е я .
Так вот откуда это деление.

Теперь, глядя на значение функции Грина , давайте предположим , что она удовлетворяет вашему дифференциальному уравнению, с дельта ( Икс с ) вместо исходного термина. В этом случае исходный термин ф ( Икс ) Е я , поэтому мы предполагаем, что

г Икс Икс Икс Икс "=" дельта ( Икс с ) .
(Чтобы получить решение какой-либо конкретной задачи, вам обычно нужно найти конкретную функцию Грина, удовлетворяющую вашему уравнению и граничным условиям. Но сейчас все, что нам нужно, это просто предположить, что она у вас уже есть.)

Тогда причина, по которой функции Грина полезны, заключается в том, что мы можем использовать решение для г найти ты путем выражения исходного члена в виде интеграла, включающего это дельта -функция, а затем заменить дельта с этой величиной функции Грина:

Икс 4 ты "=" ф ( Икс ) Е я "=" дельта ( Икс с ) ф ( с ) Е я г с "=" Икс 4 г ( Икс , с ) ф ( с ) Е я г с "=" Икс 4 г ( Икс , с ) ф ( с ) Е я г с .
Теперь, сравнивая первую строку с последней, мы видим, что одно решение для ты просто
ты ( Икс ) "=" г ( Икс , с ) ф ( с ) Е я г с .
Предполагая г также удовлетворяет вашим граничным условиям, это решение, которое вы ищете.

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


Я отмечу, что на самом деле вам не нужно делить на Е я . Вы могли относиться к своему источнику как к простому ф ( Икс ) , так что УЧП было просто

Е я ты Икс Икс Икс Икс "=" ф ( Икс ) .
Помните, что вы заменяете источник [вот это ф ( Икс ) ] к дельта , а затем предположим, что г удовлетворяет этому уравнению, поэтому вам нужно будет иметь
Е я г Икс Икс Икс Икс "=" дельта ( Икс с ) .
Результат был бы таким же, за исключением того, что то, что вы называете «функцией Грина», было бы изменено на Е я .

Вы также могли бы разделить на Е я , в таком случае

ты Икс Икс Икс Икс "=" ф ( Икс ) Е я
Говорит, что ф ( Икс ) Е я является вашим источником, поэтому ваша функция Грина должна удовлетворять
г Икс Икс Икс Икс "=" дельта ( Икс с ) .
Все дело в том, когда вы делите на любую заданную константу.

Спасибо за развернутый ответ! Это очень помогло. Сейчас я немного сбит с толку г Икс Икс Икс Икс "=" дельта ( с Икс ) . Имеет ли это значение в дельта -функция, если написать дельта ( с Икс ) или дельта ( Икс с ) ? Кроме того, как вы понимаете, что вы должны поставить отрицательный знак перед г Икс Икс Икс Икс , а не в интеграле перед г ( Икс , с ) ?
Знак внутри дельта -function на самом деле не имеет значения , потому что она "делает что-нибудь" только тогда, когда Икс с "=" 0 "=" с Икс . Теперь, что касается знака перед г Икс Икс Икс Икс , идея в том, что г должно удовлетворять точно тому же уравнению, что и ты – за исключением случаев, когда источник заменен на дельта -функция. В этом случае исходное УЧП имело отрицательный знак, поэтому оно стоит перед г Икс Икс Икс Икс , но не внутри интеграла. Я отредактировал ответ, чтобы добавить немного ясности; обратите внимание на самую последнюю часть для этого конкретного вопроса.
У меня все еще есть несколько проблем с пониманием функции Грина на абстрактном уровне, но это мне очень помогло. Я решил свою проблему, и я думаю, что получу лучшее понимание, работая с новыми проблемами еще несколько дней. Еще раз, спасибо! Я бы дал вам +1, но у меня недостаточно репутации, лол.