Большинство задач теории оптимального проектирования конструкций, в частности задачи, обсуждавшиеся в предыдущих параграфах, рассматривались в рамках детерминированного подхода, т. е. предполагались полностью известными вид прикладываемых к телу нагрузок, свойства материалов, из которых изготовлена конструкция, граничные условия. Для решения этих задач применимы методы вариационного исчисления и методы оптимизации систем с распределенными параметрами. Принципиально отличными по постановке и методам исследования оказываются задачи проектирования оптимальных конструкций при неполной информации. Обсуждая здесь различные подходы к задачам оптимизации с неполной информацией и их специфику, для определенности будем иметь в виду задачу отыскания форм упругих тел, обладающих минимальным весом и удовлетворяющих заданным ограничениям на прочность и жесткость. Формулировка и решение оптимизационных задач на основе детерминированного подхода приводит к оптимальным формам, которые, как правило, обладают тем свойством, что даже при незначительных изменениях внешних условий (например, при изменении положения точки приложения силы) конструкция данной формы уже не будет удовлетворять прочностным и геометрическим ограничениям. А так как в ряде случаев либо не имеется полной информации относительно прикладываемых нагрузок, либо известно, что на конструкцию последовательно могут действовать различные силы, то наряду с детерминированными постановками представляет интерес рассмотрение более общих задач оптимизации конструкций, в которых оптимизация проводится в расчете на целые классы сил.
Posts Tagged ‘методы’
Проектирование при неполной информации.
Пятница, июля 30, 2010Использование вариационных принципов для исключения дифференциальных связей
Среда, июля 21, 2010Вариационные принципы имеют большое значение в механике деформируемого твердого тела. Во-первых, они позволяют конпактно и в более общей форме сформулировать краевые задачи механики. Уравнения равновесия деформируемой среды и часть краевых условий (естественных) вытекают из вариационных принципов в качестве условий экстремума. Во-вторых, вариационные постановки задач о нагружении и деформировании конструкций позволяют для решения применять эффективные прямые методы вариационного исчисления. Все это справедливо не только для краевых задач механики деформируемого твердого тела, но и для многих других проблем математической физики.
Применительно к задачам оптимального проектирования вариационные принципы позволяют исключить из рассмотрения дифференциальные связи и устраняют необходимость введения сопряженных уравнений. Тем самым понижается порядок общей краевой задачи оптимизации и упрощается вывод условий оптимальности. Кроме того, вариационные принципы и вытекающие из них вариационные неравенства оказываются полезными при аналитических исследованиях оптимизационных задач и обосновании оптимальных решений.
2.4.1. Предположим, что упругое тело занимает область Q, ограниченную поверхностями Г = Ги + Га. На части поверхности Ги тело закреплено, а на остальной части Га к нему приложены внешние воздействия q. Рассматривается задача минимизации податливости
J = -L^ qudTG->minra (2.31)
за счет отыскания формы Га. Для замкнутой постановки задачи требуется сформулировать определяющие уравнения (уравнения равновесия) и указать дополнительные ограничения, а при решении задачи оптимизации ввести в рассмотрение сопряженную переменную, подчинив ее дополнительной системе уравнений. Таким способом в общем случае осуществляется учет дифференциальных связей. Однако вид функционала (2.31) и использование вариационного принципа теории упругости (принцип минимума потенциальной энергии) позволяют переформулировать задачу оптимизации таким образом, что устраняется необходимость введения сопряженной переменной. Покажем это. Согласно вариационному принципу действительное распределение вектора смещений и(х) упругой среды реализует минимум функционала [2.1, 2.8, 2.11, 2.15]
П (и, Га) = \ \ Оцец dQ — ^ uq dTG minu (2.32)
при условии (u)ru = 0. В (2.32) предполагается, что напряжения Gij и деформации б|7- выражены через перемещения при помощи кинематических условий и закона Гука. Для действительного распределения смещений и(х), т. е. для минимали функционала (2.32), имеем /= —П. Учитывая это равенство и вариационный принцип (2.32), можно выразить податливость / через функционал П:
/ = - minu П, (2.33)
следовательно, рассматриваемая задача минимизации податливости может быть сведена к последовательному выполнению операций минимума и максимума для функционала
= minra/ = niinra (— minwn) = — maxra minJI. (2.34)
Внутренний минимум по и в (2.34) вычисляется при дополнительном условии (и)ти = 0, а внешний максимум по Г0 может разыскиваться при некоторых дополнительных условиях, накладываемых на возможные вариации Га, таких, например, как изопери-метрическое условие постоянства объема области Q, занимаемой упругой средой.
Двойственные задачи оптимизации
Воскресенье, июля 4, 2010Преобразования двойственности и изучение двойственных задач приобретают значение в теории оптимального проектирования. Рассмотрение двойственных задач позволяет предложить методы построения оценок величины глобального экстремума и оценить предельные возможности оптимизации. В ряде случаев удается определить проекты, для которых значения функционала качества близки, а иногда и равны величине глобального экстремума. Следуя работе [2.10], дальнейшие рассмотрения проведем применительно к задаче оптимизации тонких пластин, совершающих свободные колебания. Заметим, что в [2.10] методы теории двойственной оптимизации применяются также к задачам оптимизации для трехмерной упругой среды (задаче минимизации массы материала при ограничениях на напряжения и задаче минимизации энергии упругой деформации при ограничениях на массу материала).
Пусть срединная поверхность пластинки занимает область Q с контуром Г в плоскости я, у. Для определения основной частоты со и собственной функции свободных колебаний имеем следующие соотношения:
h dQ = hm mes Q, hm\n
0 <^ ^min "\ hm
По определению двойственной к исходной называется задача [2.12] отыскания /ц, w%, таких, что
J w*) = inf sup J (h, w). (2.67)
wtEV ЛеЯ
Очевидно, что справедливо неравенство
sup inf / (h, w) <; inf sup J (h, w), (2.68)
Лея гиеУ w^V heH
которое может быть использовано для построения оценки сверху величины супремума в задаче (2.66). Обозначим
J0 = sup / (h, Wo), (2.69) лея
где Wo — некоторая произвольная функция из множества V. Из соотношений (2.67)—(2.69) вытекает, что величиной /0 функционал качества оценивается сверху:
sup inf / (h, w) < inf sup / (h, w) < sup / (h, w0) = Jo-
Лея №?У w^v Лея Лея
Нахождение J0 основывается на построении функции h0, такой, что / (h0, wo) = Jq.
В [2.10] доказано, что промежуточный режим hmin << h < hmax в задаче (2.69) невозможен. Укажем следствие, вытекающее из данного утверждения. Введем управляющую функцию
% = (Umax — h) I (femax — hmin). Тогда
J (Ax + В) ф (wo) dQ
Jo = sup / (x, wo), J = ° ? , (2.70)
хем J (C% + D) ф (w0) dQ
-4 = femin — uaxi -6 = hm3iy, С = femin ^max»
D = ^max-
Способы формализации ограничений
Суббота, мая 1, 2010В теории оптимального проектирования на поведение оптимизируемой конструкции накладываются как ограничения, записываемые в виде неравенств, так и ограничения типа равенств. Современные методы оптимального проектирования позволяют в принципе учитывать ограничения обоих типов. Однако наиболее просто реализуем во многих методах учет ограничений типа равенств, и поэтому в большинстве работ по оптимизации конструкций применяется сведение неравенств к ограничениям типа равенств. С этой целью применяются специальные преобразования и введение вспомогательных переменных. Укажем некоторые способы, наиболее часто используемые в теории оптимального проектирования.
Применение операции вычисления абсолютного значения величины позволяет записать ограничения типа неравенств (1.4) в эквивалентном виде:
Ь + \Ъ\=0. (2.9)
Индекс / здесь и ниже в данном параграфе принимает значения от 1 до к.
Система равенств
1 + sign ypj = 0 (2.10)
эквивалентна соотношениям (1.4).
Нетрудно заметить, что если ограничения на поведение конструкции заданы не в каноническом виде (1.4), а в виде двусторонних неравенств 0<^tyj^Cj(cj — заданные константы), то прием, аналогичный (2.9), приводит к равенствам
?/-|Ф/1 = 0, (2.11)
% - cj + | Ь — cj\= 0, а прием, сходный с (2.10),— к равенствам
1 - sign ypj = 0, (2.12)
1 + Sign (l|); — Cj) = 0.
Другая группа приемов сведения ограничений типа неравенств к ограничениям типа равенств основана на введении вспомогательных переменных. Введением вспомогательных переменных % ограничения типа неравенств (1.4) могут быть записаны в виде
% + х! = 0, (2.13)
где %j (х) — подлежащие определению неизвестные вещественные функции (— оо < Xj < °°)-
Для системы ограничений в виде двухсторонних неравенств
Ifcmin < %max (2.14)
следующий способ введения вспомогательных переменных %у позволяет осуществить преобразование к ограничениям типа равенств (—оо <^ %j <; оо):
— *jmin) (%max ~ fy) — Xi = °- (2Л5)
Можно указать и другие способы сведения ограничения типа неравенств к ограничениям типа равенств, использовавшиеся в работах по оптимальному управлению и оптимизации конструкций. Обсуждавшиеся в данном параграфе приемы сведения неравенств к равенствам позволяют применять для исследования развитые методы классического вариационного исчисления.
О замене локальных характеристик интегральными функционалами
Суббота, апреля 24, 2010Задачи оптимизации конструкций можно условно разделить по типу оптимизируемых функционалов и виду ограничений на две группы. К первой группе отнесем задачи оптимизации, для которых критерий качества и ограничения выражаются через интегралы от искомых функций. При этом ограничения имеют вид «изопериметрических» равенств и неравенств. Наиболее часто в работах по оптимальному проектированию встречаются такие интегральные характеристики, как вес, энергия деформации (податливость), сила потери устойчивости, частота собственных колебаний. Отнесем задачи оптимизации ко второй группе, если критерий качества и рассматриваемые ограничения имеют локальный (неинтегральный) характер. К этой группе отнесем смешанные задачи, когда в рассмотрение принимаются как интегральные, так и локальные характеристики конструкции. Типичными локальными функционалами, минимизируемыми при оптимальном проектировании конструкций, являются, например, максимальное смещение в деформируемом теле и максимальное значение интенсивности напряжений.
Большая часть результатов, полученных в теории оптимального проектирования, относится к задачам первой группы. Это объясняется прежде всего тем, что для решения задач с интегральными функционалами существуют известные методы классического вариационного исчисления и нелинейного программирования. Применение этих методов позволило для ряда задач выполнить аналитические и численные исследования и обнаружить интересные закономерности.
Для некоторых типов задач с интегральными функционалами исследования существенно упрощаются за счет того обстоятельства, что уравнения равновесия оказываются «естественными» для рассматриваемых функционалов и допускается исключение дифференциальных связей [2, 2.2, 3.17].
Меньшее число работ посвящено исследованию двумерных задач с локальными функционалами. Главной причиной этого является отсутствие достаточно общих эффективных методов решения для задач второй группы. Типичные трудности решения этих -задач заключаются в следующем. Например, если решается задача минимизации максимального прогиба пластинки и отыскания оптимального распределения ее толщин, то вывод необходимых условий оптимальности и их численная реализация осложнены тем, что заранее неизвестна точка максимального прогиба. Положение же этой точки существенно зависит не только от вида нагрузки, но и от искомого распределения толщин (переменной проектирования). Если же рассматриваются задачи оптимизации с ограничениями типа неравенств, наложенными на локальные характеристики, то аналогичные трудности связаны с определением положения точек или областей, для которых в рассматриваемых ограничениях реализуется знак строгого равенства. Существенные упрощения при решении задач второй группы достигаются в тех случаях, когда задано положение точек, в которых вычисляется значение локальных характеристик конструкции или положение точки экстремума функционала заранее известно, например, из условий симметрии задачи.
Заметим, что, проводя сопоставление интегральных и локаль ных функционалов, мы прежде всего имеем в виду неодномерные задачи оптимального проектирования.
Выбор расчетной схемы в теории оптимального проектирования
Пятница, марта 26, 2010В теории оптимального проектирования изучаются вопросы наилучшего выбора силовой схемы, формы, свойств материалов и условий работы конструкции, исследуются общие закономерности экстремальных решений и развиваются эффективные методы оптимизации. В результате исследований по оптимальному проектированию выясняются предельные возможности улучшения конструкций, оценивается качество традиционных (неоптимальных) сооружений и выявляются наиболее эффективные способы их совершенствования. Теорию оптимального проектирования отличает широкое разнообразие постановок задач [2, 24, 37— 41]. Это объясняется тем, что и уравнения, определяющие нагру-жение и деформирование конструкции, и требования, предъявляемые к ее механическим характеристикам, существенно отличаются при рассмотрении различных типов конструкций (балки, колонны, криволинейные стержни, пластинки, оболочки), реологических свойств (упругость, пластичность, ползучесть), внешних воздействий (поверхностные и объемные силы, статические и динамические нагрузки, „мертвые силы" и силы, зависящие от поведения самой конструкции, тепловые воздействия), видов управляющих переменных (форма конструкции), предположений о степени полноты информации об условиях работы конструкции (задачи с неполной информацией о виде внешних воздействий и способах закрепления конструкции). Точность модели и исходных данных также влияет на постановку задач.
Вопрос о выборе расчетной схемы (модели) является основным как при анализе конструкции, так и при ее оптимизации. Поэтому оптимальное проектирование невозможно без предварительной выработки представлений о существенных и несущественных аспектах поведения конструкции, схематизации условий функционирования и различного рода упрощений, сохраняющих адекватность схемы реальной конструкции. Реальный объект, освобожденный от несущественных особенностей, носит название расчетной схемы. Такое определение расчетной схемы дается в курсах сопротивления материалов [19]. Выбор расчетной схемы, по существу, неединствен.
В некоторых случаях несколько различных схем может быть предложено для одного и того же объекта. В то же время одной расчетной схеме может ставиться в соответствие много реальных объектов.
При оптимальном проектировании конструкций стремятся применять расчетные схемы, позволяющие единственным образом определить как существенные величины напряженно-деформированного состояния, так и искомые переменные проектирования. Однако этого не всегда удается достигнуть из-за отсутствия точной информации о внешних воздействиях, несовершенств изготовления изделия, разброса параметров, характеризующих материал конструкции и других факторов неполноты информации. Для адекватной схематизации в этой ситуации целесообразно смягчение требований к точности описания реального объекта и принятие либо схемы расчета конструкции на наихудший случай, либо схемы стохастического описания конструкции. Это так называемые гарантированный и вероятностный подходы.
Введение в оптимизацию конструкций
Пятница, марта 19, 2010Теория оптимального проектирования получила в последнее время значительное развитие в связи с решением стоящих перед механикой важных задач снижения материалоемкости конструкций и улучшения их механических характеристик. Расширились и сами представления о наилучших в том или ином смысле конструкциях и условиях их функционирования. Были разработаны методы численной оптимизации, позволяющие эффективно оценивать чувствительность основных характеристик конструкций к изменениям параметров проектирования и анализировать способы формирования оптимальных решений. Достигнутые результаты позволили, в частности, широко использовать методы оптимизации при разработке систем автоматизированного проектирования. Однако еще многие проблемы оптимального проектирования не получили решения и по ним в настоящее время ведутся интенсивные исследования.
В книге наряду с изложением основных понятий делается попытка отразить современное состояние теории оптимального проектирования. Книга состоит из двух разделов. Первый раздел служит введением в теорию и методы оптимального проектирования. Здесь излагаются постановки задач оптимизации конструкций и способы их преобразования, необходимые условия оптимальности, аналитические и численные методы оптимизации конструкций с распределенными параметрами, методы оптимизации дискретных систем. Обсуждаются вопросы многоцелевого проектирования конструкций, проектирования при неполной информации, основные понятия многокритериальной оптимизации. Второй раздел книги посвящен применению в оптимальном проектировании критериев прочности, жесткости, устойчивости и веса. Здесь рассматриваются полученные с использованием указанных критериев оптимальные решения для балок, криволинейных стержней, ферм, пластинок и оболочек, массивных тел.
В книге отражены результаты исследований, выполненных в Лаборатории оптимизации конструкций Института проблем механики АН СССР. Основная часть материала книги использовалась в лекциях, которые автор читал студентам Московского физико-технического института.
Тонкость помола и прочность
Среда, января 13, 2010Тонкость помола. В зависимости от степени измельчения различают гипсовые вяжущие грубого, среднего и тонкого помола с максимальным остатком на сите с размером ячеек в свету 0,2 мм не более 23, 14 и 2%, соответственно.Чем тоньше размолот гипс, тем больше его реакционная поверхность и тем быстрее протекает схватывание и дальнейшее твердение, но одновременно увеличивается водопотребность. Особое значение тонкость помола имеет для высокообжиговых вяжущих, для них максимальный остаток на сите №008 — не более 15%.
Прочность. Марку низкообжиговых гипсовых вяжущих опре деляют по прочности при сжатии образцов-балочек размеро\ 40x40x160 мм, сформованных из теста нормальной густоты в возрасте 2 часа после затворения водой. ГОСТом установлено 12 марок от Г-2 до Г-25 (цифра в обозначении марки соответствует прочности при сжатии, выраженной в МПа). Наряду с прочностью при сжатии нормируется также прочность при изгибе, которая существенно ниже.
Прочность гипса может быть повышена путем сушки. При этом растворенный двуводный гипс оседает в порах и упрочняет кристаллический сросток. Прочность гипсовых образцов, высушенных при температурах до 60°С, в 2—2,5 раза выше прочности влажных образцов после 2 часов твердения.
У высокообжиговых вяжущих свои методы оценки прочности. Марку ангидритового вяжущего определяют по прочности при сжатии в возрасте 28 суток в образцах из растворов жесткой консистенции состава 1 : 3 (вяжущее : Вольский песок). Предусмотрен выпуск четырех марок вяжущего — 5, 10, 15 и 20 МПа.
Выпускаемый высокообжиговый гипс имеет марки 10, 15 и 20 по прочности при сжатии (МПа) через 28 суток образцов-кубов из раствора пластичной консистенции 1: 0 (без песка). Кроме того, затвердевший высокопрочный гипс отличается повышенным сопротивлением истиранию.
Наряду с гипсовыми вяжущими общестроительного назначения выпускают гипс для фарфоро-фаянсовой и керамической промышленности. Он имеет сроки схватывания, установленные для нормально твердеющего гипса, а остаток на сите с ячейками размером в свету 0,2 мм (№02) не более 1%. К этим вяжущим предъявляются и особые требования: объемное расширение — не более 0,15%, содержание нерастворимых в НС1 примесей — не более 1%, водопоглощение — не менее 30%.
Гипсовые вяжущие твердеют с некоторым увеличением объема (до 1%), благодаря чему гипсовые отливки хорошо заполняют форму и передают ее очертания. При его высыхании трещин не образуется, что позволяет применять гипс без заполнителей.
Классификация измельченных материалов
Суббота, сентября 19, 2009Классификация — разделение твердых зернистых материалов на фракции ( классы ) по крупности кусков и зерен. Эта операция, с одной стороны, снижает расход энергии на измельчение, а с другой, позволяет получить продукт заданной дисперсности, не содержащий некондиционных фракций. Особенно велика значимость классификации в керамической технологии, где многие изделия получают на основе многофракционных сырьевых смесей со строго определенным соотношением размеров зерен.
Существует два основных вида классификации: механическая (разделение на ситах ) и гидравлическая, основанная на различной скорости осаждения зерен разной крупности в водной или воздушной средах.
Процесс разделения сыпучих материалов на классы по крупности путем просеивания через сита называется грохочением. Грохочение производится на грохотах, рабочий элемент которых — колосниковая решетка, состоящая из отдельных колосников, листовая штампованная решетка (при отверстиях более 3 мм) либо проволочное сито (отверстия менее 3 мм). Материал, подлежащий классификации, перемещается по поверхности решета. Если необходимо получить несколько фракций материала, он последовательно пропускается через набор сит.
При рассеве на каждом сите образуются два продукта: верхний (не прошедший через сито) и нижний (прошедший через сито).
Существует три основных схемы рассева: 1) с последовательным выделением фракций от крупных к мелким; 2) то же, от мелких к крупным; 3) комбинированный. При использовании первой схемы крупность кусков уменьшается и, следовательно, уменьшается износ полотен сит, повышается точность рассева. Однако удаление мелких фракций на последних ступенях рассева значительно увеличивает пыление. Применение второй схемы исключает этот недостаток, но увеличивает износ полотен и понижает точность рассева.
Грохочение достаточно крупных зерен применяется при подготовке фракционированного известняка для обжига в шахтных печах, а также при многостадийном дроблении твердых материалов для отделения кондиционных фракций. Грохочение на проволочных ситах широко используется в стекольной технологии и на керамических заводах.
Сита могут быть вращающимися (барабанными или призматическими), плоскими качающимися и вибрационными. На стекольных заводах часто применяют для одновременной сушки и классификации песка барабанный многогранный грохот (сито-бурат). Такой грохот представляет собой призму, грани которой являются металлическими сетками. Песок подается внутрь вращающегося бурата и просеивается через сетки. Не прошедшие через сетки зерна удаляются с другого конца бурата.
На эффективность рассева влияют влажность материала, угол наклона сит к горизонту, толщина слоя материала, лежащего на сите, форма и расположение отверстий, амплитуда колебаний сит и др. Труднее просеиваются пластинчатые зерна. С повышением влажности эффективность грохочения вначале уменьшается, а затем снова возрастает. Применяют следующие виды грохочения:
• предварительное, при котором из исходной массы выделяется негабаритный материал, либо материал, не требующий дробления;
• контрольное, применяемое для контроля крупности готового продукта и выделения отходов. Зерна, крупнее заданного размера, возвращаются на повторное дробление;
• окончательное — для разделения продукта на товарные фракции.
Эффективность грохочения резко снижается с уменьшением размера частиц. Для тонкодисперсных продуктов используют методы разделения по фракциям, основанные на различии в скорости падения частиц в воздушной или жидкой средах.
Разделение материалов в газовых средах называют воздушной сепарацией, а в жидких — гидравлической классификацией.
Оптимизация в неконсервативных задачах упругой устойчивости
Воскресенье, сентября 13, 2009В предыдущих параграфах данной главы приведены решения задач оптимального проектирования конструкций, для исследования устойчивости которых применимы статические методы. Однако при проектировании существенно неконсервативных систем анализ устойчивости должен основываться на динамических критериях. Применение динамических подходов делает задачи оптимизации более сложными, и к настоящему времени получено решение сравнительно небольшого числа задач [8.14, 8.38, 8.40, 8.73, 8.79,, 8.80, 8.86, 8.89, 8.100].
Учитывая, что исследования в данном направлении находятся в начальной стадии, ограничимся здесь лишь обсуждением некоторых результатов, полученных в работе [8.59].
Изучение устойчивости линейных систем с распределенными параметрами основывается на исследовании уравнения (8.3) для амплитудной функции,, которое с учетом обозначения X = г со записывается в виде
[С + рК + ХВ + Х2А]и = 0. (8.121)
Собственной частоте X — Хяе + iX]m соответствуют комплексные (правая и левая) собственные функции и я v. Оператор К предполагается несамосопряженным. Динамическая потеря устойчивости (флаттер) реализуется, если хотя бы одна из характеристических кривых системы (8.121) в пространстве Хце, Xim, р пересекает ПЛОСКОСТЬ ^Re = 0 При Н6К0-
торых значениях р = рп, XJm = (Х1т)п.
Задача оптимизации заключается в максимизации критического значения параметра нагрузки рц за счет соответствующего выбора переменной проектирования h при заданном значении массы конструкции. От h зависят операторы системы (8.121).
Анализ чувствительности проводится с использованием представлений об обобщенном решении. Уравнение для амплитудной функции записывается в виде скалярного произведения
(у, [С + рК + ХВ + Х2А]и) = 0, (8.122)
причем предполагается дифференцируемость (8.122) по переменной h* Далее вычисляется вариация уравнения (8.122),, обусловленная варьированием переменной проектирования, т. е. заменой h на h + 8hv где 8h — малая вариация функции h. Учитывая, что варьирование осуществляется при критических значениях параметров, будем иметь
Ы № + РпЬК + 1 Ы/i ЬВ + (К1т)Ь 8А] ип) +
(8.123)
+ (vfl,Kufl) 8рп + (vfh [В + 2i (Xlm)fl А] ип) i (8XIm)fl = 0. Вводя обозначения Ане + ibIm = (vfh [8С + pffiK + i (Xlm)fl8В + (XIm)2t 8А]), КЯе + iKim = (vn, Kun), FRe + iFlm = (vfl, [D + 2i(onA] un\ уравнение в вариациях запишем в более компактном виде: Ане + *AIm + (?Re 4- iKlm) брп + (FRe + iFlm) i (8Xlm)n = 0.