Моделирование ударно-волновых процессов в алюминии с использованием малопараметрического уравнения состояния при нешаровом тензоре деформации

Размещено на http://www.allbest.ru/

Моделирование ударно-волновых процессов в алюминии с использованием малопараметрического уравнения состояния при нешаровом тензоре деформации

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

В работе [1] на основе термодинамического подхода построено малопараметрическое уравнение состояния, т.е. связь внутренней энергии, давления, температуры и удельного объема, в котором используются только параметры материалов, приведенные в таблицах физико-механических свойств, например, [2].

В качестве основы этого уравнения состояния принят принцип Ми-Грюнайзена представления термодинамического потенциала в виде суммы холодной и тепловых составляющих:

(1.1)

где E_x - холодная (упругая) составляющая энергии вдоль нулевой изотермы, E_a - вклад в энергию от теплового движения атомной решетки, E_e - вклад в энергию от теплового возбуждения электронов. Если ограничиться интервалом сжатий не больше 2, то для описания теплового движения атомов достаточно использовать приближение Дебая, а для вклада теплового движения электронов достаточно использовать зависимость, предложенную еще в [3]. Выражение для свободной энергии при этом имеет вид:

,(1.2)

и как следствие из (1.2) получаются следующие несложные выражения для внутренней энергии и давления:

,(1.3)

,(1.4)

где c_V0 - теплоемкость материала, c_e0 - электронная теплоемкость, g - коэффициент Грюнайзена, g_e=2/3 - аналог коэффициента Грюнайзена для электронов.

Для коэффициента Грюнайзена постулируется выражение, полученное в [4] независимо от свойств конкретного материала:

,(1.5)

=> ,(1.6)

где g_s--=--bK_sV₀/c_V, K_s - адиабатический модуль объемного сжатия, P_t,0 - тепловая часть давления при нормальных условиях, Q₀ - температура Дебая при нормальных условиях.

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

,(1.7)

где a_x=a(0)=1+2/(g_s - 2/3), параметр t задает определение средней частоты ансамбля осцилляторов в соответствии с теорией Ландау-Слэтера при t=0, Дуглейла-Макдональда при t=1 и теорией свободного объема при t=2.

В итоге выражения для холодной энергии и давления имеют вид:

,(1.8)

,(1.9)

где функции H₁(V), H₂(V) являются полиномами. Константы интегрирования C₁, C₂, C₃ определяются, если известно значение g(V₀), c_V0, c_e0, модуля объемного изотермического сжатия K_t при нормальных условиях и P_t,0.

Чтобы применить описанное уравнение состояния к какому-либо материалу необходимо и достаточно знание только шести констант V₀, b, K_t, c_V0, Q₀, c_e0, соответствующих значениям величин при нормальных условиях, которые могут быть найдены в справочниках физических величин.

В работах [5-7] развивается максвелловский подход к описанию процессов необратимого деформирования твердых тел. Для замыкания системы основных уравнений в данном случае необходимо уравнение состояния, включающее зависимость удельной внутренней энергии не только от первого, но и от второго инварианта тензора деформаций. Представляется целесообразным использовать идеи построения малопараметрических уравнений состояний и для максвелловских моделей. В данной работе принципы построения таких зависимостей применены для построения уравнения состояния для алюминия.

В одномерном случае система уравнений модели имеет вид:

(1.10)

где , u, E, S, T - плотность, скорость, удельная внутренняя энергия, энтропия и температура; h_i - логарифмы коэффициентов растяжения элемента среды вдоль главных осей; _i - главные напряжения; t, r - время и пространственная переменная; F = F(h_i, T) - уравнение свободной энергии при нешаровом тензоре деформаций; t--=--t (h_i, T) - зависимость времени релаксации касательных напряжений от параметров состояния среды; - показатель симметрии: =0 - плоская, =1 - цилиндрическая, =2 - сферическая. Уравнения записаны в системе координат, совпадающей с главными осями тензоров напряжения и деформаций (в данном случае они совпадают). В уравнение для свободной энергии вводится девиаторное слагаемое:

,(1.11)

,(1.12)

, , .

Здесь c_?-----поперечная скорость звука, x₀ - интерполяционная константа.

Зависимость для времени релаксации касательных напряжений имеет вид [8]:

,(1.13)

где - интенсивность касательных напряжений, e - величина пластической деформации, t₀ -параметр, связанный с частотой колебания атомов, N₀ - число дислокаций в недеформированном состоянии, M₀ - коэффициент размножения дислокаций, s₀ - характеристическое напряжение торможения, a - деформационное упрочнение.

Методика отыскания параметров (1.13) основана на решении задачи о деформирования тонкого стержня и сравнения расчетных и экспериментальных диаграмм деформирования. Более подробное описание можно найти в [8].

На рис. 1 сплошной линией показан расчет ударной адиабаты алюминия по уравнению состояния. Пунктирной линией показан расчет по уравнению состояния [9] и приведенные там же экспериментальные данные. На рис. 2 показан расчет температуры за фронтом ударной волны в сравнении с экспериментальными данными [10].

На рис. 3, 4 показано сравнение расчетных (линии) волновых профилей с экспериментальными [11], соответствующими удару алюминиевым ударником по алюминиевой мишени, набранной из нескольких пластин, между которыми устанавливались манганиновые датчики давления на расстоянии 4, 10, 15 мм от поверхности соударения. На рис. 3 ударник толщиной 5 мм со скоростью 595 м/с, на рис. 4 ударник толщиной 4 мм со скоростью 1505 м/с.

На рис. 5 сплошной линией показан расчет скорости свободной поверхности алюминиевого образца шириной 14.9 мм, по которому произвели удар алюминиевым ударником толщиной 2 мм со скоростью 450 м/с. Символами приводятся результаты соответствующего эксперимента [12].

На рис. 6 сплошными линиями показаны результаты серии расчетов скорости свободной поверхности в образцах из алюминия. Символами показаны результаты соответствующих экспериментов для алюминиевого сплава 2139-Т8 [13]. Во всех расчетах толщиной ударника 2.95 мм, толщиной мишени 5.94 мм.

На рис. 7, 8 линиями показаны результаты расчетов задачи о распространении упруго-пластической волны с выделенным упругим предвестником . В расчетах моделировались эксперименты [14], в которых производился удар алюминиевым ударником по алюминиевой мишени. Линии на рис. 7 соответствуют толщинам мишени 12.51 и 31.58 мм. На рис. 8 графики соответствуют толщине мишени: 6.13, 12.23 и 37.91 мм.

Приведенные результаты демонстрируют хорошее соответствие расчетных и экспериментальных данных.

Список литературы

поликристаллический сжатие деформирование уравнение

1. Краус Е.И. Малопараметрическое уравнение состояния твердого вещества при высоких плотностях энергии // Вестник НГУ. Серия: Физика.-2007.-Т.2, вып. 2.-С.65-73.

2. Физические величины: Справочник. // А.П. Бабичев, Н.А. Бабушкина, А.М. Братковский и др.; Под. ред. И.С. Григорьева, Е.З. Мейлихова.- М: Энергоатомиздат,1991.-1232 с.

3. Зельдович Я.Б., Райзер Ю.П. Физика ударных волн и высокотемпературных гидродинамических явлений. М., Физматгиз, 1963, 632 с.

4. Молодец А.М. Обобщенная функция Грюнайзена для конденсированных сред // ФГВ. - 1995. - т. 1. - №5, С. 132-133.

5. Мержиевский Л.А., Реснянский А.Д. Численное моделирование ударно-волновых процессов в металлах // ФГВ. - 1984. - т. 20. - № 5, с. 114-122.

6. Мержиевский Л.А. Моделирование динамического сжатия поликристаллического Al₂O₃ // ФГВ. - 1998. - т. 34. - № 6, с. 85-94.

7. Мержиевский Л.А., Воронин М. С. Моделирование ударно-волнового деформирования полиметилметакрилата // ФГВ. - 2012. - т. 48. - №2, с. 113-123.

8. Мержиевский Л.А., Палецкий А.В. Расчет диаграмм динамического деформирования металлов и сплавов // Физ. мезомех. - 2001. - Т. 4. - № 3. - С. 85-96.

9. http://www.ficp.ac.ru/rusbank/.

10. Альтшулер Л.В., Кормер С.Б., Бражник М.И., Владимиров Л.А., Сперанский М.П., Фунтиков А.И. // ЖЭТФ, - 38, вып. 4, - 1960.

11. Дремин А.Н., Канель Г.И., Черникова О.Б. Сопротивление пластической деформации алюминия АД-1 и дюралюминия Д-16 в условия ударного сжатия // ПМТФ, №4, 1981.

12. Kanel' G. I., Razorenov A. A., Utkin A. V., Baumung K. Experimental profiles of shock waves, Preprint of Scientific Association IVTAN of RAS (1996).

13. Casem D.T., Dandekar D.P. Shock and mechanical response of 2139-T8 aluminum // J. Appl. Phys. - 2012. - v. 111, p. 063508.

14. Johnson J.N., Barker L.M. Dislocation Dynamics and Steady Plastic Wave Profiles in 6061T6 Aluminum // J. Appl. Phys. 40, 4321 (1969).

Размещено на Allbest.ru