0 0 0 0 0 0 - R 0 0 0 0 0 0 0 0 0 -1 0 0 0;
0 0 0 0 0 0 - R 0 0 0 0 0 0 0 0 0 1 0 0 0;
0 0 0 0 0 0 0 - R 0 0 0 1 0 0 0 0 0 -1 0 0;
0 0 0 0 0 0 0 - R 0 0 0 0 0 0 0 0 0 1 0 0;
0 0 0 0 0 0 0 0 - R 0 0 0 0 1 0 0 0 0 -1 0;
0 0 0 0 0 0 0 0 - R 0 0 0 0 0 0 0 0 0 1 0;
0 0 0 0 0 0 0 0 0 - R 0 0 0 0 0 0 0 0 0 -1;
0 0 0 0 0 0 0 0 0 - R 0 0 0 0 0 0 0 0 0 1];
b=zeros(20,1); Nl=-500000
Nu=500000;
% Нижняя граница решений:
x=ones(20,1); Nc=Nl
while Nc>x(i)
lb=[0 0.0 0 0 0 0.0 0 0 0 0.0 Nl Nl Nl Nl Nl Nl Nl Nl Nl Nl];
% Верхняя граница решений:
ub=[0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 Nu Nu Nu Nu Nu Nu Nu Nu Nu Nu];
% Решение задачи линейного программирования:
[x, fval, exitflag] = linprog(c, A,b, Aeq, beq, lb, ub)
% E=210e+09;
for i=11 : 20
if x(i)<0 i
d=sqrt(4*x(i-10)/pi);
J=pi*d^4/64;
Nc=-pi^3*E*J/c(i-10)^2
if Nc>x(i)
lb(i)=Nc
end;
end; %if
end; %i
end;
disp('Конец___________________________________')
В результате решения задачи в системе Matlab получены площади поперечных сечений
см2 (12)
и значения продольных сил
кН. (13)
Отсюда следует, что 6-ой и 10-ый стержни не будут нагружены. Один из них может быть оставлен по конструктивным соображениям (например, 6 для поддержания стержня 2).
Расчёты, приведённые выше, выполнены по условиям прочности (9). Для стержней, оказавшихся сжатыми, этого недостаточно, так как они могут потерять устойчивость. Поэтому им следует придать достаточную жёсткость на изгиб. Это означает, что подбирая сечение таких стержней по площадям (12), необходимо им придать такую геометрическую форму, имеющую достаточный осевой момент инерции. Определим их. С этой целью приравняем сжимающие продольные силы (13) к критическим эйлеровым силам
Отсюда имеем:
Вычисления дали
J = {1,52; 3,04; 4,30; 25,78} см4.
Трёхстержневая система
Пример 2. Заданы система (рис.1), сила P, углы наклона, длины стержней 
. Требуется определить
площади сечений 
, соответствующие наименьшему весу конструкции в условиях упруго-пластической работы материала (рис. 2).
Решение
Текучесть материала начинается при 
, R – предел текучести. Стержень будет течь при значении продольной силы 
, 
– площадь поперечного сечения. Предельное состояние данной конструкции наступит тогда, когда какие-либо два стержня будут течь. Ему соответствует значение 
. Требование задачи состоит в том, чтобы найти такие три площади сечения, что при заданном значении силы Р, два стержня будут течь.
Вырежем узел С и покажем действующие в сечениях продольные силы 
.
Уравнения равновесия данной системы сил имеют вид:
Вес конструкции можно принять пропорциональной предельным пластическим силам 
и длинам стержней 
. Тогда целевая функция примет вид
(1)
где 
.
% tri_sterjnja. m Задача линейного программирования
% Статически неопределимая трёхстежневая система,
% загруженная силой P.
% Целевая функция F=с*x --> minimum (1)
% Система ограничений: Ax<=b, lb<=x<=ub
% где f, x, b, lb, ub - векторы, A - матрица.
% x=(x1, x2, x3)- аргументы, c - коэффициенты в (1)
% Найти такой х, который доставляет
% минимум F при условии Ax<=b, lb<=x<=ub
% F(x)=l1*N10+l2*N20+l3*N20 --> inf
clear;
disp('Начало___________________________________');
P=520450; c=[1 3 2]; st=240e+06;
% Уранение равновесия: Aeq*x=beq
Aeq=[cosd(45) cosd(50) cosd(70);
- sind(45) sind(50) sind(70)];
beq=[P 0];
% Ограничения
lq=[0 0 0]; % нижние значения х
uq=[1000000 1000000 1000000]; % верхние значения х
[x, fval, exitflag] = linprog(c,[],[],Aeq, beq, lq, uq)
A=x/st*10000
disp('Конец___________________________________')
Пример. 
Н, 
, 
м, 
м, 
м, нижние значения 
, верхние значения 
Н. По этим данным составлена матлаб-программа, которая выдала результаты:
Нм, 
Н, 
, 
Н.
По этим значениям сил определены оптимальные значения площадей:
см
, 
, 
см
.
Стержень 2 не нужен.
4.Плита
Прямоугольная жёсткая плита лежит на четырёх столбах, каждый из которых способен выдержать нагрузку N (рис.1). В точке А прикладывается сосредото-ченный груз Р. Требуется найти его максимальное значение, пренебрегая дефор-мациями плиты и столбов.
Решение
В указанном случае плита является абсолютно твёрдым телом, плоскость которого поворачивается в пространстве. Целевая функция будет представляться как одно из уравнений равновесия, а именно равновесие в направлении оси z-ов
, (1)
где 
- опорные реакции, которые удовлетворяют условиям прочности
. (2)
Совокупность сил, приложенных к плите, является параллельной системой, для которой можно составить три уравнения равновесия. Одна из них использована как целевая функция (1). Два других уравнения представляют собой равенства нулю сумм моментов относительно осей x, y и служат ограничениями - равенствами для целевой функции:
(3)
|
Из за большого объема этот материал размещен на нескольких страницах:
1 2 3 4 5 |
