kill(all);

risolvo le equazioni di equilibrio nelle incognite NA, CC, NC

linsolve(
    [
        -NA    =0,
        +F/2-NC=0,

        F/2 * 0 + NA*(r+l)-CA - NC*r + CC = 0  
    ]
    ,
    [NA, NC, CC]
), globalsolve=True;

caratteristiche di sollecitazione sul tratto (curvo) AB
da equilibrio del concio compreso tra A e l'angolo theta

linsolve(
    [
        N_AB - NA * cos(theta) - F/2*sin(theta)=0,
        T_AB + F/2*cos(theta)-NA*sin(theta)=0,
        Mf_AB-CA-F/2*r*sin(theta)+NA*r*(1-cos(theta))=0
    ],

    [N_AB, T_AB, Mf_AB]
), globalsolve=True;