Я пытаюсь проверить результаты из пользовательского кода, написанного с Matlab, сравнивая результаты с Abaqus. Я делаю неявный динамический анализ с использованием прямой интеграции и, как и в Abaqus, я использую оператор, определенный Hilber, Hughes и Taylor (HHT).
Моя проблема в том, что у меня теперь есть скрипт Matlab, способный идеально соответствовать результатам Abaqus, когда я использую элементы луча или континуума (B21 и C3D8R соответственно), но не с оболочками (S4R или S4).
Я действительно застрял в этой проблеме и не понимаю, что приводит к различиям в результатах. Временная интеграция должна быть независимой от типа элемента ... Есть ли у вас какие-либо идеи или ресурсы о том, что Abaqus может сделать по-другому в случае элементов оболочки?