Программа прочность сосудов и аппаратов
Проблемы…
Расчет сосудов и аппаратов на прочность и устойчивость — работа довольно сложная и трудоемкая. Расчет отдельных элементов в различных условиях регламентирован большим количеством разрозненных нормативных документов (ГОСТ, ОСТ, РД, РТМ и др.), но в то же время на некоторые виды расчетов (например, расчет прочности штуцеров) отечественная нормативная документация или отсутствует, или ее область применения существенно ограничена. При расчетах зачастую приходится использовать множество вспомогательных величин (расчетная длина, эффективный момент инерции сечения, период собственных колебаний и т.д.), получение которых трудоемко или требует применения методов, не описанных в нормативной литературе.
Эти проблемы каждый специалист решает по-своему, используя ручной расчет, применяя собственные программы для расчета отдельных элементов или обращаясь в специализированные организации. Все это не лучшие решения: вручную — долго и очень непросто, собственные программы используются на свой страх и риск и не всегда применимы, расчет на стороне — долго и дорого… Наиболее опытные используют специальные методы (например, конечно-элементный анализ) — хорошие результаты, но дорого и трудоемко. Существует ряд зарубежных программ (Compress, Vessel, PV Elite и др.) — весьма неплохих, однако, увы, недешевых и ориентированных на зарубежную нормативную базу. В итоге проектировщик постоянно оказывается перед выбором, ищет выход… И этот выход есть!
…и решения
В 2004 году НТП «Трубопровод» представил программу ПАССАТ («Прочностной Анализ Состояния Сосудов, Аппаратов, Теплообменнико»), а в январе 2006 года была выпущена очередная версия 1.03. ПАССАТ — не только первая отечественная программа на платформе Windows, предназначенная для прочностных расчетов сосудов и аппаратов, но и, по нашему мнению, наиболее развитая и удобная для пользователя.
Рис. 1. Программа ПАССАТ
ПАССАТ позволяет рассчитывать большое число элементов аппарата в рабочих условиях, а также в условиях монтажа и испытаний, учитывает взаимное влияние элементов и требует в качестве входных данных только известные проектировщику величины: размеры и взаимное расположение элементов, материал, известные дополнительные нагрузки и т.п.
Особенностью программы является отображение трехмерной модели рассчитываемого аппарата, что позволяет визуально контролировать введенные данные.
В настоящее время ПАССАТ состоит из трех модулей:
- Базовый модуль осуществляет ввод данных, отображение модели, расчет на прочность горизонтальных и вертикальных сосудов и аппаратов с формированием отчетов на основе ГОСТ 14249–89, ГОСТ 25221–82, ГОСТ 26202–84, ГОСТ 24755–89, РД 26−15−88, РД РТМ 26−01−96- 83, РД 10−249−98, РД 26.260.09−92, РД 26−01−169−89 и др.
- Модуль «ПАССАТ-Колонны» рассчитывает аппараты колонного типа на прочность и устойчивость от внешних, ветровых и сейсмических нагрузок на основе ГОСТ Р 51273−99, ГОСТ Р 51274−99 и др.
- Модуль «ПАССАТ-Штуцер» реализует иностранные методики расчета штуцеров и аппаратных фланцев WRC-107/297, ASME VIII и др.
Еще одной особенностью программы является форма представления результатов: пользователь получает не просто заключение о работоспособности, а полный протокол расчета, включая примененные формулы, ссылки на нормативные документы и промежуточные вычисления — программа полностью имитирует расчет вручную. Этим достигаются две цели: у проектировщика есть возможность проконтролировать весь ход расчета, понять результаты; если расчет неуспешен и нормы не выполнены, сразу ясно, «кто виноват». С другой стороны, такой протокол может быть представлен в надзорные органы для согласования.
Полируем старые методики
В процессе создания программы была проведена большая работа по изучению методик, лежащих в основе нормативных документов, предложены полезные расширения, позволяющие устранить имеющиеся ошибки и неточности, уточнить методику расчета отдельных величин, а также снять ряд ограничений. Вот некоторые примеры.
При расчете обечаек от воздействия опорных нагрузок, в отличие от методики, описанной в [3], программа позволяет рассчитывать сосуды произвольной конструкции, при этом опоры можно размещать в любом месте цилиндрических обечаек.
Рис. 2. Горизонтальный аппарат на седловых опорах
На основании [10, 15−19] впервые разработана комплексная методика расчета жесткости и прочности мест соединений штуцеров, работающих под действием внутреннего давления и внешних нагрузок. Часть этой методики (прочность штуцера) реализована в модуле «ПАССАТ-Штуцер”.
Рис. 3. Расчетная схема аппарата для расчета на прочность
Для арматурных фланцевых соединений расчет проводится в соответствии с [13] – с дополнительным учетом внешних нагрузок и изгибающих моментов, а также напряжений, вызванных разницей линейных удлинений материалов фланцев и крепежа при температурном воздействии, определяемых по [5].
При расчете сейсмических и ветровых нагрузок, действующих на вертикальные аппараты колонного типа, для определения периода собственных колебаний используется (в отличие от [7]) метод Рэлея, позволяющий рассчитать конструкцию с произвольным количеством элементов.
В качестве расчетной схемы колонного аппарата принимается консольный упруго защемленный стержень. Упругое защемление позволяет учитывать свойства грунта и фундамента колонны.
Определение наиболее слабого сечения опорной обечайки аппарата колонного типа — очень трудоемкая операция. ПАССАТ выполняет этот расчет самостоятельно, используя итерационный метод. Это позволяет эффективно рассчитывать опорные обечайки с любым числом и расположением патрубков различной формы.
ПАССАТ как инструмент проектирования
Расчеты, выполняемые программой ПАССАТ, — поверочные: пользователь проверяет допустимость нагрузок на заданную им конструкцию. При этом программа содержит и элементы проектирования, упрощающие процесс конструирования аппарата.
При вводе и редактировании элемента конструкции проектировщику доступна функция Определение расчетных величин, которая рассчитывает зависимые параметры для отдельного элемента. Например, по толщине стенки можно определить расчетное давление, и наоборот, по заданному давлению — требуемую толщину стенки.
Рис. 4. Окно расчетных величин
Еще одна полезная функция программы — подготовка отчетов в формате RTF. Отчеты формируются на основе шаблонов, созданных пользователем, что позволяет быстро получить результаты расчета в требуемом виде, а также изменять их в Microsoft Word или другом текстовом редакторе. В составе программы поставляется шаблон отчета, оформленный в соответствии с ЕСКД, благодаря чему вы можете сразу же получить соответствующий раздел проекта.
Рис. 5. Отчет в формате RTF (MS Word)
Просто, быстро, удобно
Программа ПАССАТ очень проста в работе: пользователь задает конструкцию аппарата поэлементно, вводя известные ему величины (геометрические размеры, местоположение, материал, температуру, условия испытаний, нагрузки на аппарат). Правильность ввода контролируется визуально. Существуют режимы прозрачного и каркасного представления модели, а также средства навигации, позволяющие детально рассмотреть устройство аппарата — в том числе элементы, находящиеся внутри. Корректность параметров элементов (в том числе выполнение условий применимости методик) проверяется программой на этапе ввода — при обнаружении проблемы выдается соответствующее диагностическое сообщение. Диалоги ввода параметров элементов снабжены детальными схемами, иллюстрирующими назначение каждого параметра.
Рис. 6. Диалог параметров цилиндрической обечайки
При создании нового элемента часть параметров инициализируется значениями, зависящими от ранее введенных величин, что значительно сокращает время ввода.
Рис. 7. Фрагмент аппарата колонного типа
Рис. 8. Фрагмент аппарата колонного типа (прозрачный режим)
В соответствии со своими задачами пользователь может самостоятельно выбирать единицы измерения различных величин. Так, например, толщину стенки можно задавать в миллиметрах, а длину аппарата — в метрах или даже дюймах. При изменении единиц измерения ранее введенные величины будут корректно пересчитаны: значения хранятся во внутренних единицах программы, не зависящих от выбора пользователя.
ПАССАТ содержит открытую для пополнения базу данных, где собраны параметры наиболее часто используемых материалов.
Расчет выполняется нажатием на одну кнопку — через несколько секунд вы получаете отчет, причем проблемные элементы будут помечены цветом и вынесены в отдельные строки.
Рис. 9. Окно результатов расчета
Перспективы
ПАССАТ — быстро развивающаяся программа. Перечислим лишь некоторые из средств, которые будут доступны пользователю в следующих версиях:
- базы данных стандартных элементов (в первую очередь фланцев) позволят значительно сократить время ввода параметров;
- еще более гибкие средства создания отчетов обеспечат возможность получать выходные данные в произвольном виде — в том числе формировать краткие отчеты, не содержащие промежуточных вычислений, вставлять в отчет таблицы штуцеров и т.д.;
- модуль расчета теплообменных аппаратов будет выполнять расчеты по РД 26−15−88, РД 24.200.21−91 и др.;
- будут реализованы расчеты аппаратов с рубашками (ГОСТ 25867−83), расчет резервуаров для хранения нефтепродуктов (ПБ 03−381−00), расчет остаточного ресурса и др.;
- экспорт модели в другие системы позволит использовать модель аппарата при создании чертежей.
Литература
- ГОСТ 14249–89. Сосуды и аппараты. Нормы и методы расчета на прочность.
- ГОСТ 25221–82. Сосуды и аппараты. Днища и крышки сферические неотбортованные. Нормы и методы расчета на прочность.
- ГОСТ 26202–84. Сосуды и аппараты. Нормы и методы расчета на прочность обечаек и днищ от воздействия опорных нагрузок.
- ГОСТ 24755–89. Сосуды и аппараты. Нормы и методы расчета на прочность укрепления отверстий.
- РД 26−15−88. Сосуды и аппараты. Нормы и методы расчета на прочность и герметичность фланцевых соединений.
- РД РТМ 26−01−96−83. Крышки и днища плоские круглые с радиальными ребрами жесткости сосудов и аппаратов.
- ГОСТ Р 51273−99. Сосуды и аппараты. Нормы и методы расчета на прочность. Определение расчетных усилий для аппаратов колонного типа от ветровых нагрузок и сейсмических воздействий.
- ГОСТ Р 51274−99. Сосуды и аппараты. Аппараты колонного типа. Нормы и методы расчета на прочность.
- РД 26−02−62−98. Расчет на прочность элементов сосудов и аппаратов, работающих в коррозионно-активных сероводородосодержащих средах.
- Нормы расчета на прочность оборудования и трубопроводов атомных энергетических установок. — М.: Энергоатомиздат, 1989. — 525 с.
- ASME VIII, Div 1, 2002.
- WRC-107 Welding Research Council. Bulletin. «Local Stresses in Spherical and Cylindrical Shells due to External Loadings», 1979.
- WRC-297 Welding Research Council. Bulletin. «Local Stresses in Cylindrical Shells due to External Loadings on Nozzles — Supplement to WRC Bulletin № 107», 1987.
- BS-5500: 1976 Specification for Unfired fusion welded pressure vessels. British Standards Institution.
- WRC-368 Welding Research Council. Bulletin. «Stresses in Intersecting Cylinders subjected to Pressure». 1991. — 32 p.
- Bildy, Les M., 2000, «A Proposed Method for Finding Stress and Allowable Pressure in Cylinders with Radial Nozzles», PVP Vol. 399, ASME, New York, NY, pp. 77−82.
- ГОСТ 25867–83. Сосуды и аппараты. Сосуды с рубашками. Нормы и методы расчета на прочность.
- РД 26−02−63−87. Технические требования к конструированию и изготовлению сосудов, аппаратов и технологических блоков установок подготовки нефти и газа, работающих в средах, вызывающих сероводородное коррозионное растрескивание.
- РД 10−249−98. Нормы расчета на прочность стационарных котлов и трубопроводов пара и горячей воды.
- РД 26−01−169−89. Сосуды и аппараты. Нормы и методы расчета на прочность днищ в местах крепления опор-стоек.
- РД 24.200.21−91. Сосуды и аппараты. Нормы и методы расчета на прочность элементов плавающих головок кожухотрубчатых теплообменных аппаратов.
- РД 26−18−8−89. Сварные соединения приварки люков, штуцеров и муфт. Основные типы, конструктивные элементы и размеры.
- РД 26.260.09−92. Сосуды и аппараты. Нормы и методы расчета на прочность цилиндрических обечаек и выпуклых днищ в местах присоединения штуцеров при внешних статических нагрузках.
- C. E. Freese. Vibrations of Vertical Pressure Vessels. Journal of Engineering for Industry, 1959.
Источник
Программа ПАССАТ компании НТП «Трубопровод» позволяет производить расчет на прочность и устойчивость конструкций сосудов и аппаратов для оценки несущей способности в рабочих условиях, а также при испытаниях и монтаже.
ПАССАТ предназначен для организаций, проектирующих аппараты или сосуды. Название программы — аббревиатура, которая расшифровывается как «Прочностной анализ состояния сосудов, аппаратов, теплообменнико».
Пролог, в котором описываются идеи, заложенные в программу, ее назначение и используемые методы расчета
Разработчики программы ПАССАТ стремились создать эффективный, удобный и наглядный инструмент для автоматизации проектирования, поскольку на современном рынке отсутствуют отечественные программы, соответствующие нашим представлениям об инженерных инструментах.
Программу ПАССАТ выгодно отличают от зарубежных аналогов (Compress, Vessel, PV Elite и т.д.) ориентация на российскую нормативно-правовую базу, более доступная цена, интуитивно понятный русскоязычный интерфейс, наличие встроенной базы отечественных материалов.
В программе реализовано объемное графическое отображение геометрии отдельных элементов и модели в целом.
Рис. 1
Расчетная модель создается в трехмерной среде, что позволяет уже на этапе ввода исходных данных оценить габариты емкости и автоматически исключить нестыковку элементов модели из-за несовпадения размеров. Возможность вывода «каркасного» изображения (Wireless) обеспечивает полную видимость всех элементов, включая внутренние. ПАССАТ автоматически проверяет геометрию модели, позволяя пользователю выбрать способ стыковки элементов. Виртуальная расчетная модель, создаваемая в полном соответствии с заданным рассчитываемым аппаратом (сосудом), обеспечивает правильность ввода исходных размеров.
Программа создана в соответствии с основными нормативами и методиками, применяющимися в Российской Федерации. К сожалению, они не полностью регламентируют все стадии расчета, поэтому в отдельных случаях при расчете тех или иных элементов используются зарубежные методические документы.
Программа ПАССАТ производит расчет на основе:
- ГОСТ 14249–89;
- ГОСТ 25221–82;
- ГОСТ 26202–84;
- ГОСТ 24755–89;
- РД 26−15−88;
- РД РТМ 26−01−96−77;
- РД 26−02−62−98;
- норм расчета на прочность оборудования и трубопроводов атомных энергетических установок;
- ГОСТ 27772–88;
- ГОСТ 25859–83;
- ASME VIII;
- ASME II;
- WRC-107;
- WRC-297;
- BS-5500.
Для оценки прочности и устойчивости заданной конструкции сосуда или аппарата реализованы следующие функции:
- определение расчетных толщин и допускаемых значений давления, сил и моментов;
- автоматическое определение расчетных величин, таких как вес, длины и диаметры элементов, характеристики колец жесткости (в цилиндрических обечайках и в седловых опорах) и др.;
- расчет прочности места соединения штуцера с сосудом (аппаратом);
- расчет арматурных фланцевых соединений с учетом воздействия давления, внешних сил и моментов, температурных напряжений и т.д.
Кроме того, программа ПАССАТ обеспечивает выполнение ряда вспомогательных функций:
- настройка размерностей;
- ввод и анализ исходных данных. В случае, если пользователь не ввел всех данных, необходимых для выполнения расчета, или ввел их некорректно, программа выдает предупреждение до тех пор, пока все данные не будут заданы;
- автоматические изменения в смежных элементах всей модели при изменении геометрических параметров или условий нагружения в одном из элементов
- выбор используемых материалов из базы данных с возможностью ее пополнения;
- автоматическое изменение величины допускаемых напряжений, модулей упругости и т.д. при изменении материала, температуры или толщины стенки;
- формирование, просмотр и печать полного отчета по расчетам элементов модели с промежуточными результатами вычислений.
ПАССАТ содержит открытую базу данных материалов, применяемых в СНГ, доступную для пополнения и корректировки.
Программа предназначается для проектно-конструкторских бюро и отделов, которые специализируются на проектировании и реконструкции сосудов и аппаратов, работающих под давлением, а также для организаций нефтеперерабатывающей, химической, нефтехимической, газовой, нефтяной, теплоэнергетической и других отраслей промышленности.
Системные требования программы ПАССАТ весьма скромны: программа работает в среде Windows 9x/2000/XP, рекомендуется видеокарта с поддержкой OpenGL.
Глава об интерфейсе, исходных данных и результатах расчета
ПАССАТ имеет дружественный интерфейс, а также интуитивно понятную структуру создания моделей и расчета сосудов и аппаратов. Параметризованные элементы, применяемые в программе вместо геометрических примитивов, значительно упрощают работу и позволяют существенно сэкономить время.
На рисунках 2−4 в качестве примеров показано создание конической обечайки и цилиндрического элемента.
Рис. 2
Рис. 3
Рис. 4
Из иллюстраций видно, что для создания элемента, корректно пристыкованного к модели, инженеру достаточно заполнить поля в коротком диалоге, иллюстрированном наглядной схемой, и выбрать необходимую схему соединения.
Особо подчеркнем, что для работы с программой пользователю не требуется большой опыт работы с системами 3D-моделирования — достаточно владеть обычным набором знаний инженера-расчетчика.
Исходными данными в ПАССАТ являются тип, геометрические характеристики и материал элементов сосуда или аппарата, тип и расположение опор, вид испытаний, величины нагрузок. Выбор используемых материалов производится из базы данных.
Процесс завершается выдачей полного отчета по расчетам элементов модели с промежуточными результатами вычислений, который пользователь может либо экспортировать в Word, либо, при необходимости, модифицировать, внеся изменения в конструкцию и/или условия. Заметим, что экспорт расчета осуществляется поэлементно — это позволяет пользователю компоновать отчет в соответствии со своими требованиями.
В качестве примера приведен образец такого документа: в нем отражены результаты расчета эллиптического днища горизонтального сосуда в условиях гидроиспытаний.
Расчеты производятся по следующим элементам:
- цилиндрические обечайки (гладкие и подкрепленные кольцами жесткости);
- конические переходы;
- днища (сферические, эллиптические, торосферические, конические, плоские, сферические неотбортованные);
- врезки (штуцеры) в обечайки и выпуклые днища;
- седловые опоры в горизонтальных сосудах и аппаратах;
- опорные стойки и лапы в вертикальных сосудах и аппаратах.
Существующие конструкции сосудов и аппаратов, а также условия их работы часто не позволяют без значительного упрощения расчетных моделей выполнить расчеты в строгом соответствии с нормами, что приводит к искажению полученных результатов. После долгой и глубокой проработки методик нам удалось создать универсальную программу, значительно расширив сферу ее применения.
Расчет в условиях испытаний (гидроиспытания)
Условия нагружения при испытаниях:
Рабочая температура, T: 20 °C
Внутреннее давление, p: 0.5 МПа
По ГОСТ 14249–89 расчет на прочность при испытаниях не проводится, если выполнено условие:
| Pисп < 1,35 · Pрасч | [ σ ]20 [ σ ] |
| 1,35 · Pрасч | [ σ ]20 [ σ ] | = |
Допускаемые напряжения для материала ВСТ3 при температуре 20 °C (условия испытаний):
| [σ]= | Re nT | 250/1,1=227,3 МПа |
Коэффициент nT выбирается в зависимости от группы аппарата.
Модуль продольной упругости при температуре 20 °C.
E = 199000 МПа
Коэффициент линейного расширения при температуре 20°C:
α = 1.16e-0,05 C-1
Днища, нагруженные внутренним избыточным давлением (п. 3.3.1).
Расчетная толщина стенки с учетом прибавок:
| s1p+c = | p·R 2·[σ]·φ-0,5·p | +c= |
| = | 0,5·1330 2·227,3·1−0,5·0.5 | + = |
= 1,47 мм
Допускаемое давление:
| [p]= | 2·[σ]·φ·(s1-c) R+0,5·(s1-c) | = |
| = | 2·227,3·1·(10−0) 1330+0,5·(10−0) | = |
Заключение: условие прочности выполнено.
Это касается, в частности, расчета прочности и жесткости узлов врезок штуцеров в цилиндрические обечайки и выпуклые днища с учетом внутреннего давления и внешних нагрузок. В основе методики такого расчета лежат известные зарубежные разработки.
Существенно расширена область применения горизонтальных сосудов и аппаратов на седловых опорах: в отличие от методики, описанной в ГОСТ 26202–84, стало возможным определять изгибающие моменты и силы (как над опорами, так и между ними) сосудов произвольной конструкции, а также опор, расположенных в любом месте цилиндрических обечаек.
Расчет арматурных фланцев производится в соответствии с ASME VIII. При этом помимо давления учитываются внешние нагрузки и изгибающие моменты, а также напряжения, вызванные разницей линейных удлинений фланцев и шпилек (болтов) при температурном воздействии.
Эпилог, в котором описывается дальнейшая судьба программы
Внимательные читатели заметили, что мы не упомянули о теплообменниках — последней составляющей аббревиатуры ПАССАТ. Функцией расчета теплообменников программа пополнится в ближайшей перспективе.
Кроме того, планируется реализовать:
- расчет аппаратов колонного типа в соответствии с ГОСТ Р 51273−99 и ГОСТ Р 51274−99;
- расчет сейсмических и ветровых нагрузок, воздействующих на вертикальные аппараты колонного типа, в том числе — при определении периода собственных колебаний (без ограничений количества элементов);
- расчет рубашек сосудов и аппаратов по ГОСТ 25867–83;
- расчет вертикальных цилиндрических стальных резервуаров для нефти и нефтепродуктов;
- БД стандартных узлов и элементов в полном соответствии с требованиями ОСТ, ГОСТ и АТК.
К уже реализованной функции расчета арматурных фланцевых соединений будет добавлен расчет фланцевых соединений сосудов и аппаратов.
В эти дни проводится открытое бета-тестирование, участниками которого стали более 20 предприятий различных отраслей промышленности.
С 1 июня 2004 года НТП «Трубопровод» осуществляет опрос заинтересованных организаций. Опросный лист размещен по адресу: https://www.truboprovod.ru/cad/demo/PollPassat.doc
По результатам опроса и бета-тестирования мы сможем установить приоритетные направления развития программы и определить, в какие программные системы будет осуществляться импорт/экспорт расчетных моделей, созданных в программе ПАССАТ.
Литература, использованная при создании программы ПАССАТ
- ГОСТ 14249–89. Сосуды и аппараты. Нормы и методы расчета на прочность.
- ГОСТ 25221–82. Сосуды и аппараты. Днища и крышки сферические неотбортованные. Нормы и методы расчета на прочность.
- ГОСТ 26202–84. Сосуды и аппараты. Нормы и методы расчета на прочность обечаек и днищ от воздействия опорных нагрузок.
- ГОСТ 24755–89. Сосуды и аппараты. Нормы и методы расчета на прочность укрепления отверстий.
- РД 26−15−88. Сосуды и аппараты. Нормы и методы расчета на прочность и герметичность фланцевых соединений.
- РД РТМ 26−01−96−77. Плоские круглые крышки и днища с ребрами жесткости.
- РД 26−02−62−98. Расчет на прочность элементов сосудов и аппаратов, работающих в коррозионно-активных сероводородсодержащих средах.
- Нормы расчета на прочность оборудования и трубопроводов атомных энергетических установок. — М., Энергоатомиздат, 1989. — 525 с.
- ГОСТ 27772–88. Прокат для строительных стальных конструкций.
- ГОСТ 25859–83. Сосуды и аппараты стальные. Нормы и методы расчета на прочность при малоцикловых нагрузках.
- ASME VIII, Div 1, 2002, Appendix 2.
- ASME II, 1998, Appendix 2.
- WRC-107 Welding Research Council. Bulletin. — «Local Stresses in Spherical and Cylindrical Shells due to External Loadings», 1979.
- WRC-297 Welding Research Council. Bulletin. «Local Stresses in Cylindrical Shells due to External Loadings on Nozzles, — Supplement to WRC Bulletin № 107», 1987.
- BS-5500: 1976 Specification for Unfired fusion welded pressure vessels. British Standards Institution.
- WRC-368 Welding Research Council. Bulletin. «Stresses in Intersecting Cylinders subjected to Pressure», — 1991. — 32 P.
- Bildy Les M. 2000. «A Proposed Method for Finding Stress and Allowable Pressure in Cylinders with Radial Nozzles», — PVP Vol. 399, ASME, New York, NY. — pp. 77−82.
- Zick L.P. «Stresses in Large Horizontal Cylindrical Pressure Vessels on Two Saddle Supports», — Welding Research Journal Supplement, September, 1951.
Источник