ПОИСК Статьи Чертежи Таблицы Метод Гаусса из "Расчет машиностроительных конструкций методом конечных элементов " В программе, реализующей метод Гаусса, обрабатывается полная лента матрицы коэффициентов системы (см. рис. 3.1). Для осуществления t-ro шага необходимо иметь т + 1 строку (заштрихованная область на рис. 3.4). Поскольку не всегда эти строки могут разместиться в. оперативной памяти, их обработку осуществляют порциями. Размер порции ограничен доступным объемом оперативной памяти, шириной ленты матрицы и принятым алгоритмом обработки порций, сущность которого заключается в следующем. [c.31] На i-и шаге реализации процедуры прямого хода по Гауссу в массив В из рабочего файла считывается первая порция не полностью обработанных элементов (заштрихованная область на рис. 3.5). Первую половину ленты переписываем в массив А и сохраняем до конца t-ro шага. Элементы массива А с учетом свойства симметрии матрицы жесткости служат исходной информацией для обработки последующих порций на i-м шаге. [c.32] После обработки всех строк порции 1 последняя в преобразованном виде записывается в рабочий файл. Далее последовательно считываются порции 2, 3, 4, 5 (их число зависит от конкретной задачи) и преобразуются строки, расположенные в зонах 12, 13, /4, /5. На (t + 1)-м шаге перечисленные операции повторяются с той лишь разницей, что вместо порции 1 используется порция 2 -го шага. [c.32] Здесь формальные параметры те же, что и в процедуре LDLF1 (см. подразд. 3.2). Отличие в том, что в файле FL размещаются коэффициенты системы уравнений в виде массива А2 (NM, —М М + NQL). В результате выполнения процедуры GAUS1 из предварительно заполненного файла FL порциями считывается массив А2, который после преобразования вновь записывается в рабочий файл (при этом левая половина ленты оказывается обнуленной). [c.32] Формальные параметры этих процедур те же, что у процедур FA TB и FA T1. Они обеспечивают решение методом Гаусса систем уравнений соответственно с ленточной и полной матрицей коэффициентов. [c.33] Вернуться к основной статье