babkin_selivanov (550243), страница 46
Текст из файла (страница 46)
В еще более геометрически простом случаевзрыва сферического заряда, инициируемого в центре, движение обладает точечной симметрией, поэтому наиболее удобнопринять для описания движения сферическую систему координат (ж1 = г, х2 — 0, х3 = (р), которая обеспечивает зависимость параметров движения и состояния среды лишь от однойрадиальной координаты г и времени t (одномерная нестационарная задача с центральной симметрией).4.1.2. Выбор модели сплошной средыВыбор модели сплошной среды для участвующей в исследуемом процессе реальной деформируемой среды базируется294на анализе особенностей поведения этой среды в отношениисопротивления деформированию, на выделении основных факторов и игнорировании второстепенных.
Этап выбора модели заканчивается определением конкретного вида физическихсоотношений (3.7), ближе всего соответствующих особенностям физико-механического поведения реальной деформируемой среды.Например, при решении прикладной задачи прониканиятела вращения в воду с относительно небольшой начальнойскоростью взаимодействия Vo « 100 м/с в качестве модели реальной деформируемой среды (воды) вполне допустимо принять модель идеальной жидкости. Действительно, реальныежидкости обладают свойством сжимаемости и вязкости и вто же время не оказывают сопротивления непосредственно изменению формы своих частиц. При малых скоростях деформации, соответствующих малым скоростям взаимодействия,можно также пренебречь влиянием вязкости и вообще не учитывать появление касательных напряжений, используя дляописания физико-механического поведения физические соотношения (Jij — —р(Р) Т) gijy присущие модели идеальной среды.Следует отметить, что достаточно часто выбор моделисплошной среды применительно к процессам, происходящим вэкстремальных условиях (например, к взрывным и ударным),осуществляется итерационным путем, так как заранее трудно предсказать, какие именно физико-механические свойствареальных сред будут определяющими, а какими можно пренебречь.
В таких случаях последовательно используют все болеесложные модели, а критерием удовлетворительности выбораявляется соответствие получаемых расчетным путем результатов имеющимся экспериментальным данным.4.1.3. Составление системы исходных уравненийСистема исходных уравнений — это замкнутая системауравнений и соотношений, которая полностью описывает движение и состояние деформируемых сред с учетом их физикомеханических свойств.
В самом общем виде система исходныхуравнений имеет следующий вид:ю*295(4.1)(4-2)(4-3)(4-4)(4-5)(4-6)(4-7)Система исходных уравнений в обязательном порядкевключает основные общие для всех сплошных сред дифференциальные уравнения механики, выражающие фундаментальные законы сохранения массы (4.1), импульса (4.2), энергии(4.3), а также общие для всех сред кинематические соотношения (4.4) и (4.5) и геометрические соотношения (4.6). Индивидуальные особенности рассматриваемой деформируемойсреды в отношении оказания сопротивления деформированиюучитываются физическими соотношениями (4.7), обязательно включаемыми в систему исходных уравнений согласно выбранной модели сплошной среды.В зависимости от конкретного вида физических соотношений (4.7) и от характера процесса деформирования средыв систему исходных уравнений для обеспечения ее замкнутости могут быть включены дополнительные уравнения и соотношения.
Например, при отсутствии влияния температуры на физико-механическое поведение рассматриваемой среды физические соотношения имеют вид crty = сггу(£гу, eij) идля адиабатического процесса (V^1 = 0) система уравнений(4.1)—(4.7) является замкнутой и содержит 26 уравнений и296соотношений и такое же количество искомых характеристичесих функций (см. раздел 3.1). Напротив, в случаях зависимости компонент тензора напряжений от температуры илиже при учете теплообмена между частицами сплошной средыи необходимости определения температурного поля в систему исходных уравнений необходимо включать дополнительные соотношения, учитывающие закон теплопроводности Фурье (ql = qjglJ = -Xg^VjT, где А — коэффициент теплопроводности) и взаимосвязь между удельной внутренней энергиейи температурой (Е = Е(р, У)).В ряде случаев система исходных уравнений может бытьи более узкой, нежели представленная выше система(4.1) —(4.7).
Например, при постановке задачи механики идеальной жидкости, для которой компоненты тензора напряжений не зависят напрямую от компонент тензора деформаций(зависимость напряжений от деформаций имеет косвенный характер через плотность, взаимосвязанную с объемной деформацией (см. раздел 2.4.2)), не требуется включения в системуисходных уравнений кинематических соотношений (4.4) и геометрических соотношений (4.6). Однако в любом случае следует обеспечивать замкнутость системы исходных уравненийс равенством количества уравнений числу неизвестных характеристических функций, описывающих движение и состояниесплошной среды. Это является необходимым условием дляпоследующего нахождения единственного решения задачи.4.1.4.
Выбор основных неизвестныхи переход к системе разрешающих уравненийСистема разрешающих уравнений — это замкнутая система уравнений и соотношений, содержащая минимальноеколичество взаимно независимых искомых функций и получающаяся исключением остальных неизвестных функций изуравнений исходной системы.При рассмотрении системы исходных уравнений(4.1)—(4.7) нетрудно заметить, что она может быть преобразована с исключением некоторых неизвестных функций. На297пример, компонентытензора скоростей деформаций, участвующие в записи уравнения энергии (4.3) и физических соотношений (4.7), могут быть выражены через компоненты щ вектора скорости в соответствии с кинематическими соотношениями (4.5) и, таким образом, исключены из системы уравненийпри сохранении в качестве основных неизвестных компонентвектора скорости.
Аналогичным образом из уравнений движения (4.2) и энергии (4.3) могут быть исключены и компонентыajj тензора напряжений в соответствии с физическими соотношениями (4.7) и т.п. Такой процесс преобразования системы исходных уравнений сопровождается уменьшением количества искомых функций и уравнений, а по существу, соответствует получению из общих уравнений механики их частноговида и построению замкнутой системы уравнений, описывающей движение и состояние сплошной среды с определенными физико-механическими свойствами (идеальной или вязкойжидкости, упругой или упругопластической среды и т.д.).В математическом отношении система разрешающихуравнений представляет собой систему дифференциальныхуравнений в частных производных и конечных функциональных уравнений. Например, входящие в уравнения движения(4.2) абсолютные производные компонент тензора напряженийпо координатам представляются в виде суммы частных производных по координатам и дополнительных составляющих,зависящих от символов Кристоффеля выбранной системы координат и учитывающих искривленность ее координатных линий (см.
раздел 1.3.6). Входящие в эти же уравнения субстанциональные производные компонент вектора скорости повремени при использовании, например, подхода Эйлера представляются в виде суммы локальной производной (частнаяпроизводная по времени, определяемая в фиксированных точках пространства) и конвективной производной (зависящей отчастных щроизводных по координатам) и т.п. К числу конечных функциональных уравнений относятся в большинствеслучаев физические соотношения (4.7).Система разрешающих уравнений является окончательной, результирующей системой, и именно эта система требует298своего решения путем нахождения всех входящих в нее характеристических функций.
С определением зависимостей физических величин, входящих в систему разрешающих уравнений, от координат и времени получает свое логическое завершение феноменологический подход к изучению механическогодвижения деформируемых сред. В большинстве представляющих практический интерес случаев системы разрешающихуравнений, описывающие движение газов, жидкостей и твердых тел, достаточно сложны, а их решение может быть найдено при использовании методов вычислительной математики —численных методов механики сплошных сред.4.1.5. Начальные и граничные условияНеотъемлемым и важнейшим элементом постановки любой задачи механики сплошных сред является формулировканачальных и граничных условий.
Их значение определяетсятем, что та или иная система разрешающих уравнений описывает целый класс движений соответствующей деформируемой среды, и лишь задание отвечающих исследуемому процессу начальных и граничных условий позволяет выделить изэтого класса представляющий интерес частный случай, соответствующий решаемой практической задаче.Начальные условия — это условия, которыми задаютсязначения искомых характеристических функций в момент начала рассмотрения исследуемого процесса. Количество задаваемых начальных условий определяется количеством основных неизвестных функций, входящих в систему разрешающихуравнений, а также порядком входящей в эту систему высшейпроизводной по времени. Например, адиабатическое движениеидеальной жидкости или идеального газа описывается системой шести уравнений с шестью основными неизвестными —тремя компонентами вектора скорости = ^(х1, /), давлением р —/), плотностью р = р(х') t) и удельной внутреннейэнергией Е = Е(х', £), при этом порядок производных этих физических величин по времени не превышает первый порядок299(см.
раздел 4.2). Соответственно этому в качестве начальных условий должны быть заданы начальные поля этих шести физических величин: при t = 0 v, = vfa', 0), р = р(х\ 0),р = р(ж*, 0), Е = Е(хг, 0). В некоторых случаях (например,в динамической теории упругости; см. раздел 4.4) в качестве основных неизвестных в системе разрешающих уравнений используются не компонентывектора скорости, а компоненты щ вектора перемещения, а уравнение движения содержит производные второго порядка компонент перемещениядРщ/dt2, что требует задания двух начальных условий дляискомой функции щ — щ(х\ t): при t = 0 щ = щ(х', 0) иdui/dt — ^(ж*, 0).Более сложным и разнообразным образом при постановкезадач механики сплошных сред задаются граничные условия.Граничные условия— это условия, которыми задаются значения искомых функций (или их производных по координатам ивремени) на поверхности S области, занимаемой деформируемой средой.