IOP Conference Series: Materials Science and Engineering
PAPER • OPEN ACCESS
Using the Program “MatLab” to Solve the Problems of Mechanics:
Determining the Calculated Values of Internal Forces in a Statically
Indeterminate Continuous Beam from the Actions of Constant Stationary
and Temporary Mobile Loads by the Matrix Method
To cite this article: N K Baishev et al 2018 IOP Conf. Ser.: Mater. Sci. Eng. 463 042023
View the article online for updates and enhancements.
This content was downloaded from IP address 179.61.186.83 on 31/12/2018 at 17:44
International Multi-Conference on Industrial Engineering and Modern technologies IOP Publishing
IOP Conf. Series: Materials Science and Engineering 463 (2018) 042023 doi:10.1088/1757-899X/463/4/042023
Using the Program "MatLab" to Solve the Problems of
Mechanics: Determining the Calculated Values of Internal
Forces in a Statically Indeterminate Continuous Beam from
the Actions of Constant Stationary and Temporary Mobile
Loads by the Matrix Method
N K Baishev1, V N Kononov1, M N Tolstyakova1,*
1
Engineering-technical institute, North-Eastern Federal University named after M K
Ammosov, Kulakovskii str 50, Yakutsk 677000, Russia
*
E-mail:
[email protected] Abstract. This article discusses the using of the MATLAB program to solve problems of
mechanics. Analyzing the calculated values of internal forces in a statically indeterminate
continuous beam from the actions of constant quiescent and temporary moving loads by the
matrix method. The main algorithms for solving problems using MatLab program and
graphical design of tasks are given.
1. Introduction
The calculated values of internal forces in the cross-sections of the elements under the action of a
constant quiescent and temporary moving loads are determined by the formula:
and from these calculated values, their enveloping diagrams are constructed.
When calculating a statically indeterminate continuous beam, these formulas will have the form:
In the analytical calculation of a continuous beam, the values of the support moments from the
action of a constant quiescent load are determined from the solution of the system of equations of
three moments:
The values of the bending moments in any cross-section of each span of the beam's main system
are determined by:
Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution
of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
Published under licence by IOP Publishing Ltd 1
International Multi-Conference on Industrial Engineering and Modern technologies IOP Publishing
IOP Conf. Series: Materials Science and Engineering 463 (2018) 042023 doi:10.1088/1757-899X/463/4/042023
and from these calculated values the diagram of the bending moments from the action of a constant
quiescent load is constructed.
The calculation from the action of a temporary quiescent load, which can have any position with
any discontinuities, is made by the method of focal relations. The focal relations of the beams on each
span is determined by:
– left focal relations
– right focal relations
The support moments in each loaded span of the beam from the action of the temporary moving
load are determined by:
The support moments in the unloaded spans are determined by:
The calculation of continuous beams using this method is quite laborious. Therefore, the solution of
such problems is proposed by the matrix method, using the Maxwell-Mohr’s integral in solving the
Simpson method.
2. Methodology
When solving a statically indeterminate continuous beam, the basic system and the resolving equations
of calculation are assumed to be the same as the method of equations of three moments. The basic
system is obtained by introducing complete hinges in all the support sections of the beam and
replacing them with the support moments .
In this case, the resolving equations will have the form:
The solution of these equations is carried out by the Maxwell-Mohr’s integral
When composing matrices, these integral equations are solved by Simpson's numerical method
Then we obtain the matrix equation
Where: – matrix obtained from vectors of bending moments
2
International Multi-Conference on Industrial Engineering and Modern technologies IOP Publishing
IOP Conf. Series: Materials Science and Engineering 463 (2018) 042023 doi:10.1088/1757-899X/463/4/042023
from actions of constant and temporary loads, which occupies each span of the main beam system
separately;
– matrix composed of vectors of bending moments obtained from the action
of single moments in the main beam system;
– quasi-diagonal stiffness matrix of a beam;
– matrix from the vectors of the final diagrams of the
bending moments from the actions of constant quiescent and temporary moving loads;
matrix from the vectors of the calculated values of the bending moments, on which
the enveloping (calculated) diagram Мх is constructed.
From the values of the support moments for each individual loading, the values of the shearing
forces in any cross-section of the beam of the main system are obtained from
and we obtain an analogous matrix of the calculated values of shearing forces
The calculated values of the support reactions will be from the excision of the support nodes
Example: When the main loads on the beam are represented: a constant load with intensity qconst=4
kN/m, temporary moving loads with intensity qtemp=qeq=6 kN/m and a beam with spans l1=4m,
l2=10m, l3=6m, l4=8m, stiffness EIх=constant.
Figure 1. Initial scheme
Figure 2. Diagrams from a single loading
3
International Multi-Conference on Industrial Engineering and Modern technologies IOP Publishing
IOP Conf. Series: Materials Science and Engineering 463 (2018) 042023 doi:10.1088/1757-899X/463/4/042023
Figure 3. Diagrams of bending moments from the actions of constant and moving loads when each
span is loaded separately of the main beam system.
Create m fail:
l1=4; l2=10; l3=6; l4=8;
p=4; q=6; E=2*10^8; Ix=2000*10^-8;
Lm=[0 0.5 1 1 0.5 0 0 0 0 0 0 0; 0 0 0 0 0.5 1 1 0.5 0 0 0 0;0 0 0 0 0 0 0 0.5 1 1 0.5 0];
Lmp=[0 p*l1^2/8 0 0 p*l2^2/8 0 0 p*l3^2/8 0 0 p*l4^2/8 0; 0 q*l1^2/8 0 0 0 0 0 0 0 0 0 0; 0 0 0 0
q*l2^2/8 0 0 0 0 0 0 0; 0 0 0 0 0 0 0 q*l3^2/8 0 0 0 0; 0 0 0 0 0 0 0 0 0 0 q*l4^2/8 0];
D=[1 0 0; 0 4 0; 0 0 1];
O=[0 0 0; 0 0 0; 0 0 0];
Dg1=(l1/(6*E*Ix))*D;
Dg2=(l2/(6*E*Ix))*D;
Dg3=(l3/(6*E*Ix))*D;
Dg4=(l4/(6*E*Ix))*D;
Dg=[Dg1 O O O; O Dg2 O O; O O Dg3 O; O O O Dg4];
Lmx=Lmp+Lm*(inv(-1*(Lm')*Dg*Lm)*((Lm')*Dg*Lmp))
Lmax=zeros(12,2);
for i=1:12
Z=0; Y=0;
for j=2:5
if (Lmx(i,j)>0) Z=Z+Lmx(i,j);
else Y=Y+Lmx(i,j);
end
end
Lmax(i,1)=Lmax(i,1)+Lmx(i,1)+Z;
Lmax(i,2)=Lmax(i,2)+Lmx(i,1)+Y;
end
Lmax
Get
Lmx =
0 0 0 0 0
-101.9211 10.0602 -20.4417 1.6748 -1.0827
-219.8421 -3.8797 -40.8835 3.3496 -2.1654
4
International Multi-Conference on Industrial Engineering and Modern technologies IOP Publishing
IOP Conf. Series: Materials Science and Engineering 463 (2018) 042023 doi:10.1088/1757-899X/463/4/042023
-219.8421 -3.8797 -40.8835 3.3496 -2.1654
194.6579 -1.3083 36.7951 -3.0147 1.9489
-190.8421 1.2632 -35.5263 -9.3789 6.0632
-190.8421 1.2632 -35.5263 -9.3789 6.0632
-69.9737 0.4962 -13.9568 17.5297 -11.3323
14.8947 -0.2707 7.6128 -9.5617 -28.7278
14.8947 -0.2707 7.6128 -9.5617 -28.7278
39.4474 -0.1353 3.8064 -4.7808 33.6361
0 0 0 0 0
Lmax =
0 0
-90.1861 -123.4455
-216.4925 -266.7707
-216.4925 -266.7707
233.4019 190.3350
-183.5158 -235.7474
-183.5158 -235.7474
-51.9477 -95.2628
22.5075 -23.6654
22.5075 -23.6654
76.8898 34.5312
0 0
Figure 4. Enveloping (calculated) diagram Мх
3. Conclusion
From the example you can see the solution of the technical problem with the help of the MatLab
program are simplified. The practical significance of the research is that the application of the program
"MATLAB" for the solution of problems of mechanics ", with the use of information technologies, is
realized in the practice of professional training of bachelors of a technical university (engineering-
technical institute of the North-Eastern Federal University).
References
[1] Sargsyan A E 2004 Structural mechanics (Moscow: Vusshaya shkola)
[2] Shaposhnikov N N, Kristalinskii R E, Darkov A V 2017 Structural mechanics (Saint-
Petersburg: Lan)
[3] Smirnov V A, Gorodetskii A S 2016 Structural mechanics (Moscow: Urait)
