Форум программистов, компьютерный форум, киберфорум
Matlab
Войти
Регистрация
Восстановить пароль
 
Рейтинг 4.69/13: Рейтинг темы: голосов - 13, средняя оценка - 4.69
0 / 0 / 0
Регистрация: 09.01.2016
Сообщений: 4
1

Построение эпюр

22.12.2016, 00:21. Просмотров 2326. Ответов 1
Метки нет (Все метки)

В матлабе нужно построить эпюры для 3-х участков. все выходит идеально, кроме эпюры Q, причем если не делать 3-й участок, то вышеуказанная эпюра совпадает с рассчитанной. На третьем участке значение Q должно быть равно 0 на всем отрезке. Подскажите, где ошибка, возможно в задании глобальных координат.

Вот код
Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
clear all;
%Построение упругой линии балки при прямом изгибе методом начальных параметров.
%Шарнирно-опертая балка под действием центральной сосредоточенной силы.
 
%Исходные данные:
l = 2.4;                % м, длина первого и второго участка балки
a = 0.8;                % м, длина третьего участка балки
I = 17.5*10^-5;        % м^4, момент инерции сечения
E = 2*10^11;            % Па, модуль Юнга для стали
EI = E*I;               % жесткость при изгибе
q = 2300;               % H/м, максимальная допускаемая нагрузка
 
%Начальные параметры для первого участка:
v_0 = 0;                % м, прогиб               
fi_0 = 0;               % рад,  угол поворота 
M_0 = q*(l^2+a^2);      % Нм, момент
Q_0 = (q*(l^2+a^2))/l;  % Н, сила
q_0 = 2300;             % Н/м, интенсивность распределенной нагрузки
 
%Первый участок
N1 = 20;                    % число точек  
z1 = linspace(0, l, N1);     % координата 
 
for i=1:N1
    
    % координата
    z(i) = z1(i);
    
    %поперечная сила
    Q(i) = Q_0 - q_0*z1(i);  
    
    %изгибающий момент
    M(i) = M_0 - Q_0*z1(i) + q*(z1(i))^2/2; 
    
    %угол поворота сечения
    fi(i)= fi_0 + M_0*z1(i)/EI + Q_0*(z1(i))^2/(2*EI) + q_0*(z1(i))^3/(6*EI);  
    
    %прогиб
    v(i) = v_0 + fi_0*z1(i) + M_0*z1(i)*z1(i)/(2*EI) + Q_0*(z1(i))^3/(6*EI) + q_0*(z1(i))^4/(24*EI);  
    
end;
 
%Начальные параметры для второго участка:
v_0 = v(N1);            % м, прогиб               
fi_0 = fi(N1);          % рад,  угол поворота 
M_0 = M(N1);            % Нм, момент
Q_0 = q*a^2/l;        % Н, сила
q_0 = 2300;                % Н/м, интенсивность распределенной нагрузки
 
%Второй участок
N2 = 20;                     % число точек  
z2 = linspace(l, 2*l, N2);   % координата 
 
for i=1:N2
    
    ii=i + N1;
    
    %глобальная координата
    z(ii) = z2(i);
    
    %поперечная сила
    Q(ii) = Q_0 + q_0*(z2(i));  
    
    %изгибающий момент
    M(ii) = q*l^2/2 -q*l^2*(z2(i))/l  + q_0*(z2(i))^2/2; 
    
    %угол поворота сечения
    fi(ii)= fi_0 + M_0*z2(i)/EI + Q_0*(z(i))^2/(2*EI) + q_0*(z2(i))^3/(6*EI);  
    
    %прогиб
    v(ii) = v_0 + fi_0*z2(i) + M_0*z2(i)*z2(i)/(2*EI) + Q_0*(z2(i))^3/(6*EI) + q_0*(z2(i))^4/(24*EI);  
    
end;
 
%Начальные параметры для третьего участка:
v_0 = v(N2);            % м, прогиб               
fi_0 = fi(N2);          % рад,  угол поворота 
M_0 = q*a^2;            % Нм, момент
Q_0 = (q*(l^2+a^2))/l ;        % Н, сила
q_0 = 0;                % Н/м, интенсивность распределенной нагрузки
 
%Третий участок
N3 = 20;                     % число точек  
z3 = linspace(2*l, 2*l+a, N3);   % координата 
 
 for i=1 : N3
    iii = i + N2;
    %глобальная координата
    z(iii) = z3(i);
    
    %поперечная сила
    Q(iii) = q_0*(z3(i))^2/l;  
    
    %изгибающий момент
    M(iii) = M_0; 
    
    %угол поворота сечения
    fi(iii)= fi_0 + M_0*z3(i)/EI + Q_0*(z3(i))^2/(2*EI) + q_0*(z3(i))^3/(6*EI);  
    
    %прогиб
    v(iii) = v_0 + fi_0*z3(i) + M_0*z3(i)*z(i)/(2*EI) + Q_0*(z3(i))^3/(6*EI) + q_0*(z3(i))^4/(24*EI);  
end;
 
subplot(4,1,1), area(z, Q, 'FaceColor', [0.75 1 0.75]), xlabel('z, m'); ylabel('Q, H'), title([' Поперечная сила Q ']), grid on,
subplot(4,1,2), area(z, M, 'FaceColor', [0.75 1 0.75]); xlabel('z, m'); ylabel('M, Hm'), title([' Изгибающий момент М ']),grid on,
subplot(4,1,3), area(z, fi, 'FaceColor', [0.75 1 0.75]); xlabel('z, m'); ylabel('\phi, rad'), title([' Угол поворота \phi ']),grid on,
subplot(4,1,4), area(z, v, 'FaceColor', [0.75 1 0.75]); xlabel('z, m'); ylabel('v, m'), title([' Прогиб v ']),grid on
0
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
22.12.2016, 00:21
Ответы с готовыми решениями:

Построение эпюр изгибающих моментов и поперечных сил
Добрый день. Пытаюсь построить эпюры изгибающих моментов и поперечных сил. Но на графике выходит...

Построение эпюр
Здравствуйте. Как можно построить данные эпюры, сложно ли это? На что следует обратить внимание?...

Построение эпюр
проверьте , пожалуйста начало задачи двухступенчатый стальной брус , длины степеней которого...

Построение эпюр в пространственных рамах
Вопрос такой: как в Autocad 2020 нарисовать раму, приложить к ней силы и затем построить эпюры...

1
Эксперт по математике/физике
3360 / 1886 / 568
Регистрация: 09.04.2015
Сообщений: 5,296
22.12.2016, 09:07 2
Веселая задача для програмистов, найти ошибку в задаче по сопромату при отсутствии схемы нагружения.
И это при том что программисты о сопромате слышали только страшные былинные расказы собратьев студентов с других факультетов.

При запуске Вашей программы на третьем участке Q=0.
Но в формулах у Вас похоже наврано.
Например второй участок.
Судя по Q там есть распределенная нагрузка, а M на этом участке прямая линия, а должна быть парабола.
0
Миниатюры
Построение эпюр  
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
22.12.2016, 09:07

Заказываю контрольные, курсовые, дипломные и любые другие студенческие работы здесь.

Сопротивление балки, построение эпюр Q и N. проверьте решение
построить эпюру Q и N.( построение эпюры- они для формул Q с +,а расчеты приведены когда в формулы...

Построение эпюр изгибающих моментов статически неопределимой системы
Разработать алгоритм и программу решения задачи построения эпюр изгибающих моментов статически...

Эпюр Монжа
Помогите растолковать эпюр Монжа. Поправьте меня, если я ошибаюсь. Рис.12 можно превратить в...

Как из эпюр напряжения и тока ДЛ определить значение индуктивности и емкости?
Имеются эпюры распределения напряжения в линии, нагруженной на различные элементы: а)...

Построение группы на примере вращений правильного треугольника и построение групп вообще
Вот есть группа, образованная в результате вращений правильного треугольника. Эта группа замкнутая,...

Построение
Здравствуйте вот нужно отредактировать программу т.к как не работает. Заранее спасибо все что не...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
2
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.