Исследования прочности коры нейтронных звезд имеют фундаментальное \r\nзначение для моделирования магнитарных вспышек, глитчей радио-пульсаров \r\nи излучения нейтронными звездами гравитационных волн. В данной работе \r\nпредпринята попытка пролить свет на эту проблему путем изучения \r\nдеформаций растяжения и сжатия идеальных объемно-центрированных \r\nкубических кулоновских кристаллов. Наш анализ основан на полуаналитическом подходе Байко и Кожберова (2017), который позволяет для любой однородной \r\nдеформации рассчитать в полностью нелинейном режиме критические \r\nпараметры деформации, при превышении которых решетка теряет \r\nдинамическую устойчивость. Мы определили предельные величины деформации\r\nи анизотропии давления для любого \r\nнаправления растяжения или сжатия относительно кристаллографических \r\nосей. Эти величины оказались сильно анизотропными: они могут \r\nизменяться почти в 10 раз в зависимости от ориентации оси деформации.\r\nДля поликристаллической коры максимальная величина упругих \r\nрастяжений и сжатий в нашей модели равна 0.04. Это меньше, чем \r\nпредельная деформация 0.1, полученная с помощью МД \r\nсимуляций сдвиговых деформаций Хоровицем и Кадау (2009). Максимальная \r\nанизотропия давления поликристаллического вещества, согласно нашим \r\nоценкам, лежит в диапазоне от 0.005 до 0.014 $n Z^2 e^2/a$. Мы обсуждаем механизмы пластического движения, \r\nобразование больших кристаллитов в коре нейтронной звезды, а также \r\nанализируем энерговыделение в результате разрушения таких кристаллитов \r\nв контексте эволюции магнитного поля и вспышечной активности магнитаров.