clear
//Given
load_distributed = 20 //kN/sq.m, This is the load distributed over the pier
H = 2          // m, Total height
h = 1          //m , point of investigation
base = 1.5     //m The length of crossection in side view
top = 0.5      //m ,The length where load is distributed on top
base_inv = 1   //m , the length at the point of investigation
area = 0.5*1   //m ,The length at a-a cross-section
density_conc = 25 //kN/sq.m
//calculation of total weight

v_total = ((top+base)/2)*top*H       //sq.m ,The total volume
w_total = v_total* density_conc  //kN , The total weight
R_top = (top**2)*load_distributed    //kN , THe reaction force due to load distribution
reaction_net = w_total + R_top

//calculation of State of stress at 1m
v_inv = ((top+base_inv)/2)*top*h    //sq.m ,The total volume from 1m to top
w_inv = v_inv*density_conc          //kN , The total weight from 1m to top
reaction_net = w_inv + R_top        //kN
Stress = reaction_net/area           //kN/sq.m
printf("\n The total weight of pier is %0.3f kN",w_total)
printf("\n The stress at 1 m above is %0.1f kN/sq.m",Stress)