Физический закон: сила, необходимая для деформации пружины, линейно зависит от расстояния.
В физике закон Гука — это эмпирический закон , который гласит, что сила ( F ) , необходимая для растяжения или сжатия пружины на некоторое расстояние ( x ), линейно масштабируется относительно этого расстояния, то есть F s = kx , где k — постоянный фактор, характерный для пружины (т. е. ее жесткость ), а x мало по сравнению с полной возможной деформацией пружины. Закон назван в честь британского физика 17-го века Роберта Гука . Он впервые сформулировал закон в 1676 году в виде латинской анаграммы . [1] [2] Он опубликовал решение своей анаграммы в 1678 году [3] как: ut tensio, sic vis («как растяжение, так и сила» или «растяжение пропорционально силе»). Гук утверждает в работе 1678 года, что он знал об этом законе с 1660 года.
Уравнение Гука справедливо (в некоторой степени) во многих других ситуациях, когда упругое тело деформируется , например, когда ветер дует на высокое здание, или когда музыкант дергает струну гитары . Упругое тело или материал, для которого можно предположить это уравнение, называется линейно-упругим или гуковским .
Закон Гука — это всего лишь линейное приближение первого порядка к реальному отклику пружин и других упругих тел на приложенные силы. Он должен в конечном итоге перестать работать, как только силы превысят некоторый предел, поскольку ни один материал не может быть сжат сверх определенного минимального размера или растянут сверх максимального размера без некоторой постоянной деформации или изменения состояния. Многие материалы будут заметно отклоняться от закона Гука задолго до того, как будут достигнуты эти пределы упругости .
С другой стороны, закон Гука является точным приближением для большинства твердых тел, пока силы и деформации достаточно малы. По этой причине закон Гука широко используется во всех отраслях науки и техники и является основой многих дисциплин, таких как сейсмология , молекулярная механика и акустика . Он также является фундаментальным принципом, лежащим в основе пружинных весов , манометра , гальванометра и балансира механических часов .
Современная теория упругости обобщает закон Гука, утверждая, что напряжение (деформация) упругого объекта или материала пропорционально приложенному к нему напряжению . Однако, поскольку общие напряжения и деформации могут иметь несколько независимых компонентов, «коэффициент пропорциональности» может больше не быть просто одним действительным числом, а скорее линейным отображением ( тензором ), которое может быть представлено матрицей действительных чисел.
В этой общей форме закон Гука позволяет вывести соотношение между деформацией и напряжением для сложных объектов с точки зрения внутренних свойств материалов, из которых они сделаны. Например, можно вывести, что однородный стержень с равномерным поперечным сечением будет вести себя как простая пружина при растяжении, с жесткостью k, прямо пропорциональной площади его поперечного сечения и обратно пропорциональной его длине.
Формальное определение
Линейные пружины
Рассмотрим простую винтовую пружину, один конец которой прикреплен к некоторому фиксированному объекту, в то время как свободный конец тянется силой, величина которой равна F s . Предположим, что пружина достигла состояния равновесия , в котором ее длина больше не меняется. Пусть x будет величиной, на которую свободный конец пружины был смещен из своего «расслабленного» положения (когда он не растягивается). Закон Гука гласит, что F s = k x {\displaystyle F_{s}=kx} или, что то же самое,
где k — положительное действительное число, характерное для пружины. Пружина с промежутками между витками может быть сжата, и та же формула справедлива для сжатия, причем F s и x в этом случае отрицательны. [4]
Согласно этой формуле график приложенной силы F s как функции смещения x будет представлять собой прямую линию, проходящую через начало координат , наклон которой равен k .
Закон Гука для пружины также формулируется в соответствии с соглашением, что F s — это восстанавливающая сила, оказываемая пружиной на то, что тянет ее свободный конец. В этом случае уравнение становится таким, поскольку направление восстанавливающей силы противоположно направлению смещения.
Пружины кручения
Торсионный аналог закона Гука применим к торсионным пружинам . Он гласит, что крутящий момент (τ), необходимый для вращения объекта, прямо пропорционален угловому смещению (θ) от положения равновесия. Он описывает связь между крутящим моментом, приложенным к объекту, и результирующей угловой деформацией из -за кручения. Математически это можно выразить как:
k — крутильная постоянная (измеряется в Н·м/радиан), характеризующая жесткость торсионной пружины или сопротивление угловому смещению.
θ — угловое смещение (измеряется в радианах) от положения равновесия.
Как и в линейном случае, этот закон показывает, что крутящий момент пропорционален угловому смещению, а отрицательный знак указывает на то, что крутящий момент действует в направлении, противоположном угловому смещению, обеспечивая возвращающую силу для возвращения системы в состояние равновесия.
Общие "скалярные" пружины
Закон упругости Гука обычно применим к любому упругому объекту произвольной сложности, если и деформация, и напряжение могут быть выражены одним числом, которое может быть как положительным, так и отрицательным.
Например, когда блок резины, прикрепленный к двум параллельным пластинам, деформируется сдвигом , а не растяжением или сжатием, сдвигающая сила Fs и боковое смещение пластин x подчиняются закону Гука (для достаточно малых деформаций).
Закон Гука также применим, когда прямой стальной стержень или бетонная балка (например, те, что используются в зданиях), поддерживаемые с обоих концов, изгибаются под действием веса F, помещенного в некоторой промежуточной точке. Смещение x в этом случае представляет собой отклонение балки, измеренное в поперечном направлении, относительно ее ненагруженной формы.
Векторная формулировка
В случае спиральной пружины, которая растягивается или сжимается вдоль своей оси , приложенная (или восстанавливающая) сила и результирующее удлинение или сжатие имеют одно и то же направление (которое является направлением указанной оси). Поэтому, если F s и x определены как векторы , уравнение Гука по-прежнему справедливо и гласит, что вектор силы является вектором удлинения, умноженным на фиксированный скаляр .
Общая тензорная форма
Некоторые упругие тела будут деформироваться в одном направлении, если подвергнуть их воздействию силы с другим направлением. Одним из примеров является горизонтальная деревянная балка с неквадратным прямоугольным поперечным сечением, которая изгибается поперечной нагрузкой, которая не является ни вертикальной, ни горизонтальной. В таких случаях величина смещения x будет пропорциональна величине силы F s , пока направление последней остается тем же (и ее значение не слишком велико); поэтому скалярная версия закона Гука F s = − kx будет верна. Однако векторы силы и смещения не будут скалярными кратными друг другу, поскольку они имеют разные направления. Более того, отношение k между их величинами будет зависеть от направления вектора F s .
Тем не менее, в таких случаях часто существует фиксированная линейная связь между векторами силы и деформации, пока они достаточно малы. А именно, существует функция κ от векторов к векторам, такая, что F = κ ( X ) , и κ ( α X 1 + β X 2 ) = α κ ( X 1 ) + β κ ( X 2 ) для любых действительных чисел α , β и любых векторов смещения X 1 , X 2 . Такая функция называется тензором (второго порядка) .
Относительно произвольной декартовой системы координат векторы силы и смещения могут быть представлены матрицами действительных чисел размером 3 × 1. Тогда тензор κ, соединяющий их, может быть представлен матрицей действительных коэффициентов размером 3 × 3 κ , которая при умножении на вектор смещения дает вектор силы:
То есть для i = 1, 2, 3. Поэтому можно сказать, что закон Гука F = κ X выполняется также и тогда, когда X и F являются векторами с переменными направлениями, за исключением того, что жесткость объекта является тензором κ , а не одним действительным числом k .
Закон Гука для сплошных сред
Напряжения и деформации материала внутри сплошного упругого материала (такого как блок резины, стенка котла или стальной стержень) связаны линейной зависимостью, которая математически аналогична закону упругости Гука и часто называется этим именем.
Однако состояние деформации в твердой среде вокруг некоторой точки не может быть описано одним вектором. Один и тот же кусок материала, каким бы малым он ни был, может быть сжат, растянут и сдвинут одновременно в разных направлениях. Аналогично, напряжения в этом куске могут быть одновременно толкающими, тянущими и сдвигающими.
Чтобы охватить эту сложность, соответствующее состояние среды вокруг точки должно быть представлено двумя тензорами второго порядка, тензором деформации ε (вместо смещения X ) и тензором напряжения σ (заменяющим восстанавливающую силу F ). Аналог закона пружины Гука для сплошных сред тогда имеет вид где c — тензор четвертого порядка (то есть линейное отображение между тензорами второго порядка), обычно называемое тензором жесткости или тензором упругости . Его также можно записать как где тензор s , называемый тензором податливости , представляет собой обратное указанному линейному отображению.
В декартовой системе координат тензоры напряжений и деформаций могут быть представлены матрицами 3 × 3.
Будучи линейным отображением между девятью числами σ ij и девятью числами ε kl , тензор жесткости c представлен матрицей 3 × 3 × 3 × 3 = 81 действительных чисел c ijkl . Закон Гука тогда гласит, что
где i , j = 1,2,3 .
Все три тензора, как правило, изменяются от точки к точке внутри среды и могут также меняться со временем. Тензор деформации ε просто определяет смещение частиц среды в окрестности точки, в то время как тензор напряжения σ определяет силы, которые соседние участки среды оказывают друг на друга. Поэтому они не зависят от состава и физического состояния материала. Тензор жесткости c , с другой стороны, является свойством материала и часто зависит от переменных физического состояния, таких как температура, давление и микроструктура .
Из-за присущих симметрий σ , ε и c только 21 упругий коэффициент последнего являются независимыми. [6] Это число может быть дополнительно уменьшено с помощью симметрии материала: 9 для орторомбического кристалла, 5 для гексагональной структуры и 3 для кубической симметрии. [7] Для изотропных сред (которые имеют одинаковые физические свойства в любом направлении) c может быть уменьшен только до двух независимых чисел, объемного модуля K и модуля сдвига G , которые количественно определяют сопротивление материала изменениям объема и сдвиговым деформациям соответственно.
В частности, тензорное уравнение σ = cε, связывающее упругие напряжения с деформациями, полностью аналогично уравнению τ = με̇, связывающему тензор вязких напряжений τ и тензор скорости деформации ε̇ в потоках вязких жидкостей; хотя первое относится к статическим напряжениям (связанным с величиной деформации), а второе — к динамическим напряжениям (связанным со скоростью деформации).
Единицы измерения
В единицах СИ смещения измеряются в метрах (м), а силы в ньютонах (Н или кг·м/с 2 ). Таким образом, жесткость пружины k и каждый элемент тензора κ измеряются в ньютонах на метр (Н/м) или килограммах на секунду в квадрате (кг/с 2 ).
Для сплошных сред каждый элемент тензора напряжений σ представляет собой силу, деленную на площадь; поэтому он измеряется в единицах давления, а именно паскалях (Па, или Н/м 2 , или кг/(м·с 2 ). Элементы тензора деформаций ε безразмерны (перемещения , деленные на расстояния). Поэтому записи c ijkl также выражаются в единицах давления.
Общее применение для эластичных материалов
Объекты, которые быстро восстанавливают свою первоначальную форму после деформации под действием силы, при этом молекулы или атомы их материала возвращаются в исходное состояние устойчивого равновесия, часто подчиняются закону Гука.
Закон Гука справедлив только для некоторых материалов при определенных условиях нагрузки. Сталь проявляет линейно-упругое поведение в большинстве инженерных приложений; закон Гука справедлив для нее во всем ее упругом диапазоне (т. е. для напряжений ниже предела текучести ). Для некоторых других материалов, таких как алюминий, закон Гука справедлив только для части упругого диапазона. Для этих материалов определено пропорциональное предельное напряжение, ниже которого ошибки, связанные с линейным приближением, пренебрежимо малы.
Резина обычно считается «негуковским» материалом, поскольку ее эластичность зависит от напряжения и чувствительна к температуре и скорости нагрузки.
Стержень из любого упругого материала можно рассматривать как линейную пружину . Стержень имеет длину L и площадь поперечного сечения A. Его растягивающее напряжение σ линейно пропорционально его относительному удлинению или деформации ε по модулю упругости E :
Модуль упругости часто можно считать постоянным. В свою очередь,
(то есть дробное изменение длины), и поскольку
следует, что:
Изменение длины можно выразить как
Весенняя энергия
Потенциальная энергия U el ( x ), хранящаяся в пружине, определяется как которая получается путем сложения энергии, необходимой для постепенного сжатия пружины. То есть интеграла силы по смещению. Поскольку внешняя сила имеет то же общее направление, что и смещение, потенциальная энергия пружины всегда неотрицательна. Подстановка дает
Этот потенциал U el можно визуализировать как параболу на плоскости Ux, такую, что U el ( x ) = 1/2 kx 2 . При растяжении пружины в положительном направлении x потенциальная энергия увеличивается параболически (то же самое происходит при сжатии пружины). Поскольку изменение потенциальной энергии изменяется с постоянной скоростью:
Обратите внимание, что изменение изменения U постоянно, даже когда смещение и ускорение равны нулю.
Релаксированные силовые константы (обратные обобщенным константам податливости ) однозначно определены для молекулярных систем, в отличие от обычных «жестких» силовых констант, и, таким образом, их использование позволяет устанавливать значимые корреляции между силовыми полями, рассчитанными для реагентов , переходных состояний и продуктов химической реакции . Так же, как потенциальная энергия может быть записана в виде квадратичной формы во внутренних координатах, ее также можно записать в терминах обобщенных сил. Полученные коэффициенты называются константами податливости . Существует прямой метод расчета константы податливости для любой внутренней координаты молекулы без необходимости проведения анализа нормального режима. [8] Пригодность расслабленных силовых констант (обратных констант податливости) в качестве дескрипторов прочности ковалентных связей была продемонстрирована еще в 1980 году. Недавно была продемонстрирована также пригодность в качестве дескрипторов прочности нековалентных связей. [9]
Гармонический осциллятор
Масса m, прикрепленная к концу пружины, является классическим примером гармонического осциллятора . Если слегка потянуть за массу и затем отпустить ее, система будет приведена в синусоидальное колебательное движение около положения равновесия. В той степени, в которой пружина подчиняется закону Гука, и можно пренебречь трением и массой пружины, амплитуда колебания останется постоянной; а ее частота f будет независима от ее амплитуды, определяемой только массой и жесткостью пружины:
Это явление сделало возможным создание точных механических часов , которые можно было носить на кораблях и в карманах людей.
Вращение в пространстве без гравитации
Если бы масса m была прикреплена к пружине с силой постоянной k и вращалась в свободном пространстве, то натяжение пружины ( F t ) обеспечивало бы требуемую центростремительную силу ( F c ):
Поскольку F t = F c и x = r , то:
Учитывая, что ω = 2π f , это приводит к тому же уравнению частоты, что и выше:
Линейная теория упругости для сплошных сред
Изотропные материалы
Изотропные материалы характеризуются свойствами, которые не зависят от направления в пространстве. Физические уравнения, включающие изотропные материалы, должны, следовательно, быть независимыми от системы координат, выбранной для их представления. Тензор деформации является симметричным тензором. Поскольку след любого тензора не зависит от любой системы координат, наиболее полное разложение симметричного тензора без учета координат состоит в его представлении в виде суммы постоянного тензора и бесследового симметричного тензора. [10] Таким образом, в индексной записи :
где δ ij — дельта Кронекера . В прямой тензорной записи:
где I — тензор тождественности второго порядка.
Первый член справа — постоянный тензор, также известный как тензор объемной деформации , а второй член — бесследовый симметричный тензор, также известный как тензор девиаторной деформации или тензор сдвига.
Наиболее общую форму закона Гука для изотропных материалов теперь можно записать в виде линейной комбинации этих двух тензоров:
где K — модуль объемной упругости , а G — модуль сдвига .
Используя соотношения между упругими модулями , эти уравнения могут быть выражены и другими способами. Распространенная форма закона Гука для изотропных материалов, выраженная в прямой тензорной записи, есть [11],
где λ = K − 2/3 G = c 1111 − 2 c 1212 и μ = G = c 1212 — константы Ламе , I — тензор идентичности второго ранга, а I — симметричная часть тензора идентичности четвертого ранга. В индексной записи:
Обратное соотношение равно [12]
Следовательно, тензор податливости в соотношении ε = s : σ равен
В терминах модуля Юнга и коэффициента Пуассона закон Гука для изотропных материалов может быть выражен как
Это форма, в которой деформация выражается через тензор напряжений в машиностроении. Выражение в развернутом виде равно
где E — модуль Юнга , а ν — коэффициент Пуассона . (См. 3-D эластичность ).
Вывод закона Гука в трех измерениях
Трехмерную форму закона Гука можно вывести с помощью коэффициента Пуассона и одномерной формы закона Гука следующим образом. Рассмотрим соотношение деформации и напряжения как суперпозицию двух эффектов: растяжения в направлении нагрузки (1) и сжатия (вызванного нагрузкой) в перпендикулярных направлениях (2 и 3),
где ν — коэффициент Пуассона, а E — модуль Юнга.
Аналогичные уравнения мы получаем для нагрузок в направлениях 2 и 3,
и
Суммируя три случая вместе ( ε i = ε i ′ + ε i ″ + ε i ‴ ), мы получаем
или путем прибавления и вычитания одного νσ и далее, решая σ 1 ,
получаем
Вычислив сумму
и подставив ее в уравнение, решенное относительно σ 1 , получим
, где μ и λ — параметры Ламе .
Аналогичная обработка направлений 2 и 3 дает закон Гука в трех измерениях.
В матричной форме закон Гука для изотропных материалов можно записать как
где γ ij = 2 ε ij — инженерная деформация сдвига . Обратное соотношение можно записать как
что можно упростить благодаря константам Ламе:
В векторной записи это становится
где I — единичный тензор.
Обратное отношение обычно записывается в сокращенной форме
Плоская деформация
В условиях плоской деформации ε 31 = ε 13 = ε 32 = ε 23 = ε 33 = 0. В этом случае закон Гука принимает вид
Анизотропные материалы
Симметрия тензора напряжений Коши ( σ ij = σ ji ) и обобщенных законов Гука ( σ ij = c ijkl ε kl ) подразумевает, что c ijkl = c jikl . Аналогично, симметрия тензора бесконечно малых деформаций подразумевает, что c ijkl = c ijlk . Эти симметрии называются малыми симметриями тензора жесткости c . Это уменьшает число упругих констант с 81 до 36.
Если в дополнение, поскольку градиент смещения и напряжение Коши сопряжены по работе, соотношение напряжение-деформация может быть выведено из функционала плотности энергии деформации ( U ), то
Произвольность порядка дифференциации подразумевает, что c ijkl = c klij . Они называются главными симметриями тензора жесткости. Это уменьшает число упругих констант с 36 до 21. Главная и второстепенная симметрии указывают, что тензор жесткости имеет только 21 независимую компоненту.
Матричное представление (тензор жесткости)
Часто бывает полезно выразить анизотропную форму закона Гука в матричной записи, также называемой записью Фойгта . Для этого мы воспользуемся симметрией тензоров напряжений и деформаций и выразим их как шестимерные векторы в ортонормированной системе координат ( e 1 , e 2 , e 3 ) как
Тогда тензор жесткости ( c ) можно выразить как
и закон Гука запишется как
Аналогично тензор податливости ( s ) можно записать как
Изменение системы координат
Если линейный упругий материал поворачивается из одной исходной конфигурации в другую, то материал симметричен относительно поворота, если компоненты тензора жесткости в повернутой конфигурации связаны с компонентами в исходной конфигурации соотношением [13]
, где l ab — компоненты ортогональной матрицы поворота [ L ] . Такое же соотношение справедливо и для инверсий.
В матричной записи, если преобразованный базис (повернутый или инвертированный) связан с опорным базисом соотношением
, то
Кроме того, если материал симметричен относительно преобразования [ L ] , то
Ортотропные материалы
Ортотропные материалы имеют три ортогональные плоскости симметрии . Если базисные векторы ( e 1 , e 2 , e 3 ) являются нормалями к плоскостям симметрии, то соотношения преобразования координат подразумевают, что
Обратное отношение этого соотношения обычно записывается как [14] [ нужна страница ]
где
G ij — модуль сдвига в направлении j на плоскости, нормаль которой направлена в направлении i.
ν ij — коэффициент Пуассона , соответствующий сжатию в направлении j при приложении растяжения в направлении i .
В условиях плоского напряженного состояния σ zz = σ zx = σ yz = 0 закон Гука для ортотропного материала принимает вид
Обратное соотношение имеет
вид Часто также используется транспонированная форма приведенной выше матрицы жесткости.
Поперечно-изотропные материалы
Трансверсально изотропный материал симметричен относительно вращения вокруг оси симметрии . Для такого материала, если e 3 является осью симметрии, закон Гука можно выразить как
Чаще всего ось x ≡ e 1 принимают за ось симметрии, а обратный закон Гука записывают в виде [15]
Универсальный индекс упругой анизотропии
Для понимания степени анизотропии любого класса был сформулирован универсальный индекс упругой анизотропии (AU) [16] . Он заменяет отношение Ценера , которое подходит для кубических кристаллов .
Термодинамическая основа
Линейные деформации упругих материалов можно аппроксимировать как адиабатические . При этих условиях и для квазистатических процессов первый закон термодинамики для деформированного тела можно выразить как
где δU - увеличение внутренней энергии , а δW - работа, совершаемая внешними силами. Работу можно разделить на два члена
, где δW s - работа, совершаемая поверхностными силами , а δW b - работа, совершаемая объемными силами . Если δ u - это изменение поля смещения u в теле, то два члена внешней работы можно выразить как
где t - вектор поверхностного натяжения , b - вектор объемной силы, Ω представляет тело, а ∂ Ω представляет его поверхность. Используя соотношение между напряжением Коши и поверхностным натяжением, t = n · σ (где n — единичная внешняя нормаль к ∂ Ω ), мы имеем
Преобразование интеграла поверхности в объемный интеграл с помощью теоремы о расходимости дает
Используя симметрию напряжения Коши и тождество,
мы имеем следующее
Из определения деформации и из уравнений равновесия имеем
Следовательно, мы можем записать
и, следовательно, изменение плотности внутренней энергии определяется как
Упругий материал определяется как материал, в котором полная внутренняя энергия равна потенциальной энергии внутренних сил (также называемой энергией упругой деформации ). Следовательно, плотность внутренней энергии является функцией деформаций, U 0 = U 0 ( ε ) , а изменение внутренней энергии можно выразить как
Поскольку изменение деформации произвольно, соотношение напряжение-деформация упругого материала определяется как
Для линейного упругого материала величина ∂ U0/∂ ε является линейной функцией ε и, следовательно, может быть выражена как ,
где c — тензор четвертого ранга материальных констант, также называемый тензором жесткости . Мы можем понять, почему c должен быть тензором четвертого ранга, заметив, что для линейного упругого материала
В индексной нотации
Константа правой стороны требует четырех индексов и является величиной четвертого ранга. Мы также можем увидеть, что эта величина должна быть тензором, поскольку она является линейным преобразованием, которое переводит тензор деформации в тензор напряжения. Мы также можем показать, что константа подчиняется правилам преобразования тензора для тензоров четвертого ранга.
^ Анаграмма была дана в алфавитном порядке, ceiiinosssttuv , представляющая Ut tensio, sic vis – «Как расширение, так и сила»: Петроски, Генри (1996). Изобретение по замыслу: как инженеры переходят от мысли к вещи . Кембридж, Массачусетс: Издательство Гарвардского университета. стр. 11. ISBN 978-0674463684.
^ См. http://civil.lindahall.org/design.shtml, где также можно найти анаграмму для слова «цепная линия» .
^ Роберт Гук , De Potentia Restitutiva, или о пружине. Объяснение силы пружинящих тел , Лондон, 1678.
^ Янг, Хью Д.; Фридман, Роджер А.; Форд, А. Льюис (2016). Университетская физика Сирса и Земанского: с современной физикой (14-е изд.). Пирсон. стр. 209.
^ Беленький; Салаев (1988). «Эффекты деформации в слоистых кристаллах». Успехи физических наук . 155 (5): 89. doi : 10.3367/УФНр.0155.198805c.0089 .
^ Mouhat, Félix; Coudert, François-Xavier (5 декабря 2014 г.). «Необходимые и достаточные условия упругой устойчивости в различных кристаллических системах». Physical Review B. 90 ( 22): 224104. arXiv : 1410.0065 . Bibcode : 2014PhRvB..90v4104M. doi : 10.1103/PhysRevB.90.224104. ISSN 1098-0121. S2CID 54058316.
^ Виджай Мадхав, М.; Маногаран, С. (2009). «Пересмотр констант соответствия в избыточных внутренних координатах и некоторые новые идеи». J. Chem. Phys . 131 (17): 174112–174116. Bibcode : 2009JChPh.131q4112V. doi : 10.1063/1.3259834. PMID 19895003.
^ Пономарева, Алла; Юренко, Евгений; Журакивский, Роман; Ван Мурик, Таня; Говорун, Дмитрий (2012). «Полное конформационное пространство потенциальных ингибиторов обратной транскриптазы ВИЧ-1 d4U и d4C. Квантово-химическое исследование». Phys. Chem. Chem. Phys . 14 (19): 6787–6795. Bibcode :2012PCCP...14.6787P. doi :10.1039/C2CP40290D. PMID 22461011.
^ Symon, Keith R. (1971). "Глава 10". Механика . Рединг, Массачусетс: Addison-Wesley. ISBN9780201073928.
^ Симо, Дж. К.; Хьюз, Т. Дж. Р. (1998). Вычислительная неэластичность . Springer. ISBN9780387975207.
^ Милтон, Грэм В. (2002). Теория композитов . Кембриджские монографии по прикладной и вычислительной математике. Cambridge University Press. ISBN9780521781251.
^ Слотер, Уильям С. (2001). Линеаризованная теория упругости . Биркхойзер. ISBN978-0817641177.
^ Бореси, А. П.; Шмидт, Р. Дж.; Сайдботтом, О. М. (1993). Advanced Mechanics of Materials (5-е изд.). Wiley. ISBN9780471600091.