Термодинамические основы термоупругости. I. Граничные условия первого рода Граничные и начальные условия дифференциального уравнения

), задающее его поведение в начальный момент времени или на границе рассматриваемой области соответственно.

Обычно дифференциальное уравнение имеет не одно решение, а целое их семейство. Начальные и граничные условия позволяют выбрать из него одно, соответствующее реальному физическому процессу или явлению. В теории обыкновенных дифференциальных уравнений доказана теорема существования и единственности решения задачи с начальным условием (т. н. задачи Коши). Для уравнений в частных производных получены некоторые теоремы существования и единственности решений для определённых классов начальных и краевых задач.

Терминология

Иногда к граничным относят и начальные условия в нестационарных задачах, таких как решение гиперболических или параболических уравнений .

Для стационарных задач существует разделение граничных условий на главные и естественные .

Главные условия обычно имеют вид u (∂ Ω) = g {\displaystyle u(\partial \Omega)=g} , где ∂ Ω {\displaystyle \partial \Omega } - граница области Ω {\displaystyle \Omega } .

Естественные условия содержат также и производную решения по нормали к границе.

Пример

Уравнение d 2 y d t 2 = − g {\displaystyle {\frac {d^{2}y}{dt^{2}}}=-g} описывает движение тела в поле земного тяготения . Ему удовлетворяет любая квадратичная функция вида y (t) = − g t 2 / 2 + a t + b , {\displaystyle y(t)=-gt^{2}/2+at+b,} где a , b {\displaystyle a,b} - произвольные числа. Для выделения конкретного закона движения необходимо указать начальную координату тела и его скорость, то есть начальные условия .

Корректность постановки граничных условий

Задачи математической физики описывают реальные физические процессы, а потому их постановка должна удовлетворять следующим естественным требованиям:

  1. Решение должно существовать в каком-либо классе функций;
  2. Решение должно быть единственным в каком-либо классе функций;
  3. Решение должно непрерывно зависеть от данных (начальных и граничных условий, свободного члена, коэффициентов и т. д.).

Требование непрерывной зависимости решения обусловливается тем обстоятельством, что физические данные, как правило, определяются из эксперимента приближённо, и поэтому нужно быть уверенным в том, что решение задачи в рамках выбранной математической модели не будет существенно зависеть от погрешности измерений. Математически это требование можно записать, например, так (для независимости от свободного члена):

Пусть задано два дифференциальных уравнения: L u = F 1 , L u = F 2 {\displaystyle Lu=F_{1},~Lu=F_{2}} с одинаковыми дифференциальными операторами и одинаковыми граничными условиями, тогда их решения будут непрерывно зависеть от свободного члена, если:

∀ ε > 0 ∃ δ > 0: (‖ F 1 − F 2 ‖ < δ) ⇒ (‖ u 1 − u 2 ‖ < ε) {\displaystyle \forall \varepsilon >0~\exists \delta >0:~\left(\|F_{1}-F_{2}\|<\delta \right)\Rightarrow \left(\|u_{1}-u_{2}\|<\varepsilon \right)} , где u 1 {\displaystyle u_{1}} , u 2 {\displaystyle u_{2}} - решения соответствующих уравнений.

Множество функций, для которых выполняются перечисленные требования, называется классом корректности . Некорректную постановку граничных условий хорошо иллюстрирует

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

Начальные условия -- это условия, которыми задаются значения искомых характеристических функций в момент начала рассмотрения исследуемого процесса. Количество задаваемых начальных условий определяется количеством основных неизвестных функций, входящих в систему разрешающих уравнений, а также порядком входящей в эту систему высшей производной по времени. Например, адиабатическое движение идеальной жидкости или идеального газа описывается системой шести уравнений с шестью основными неизвестными -- тремя компонентами вектора скорости,давлением,плотностью и удельной внутренней энергией, при этом порядок производных этих физических величин по времени не превышает первый порядок. Соответственно этому в качестве начальных условий должны быть заданы начальные поля этих шести физических величин: при t =0 ,. В некоторых случаях (например, в динамической теории упругости) в качестве основных неизвестных в системе разрешающих уравнений используются не компоненты вектора скорости, а компоненты вектора перемещения, а уравнение движения содержит производные второго порядка компонент перемещения, что требует задания двух начальных условий для искомой функции: при t = 0

Более сложным и разнообразным образом при постановке задач механики сплошных сред задаются граничные условия. Граничные условия -- это условия, которыми задаются значения искомых функций (или их производных по координатам и времени) на поверхности S области, занимаемой деформируемой средой. Различают граничные условия нескольких типов: кинематические, динамические, смешанные и температурные.

Кинематические граничные условия соответствуют случаю, когда на поверхности S тела (или ее части) задаются перемещения или скорости где -- координаты точек поверхности S, изменяющиеся в общем случае в зависимости от времени.

Динамические граничные условия (или граничные условия в напряжениях) задаются, когда на поверхности S действуют поверхностные силы р. Как следует из теории напряжений, в этом случае на любой элементарной площадке поверхности с единичным вектором нормали п вектор удельных поверхностных сил рп принудительно задает вектор полного напряжения?п = рn, действующий в сплошной среде в точке на данном участке поверхности, что приводит к взаимосвязи тензора напряжений (?) в этой точке с поверхностной силой и ориентацией вектора п соответствующего участка поверхности: (?) · п = рп или.

Смешанные граничные условия соответствуют случаю, когда на поверхности S задаются значения и кинематических, динамических величин или устанавливаются взаимосвязи между ними.

Температурные граничные условия подразделяются на несколько групп (родов). Граничные условия первого рода задают на поверхности S деформируемой среды определенные значения температуры Т. Граничные условия второго рода задают на границе вектор теплового потока q, что с учетом закона теплопроводности Фурье q = -- ? grad T, по существу, накладывает ограничения на характер температурного распределения в окрестности граничной точки. Граничные условия третьего рода устанавливают зависимость между вектором теплового потока q, направленным к данной среде со стороны окружающей среды, и температурным перепадом между этими средами и т.д.

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

На рис. 3 схематично представлен процесс взаимодействия при проникании деформируемого тела I в деформируемую преграду II. Тело I ограничено поверхностями S1 и S5, а тело II -- поверхностями S2, S3, S4, S5. По -верхность S5 является границей раздела взаимодействующих деформируемых тел. Будем полагать, что движение тела I до начала взаимодействия, а также в его процессе происходит в жидкости, создающей определенное гидростатическое давление

Рисунок 3

и задающей внешние по отношению к обоим телам поверхностные силы рп = -- рп= -- рni ri, действующие на любой из элементарных площадок поверхностей S1 тела I и S2 преграды II, граничащих с жидкостью. Будем также считать, что поверхность Sз преграды жестко закреплена, а поверхность S4 свободна от действия поверхностных сил (рп = 0).

Для приведенного примера на различных поверхностях, ограничивающих деформируемые среды I и II, должны задаваться граничные условия всех трех основных типов. Очевидно, что на жестко закрепленной поверхности Sз следует задать кинематические граничные условия?(S3) = ?(, t) = 0. Граничные условия на поверхностях S1 и S2 однотипны и относятся к динамическим условиям, накладывающим ограничения на компоненты тензора напряжений в граничных точках соответствующих тел: или Компоненты тензора напряжений на поверхности S4 преграды также не могут быть произвольными, а взаимосвязаны с ориентацией ее элементарных площадок как.

Граничные условия на границе раздела (поверхность S5) взаимодействующих деформируемых сред являются наиболее сложными и относятся к условиям смешанного типа, включающим, в свою очередь, кинематическую и динамическую части (см. рис. 3). Кинематическая часть смешанных граничных условий накладывает ограничения на скорости движения индивидуальных точек обеих сред, находящихся в контакте в каждой пространственной точке поверхности S5. Возможны два варианта задания этих ограничений, проиллюстрированные на рис. 4, а и б. По наиболее простому первому варианту предполагается, что скорости движения любых двух находящихся в контакте индивидуальных точек одинаковы (? = ?) -- это так называемое условие "прилипания", или условие "сварки" (см. рис. 4, а). Более сложным и в то же время более адекватным для рассматриваемого процесса является задание условия "непроницаемости", или условия "непротекания" (? · n= ? · n; см. рис. 4, б), которое соответствует экспериментально подтверждающемуся факту: взаимодействующие деформируемые среды не могут проникать


Рисунок 4

друг в друга или отставать друг от друга, а могут проскальзывать одна относительно другой со скоростью? - ?, направленной по касательной к границе раздела ((?I - ?II) · n = 0). Динамическая часть смешанных граничных условий на границе раздела двух сред формулируется на основе третьего закона Ньютона с использованием соотношений теории напряжений (рис. 4, в). Так, в каждой из двух находящихся в контакте индивидуальных частиц деформируемых сред I и II реализуется свое напряженное состояние, характеризуемое тензорами напряжений (?)I и (?) II.При этом в среде I на каждой элементарной площадке границы раздела с единичным вектором нормали nII, внешней по отношению к данной среде, действует вектор полного напряжения?nI = (?)·nI. В среде II на той же площадке, но с единичным вектором нормали nII , внешней по отношению к этой среде, действует вектор полного напряжения?nII =(?)II · пII. С учетом взаимности действия и противодействия?nI = - ? n II , а также очевидного условия nI = --nII = n устанавливается взаимосвязь между тензорами напряжений в обеих взаимодействующих средах на границе их раздела: (?)I · п = (?) II ·п или же (?ijI - ?ijII) nj = 0.Возможные варианты задания граничных условий не исчерпываются рассмотренным частным примером. Вариантов задания начальных и граничных условий столь же много, сколь много существует в природе и технике процессов взаимодействия деформируемых тел или сред. Они определяются особенностями решаемой практической задачи и задаются в соответствии с приведенными выше общими принципами.

Одного уравнения движения (1.116) при математическом описании физического процесса недостаточно. Надо сформулировать условия, достаточные для однозначного определения процесса. При рассмотрении задачи о колебании струны дополнительные условия могут быть двух видов: начальные и граничные (краевые).

Сформулируем дополнительные условия для струны с закрепленными концами. Так как концы струны длины закреплены, то их отклонения в точках и должны быть равны нулю при любых :

, . (1.119)

Условия (1.119) называются граничными условиями; они показывают, что происходит на концах струны на протяжении процесса колебания.

Очевидно, процесс колебаний будет зависеть от того, каким способом струна выводится из состояния равновесия. Удобнее считать, что струна начала колебаться в момент времени . В начальный момент времени всем точкам струны сообщаются некоторые смещения и скорости:

,

, , (1.120)

где и - заданные функции.

Условия (1.120) называются начальными условиями.

Итак, физическая задача о колебаниях струны свелась к следующей математической задаче: найти такое решение уравнения (1.116) (или (1.117) или (1.118)), которое удовлетворяло бы граничным условиям (1.119) и начальным условиям (1.120). Эта задача называется смешанной краевой задачей, так как включает в себя и граничные и начальные условия. Доказано, что при некоторых ограничениях, наложенных на функции и , смешанная задача имеет единственное решение.

Оказывается, что к задаче (1.116), (1.119), (1.120), помимо задачи о колебаниях струны, сводятся многие другие физические задачи: продольные колебания упругого стержня, крутильные колебания вала, колебания жидкостей и газа в трубе и др.

Помимо граничных условий (1.119) возможны граничные условия других типов. Наиболее распространенными являются следующие:

I. , ;

II. , ;

III. , ,

где , - известные функции, а , - известные постоянные.

Приведенные граничные условия называют соответственно граничными условиями первого, второго, третьего рода. Условия I имеют место в том случае, если концы объекта (струна, стержень и т.д.) перемещаются по заданному закону; условия II – в случае, если к концам приложены заданные силы; условия III – в случае упругого закрепления концов.

Если функции, заданные в правой части равенств, равны нулю, то граничные условия называются однородными. Так, граничные условия (1.119) – однородные.

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

Для уравнения (1.116) может быть поставлена и другая задача. Пусть струна достаточно длинная и нас интересует колебание ее точек, достаточно удаленных от концов, причем в течение малого промежутка времени. В этом случае режим на концах не будет оказывать существенного влияния и поэтому его не учитывают; струну же при этом считают бесконечной. Вместо полной задачи ставят предельную задачу с начальными условиями для неограниченной области: найти решение уравнения (1.116) для при , удовлетворяющее начальным условиям:

, .

U| x=0 = g 1 (t), U| x=l = g 2 (t)

Эти условия физически означают, что на концах заданы режимы колебаний.

II. Граничные условия второго рода

U x | x=0 = g 1 (t), U x | x=l = g 2 (t)

Такие условия соответствуют тому, что на концах заданы силы.

III. Граничные условия третьего рода

(U x 1 U)| x=0 = g 1 (t) , (U x –σ 2 U)| x=l = g 2 (t)

Эти условия соответствуют упругому закреплению концов.

Граничные условия (5), (6) и (7) называются однородными, если правые части g 1 (t) и g 2 (t) тождественно равны нулю при всех значениях t. Если хотя бы одна из функций в правых частях не равна нулю, то граничные условия называются неоднородными.

Аналогично формулируются граничные условия и в случае трех или четырех переменных при условии, что одна из этих переменных - время. Границей в этих случаях будет или замкнутая кривая Г, ограничивающая некоторую плоскую область, или замкнутая поверхность Ω, ограничивающая область в пространстве. Соответственно изменится и производная от функции, фигурирующая в граничных условиях второго и третьего рода. Это будет производная по нормали n к кривой Г на плоскости или к поверхности Ω в пространстве, причем, как правило, рассматривают нормаль, внешнюю по отношению к области(см.рис.5).

К примеру, граничное условие (однородное) первого рода на плоскости записывается в виде U| Γ =О, в пространствеU| Ω =0. Граничное условие второго рода на плоскости имеет вид ,а в пространстве . Конечно, физический смысл этих условий разный для различных задач.

При постановке начальных и граничных условий возникает задача об отыскании решения дифференциального уравнения, удолетворяющего заданным начальным и граничным (краевым) условиям. Для волнового уравнения (3) или (4), начальных условий U(x,0)=φ(x), U t (x,0)=ψ(x) и в случае граничных условий первого рода (5), задача называетсяпервой начально-краевой задачей для волнового уравнения . Если вместо граничных условий первого рода задавать условия второго рода (6) или третьего рода (7), то задача будет называться, соответственно, второй и третьей начально-краевой задачей . Если граничные условия на разных участках границы имеют различные типы, то такие начально-краевые задачи называют смешанными .

Рассмотрим две типичных электростатических задачи :

1) Найти потенциал электрического поля при неизвестном местоположении исходных зарядов, но заданном электрическом потенциале на границах области. (Например, задача о распределении потенциала электрического поля, создаваемого системой неподвижных проводников, помещенных в вакуум и подключенных к батареям. Здесь можно измерить потенциал каждого проводника, но задать распределение электрических зарядов на проводниках, зависящее от их формы, весьма сложно.)

2) Найти потенциал электрического поля, создаваемого заданным распределением в пространстве электрических зарядов .

Хорошо известно, что прямой метод вычисления потенциала электрического поля в этих задачах состоит в решении уравнения Лапласа (задача 1)

(1)

и уравнения Пуассона (задача 2)

. (2)

Уравнения (1), (2) относится к классу дифференциальных уравнений в частных производных эллиптического типа .

Далее мы будем рассматривать только частный случай эллиптических уравнений для поля  , зависящего от двух пространственных переменных. Совершенно очевидно, что для полного решения задачи уравнения (1), (2) необходимо дополнить граничными условиями. Различают три типа граничных условий:

1) граничные условия Дирихле (значения  задаются на некоторой замкнутой кривой в плоскости (х,у) и, возможно, на некоторых дополнительных кривых, расположенных внутри области (рис. 1));

2) граничные условия Неймана (на границе задается нормальная производная потенциала );

3) смешанная краевая задача (на границе задается линейная комбинация потенциала  и его нормальной производной).

Начальные условия отвечают на вопрос о том, каково было температурное поле в момент вре­мени, принятый за начало отсчета. Они описываются выраже­нием . Очень часто температура компонентов технологических подсистем в начальный момент времени может быть принята равной температуре окружающей среды, т. е. . В этом случае удобно, как отмечалось выше, вести расчет в так называемых избыточных температурах, условно считая, что , а затем по окончании расчета к ре­зультату прибавляя . Граничными называются условия взаимодействия поверх­ностей тел с окружающей средой или другими телами. Различают несколько разновидностей граничных условий. При граничных условиях первого рода (ГУ1) предполагают, что известен закон распределения температур на граничных поверхностях тела . Пусть, например, требуется определить темпе­ратурное поле внутри какой-либо детали или инструмента. Сде­лать это экспериментальным путем, не разрушая объект измерения, довольно трудно, измерить же температуру на поверхности де­тали, инструмента или другого твердого тела экспериментальным путем значительно проще, это может быть выполнено без повреж­дения объекта. Если мы знаем ГУ1 в виде закона распределения температур на поверхностях тела, то, решая дифференциальное уравнение теплопроводности, можем рассчитать поле температур внутри детали, инструмента и т. д. Частным случаем ГУ1 яв­ляется условие изотермичности поверхностей тела, т. е. .

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

Граничные условия третьего рода (ГУЗ) используют в том случае, когда теплообменом поверхности с окружающей средой пренебречь нельзя. В этом случае должны быть заданы температура среды, с которой соприкасается данное тело, и так называемый коэффициент теплоотдачи , Вт/(м 2 × °С), характеризующий теплообмен между средой и поверхностью.

Согласно закону Ньютона-Рихмана плотность теплового потока пропорциональна разности температур поверхности и окружающей ее среды , т.


е.

Формула (2.1) дает возможность определить количество теп­лоты , Вт/м 2 , которое в единицу времени с единицы поверхности отводится в окружающую среду. Как следует из закона Фурье, к поверхности тела подводится поток

.

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

или . (2.2)

Выражение (2.2) представляет собой математическое описание граничных условий третьего рода.

Граничные условия четвертого рода (ГУ4) возникают тогда, когда рассматриваемое твердое тело находится в беззазорном контакте с другим твердым телом и между ними происходит теплообмен. Этот вариант граничных условий весьма часто встре­чается в теплофизике технологических процессов. Например, при обработке давлением детали штампа практически беззазорно соприкасаются с обрабатываемой заготовкой; при резании ме­талла поверхности инструмента на определенных участках сопри­касаются со стружкой и заготовкой. При граничных условиях четвертого рода, когда контакт между телами идеален, темпе­ратура в любой точке поверхности соприкосновения как со сто­роны одного, так и со стороны другого тела одна и та же, т. е.

С целью упрощения расчетов часто вместо равенства темпе­ратур в каждой точке контакта в качестве ГУ4 принимают равенство средних температур на поверхности контакта, т. е. вместо формулы (2.3) полагают

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

Заканчивая рассмотрение вопроса о граничных условиях, отметим, что на разных участках реальных тел могут иметь место различные граничные условия. Рассмотрим, например, процесс плоского шлифования заготовки торцом чашечного круга (см. рис. 2.5). Если решена задача о распределении теплоты шлифо­вания между кругом и заготовкой, то по отношению к заготовке имеем следующие граничные условия: ГУ3 - на поверхности соприкосновения с жидкостью; ГУ2 - на контактной поверхности с кругом, где известна плотность теплового потока, и на торце заготовки, который можно считать адиабатическим, если пре­небречь его теплоотдачей в воздух; ГУ4 - на поверхности, где заготовка соприкасается с магнитным столом станка.