Scribd
Upload
302 views3 pages

Chapter 8 MultiFreedom Constraints Solutions

This document contains solutions to homework exercises from Chapter 8 on multi-freedom constraints. It provides the results of running Mathematica programs to solve Exercises 8.1 and 8.2 numerically. For Exercise 8.1, it shows the stiffness matrix, applied forces, transformation matrix, constraint gap vector, modified stiffness matrix, and solution for a 6-element bar problem. For Exercise 8.2, it shows the code used to solve a similar problem with multiple constraints using a master-slave approach. The remaining exercises are not assigned or suggested as a challenging problem for employment.

Uploaded by

halcyon304
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
302 views3 pages

Chapter 8 MultiFreedom Constraints Solutions

This document contains solutions to homework exercises from Chapter 8 on multi-freedom constraints. It provides the results of running Mathematica programs to solve Exercises 8.1 and 8.2 numerically. For Exercise 8.1, it shows the stiffness matrix, applied forces, transformation matrix, constraint gap vector, modified stiffness matrix, and solution for a 6-element bar problem. For Exercise 8.2, it shows the code used to solve a similar problem with multiple constraints using a master-slave approach. The remaining exercises are not assigned or suggested as a challenging problem for employment.

Uploaded by

halcyon304
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

817

Solutions to Exercises

Homework Exercises for Chapter 8. - Solutions


MultiFreedom Constraints I

EXERCISE 8.1 Here are the results of running the Mathematica program under version 4.2:

0
0
0
0
0 
 100 100







100
200
100
0
0
0
0 







 0


100
200
100
0
0
0










0
0
100
200
100
0
0
Stiffness K








0
0
0
100 200 100
0 












0
0
0
0
100
200
100






0
0
0
0
100 100 
 0
Applied forces1, 2, 3, 4, 5, 6, 7
0



0




0




0





0




0



1
1
Constraint gap vector g 0, 0, 0, 0, 0,   , 0

5
1
0
0
0
0
0 









0
400
100
0
100
100










0
100
200
100
0
0






Modified Stiffness upon fixing node 1:



0
0
100
200
100
0









0 100
0
100 200
0 






0
100
0
0
0
100


1





0





0




0
Transformation matrix T





0






0



0

0
1
0
0
0
1
0

0
0
1
0
0
0
0

0
0
0
1
0
0
0

0
0
0
0
1
0
0

Modified RHS upon fixing node 1:0, 48, 3, 4, 15, 13


27
11
1
37
7
Computed umod lacks slave u6  0,  ,  ,  ,  , 

100
40
4
200
50
27
11
1
37
7
7
Complete solution u 0,  ,  ,  ,  ,  , 

100
40
4
200
100
50
Numerical u0., 0.27, 0.275, 0.25, 0.185, 0.07, 0.14
53
37
Recovered forces K.u with reactions 27,  , 3, 4, 5,   , 7

2
2
Numerical K.u27., 26.5, 3., 4., 5., 18.5, 7.

EXERCISE 8.2

The following Mathematica script solves this Exercise:

817

818

Chapter 8: MULTIFREEDOM CONSTRAINTS I

(* Exercise 8.2 - Master-Slave Method *)


(* MFCs: u2-u6=1/5, u3+2u4=-2/3, 2u3-u4+u5=0 - slaves: u4,u5,u6 *)
K=MasterStiffnessOfSixElementBar[100];
Print["Stiffness K=",K//MatrixForm];
f={1,2,3,4,5,6,7}; Print["Applied forces=",f];
T={{1, 0, 0, 0}, {0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, -1/2, 0},
{0, 0, -5/2, 0}, {0, 1, 0, 0}, {0, 0, 0, 1}};
Print["Transf matrix T (transposed)=",Transpose[T]//MatrixForm];
g={0, 0, 0, -1/3, -1/3, -1/5, 0};
Print["Constraint gap vector g=",g];
Khat=Simplify[Transpose[T].K.T]; fhat=Simplify[Transpose[T].(f-K.g)];
{Kmod,fmod}=FixLeftEndOfSixElementBar[Khat,fhat]; (* fix left end *)
Print["Modified Stiffness upon fixing node 1:",Kmod//MatrixForm];
Print["Modified RHS upon fixing node 1:",fmod];
umod=LinearSolve[Kmod,fmod];
Print["Computed umod (lacks slaves)=",umod];
u=T.umod+g; Print["Complete solution u=",u];
Print["Numerical u=",SetPrecision[N[u],5]];
fu=K.u; Print["Recovered forces K.u with reactions=",fu];
Print["Numerical K.u=",SetPrecision[N[fu],5]];

Modules MasterStiffnessOfSixElementBar and FixLeftEndOfSixElementBar are listed in the statement of Exercise 8.1. T and g were actually computed by the commands shown in the statement of the Exercise,
printed in InputForm and fed into the script. But they could have also been computed by hand. Running this
code gives
0
0
0
0
0 
 100 100








100
200
100
0
0
0
0







 0


100
200
100
0
0
0








 0
0
100 200 100
0
0 
Stiffness K










0
0
0
100
200
100
0










 0
0
0
0
100
200
100






0
0
0
0
0
100
100


Applied forces1, 2, 3, 4, 5, 6, 7
1




0



Transf matrix T transposed 


0




0

0 0
1 0

0
0

1
0 1  

2

0
0
5
 

2

0 0



1 0





0 0




0 1

0 0 0
0
1
1
1
Constraint gap vector g 0, 0, 0,   ,   ,   , 0

3
3
5
1
0
0
0 






 0 400 150 100 





Modified Stiffness upon fixing node 1:





0
150
1350
0






0
100 
 0 100
44
569
Modified RHS upon fixing node 1: 0,  ,   , 13

3
6

659
287
133
Computed umod lacks slaves  0,  ,   ,  

15300
3825
1530
659
287
2263
223
2401
133
Complete solution u 0,  ,   ,   ,   ,   ,  

15300
3825
7650
1530
15300
1530
Numerical u0, 0.043072, 0.075033, 0.29582, 0.14575, 0.15693, 0.086928
659
274
1571
5674
2467
138
Recovered forces K.u with reactions   ,  ,  ,   ,  ,   , 7

153
17
153
153
153
17
Numerical K.u4.3072, 16.118, 10.268, 37.085, 16.124, 8.1176, 7.0000

818

819

Solutions to Exercises

EXERCISE 8.3 Never assigned.


EXERCISE 8.4 Never assigned.
EXERCISE 8.5 Never assigned.
EXERCISE 8.6 Never assigned.
EXERCISE 8.7 If you can solve this one as a hobby, apply for a job with MSC Software.

819

You might also like