Преобразование скорости удара в величину давления

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

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

Поскольку для моделирования этого я использую код Эйлера, мне нужно знать давление результирующей ударной волны в материале. Есть ли способ определить это давление?

Предположим, что я знаю все свойства материала (плотность, скорость звука, модуль Юнга и т. д.), что ударник абсолютно жесткий и что удар находится в пределах упругости образца материала.

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

Соответствующий URL: lacan.upc.edu/sites/default/files/ppt_cdc_1.pdf

Ответы (1)

Предыстория
Существует хороший справочник 1 по физике звуковых/ударных волн в твердых телах (см. главу XI). Я нашел следующее (на странице 688) очень интересным и относящимся к вашему вопросу:

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

Давайте определим п как давление и ε как внутренняя энергия твердого тела. Их можно разделить на две части: эластичная (индекс с ) и термическая часть. п с а также ε с зависит только от плотности материала, р , или конкретный объем, В знак равно 1 / р . Они равны полному давлению и удельной внутренней энергии при абсолютном нуле или Т знак равно 0   К . Предположим, что удельный объем при Т знак равно 0 а также п знак равно 0 дан кем-то В о с , что всего на ~1-2% меньше, чем удельный объем на СТП , В о , для большинства металлов.

Кривая потенциальной энергии, или кривая, определяющая ε с , качественно похожа на кривую потенциальной энергии, описывающую взаимодействие между двумя атомами в зависимости от внутриядерного расстояния, Δ Икс н . Когда В > В о с , силы притяжения преобладают, но быстро падают по мере увеличения внутриядерных расстояний (например, по мере Т увеличивается). Другими словами, когда атомы расходятся дальше друг от друга ε с будет асимптотически возрастать до некоторого значения U , что примерно соответствует энергии связи атомов в теле. Таким образом, U представляет собой энергию, необходимую для удаления всех атомов из объекта в бесконечность, которая примерно равна теплоте парообразования для материала (более подробно я написал о теплоте парообразования и предоставил несколько полезных ссылок в этом ответе ). Например, теплота атомизации (аналогично испарению) железа примерно равна 415 к Дж м о л 1 или ~4,3 эВ/атом. Таким образом, ε с ( В ) U в качестве Δ Икс н →∼ 2 .

И наоборот, силы отталкивания преобладают, если В < В о с . Мы можем определить это количественно, считая, что работа, совершаемая при сжатии материала, будет равна увеличению внутренней энергии. Другими словами:

п с знак равно ( д ε с д В ) Т знак равно 0
что эквивалентно утверждению, что это изотермическое / изоэнтропическое уравнение для холодного сжатия. Знак минус показывает, что если бы к телу приложили растягивающую силу , силы связи между атомами действовали бы как восстанавливающая сила. Наклон п с кривая в п знак равно 0 (или 1 атм) определяет сжимаемость материала при нормальных условиях (т.е. Т знак равно Т о 300   К ). Это дается:
κ о знак равно 1 В о ( В п ) Т о
Обратите внимание, что наклон κ о определяет скорость упругих волн внутри объекта. Таким образом, давайте определим скорость звука в твердом теле как эту скорость, определяемую выражением:
С о знак равно В 2 ( п В ) С
где индекс С указывает изоэнтропическую производную, а частная производная будет отрицательной, чтобы избежать мнимых скоростей звука.

Простая аппроксимация нулевого порядка
Мое рефлекторное предположение состоит в том, что простейший подход, учитывая, что вы допускаете упругие отношения столкновения, состоит в том, чтобы просто аппроксимировать Δ ε с по конечной кинетической энергии вашего ударяющего объекта, предполагая, что ударник (?) не двигается после удара.

Аппроксимация первого
порядка

Ниже мы рассмотрим эффекты на цилиндрическом стержне (используемом для симметрии и простоты).

При малых деформациях относительное изменение объема Δ В / В , дан кем-то:

Δ В В знак равно κ   п знак равно п К
куда К знак равно 1 / κ объемный модуль .

Давайте определим С 1 как скорость волны сжатия в материале из-за приложения постоянного давления, п , приложенный к одному концу стержня в некоторый начальный момент времени. Материал между фронтом волны и концом стержня сжимается с постоянной скоростью, ты . В этих условиях мы можем использовать закон Гука и показать, что для малых нагрузок и деформаций имеем:

ты С 1 знак равно п Е
куда Е модуль Юнга . Через некоторое время, т , масса материала, охваченная волной, приобретет импульс р   С 1   т   ты , который должен быть равен п   т из закона Ньютона, который дает нам:
п знак равно р   ты   С 1

Более подробный ответ
К сожалению, у меня нет времени на полный вывод, но я предлагаю главу XI в Справке 1 и использовать Ссылку 2 для вспомогательной информации. Зельдович и Райзер в основном тратят всю главу XI на обсуждение этой темы и вникают во все нюансы, которые применимы к вашей проблеме (например, волна давления, вызванная сжатием ударной волны). Я предполагаю, что многое из этого потребует численного анализа, но есть отправные точки и аналитические приближения, которые, вероятно, сэкономят вам много времени.

использованная литература

  1. Зельдович Я.Б., Ю.П. Райзер (2002) Физика ударных волн и высокотемпературных гидродинамических явлений , под ред. WD Hayes и RF Probstein, Mineola, NY, Dover Publications, Inc., The Dover Edition; ISBN-13: 978-0486420028.
  2. Уитэм, Великобритания (1999), Линейные и нелинейные волны , Нью-Йорк, штат Нью-Йорк: John Wiley & Sons, Inc.; ISBN: 0-471-35942-4.