This question is perfectly described in the 7th edition of this book.

Zienkiewicz OC, Taylor RL, Fox DD (2014) The finite element method for solid and structural mechanics, 7th edn. Butterworth-Heinemann, Oxford.

Hi Pratik Bhala!

Here is a simple code for curve element modeling

CSYS, 5 ! Activate cylindrical coordinate with Y as axis of rotation

Rbr=50;

PI=3.1416

Lapp1=20

ELE_SIZE=0.5

K, (KPmax+115), (Rbr-3.25), ((Lapp1*180)/(Rbr*PI)), -0.3

K, (KPmax+116), (Rbr-2.75), ((Lapp1*180)/(Rbr*PI)), -0.3

K, (KPmax+215), (Rbr-3.25), (((Lapp1+Lbr)*180)/(Rbr*PI)/2), -0.3

K, (KPmax+216), (Rbr-2.75), (((Lapp1+Lbr)*180)/(Rbr*PI)/2), -0.3

K, (KPmax+315), (Rbr-3.25), (((Lapp1+Lbr)*180)/(Rbr*PI)), -0.3

K, (KPmax+316), (Rbr-2.75), (((Lapp1+Lbr)*180)/(Rbr*PI)), -0.3

ASEL, NONE $ LSEL, NONE $ ESEL, NONE

LARC, (KPmax+115), (KPmax+315), (KPmax+215),

LARC, (KPmax+116), (KPmax+316), (KPmax+216),

LESIZE, ALL, , , (Lbr/ELE_SIZE), , , , , 0

L, (KPmax+115), (KPmax+116), 2

L, (KPmax+315), (KPmax+316), 2

AL, ALL

Analysis of Distortion-Induced Stress and Retrofitting Technique of Curved Twin I-Girder Composite Bridge | SpringerLink