Beam Element (3D)

Stiffness matrix of beam element in 3D

This function finds the stiffness matrix of beam element (frame element) in 3D.

function k=beam_stiffness(E,Iy,Iz,A,L,G,J,Cx,Cy,Cz,omega)

Input:

E: Modulus of Elasticity

Iy, Iz: Moment of Inertia around y and z axes

A: Section Area

L: Length

G: Shear Modulus

J: Polar Moment of Inertia

Cx, Cy, and Cz: Direction cosines of the beam with respect to x, y, and z axes

omega: Rotation of the beam around the beam major axis (x-axis)

Output:

k: Stiffness Matrix