diff options
Diffstat (limited to '3588/CH9/EX9.2/EX9_2.sce')
-rw-r--r-- | 3588/CH9/EX9.2/EX9_2.sce | 25 |
1 files changed, 25 insertions, 0 deletions
diff --git a/3588/CH9/EX9.2/EX9_2.sce b/3588/CH9/EX9.2/EX9_2.sce new file mode 100644 index 000000000..659f54494 --- /dev/null +++ b/3588/CH9/EX9.2/EX9_2.sce @@ -0,0 +1,25 @@ +//Clearing console +clc +clear + +//Intializing Variables +t = 0.2 +p = 7.4*(10^-4) +FBx = 0 +FBy = -386.4 + +//As the x component of the body force is zero, the x components of the nodal force vector will be zero +f1x = 0 +f2x = 0 +f3x = 0 + +//The y components of the nodal force vector +f1y = p*t*FBy*integrate('2*x*(1-x)','x',0,1) +f2y = p*t*FBy*integrate('2*x*(x)','x',0,1)/2 +f3y = p*t*FBy*integrate('2*x*(x)','x',0,1)/2 + +printf('\nResults\n') +printf('\nThe x components of the nodal force vector\nf1x=%flb \nf2x=%flb \nf3x=%flb',f1x,f2x,f3x) +printf('\nThe y components of the nodal force vector\nf1y=%flb \nf2y=%flb \nf3y=%flb',f1y,f2y,f3y) +printf('\nbody force is equally distributed to the element nodes') + |