1. Solve the deflection, and stresses in the cantilever beam manually by FEM.

The section
has a width of 40 cm, height of 20. The modulus elasticity of the material is given as
E=2.1x 1011 Pa. Compare the result with that of ANSYS.

Method 1: By manual calculation

Given Required
w=40cm Deflections, V2, V3, Ө2, Ө3
h=20cm Stresses, 𝝈𝒙
11
E=2.1*10 pa

Solution
 The beam has two elements and three nodes shown as follow.

1 2 3
1 2

The stiffness matrix[𝐾]𝑒 for a beam element is given as

𝟏𝟐 𝟔𝑳 −𝟏𝟐 𝟔𝑳
𝐸𝐼 𝟒𝑳𝟐 𝟐𝑳𝟐
[𝐾]𝑒 = 3 [ 𝟔𝑳 −𝟔𝑳
]
𝐿 −𝟏𝟐 −𝟔𝑳 𝟏𝟐 −𝟔𝑳
𝟔𝑳 𝟐𝑳𝟐 −𝟔𝑳 𝟒𝑳𝟐

For element (1) For element (2)

 12 6 −12 6 12 6 −12 6
3 2 3 2 3 2 3
𝐿1 𝐿1 𝐿1 𝐿1 𝐿2 𝐿2 𝐿2 𝐿2 2
6 4 6 2 6 4 6 2
−𝐿 −𝐿
𝐿1 2 𝐿1 2
𝐿1 𝐿2 2 𝐿2 2
𝐿2
[𝐾]1 =EI/L^3 −12 −6 12
1
−6
[𝐾]2 = 𝐸𝐼 −12 −6 12
2
−6
𝐿1 3 𝐿1 2 𝐿1 3 𝐿1 2 𝐿2 3 𝐿2 2 𝐿2 3 𝐿2 2
6 2 6 4 6 2 6 4
[ 𝐿1 2 −𝐿 2 ] [ −𝐿 ]
𝐿1 1 𝐿1 𝐿2 2 𝐿2 2
2
𝐿2

 When the two element matrix are assembled and we get the global matrix
12 6 −12 6
3 2 3 𝟎 𝟎
𝐿1 𝐿1 𝐿1 𝐿2 2
6 4 6 2
2 − 𝟎 𝟎
𝐿1 𝐿1 𝐿1 2 𝐿1
−12 −6 12 12 −6 6 −12 6
3 2 3 + 3 2 + 2 3
𝐿1 𝐿1 𝐿1 𝐿2 𝐿1 𝐿2 𝐿2 𝐿2 2
[𝑲](𝑮) = 𝑬𝑰 6 2 −6 6 4 4 6 2 …… (1)
2 + + −
𝐿1 2 𝐿1 𝐿1 𝐿2 2 𝐿1 𝐿2 𝐿2 2 𝐿2
−12 −6 12 −6
𝟎 𝟎
𝐿2 3 𝐿2 2 𝐿2 3 𝐿2 2
6 2 6 4
[ 𝟎 𝟎
𝐿2 2 𝐿2
−
𝐿2 2 𝐿2 ]
 Then the force vector is defined as follows.
{𝑭}=[K]G {𝒅}
12 6 −12 6
𝑭𝟏𝒚 3 2
𝐿1 3
𝟎 𝟎
𝐿1 𝐿1 𝐿2 2
6 4 6 2 𝑽𝟏
𝑴𝟏 2 𝐿1
− 2 𝐿1
𝟎 𝟎
𝐿1 𝐿1
−12 −6 12 12 −6
+
6 −12 6 𝜽𝟏
3 2 3 + 3 2 2 3 2
𝑭𝟐𝒚 𝐿1 𝐿1 𝐿1 𝐿2 𝐿1 𝐿2 𝐿2 𝐿2 𝑽𝟐
= 𝑬𝑰 6 2 −6 6 4 4 6 2 ……………..(2)
+ +𝐿 − 𝜽𝟐
𝑴𝟐 𝐿1 2 𝐿1 𝐿1 2 𝐿2 2 𝐿1 2 𝐿2 2 𝐿2
−12 −6 12 −6 𝑽𝟑
𝑭𝟑𝒚 𝟎 𝟎
𝐿2 3 𝐿2 2 𝐿2 3 𝐿2 2 {𝜽𝟑 }
6 2 6 4
{ 𝑴𝟑 } 𝟎 𝟎 −
[ 𝐿2 2 𝐿2 𝐿2 2 𝐿2 ]

Boundary conditions
 At node-1 V1=Ө1=0, and at node-2 and 3 F2y = -P, M2=0, M3=0 F3y=0
 Substituting these values in the above equation (1) and we will have:
 12 6 −12 6
𝐿1 3 𝐿1 2 𝐿1 3
𝟎 𝟎
𝐿2 2
6 4 6 2
𝑭𝟏𝒚 − 𝟎 𝟎
𝐿1 2 𝐿1 𝐿1 2 𝐿1 𝟎
𝑴𝟏 −12 −6 12 12 −6 6 −12 6 𝟎
+ +
−𝑷 𝐿1 3 𝐿1 2 𝐿1 3 𝐿2 3 𝐿1 2 𝐿2 2 𝐿2 3 𝐿2 2 𝑽𝟐
𝟎 = 𝑬𝑰 6 2 −6 6 4 4 6 2 𝜽𝟐 …………… (3)
+ +𝐿 −
𝟎 𝐿1 2 𝐿1 𝐿1 2 𝐿2 2 𝐿1 2 𝐿2 2 𝐿2 𝑽𝟑
−12 −6 12 −6 {𝜽𝟑 }
{ 𝟎 } 𝟎 𝟎
𝐿2 3 𝐿2 2 𝐿2 3 𝐿2 2
6 2 6 4
𝟎 𝟎 −
[ 𝐿2 2 𝐿2 𝐿2 2 𝐿2 ]

 Equation (3) will be reduced to

12 12 −6 6 −12 6
3 +𝐿 3 2 +𝐿 2 3
𝐿1 2 𝐿1 2 𝐿2 𝐿2 2
−𝑷 −6 6 4 4 6 2 𝑽𝟐
2 +𝐿 2 +𝐿 −𝐿 2 𝜽
{ 𝟎 } = 𝐸𝐼
𝐿1 𝐿1 𝐿2
2 2 2
{ 𝟐 })
𝟎 −12 −6 12 −6 𝑽𝟑
𝟎 𝐿2 3 𝐿2 2 𝐿2 3 𝐿2 2 𝜽𝟑
6 2 6 4
[ −𝐿 ]
𝐿2 2 𝐿2 2
2
𝐿2

 For P=10,000N,I=2.66*10-4m4,E=2.1*1011 Pa, L1=3m and L2=1.5m, then the above

equation can be written as:

−𝟏𝟎, 𝟎𝟎𝟎 4 2 −3.5556 2.6667 𝑽𝟐

𝟎 2 4 −2.6667 1.3333 𝜽
{ } = 5.58 ∗e7[ ] { 𝟐 } ……….(4)
𝟎 −3.5556 −2.6667 3.5556 −2.6667 𝑽 𝟑
𝟎 2.6667 1.3333 −2.6667 2.6667 𝜽𝟑
 Then from equation (4)
4V2+2𝜃2 -3.5556V3+2.6667𝜃3 = −1.78 ∗ 𝑒 − 4………………………... (a)
2V2+4𝜃2 -2.6667V3+1.3333𝜃3 = 0……………………………….............(b)
-3.5556V2-2.6667𝜃2 +3.5556V3-2.6667𝜃3 = 0…………………………...(c)
2.6667V2+1.3333𝜃2 -2.6667V3+2.6667𝜃3 =0 …………………………… (d)
 Thus solving the above equations and we get the following solutions
V2=-0.0016m V3=-0.0028m
𝜽𝟐 = −0.0008 𝒓𝒂𝒅 𝜽𝟑 =-0.0008rad

 To find the reaction force and moment at point -1, we can use the static equilibrium
equations.
∑ 𝐹𝑦=0 , F1y=P
∑ 𝑀=0 , M1=3m*10*103N=30000N.m
 Therefore the stress at point -1 is
𝑀𝑦
𝜎𝑥 = Where M=30000N.m y=0.2/2=0.1m and I=2.6667*10-4m4
𝐼
𝟑𝟎𝟎𝟎𝟎∗𝟎.𝟏
𝝈𝒙 = 𝟐.𝟔𝟔𝟔𝟕∗𝟏𝟎−𝟒 = 11.248Mpa

 Therefore the general solution of this beam is summarized as:

V1=0 𝜽𝟏 = 𝟎
V2= -0.0016m V3= -0.0028m
𝜽𝟐 = −𝟎. 𝟎𝟎𝟎𝟖𝒓𝒂𝒅 𝜽𝟑 = −𝟎. 𝟎𝟎𝟎𝟖𝒓𝒂𝒅
𝝈𝒙 = 11.248Mpa

Method 2: By using ANSYS software [Mechanical APDL]

 Deformed and undeformed shape

 Displacement vector sum

  Rotation vector sum

Nodal solutions
a) Displacement vector sum
PRINT U NODAL SOLUTION PER NODE

***** POST1 NODAL DEGREE OF FREEDOM LISTING *****

LOAD STEP= 1 SUBSTEP= 1

TIME= 1.0000 LOAD CASE= 0

THE FOLLOWING DEGREE OF FREEDOM RESULTS ARE IN THE GLOBAL COORDINATE

SYSTEM

NODE UX UY UZ USUM
1 0.0000 0.0000 0.0000 0.0000
2 0.0000 -0.32687E-03 0.0000 0.32687E-03
3 0.0000 -0.62821E-03 0.0000 0.62821E-03

MAXIMUM ABSOLUTE VALUES

NODE 0 3 0 3
VALUE 0.0000 -0.62821E-03 0.0000 0.62821E-03

b) Rotational deflection
PRINT ROT NODAL SOLUTION PER NODE

***** POST1 NODAL DEGREE OF FREEDOM LISTING *****

LOAD STEP= 1 SUBSTEP= 1

TIME= 1.0000 LOAD CASE= 0

THE FOLLOWING DEGREE OF FREEDOM RESULTS ARE IN THE GLOBAL COORDINATE

SYSTEM

NODE ROTX ROTY ROTZ RSUM

1 0.0000 0.0000 0.0000 0.0000
2 0.0000 0.0000 -0.20089E-03 0.20089E-03
3 0.0000 0.0000 -0.20089E-03 0.20089E-03

MAXIMUM ABSOLUTE VALUES

NODE 0 0 3 3
VALUE 0.0000 0.0000 -0.20089E-03 0.20089E-03
c) All structural forces
PRINT F SUMMED NODAL LOADS

***** POST1 SUMMED TOTAL NODAL LOADS LISTING *****

LOAD STEP= 1 SUBSTEP= 1

TIME= 1.0000 LOAD CASE= 0

THE FOLLOWING X,Y,Z SOLUTIONS ARE IN THE GLOBAL COORDINATE SYSTEM

NODE FX FY FZ
1 -10000.
2 10000.

TOTAL VALUES
VALUE 0.0000 0.50932E-10 0.0000

d) All structural moment

PRINT M SUMMED NODAL LOADS

***** POST1 SUMMED TOTAL NODAL LOADS LISTING *****

LOAD STEP= 1 SUBSTEP= 1

TIME= 1.0000 LOAD CASE= 0

THE FOLLOWING X,Y,Z SOLUTIONS ARE IN THE GLOBAL COORDINATE SYSTEM

NODE MX MY MZ
1 -0.61126E-12-0.19058E-12 -30000.
2 0.61126E-12 0.19058E-12

TOTAL VALUES
VALUE 0.31302E-26 0.53011E-27 -30000.
Elemental solutions
a) 1st principal stress
PRINT S PRIN ELEMENT SOLUTION PER ELEMENT

STRESSES AT BEAM SECTION NODAL POINTS

ELEMENT = 1 SECTION ID = 1

ELEMENT NODE = 1

SEC NODE S1 S2 S3 SINT SEQV

1 0.0000 0.0000 -0.28125E+07 0.28125E+07 0.28125E+07
3 0.0000 0.0000 0.0000 0.0000 0.0000
13 0.14552E-10 0.0000 0.0000 0.14552E-10 0.14552E-10
11 0.0000 0.0000 -0.28125E+07 0.28125E+07 0.28125E+07
5 0.28125E+07 0.0000 0.0000 0.28125E+07 0.28125E+07
15 0.28125E+07 0.0000 0.0000 0.28125E+07 0.28125E+07
23 0.0000 0.0000 -0.23283E-09 0.23283E-09 0.23283E-09
21 0.0000 0.0000 -0.28125E+07 0.28125E+07 0.28125E+07
25 0.28125E+07 0.0000 0.0000 0.28125E+07 0.28125E+07

Max= 0.28125E+07 0.0000 0.0000 0.28125E+07 0.28125E+07

Min= 0.0000 0.0000 -0.28125E+07 0.0000 0.0000

ELEMENT NODE = 2

SEC NODE S1 S2 S3 SINT SEQV

1 0.0000 0.0000 -0.28125E+07 0.28125E+07 0.28125E+07
3 0.0000 0.0000 0.0000 0.0000 0.0000
13 0.14552E-10 0.0000 0.0000 0.14552E-10 0.14552E-10
11 0.0000 0.0000 -0.28125E+07 0.28125E+07 0.28125E+07
5 0.28125E+07 0.0000 0.0000 0.28125E+07 0.28125E+07
15 0.28125E+07 0.0000 0.0000 0.28125E+07 0.28125E+07
23 0.0000 0.0000 -0.23283E-09 0.23283E-09 0.23283E-09
21 0.0000 0.0000 -0.28125E+07 0.28125E+07 0.28125E+07
25 0.28125E+07 0.0000 0.0000 0.28125E+07 0.28125E+07

Max= 0.28125E+07 0.0000 0.0000 0.28125E+07 0.28125E+07

Min= 0.0000 0.0000 -0.28125E+07 0.0000 0.0000

STRESSES AT BEAM SECTION NODAL POINTS

ELEMENT = 2 SECTION ID = 1

ELEMENT NODE = 2

SEC NODE S1 S2 S3 SINT SEQV

1 0.79689E-08 0.0000 0.0000 0.79689E-08 0.79689E-08
3 0.0000 0.0000 0.0000 0.0000 0.0000
13 0.0000 0.0000 0.0000 0.0000 0.0000
11 0.79689E-08 0.0000 0.0000 0.79689E-08 0.79689E-08
5 0.0000 0.0000 -0.79689E-08 0.79689E-08 0.79689E-08
15 0.0000 0.0000 -0.79689E-08 0.79689E-08 0.79689E-08
 23 0.0000 0.0000 0.0000 0.0000 0.0000
21 0.79689E-08 0.0000 0.0000 0.79689E-08 0.79689E-08
25 0.0000 0.0000 -0.79689E-08 0.79689E-08 0.79689E-08

Max= 0.79689E-08 0.0000 0.0000 0.79689E-08 0.79689E-08

Min= 0.0000 0.0000 -0.79689E-08 0.0000 0.0000
ELEMENT NODE = 3

SEC NODE S1 S2 S3 SINT SEQV

1 0.79689E-08 0.0000 0.0000 0.79689E-08 0.79689E-08
3 0.0000 0.0000 0.0000 0.0000 0.0000
13 0.0000 0.0000 0.0000 0.0000 0.0000
11 0.79689E-08 0.0000 0.0000 0.79689E-08 0.79689E-08
5 0.0000 0.0000 -0.79689E-08 0.79689E-08 0.79689E-08
15 0.0000 0.0000 -0.79689E-08 0.79689E-08 0.79689E-08
23 0.0000 0.0000 0.0000 0.0000 0.0000
21 0.79689E-08 0.0000 0.0000 0.79689E-08 0.79689E-08
25 0.0000 0.0000 -0.79689E-08 0.79689E-08 0.79689E-08

Max= 0.79689E-08 0.0000 0.0000 0.79689E-08 0.79689E-08

Min= 0.0000 0.0000 -0.79689E-08 0.0000 0.0000

Conclusion [comparison the two solutions]

 The ANSYS solutions are lower than the solutions obtained by manual hand calculation.
 Thus by manual hand calculation there may be some errors like human error and
numbering error. Therefore the answer obtained by hand calculation is somewhat
different from ANSYS solutions.
2. Solve the given frame structure using MATLAB program and compare the
results with that of Ansys.
F1 F2 F2 F
F33

F1

Where:
A = 0.040 m2 I= 5 x10-4 m2
Material is steel with E=2.1x 1011 Pa.
𝐹1 = 35𝐾𝑁, 𝐹2 = 50𝐾𝑁, 𝐹3 = 75𝐾𝑁 𝑎𝑛𝑑 𝐿 = 2𝑚 is given for our group
Method 1: By using MATLAB Software
%Q#1 (For Group 1 members)
%_________________________________
% Problem description
% Find the deflection of a frame structure which is made of three beams.
% The length of beams: 10m, 10m and 6+6+4+5+3=24m(if L=4).
% The cross-section of beams equals 0.040m^2
% The elastic modulus is 2.1x10^11 pa.
% The frame is subjected to a concenterated load of
% 40, 60 & 60kN at the longest beam and the two short members are fixed.
% SOLUTION:
%_______________________________
% Variable descriptions
% x and y = global x and y coordiates of each node
% kk = system stiffness matrix
% ff = system force vector
% index = a vector containing system dofs associated with each element
% bcdof = a vector containing dofs associated with boundary conditions
% bcval = a vector containing boundary condition values associated
% with the dofs in 'bcdof'
%_______________________________________________________________________
clear all,clc
nel=7; % number of elements
nnel=2; % number of nodes per element
ndof=3; % number of dofs per node
nnode=(nnel-1)*nel+1; % total number of nodes in system
sdof=nnode*ndof; % total system dofs
x(1)=0; y(1)=0; % x, y coord. values of node 1 in terms of the global axis
n1=[x(1);y(1)];
x(2)=0; y(2)=8; % x, y coord. values of node 2 in terms of the global axis
n2=[x(2);y(2)];
x(3)=6; y(3)=8; % x, y coord. values of node 3 in terms of the global axis
n3=[x(3);y(3)];
x(4)=10; y(4)=8; % x, y coord. values of node 4 in terms of the global axis
n4=[x(4); y(4)];
x(5)=13; y(5)=8; % x, y coord. values of node 5 in terms of the global axis
n5=[x(5); y(5)];
x(6)=18; y(6)=8; % x, y coord. values of node 6 in terms of the global axis
n6=[x(6); y(6)];
x(7)=24; y(7)=8; % x, y coord. values of node 7 in terms of the global axis
n7=[x(7); y(7)];
x(8)=24; y(8)=0; % x, y coord. values of node 8 in terms of the global axis
n8=[x(8); y(8)];
el=2.1*(10^11); % elastic modulus
area=0.04; % cross-sectional area
xi=5*(10^-4); % moment of inertia of cross-section
bcdof(1)=1; % transverse deflection at node 1 is constrained
bcval(1)=0; % whose described value is 0
bcdof(2)=2; % axial displacement at node 1 is constrained
bcval(2)=0; % whose described value is 0
bcdof(3)=3; % slope at node 1 is constrained
bcval(3)=0; % whose described value is 0
bcdof(4)=22; % transverse deflection at node 8 is constrained
bcval(22)=0; % whose described value is 0
bcdof(5)=23; % axial displacement at node 8 is constrained
bcval(23)=0; % whose described value is 0
bcdof(6)=24; % slope at node 8 is constrained
bcval(24)=0; % whose described value is 0
ff=zeros(sdof,1); % initialization of system force vector
kk=zeros(sdof,sdof); % initialization of system matrix
index=zeros(nel*ndof,1); % initialization of index vector
ff(10)=-40000; % load applied at node 4 in the negative y direction
ff(13)=-60000; % load applied at node 5 in the negative y direction
ff(16)=-60000; % load applied at node 6 in the negative y direction
nodes=[1 3;2 3;3 4;4 5;5 6;6 7;6 8;];
for iel=1:nel % loop for the total number of elements
nd(1)=nodes(iel,1); % starting node number for element 'iel'
nd(2)=nodes(iel,2); % ending node number for element 'iel'
index=feeldof(nd,nnel,ndof); % extract system dofs for the element
x1=x(nd(1)); y1=y(nd(1)); % x and y coordinate values of 'node1'
x2=x(nd(2)); y2=y(nd(2)); % x and y coordinate values of 'node2'
leng=sqrt((x2-x1)^2+(y2-y1)^2); % length of element 'iel'
if (x2-x1)==0 % compute the angle between the local and global axes
if y2>y1
beta=pi/2;
else
beta=-pi/2;
end
else
beta=atan((y2-y1)/(x2-x1));
end
k=frameel(el,xi,leng,area,beta); % compute element stiffness matrix
kk=feasmbl1(kk,k,index); % assemble each element matrix into system matrix
end
[kk,ff]=feaplyc2(kk,ff,bcdof,bcval); % apply the boundary conditions
fsol=kk\ff; % solve the matrix equation and print
% nodal connectivity
% solution display loop:
% in the result desplayed we have four columns
% first column is for node numbers
% second column is for deflection u (x - direction)
% third column is for deflection v (y - direction)
% fourth column is for rotation theta (z - axis rotation)
% so, each node will have its own coresponding deflection in both x and y
% axis in meter and rotation about z-axis in radians.
store=zeros(nnode,4); % intiation of result desplay matirix
for i=1:nnode % loop that starts from 1 to total number of node
store(i,1)=i;
store(i,2)=fsol(3*i-2);
store(i,3)=fsol(3*i-1);
store(i,4)=fsol(3*i);
end
num=1:1:sdof;
disp ('deformation results are:')
store=[num' fsol]
function [index]=feeldof(nd,nnel,ndof)
%----------------------------------------------------------
% Purpose:
% Compute system dofs associated with each element
% Synopsis:
% [index]=feeldof(nd,nnel,ndof)
% Variable Description:
% index - system dof vector associated with element "iel"
% iel - element number whose system dofs are to be determined
% nnel - number of nodes per element
% ndof - number of dofs per node
%-----------------------------------------------------------
edof = nnel*ndof;
k=0;
for i=1:nnel
start = (nd(i)-1)*ndof;
for j=1:ndof
k=k+1;
index(k)=start+j;
end
end
function [k]=frameel(el,xi,leng,area,beta)
%--------------------------------------------------------------
% Purpose:
% Stiffness and mass matrices for the 2-d frame element
% nodal dof {u_1 v_1 theta_1 u_2 v_2 theta_2}
% Synopsis:
% [k,m]=feframe2(el,xi,leng,area,rho,beta,ipt)
% Variable Description:
% k - element stiffness matrix (size of 6x6)
% m - element mass matrix (size of 6x6)
% el - elastic modulus
% xi - second moment of inertia of cross-section
% leng - element length
% area - area of beam cross-section
% beta - angle between the local and global axes ipt = 1: consistent mass
matrix
% is positive if the local axis is in the ccw direction from
% the global axis
%---------------------------------------------------------------------
% stiffness matrix at the local axis
a=el*area/leng;
c=el*xi/(leng^3);
kl=[a 0 0 -a 0 0;...
0 12*c 6*leng*c 0 -12* c 6*leng*c;...
0 6*leng*c 4*leng^2*c 0 -6*leng*c 2*leng^2*c; ...
-a 0 0 a 0 0;...
0 -12*c -6*leng*c 0 12*c -6*leng*c;...
0 6*leng*c 2*leng^2*c 0 -6*leng*c 4*leng^2*c];
% coordinate transformation matrix
r=[ cos(beta) sin(beta) 0 0 0 0;...
-sin(beta) cos(beta) 0 0 0 0;...
0 0 1 0 0 0;...
0 0 0 cos(beta) sin(beta) 0;...
0 0 0 -sin(beta) cos(beta) 0;...
0 0 0 0 0 1];
% stiffness matrix at the global axis
k=r'*kl*r;
function [kk]=feasmbl1(kk,k,index)
%----------------------------------------------------------
% Purpose:
% Assembly of element matrices into the system matrix
% Synopsis:
% [kk]=feasmbl1(kk,k,index)
% Variable Description:
% kk - system matrix
% k - element matri
% index - d.o.f. vector associated with an element
%-----------------------------------------------------------
edof = length(index);
for i=1:edof
ii=index(i);
for j=1:edof
jj=index(j);
kk(ii,jj)=kk(ii,jj)+k(i,j);
end
end
function [kk,ff]=feaplyc2(kk,ff,bcdof,bcval)
%----------------------------------------------------------
% Purpose:
% Apply constraints to matrix equation [kk]{x}={ff}
% Synopsis:
% [kk,ff]=feaplybc(kk,ff,bcdof,bcval)
% Variable Description:
% kk - system matrix before applying constraints
% ff - system vector before applying constraints
% bcdof - a vector containging constrained d.o.f
% bcval - a vector containing contained value
% For example, there are constraints at d.o.f=2 and 10
% and their constrained values are 0.0 and 2.5,
% respectively. Then, bcdof(1)=2 and bcdof(2)=10; and
% bcval(1)=1.0 and bcval(2)=2.5.
%-----------------------------------------------------------
n=length(bcdof);
sdof=size(kk);
for i=1:n
c=bcdof(i);
for j=1:sdof
kk(c,j)=0;
end
kk(c,c)=1;
ff(c)=bcval(i);
disp=zeros(sdof,1);
stress=zeros(nel,3);
eldisp=zeros(ndof,1);
for i=1:(nnel*ndof)
eldisp(i)=eldisp(index(i));
end
elforce=k*eldisp;
stress(iel)=sqrt(elforce(1)^2+elforce(2)^2)/area;
num=1:1:sdof;
eldisp1=[num;disp]
numm=1:1:nel;
stresses=[numm stress]
store=[num eldisp]
end
end
end
end
end

RESULTS DISPLAYED ON THE COMMAND WINDOW

deformation results are:
store =
1.0000 0.0000
2.0000 -0.0000
3.0000 0.0000
4.0000 -0.0289
5.0000 0.0192
6.0000 0.0004
7.0000 -0.0289
8.0000 0.0216
9.0000 0.0004
10.0000 -0.0290
11.0000 0.0107
12.0000 -0.0049
13.0000 -0.0290
14.0000 -0.0056
15.0000 -0.0054
16.0000 -0.0290
17.0000 -0.0216
18.0000 0.0004
19.0000 -0.0290
20.0000 -0.0192
21.0000 0.0004
22.0000 0
23.0000 0
24.0000 0
>>
 Method 2: By using ANSYS software [mechanical APDL]
Naming of elements and key points

 Deformed and undeformed shape

 Displacement vector sum

 Rotation vector sum

Nodal solutions
a) Displacement vector sum

PRINT U NODAL SOLUTION PER NODE

***** POST1 NODAL DEGREE OF FREEDOM LISTING *****

LOAD STEP= 1 SUBSTEP= 1

TIME= 1.0000 LOAD CASE= 0

THE FOLLOWING DEGREE OF FREEDOM RESULTS ARE IN THE GLOBAL COORDINATE

SYSTEM

NODE UX UY UZ USUM
1 0.0000 0.0000 0.0000 0.0000
2 -0.32451E-01 0.24193E-01 0.0000 0.40477E-01
3 -0.32451E-01 0.41619E-01 0.0000 0.52775E-01
4 -0.32470E-01-0.64584E-03 0.0000 0.32477E-01
5 -0.32498E-01-0.42023E-01 0.0000 0.53123E-01
6 -0.32545E-01-0.24595E-01 0.0000 0.40794E-01
7 -0.32545E-01 0.36853E-01 0.0000 0.49167E-01
8 0.0000 0.0000 0.0000 0.0000

MAXIMUM ABSOLUTE VALUES

NODE 6 5 0 5
VALUE -0.32545E-01-0.42023E-01 0.0000 0.53123E-01

b) Rotation vector sum

PRINT ROT NODAL SOLUTION PER NODE

***** POST1 NODAL DEGREE OF FREEDOM LISTING *****

LOAD STEP= 1 SUBSTEP= 1

TIME= 1.0000 LOAD CASE= 0

THE FOLLOWING DEGREE OF FREEDOM RESULTS ARE IN THE GLOBAL COORDINATE

SYSTEM

NODE ROTX ROTY ROTZ RSUM

1 0.0000 0.0000 0.0000 0.0000
2 0.0000 0.0000 -0.29042E-02 0.29042E-02
3 0.0000 0.0000 -0.29042E-02 0.29042E-02
4 0.0000 0.0000 -0.20807E-01 0.20807E-01
5 0.0000 0.0000 -0.56800E-02 0.56800E-02
6 0.0000 0.0000 0.10241E-01 0.10241E-01
7 0.0000 0.0000 0.10241E-01 0.10241E-01
8 0.0000 0.0000 0.0000 0.0000
 MAXIMUM ABSOLUTE VALUES
NODE 0 0 4 4
VALUE 0.0000 0.0000 -0.20807E-01 0.20807E-01

Elemental solutions
i. Moments and loads

Structural force at reactions

PRINT F REACTION SOLUTIONS PER NODE

***** POST1 TOTAL REACTION SOLUTION LISTING *****

LOAD STEP= 1 SUBSTEP= 1

TIME= 1.0000 LOAD CASE= 0

THE FOLLOWING X,Y,Z SOLUTIONS ARE IN THE GLOBAL COORDINATE SYSTEM

NODE FX FY FZ
1 78945. 62783. 0.0000
8 -78945. 97217. 0.0000

TOTAL VALUES
VALUE 0.24884E-08 0.16000E+06 0.0000

Nodal loads
PRINT F SUMMED NODAL LOADS

***** POST1 SUMMED TOTAL NODAL LOADS LISTING *****

LOAD STEP= 1 SUBSTEP= 1

TIME= 1.0000 LOAD CASE= 0

THE FOLLOWING X,Y,Z SOLUTIONS ARE IN THE GLOBAL COORDINATE SYSTEM

NODE FX FY FZ
1 -78945. -62783.
4 35000.
5 50000.
6 75000.
8 78945. -97217.

TOTAL VALUES
VALUE -0.24884E-08-0.12224E-08 0.0000
 Structural moment solutions at Reaction
PRINT M REACTION SOLUTIONS PER NODE

***** POST1 TOTAL REACTION SOLUTION LISTING *****

LOAD STEP= 1 SUBSTEP= 1

TIME= 1.0000 LOAD CASE= 0

THE FOLLOWING X,Y,Z SOLUTIONS ARE IN THE GLOBAL COORDINATE SYSTEM

NODE MX MY MZ
1 -0.42579E-12-0.73991E-12-0.12526E+06
8 0.19520E-12 0.34697E-12 16484.

TOTAL VALUES
VALUE -0.23058E-12-0.39294E-12-0.10878E+06

Nodal moments
PRINT M SUMMED NODAL LOADS

***** POST1 SUMMED TOTAL NODAL LOADS LISTING *****

LOAD STEP= 1 SUBSTEP= 1

TIME= 1.0000 LOAD CASE= 0

THE FOLLOWING X,Y,Z SOLUTIONS ARE IN THE GLOBAL COORDINATE SYSTEM

NODE MX MY MZ
1 0.42579E-12 0.73991E-12 0.12526E+06
2 -0.25132E-11 0.24445E-11
4 0.11637E-11-0.49783E-11
5 0.16624E-11 0.66099E-12
6 -0.54350E-12 0.14798E-11
8 -0.19520E-12-0.34697E-12 -16484.

TOTAL VALUES
VALUE -0.88352E-27 0.20195E-27 0.10878E+06

ii. Strain intensity

PRINT EPEL ELEMENT SOLUTION PER ELEMENT

ELASTIC STRAINS AT BEAM SECTION NODAL POINTS

ELEMENT = 1 SECTION ID = 1

ELEMENT NODE = 1

SEC NODE E1 E2 E3 EINT EEQV

1 0.0000 0.0000 -0.26615E-04 0.26615E-04 0.26615E-04
3 0.0000 0.0000 -0.11619E-04 0.11619E-04 0.11619E-04
13 0.0000 0.0000 -0.11619E-04 0.11619E-04 0.11619E-04
 11 0.0000 0.0000 -0.26615E-04 0.26615E-04 0.26615E-04
5 0.33767E-05 0.0000 0.0000 0.33767E-05 0.33767E-05
15 0.33767E-05 0.0000 0.0000 0.33767E-05 0.33767E-05
23 0.0000 0.0000 -0.11619E-04 0.11619E-04 0.11619E-04
21 0.0000 0.0000 -0.26615E-04 0.26615E-04 0.26615E-04
25 0.33767E-05 0.0000 0.0000 0.33767E-05 0.33767E-05

Max= 0.33767E-05

Min= -0.26615E-04

ELEMENT NODE = 2

SEC NODE E1 E2 E3 EINT EEQV

1 0.0000 0.0000 -0.26615E-04 0.26615E-04 0.26615E-04
3 0.0000 0.0000 -0.11619E-04 0.11619E-04 0.11619E-04
13 0.0000 0.0000 -0.11619E-04 0.11619E-04 0.11619E-04
11 0.0000 0.0000 -0.26615E-04 0.26615E-04 0.26615E-04
5 0.33767E-05 0.0000 0.0000 0.33767E-05 0.33767E-05
15 0.33767E-05 0.0000 0.0000 0.33767E-05 0.33767E-05
23 0.0000 0.0000 -0.11619E-04 0.11619E-04 0.11619E-04
21 0.0000 0.0000 -0.26615E-04 0.26615E-04 0.26615E-04
25 0.33767E-05 0.0000 0.0000 0.33767E-05 0.33767E-05

Max= 0.33767E-05

Min= -0.26615E-04

ELASTIC STRAINS AT BEAM SECTION NODAL POINTS

ELEMENT = 2 SECTION ID = 1

ELEMENT NODE = 2

SEC NODE E1 E2 E3 EINT EEQV

1 0.0000 0.0000 0.0000 0.0000 0.0000
3 0.0000 0.0000 0.0000 0.0000 0.0000
13 0.0000 0.0000 0.0000 0.0000 0.0000
11 0.0000 0.0000 0.0000 0.0000 0.0000
5 0.0000 0.0000 0.0000 0.0000 0.0000
15 0.0000 0.0000 0.0000 0.0000 0.0000
23 0.0000 0.0000 0.0000 0.0000 0.0000
21 0.0000 0.0000 0.0000 0.0000 0.0000
25 0.0000 0.0000 0.0000 0.0000 0.0000

Max= 0.55983E-20

Min= -0.55983E-20

ELEMENT NODE = 3

SEC NODE E1 E2 E3 EINT EEQV

1 0.0000 0.0000 0.0000 0.0000 0.0000
3 0.0000 0.0000 0.0000 0.0000 0.0000
13 0.0000 0.0000 0.0000 0.0000 0.0000
11 0.0000 0.0000 0.0000 0.0000 0.0000
5 0.0000 0.0000 0.0000 0.0000 0.0000
15 0.0000 0.0000 0.0000 0.0000 0.0000
23 0.0000 0.0000 0.0000 0.0000 0.0000
21 0.0000 0.0000 0.0000 0.0000 0.0000
25 0.0000 0.0000 0.0000 0.0000 0.0000

Max= 0.55983E-20
 Min= -0.55983E-20

ELASTIC STRAINS AT BEAM SECTION NODAL POINTS

ELEMENT = 3 SECTION ID = 1

ELEMENT NODE = 2

SEC NODE E1 E2 E3 EINT EEQV

1 0.0000 0.0000 -0.47160E-03 0.47160E-03 0.47160E-03
3 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
13 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
11 0.0000 0.0000 -0.47160E-03 0.47160E-03 0.47160E-03
5 0.45280E-03 0.0000 0.0000 0.45280E-03 0.45280E-03
15 0.45280E-03 0.0000 0.0000 0.45280E-03 0.45280E-03
23 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
21 0.0000 0.0000 -0.47160E-03 0.47160E-03 0.47160E-03
25 0.45280E-03 0.0000 0.0000 0.45280E-03 0.45280E-03

Max= 0.45280E-03

Min= -0.47160E-03

ELEMENT NODE = 4

SEC NODE E1 E2 E3 EINT EEQV

1 0.0000 0.0000 -0.47160E-03 0.47160E-03 0.47160E-03
3 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
13 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
11 0.0000 0.0000 -0.47160E-03 0.47160E-03 0.47160E-03
5 0.45280E-03 0.0000 0.0000 0.45280E-03 0.45280E-03
15 0.45280E-03 0.0000 0.0000 0.45280E-03 0.45280E-03
23 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
21 0.0000 0.0000 -0.47160E-03 0.47160E-03 0.47160E-03
25 0.45280E-03 0.0000 0.0000 0.45280E-03 0.45280E-03

Max= 0.45280E-03

Min= -0.47160E-03

ELASTIC STRAINS AT BEAM SECTION NODAL POINTS

ELEMENT = 4 SECTION ID = 1

ELEMENT NODE = 4

SEC NODE E1 E2 E3 EINT EEQV

1 0.25096E-03 0.0000 0.0000 0.25096E-03 0.25096E-03
3 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
13 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
11 0.25096E-03 0.0000 0.0000 0.25096E-03 0.25096E-03
5 0.0000 0.0000 -0.26976E-03 0.26976E-03 0.26976E-03
15 0.0000 0.0000 -0.26976E-03 0.26976E-03 0.26976E-03
23 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
21 0.25096E-03 0.0000 0.0000 0.25096E-03 0.25096E-03
25 0.0000 0.0000 -0.26976E-03 0.26976E-03 0.26976E-03

Max= 0.25096E-03

Min= -0.26976E-03

ELEMENT NODE = 5
 SEC NODE E1 E2 E3 EINT EEQV
1 0.25096E-03 0.0000 0.0000 0.25096E-03 0.25096E-03
3 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
13 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
11 0.25096E-03 0.0000 0.0000 0.25096E-03 0.25096E-03
5 0.0000 0.0000 -0.26976E-03 0.26976E-03 0.26976E-03
15 0.0000 0.0000 -0.26976E-03 0.26976E-03 0.26976E-03
23 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
21 0.25096E-03 0.0000 0.0000 0.25096E-03 0.25096E-03
25 0.0000 0.0000 -0.26976E-03 0.26976E-03 0.26976E-03

Max= 0.25096E-03

Min= -0.26976E-03

ELASTIC STRAINS AT BEAM SECTION NODAL POINTS

ELEMENT = 5 SECTION ID = 1

ELEMENT NODE = 5

SEC NODE E1 E2 E3 EINT EEQV

1 0.15502E-03 0.0000 0.0000 0.15502E-03 0.15502E-03
3 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
13 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
11 0.15502E-03 0.0000 0.0000 0.15502E-03 0.15502E-03
5 0.0000 0.0000 -0.17382E-03 0.17382E-03 0.17382E-03
15 0.0000 0.0000 -0.17382E-03 0.17382E-03 0.17382E-03
23 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
21 0.15502E-03 0.0000 0.0000 0.15502E-03 0.15502E-03
25 0.0000 0.0000 -0.17382E-03 0.17382E-03 0.17382E-03

Max= 0.15502E-03

Min= -0.17382E-03

ELEMENT NODE = 6

SEC NODE E1 E2 E3 EINT EEQV

1 0.15502E-03 0.0000 0.0000 0.15502E-03 0.15502E-03
3 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
13 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
11 0.15502E-03 0.0000 0.0000 0.15502E-03 0.15502E-03
5 0.0000 0.0000 -0.17382E-03 0.17382E-03 0.17382E-03
15 0.0000 0.0000 -0.17382E-03 0.17382E-03 0.17382E-03
23 0.0000 0.0000 -0.93991E-05 0.93991E-05 0.93991E-05
21 0.15502E-03 0.0000 0.0000 0.15502E-03 0.15502E-03
25 0.0000 0.0000 -0.17382E-03 0.17382E-03 0.17382E-03

Max= 0.15502E-03

Min= -0.17382E-03

ELASTIC STRAINS AT BEAM SECTION NODAL POINTS

ELEMENT = 6 SECTION ID = 1

ELEMENT NODE = 6

SEC NODE E1 E2 E3 EINT EEQV

1 0.0000 0.0000 -0.11197E-19 0.11197E-19 0.11197E-19
3 0.0000 0.0000 0.0000 0.0000 0.0000
 13 0.0000 0.0000 0.0000 0.0000 0.0000
11 0.0000 0.0000 -0.11197E-19 0.11197E-19 0.11197E-19
5 0.11197E-19 0.0000 0.0000 0.11197E-19 0.11197E-19
15 0.11197E-19 0.0000 0.0000 0.11197E-19 0.11197E-19
23 0.0000 0.0000 0.0000 0.0000 0.0000
21 0.0000 0.0000 -0.11197E-19 0.11197E-19 0.11197E-19
25 0.11197E-19 0.0000 0.0000 0.11197E-19 0.11197E-19

Max= 0.11197E-19

Min= -0.11197E-19

ELEMENT NODE = 7

SEC NODE E1 E2 E3 EINT EEQV

1 0.0000 0.0000 -0.11197E-19 0.11197E-19 0.11197E-19
3 0.0000 0.0000 0.0000 0.0000 0.0000
13 0.0000 0.0000 0.0000 0.0000 0.0000
11 0.0000 0.0000 -0.11197E-19 0.11197E-19 0.11197E-19
5 0.11197E-19 0.0000 0.0000 0.11197E-19 0.11197E-19
15 0.11197E-19 0.0000 0.0000 0.11197E-19 0.11197E-19
23 0.0000 0.0000 0.0000 0.0000 0.0000
21 0.0000 0.0000 -0.11197E-19 0.11197E-19 0.11197E-19
25 0.11197E-19 0.0000 0.0000 0.11197E-19 0.11197E-19

Max= 0.11197E-19

Min= -0.11197E-19

ELASTIC STRAINS AT BEAM SECTION NODAL POINTS

ELEMENT = 7 SECTION ID = 1

ELEMENT NODE = 6

SEC NODE E1 E2 E3 EINT EEQV

1 0.0000 0.0000 -0.67781E-04 0.67781E-04 0.67781E-04
3 0.0000 0.0000 -0.14899E-04 0.14899E-04 0.14899E-04
13 0.0000 0.0000 -0.14899E-04 0.14899E-04 0.14899E-04
11 0.0000 0.0000 -0.67781E-04 0.67781E-04 0.67781E-04
5 0.37983E-04 0.0000 0.0000 0.37983E-04 0.37983E-04
15 0.37983E-04 0.0000 0.0000 0.37983E-04 0.37983E-04
23 0.0000 0.0000 -0.14899E-04 0.14899E-04 0.14899E-04
21 0.0000 0.0000 -0.67781E-04 0.67781E-04 0.67781E-04
25 0.37983E-04 0.0000 0.0000 0.37983E-04 0.37983E-04

Max= 0.37983E-04

Min= -0.67781E-04

ELEMENT NODE = 8

SEC NODE E1 E2 E3 EINT EEQV

1 0.0000 0.0000 -0.67781E-04 0.67781E-04 0.67781E-04
3 0.0000 0.0000 -0.14899E-04 0.14899E-04 0.14899E-04
13 0.0000 0.0000 -0.14899E-04 0.14899E-04 0.14899E-04
11 0.0000 0.0000 -0.67781E-04 0.67781E-04 0.67781E-04
5 0.37983E-04 0.0000 0.0000 0.37983E-04 0.37983E-04
15 0.37983E-04 0.0000 0.0000 0.37983E-04 0.37983E-04
23 0.0000 0.0000 -0.14899E-04 0.14899E-04 0.14899E-04
21 0.0000 0.0000 -0.67781E-04 0.67781E-04 0.67781E-04
25 0.37983E-04 0.0000 0.0000 0.37983E-04 0.37983E-04
 Max= 0.37983E-04

Min= -0.67781E-04

iii. Stress intensity

PRINT S PRIN ELEMENT SOLUTION PER ELEMENT

STRESSES AT BEAM SECTION NODAL POINTS

ELEMENT = 1 SECTION ID = 1

ELEMENT NODE = 1

SEC NODE S1 S2 S3 SINT SEQV

1 0.0000 0.0000 -0.55892E+07 0.55892E+07 0.55892E+07
3 0.0000 0.0000 -0.24400E+07 0.24400E+07 0.24400E+07
13 0.0000 0.0000 -0.24400E+07 0.24400E+07 0.24400E+07
11 0.0000 0.0000 -0.55892E+07 0.55892E+07 0.55892E+07
5 0.70911E+06 0.0000 0.0000 0.70911E+06 0.70911E+06
15 0.70911E+06 0.0000 0.0000 0.70911E+06 0.70911E+06
23 0.0000 0.0000 -0.24400E+07 0.24400E+07 0.24400E+07
21 0.0000 0.0000 -0.55892E+07 0.55892E+07 0.55892E+07
25 0.70911E+06 0.0000 0.0000 0.70911E+06 0.70911E+06

Max= 0.70911E+06 0.0000 0.0000 0.55892E+07 0.55892E+07

Min= 0.0000 0.0000 -0.55892E+07 0.70911E+06 0.70911E+06

ELEMENT NODE = 2

SEC NODE S1 S2 S3 SINT SEQV

1 0.0000 0.0000 -0.55892E+07 0.55892E+07 0.55892E+07
3 0.0000 0.0000 -0.24400E+07 0.24400E+07 0.24400E+07
13 0.0000 0.0000 -0.24400E+07 0.24400E+07 0.24400E+07
11 0.0000 0.0000 -0.55892E+07 0.55892E+07 0.55892E+07
5 0.70911E+06 0.0000 0.0000 0.70911E+06 0.70911E+06
15 0.70911E+06 0.0000 0.0000 0.70911E+06 0.70911E+06
23 0.0000 0.0000 -0.24400E+07 0.24400E+07 0.24400E+07
21 0.0000 0.0000 -0.55892E+07 0.55892E+07 0.55892E+07
25 0.70911E+06 0.0000 0.0000 0.70911E+06 0.70911E+06

Max= 0.70911E+06 0.0000 0.0000 0.55892E+07 0.55892E+07

Min= 0.0000 0.0000 -0.55892E+07 0.70911E+06 0.70911E+06

STRESSES AT BEAM SECTION NODAL POINTS

ELEMENT = 2 SECTION ID = 1

ELEMENT NODE = 2

SEC NODE S1 S2 S3 SINT SEQV

1 0.11756E-08 0.0000 0.0000 0.11756E-08 0.11756E-08
3 0.0000 0.0000 0.0000 0.0000 0.0000
13 0.0000 0.0000 0.0000 0.0000 0.0000
11 0.11756E-08 0.0000 0.0000 0.11756E-08 0.11756E-08
5 0.0000 0.0000 -0.11756E-08 0.11756E-08 0.11756E-08
15 0.0000 0.0000 -0.11756E-08 0.11756E-08 0.11756E-08
23 0.0000 0.0000 0.0000 0.0000 0.0000
21 0.11756E-08 0.0000 0.0000 0.11756E-08 0.11756E-08
 25 0.0000 0.0000 -0.11756E-08 0.11756E-08 0.11756E-08

Max= 0.11756E-08 0.0000 0.0000 0.11756E-08 0.11756E-08

Min= 0.0000 0.0000 -0.11756E-08 0.0000 0.0000

ELEMENT NODE = 3

SEC NODE S1 S2 S3 SINT SEQV

1 0.11756E-08 0.0000 0.0000 0.11756E-08 0.11756E-08
3 0.0000 0.0000 0.0000 0.0000 0.0000
13 0.0000 0.0000 0.0000 0.0000 0.0000
11 0.11756E-08 0.0000 0.0000 0.11756E-08 0.11756E-08
5 0.0000 0.0000 -0.11756E-08 0.11756E-08 0.11756E-08
15 0.0000 0.0000 -0.11756E-08 0.11756E-08 0.11756E-08
23 0.0000 0.0000 0.0000 0.0000 0.0000
21 0.11756E-08 0.0000 0.0000 0.11756E-08 0.11756E-08
25 0.0000 0.0000 -0.11756E-08 0.11756E-08 0.11756E-08

Max= 0.11756E-08 0.0000 0.0000 0.11756E-08 0.11756E-08

Min= 0.0000 0.0000 -0.11756E-08 0.0000 0.0000

STRESSES AT BEAM SECTION NODAL POINTS

ELEMENT = 3 SECTION ID = 1

ELEMENT NODE = 2

SEC NODE S1 S2 S3 SINT SEQV

1 0.0000 0.0000 -0.99036E+08 0.99036E+08 0.99036E+08
3 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
13 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
11 0.0000 0.0000 -0.99036E+08 0.99036E+08 0.99036E+08
5 0.95088E+08 0.0000 0.0000 0.95088E+08 0.95088E+08
15 0.95088E+08 0.0000 0.0000 0.95088E+08 0.95088E+08
23 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
21 0.0000 0.0000 -0.99036E+08 0.99036E+08 0.99036E+08
25 0.95088E+08 0.0000 0.0000 0.95088E+08 0.95088E+08

Max= 0.95088E+08 0.0000 0.0000 0.99036E+08 0.99036E+08

Min= 0.0000 0.0000 -0.99036E+08 0.19738E+07 0.19738E+07

ELEMENT NODE = 4

SEC NODE S1 S2 S3 SINT SEQV

1 0.0000 0.0000 -0.99036E+08 0.99036E+08 0.99036E+08
3 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
13 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
11 0.0000 0.0000 -0.99036E+08 0.99036E+08 0.99036E+08
5 0.95088E+08 0.0000 0.0000 0.95088E+08 0.95088E+08
15 0.95088E+08 0.0000 0.0000 0.95088E+08 0.95088E+08
23 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
21 0.0000 0.0000 -0.99036E+08 0.99036E+08 0.99036E+08
25 0.95088E+08 0.0000 0.0000 0.95088E+08 0.95088E+08

Max= 0.95088E+08 0.0000 0.0000 0.99036E+08 0.99036E+08

Min= 0.0000 0.0000 -0.99036E+08 0.19738E+07 0.19738E+07

STRESSES AT BEAM SECTION NODAL POINTS

ELEMENT = 4 SECTION ID = 1

ELEMENT NODE = 4

SEC NODE S1 S2 S3 SINT SEQV

1 0.52701E+08 0.0000 0.0000 0.52701E+08 0.52701E+08
3 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
13 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
11 0.52701E+08 0.0000 0.0000 0.52701E+08 0.52701E+08
5 0.0000 0.0000 -0.56649E+08 0.56649E+08 0.56649E+08
15 0.0000 0.0000 -0.56649E+08 0.56649E+08 0.56649E+08
23 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
21 0.52701E+08 0.0000 0.0000 0.52701E+08 0.52701E+08
25 0.0000 0.0000 -0.56649E+08 0.56649E+08 0.56649E+08

Max= 0.52701E+08 0.0000 0.0000 0.56649E+08 0.56649E+08

Min= 0.0000 0.0000 -0.56649E+08 0.19738E+07 0.19738E+07

ELEMENT NODE = 5

SEC NODE S1 S2 S3 SINT SEQV

1 0.52701E+08 0.0000 0.0000 0.52701E+08 0.52701E+08
3 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
13 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
11 0.52701E+08 0.0000 0.0000 0.52701E+08 0.52701E+08
5 0.0000 0.0000 -0.56649E+08 0.56649E+08 0.56649E+08
15 0.0000 0.0000 -0.56649E+08 0.56649E+08 0.56649E+08
23 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
21 0.52701E+08 0.0000 0.0000 0.52701E+08 0.52701E+08
25 0.0000 0.0000 -0.56649E+08 0.56649E+08 0.56649E+08

Max= 0.52701E+08 0.0000 0.0000 0.56649E+08 0.56649E+08

Min= 0.0000 0.0000 -0.56649E+08 0.19738E+07 0.19738E+07

STRESSES AT BEAM SECTION NODAL POINTS

ELEMENT = 5 SECTION ID = 1

ELEMENT NODE = 5

SEC NODE S1 S2 S3 SINT SEQV

1 0.32554E+08 0.0000 0.0000 0.32554E+08 0.32554E+08
3 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
13 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
11 0.32554E+08 0.0000 0.0000 0.32554E+08 0.32554E+08
5 0.0000 0.0000 -0.36502E+08 0.36502E+08 0.36502E+08
15 0.0000 0.0000 -0.36502E+08 0.36502E+08 0.36502E+08
23 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
21 0.32554E+08 0.0000 0.0000 0.32554E+08 0.32554E+08
25 0.0000 0.0000 -0.36502E+08 0.36502E+08 0.36502E+08

Max= 0.32554E+08 0.0000 0.0000 0.36502E+08 0.36502E+08

Min= 0.0000 0.0000 -0.36502E+08 0.19738E+07 0.19738E+07

ELEMENT NODE = 6

SEC NODE S1 S2 S3 SINT SEQV

1 0.32554E+08 0.0000 0.0000 0.32554E+08 0.32554E+08
3 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
13 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
 11 0.32554E+08 0.0000 0.0000 0.32554E+08 0.32554E+08
5 0.0000 0.0000 -0.36502E+08 0.36502E+08 0.36502E+08
15 0.0000 0.0000 -0.36502E+08 0.36502E+08 0.36502E+08
23 0.0000 0.0000 -0.19738E+07 0.19738E+07 0.19738E+07
21 0.32554E+08 0.0000 0.0000 0.32554E+08 0.32554E+08
25 0.0000 0.0000 -0.36502E+08 0.36502E+08 0.36502E+08

Max= 0.32554E+08 0.0000 0.0000 0.36502E+08 0.36502E+08

Min= 0.0000 0.0000 -0.36502E+08 0.19738E+07 0.19738E+07

STRESSES AT BEAM SECTION NODAL POINTS

ELEMENT = 6 SECTION ID = 1

ELEMENT NODE = 6

SEC NODE S1 S2 S3 SINT SEQV

1 0.0000 0.0000 -0.23513E-08 0.23513E-08 0.23513E-08
3 0.0000 0.0000 0.0000 0.0000 0.0000
13 0.0000 0.0000 0.0000 0.0000 0.0000
11 0.0000 0.0000 -0.23513E-08 0.23513E-08 0.23513E-08
5 0.23513E-08 0.0000 0.0000 0.23513E-08 0.23513E-08
15 0.23513E-08 0.0000 0.0000 0.23513E-08 0.23513E-08
23 0.0000 0.0000 0.0000 0.0000 0.0000
21 0.0000 0.0000 -0.23513E-08 0.23513E-08 0.23513E-08
25 0.23513E-08 0.0000 0.0000 0.23513E-08 0.23513E-08

Max= 0.23513E-08 0.0000 0.0000 0.23513E-08 0.23513E-08

Min= 0.0000 0.0000 -0.23513E-08 0.0000 0.0000

ELEMENT NODE = 7

SEC NODE S1 S2 S3 SINT SEQV

1 0.0000 0.0000 -0.23513E-08 0.23513E-08 0.23513E-08
3 0.0000 0.0000 0.0000 0.0000 0.0000
13 0.0000 0.0000 0.0000 0.0000 0.0000
11 0.0000 0.0000 -0.23513E-08 0.23513E-08 0.23513E-08
5 0.23513E-08 0.0000 0.0000 0.23513E-08 0.23513E-08
15 0.23513E-08 0.0000 0.0000 0.23513E-08 0.23513E-08
23 0.0000 0.0000 0.0000 0.0000 0.0000
21 0.0000 0.0000 -0.23513E-08 0.23513E-08 0.23513E-08
25 0.23513E-08 0.0000 0.0000 0.23513E-08 0.23513E-08

Max= 0.23513E-08 0.0000 0.0000 0.23513E-08 0.23513E-08

Min= 0.0000 0.0000 -0.23513E-08 0.0000 0.0000

STRESSES AT BEAM SECTION NODAL POINTS

ELEMENT = 7 SECTION ID = 1

ELEMENT NODE = 6

SEC NODE S1 S2 S3 SINT SEQV

1 0.0000 0.0000 -0.14234E+08 0.14234E+08 0.14234E+08
3 0.0000 0.0000 -0.31288E+07 0.31288E+07 0.31288E+07
13 0.0000 0.0000 -0.31288E+07 0.31288E+07 0.31288E+07
11 0.0000 0.0000 -0.14234E+08 0.14234E+08 0.14234E+08
5 0.79764E+07 0.0000 0.0000 0.79764E+07 0.79764E+07
15 0.79764E+07 0.0000 0.0000 0.79764E+07 0.79764E+07
23 0.0000 0.0000 -0.31288E+07 0.31288E+07 0.31288E+07
 21 0.0000 0.0000 -0.14234E+08 0.14234E+08 0.14234E+08
25 0.79764E+07 0.0000 0.0000 0.79764E+07 0.79764E+07

Max= 0.79764E+07 0.0000 0.0000 0.14234E+08 0.14234E+08

Min= 0.0000 0.0000 -0.14234E+08 0.31288E+07 0.31288E+07

ELEMENT NODE = 8

SEC NODE S1 S2 S3 SINT SEQV

1 0.0000 0.0000 -0.14234E+08 0.14234E+08 0.14234E+08
3 0.0000 0.0000 -0.31288E+07 0.31288E+07 0.31288E+07
13 0.0000 0.0000 -0.31288E+07 0.31288E+07 0.31288E+07
11 0.0000 0.0000 -0.14234E+08 0.14234E+08 0.14234E+08
5 0.79764E+07 0.0000 0.0000 0.79764E+07 0.79764E+07
15 0.79764E+07 0.0000 0.0000 0.79764E+07 0.79764E+07
23 0.0000 0.0000 -0.31288E+07 0.31288E+07 0.31288E+07
21 0.0000 0.0000 -0.14234E+08 0.14234E+08 0.14234E+08
25 0.79764E+07 0.0000 0.0000 0.79764E+07 0.79764E+07

Max= 0.79764E+07 0.0000 0.0000 0.14234E+08 0.14234E+08

Min= 0.0000 0.0000 -0.14234E+08 0.31288E+07 0.31288E+07

Conclusion
The comparison of the results by ANSYS and MATLAB software according to the above results
can concluded as follow

The solutions obtained by MATLAB are a little bit different from that of ANSYS
solution. Thus, the MATLAB solutions are lower than ANSYS solutions. This indicates
that if we increase mesh number on the ANSYS, the solution will approach to that of
MATLAB solution.
Therefore, the variation between MATLAB and ANSYS solutions are due to
discretization error. The fact that number of discretization element (mesh number)
increases, the answer will approach to the accurate solution.
Fine discretization gives more accurate solutions and the reverse is true that of coarse
discretization.