Форум программистов, компьютерный форум, киберфорум
Matlab
Войти
Регистрация
Восстановить пароль
Другие темы раздела
Matlab Решение сложной системы нелинейных уравнений https://www.cyberforum.ru/ matlab/ thread1885419.html
Здравствуйте. Необходимо решить систему уравнений. Система из 10 уравнений, 9 неизвестных: F = ; Дано: h=0.0001; f=50;
Matlab Визуализация трехмерной матрицы по сечениям (по страницам)
Доброго вечера всем. Такой вопрос: Есть трехмерная матрица K с размерностями (N,N,3). Как визуализировать эту матрицу на одном графике, чтобы было три сечения? я использовал imagesc(K(:,:,1)) для...
От программы требуется перемножить все элементы массива Matlab
От программы требуется перемножить все элементы массива. Я пишу программу так, чтобы она делила массив на потоки, а затем в каждом потоке находила произведение. Затем собирала произведения с...
Matlab Сколько раз находится одна и таже переменная в строчке? Привет. У меня возникла проблема. В моём Матлаб коде я генерирую себе из файла код ассемблера, который хочу анализировать. Этот код Матлабу представляетсе как cell array с именем proc_asm. У меня... https://www.cyberforum.ru/ matlab/ thread1885032.html
Matlab Задать случайные значение переменной https://www.cyberforum.ru/ matlab/ thread1884955.html
Помогите, пожалуйста psi = e^(i*gamma)*(cos(tetta/2) + e^(i*phi)*sin(tetta/2)) Нужно, чтобы параметр tetta задавался случайным образом в диапазоне 0:pi, а параметр phi - 0:2*pi Как это...
С использованием аналитико-численного метода промоделировать работу терминального комплекса с буферизацией Matlab
Есть такое задание: С использованием аналитико-численного метода промоделировать работу терминального комплекса с буферизацией запросов к однопроцессорной ВС при λ=0,2сек, µ=5...
Matlab Ошибка в программе расчета изотермы сварочного процесса
Здравствуйте. Подскажите пожалуйста, что не так в участке кода с 79 по 88 строки? Массив чисел yma не фиксируется в поле worckspace. И ругается на 88-ую строку. function...
Matlab Анализ котировок Здравствуйте. Подскажите пожалуйста, как получить котировки какого-либо финансового инструмента для анализа MATLAB. Мне нужна тиковая детализация за весь доступный в электронном виде период... https://www.cyberforum.ru/ matlab/ thread1884537.html
Matlab Найти область определения функции одной переменной https://www.cyberforum.ru/ matlab/ thread1884304.html
Нужно найти область определения функции одной переменной. Кто-нибудь занимался таким? Подскажите, что можно использовать или алгоритм действий. Функция вводится таким образом: y=input('y=','s');
Matlab Subs - объясните мне, почему так работает y=input('y=','s'); %ввожу простейшее 2*x syms x; subs(y,x,-1); %в итоге ответ 3, а не -2 Объясните мне недоразвитому, почему так работает. Я подозреваю, что это происходит из-за задания... https://www.cyberforum.ru/ matlab/ thread1883619.html
0 / 0 / 0
Регистрация: 09.01.2016
Сообщений: 4
0

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

22.12.2016, 00:21. Просмотров 2480. Ответов 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
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.