Static¶
Static linear solver
Assemble global matrices
Generate profiling report for assembly and save it in self.prof_out.
- To read the report:
import pstats p=pstats.Stats(self.prof_out)
Gains to reproduce rigid-body motion such that grid displacements and velocities are given by:
dzeta = Ktra*u_tra + Krot*u_rot
dzeta_dot = Ktra_vel*u_tra_dot + Krot*u_rot_dot
Rotations are assumed to happen independently with respect to the zeta_rotation point and about the x,y and z axes of the inertial frame.
Gains to computes sectional forces. Moments are computed w.r.t. mid-vertex (chord-wise index M/2) of each section.
Calculates gain matrices to calculate the total force (Kftot) and moment (Kmtot, Kmtot_disp) about the pole zeta_pole.
Being \(f\) and \(\zeta\) the force and position at the vertex (m,n) of the lattice these are produced as:
ftot=sum(f) -> dftot += df
mtot-sum((zeta-zeta_pole) x f) -> dmtot += cross(zeta0-zeta_pole) df - cross(f0) dzeta
Reshapes state/output according to SHARPy format
Solve for bound \(\\Gamma\) using the equation;
\[\begin{split}\\mathcal{A}(\\Gamma^n) = u^n\end{split}\]# … at constant rotation speed
self.Dfqsdzeta+=scalg.block_diag(*ass.dfqsdzeta_omega(MS.Surfs,MS.Surfs_star))
Calculates total force (Ftot) and moment (Mtot) (about pole zeta_pole).