Skip to main content

SiMP

oxiphysics_core::topology_optimization

Struct SiMP

pub struct SiMP {
    pub penalty: f64,
    pub material: MaterialParams,
    pub volume_fraction: f64,
    pub filter_radius: f64,
    pub grid: GridSize,
    pub density: Vec<f64>,
}

Expand description

Solid Isotropic Material with Penalization (SIMP) model.

SIMP interpolates material stiffness as E(ρ) = E_void + ρ^p * (E_solid - E_void). The penalization parameter p discourages intermediate densities.

Fields§

§penalty: f64

Penalization exponent.

§material: MaterialParams

Material parameters.

§volume_fraction: f64

Volume fraction target.

§filter_radius: f64

Filter radius for density filtering.

§grid: GridSize

Grid size.

§density: Vec<f64>

Current density field (length = grid.n_elem()).

Implementations§

impl SiMP

pub fn new( grid: GridSize, volume_fraction: f64, penalty: f64, filter_radius: f64, material: MaterialParams, ) -> Self

Creates a new SIMP model with uniform initial density equal to volume_fraction.

pub fn interpolated_stiffness(&self, rho: f64) -> f64

Computes the interpolated Young’s modulus for a given density.

E(ρ) = E_void + ρ^p * (E_solid - E_void)

pub fn stiffness_sensitivity(&self, rho: f64) -> f64

Computes the derivative of interpolated stiffness with respect to density.

dE/dρ = p * ρ^(p-1) * (E_solid - E_void)

pub fn stiffness_field(&self) -> Vec<f64>

Returns vector of interpolated stiffnesses for current density field.

pub fn sensitivity_field(&self) -> Vec<f64>

Returns sensitivities dE/dρ for current density field.

pub fn compliance_sensitivity(&self, element_strain_energy: &[f64]) -> Vec<f64>

Computes element-wise compliance sensitivity given element displacements.

dc/dρ_e = -p * ρ_e^(p-1) * (E_solid - E_void) * u_e^T * K_e_unit * u_e

pub fn current_volume(&self) -> f64

Computes current volume fraction.

pub fn is_valid(&self) -> bool

Returns whether the density field is valid (all values in [0,1]).

Trait Implementations§

impl Clone for SiMP

fn clone(&self) -> SiMP

Returns a duplicate of the value. Read more

1.0.0 (const: unstable) · Source§

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source. Read more

impl Debug for SiMP

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter. Read more

Auto Trait Implementations§

impl Freeze for SiMP

impl RefUnwindSafe for SiMP

impl Send for SiMP

impl Sync for SiMP

impl Unpin for SiMP

impl UnsafeUnpin for SiMP

impl UnwindSafe for SiMP

Blanket Implementations§

impl<T> Any for T
where T: 'static + ?Sized,

fn type_id(&self) -> TypeId

Gets the TypeId of self. Read more

impl<T> Borrow<T> for T
where T: ?Sized,

fn borrow(&self) -> &T

Immutably borrows from an owned value. Read more

impl<T> BorrowMut<T> for T
where T: ?Sized,

fn borrow_mut(&mut self) -> &mut T

Mutably borrows from an owned value. Read more

impl<T> CloneToUninit for T
where T: Clone,

unsafe fn clone_to_uninit(&self, dest: *mut u8)

🔬This is a nightly-only experimental API. (clone_to_uninit)

Performs copy-assignment from self to dest. Read more

impl<T> From<T> for T

fn from(t: T) -> T

Returns the argument unchanged.

impl<T, U> Into<U> for T
where U: From<T>,

fn into(self) -> U

Calls U::from(self).

That is, this conversion is whatever the implementation of From<T> for U chooses to do.

impl<T> Same for T

type Output = T

Should always be Self

impl<SS, SP> SupersetOf<SS> for SP
where SS: SubsetOf<SP>,

fn to_subset(&self) -> Option<SS>

The inverse inclusion map: attempts to construct self from the equivalent element of its superset. Read more

fn is_in_subset(&self) -> bool

Checks if self is actually part of its subset T (and can be converted to it).

fn to_subset_unchecked(&self) -> SS

Use with care! Same as self.to_subset but without any property checks. Always succeeds.

fn from_subset(element: &SS) -> SP

The inclusion map: converts self to the equivalent element of its superset.

impl<T> ToOwned for T
where T: Clone,

type Owned = T

The resulting type after obtaining ownership.

fn to_owned(&self) -> T

Creates owned data from borrowed data, usually by cloning. Read more

fn clone_into(&self, target: &mut T)

Uses borrowed data to replace owned data, usually by cloning. Read more

impl<T, U> TryFrom<U> for T
where U: Into<T>,

type Error = Infallible

The type returned in the event of a conversion error.

fn try_from(value: U) -> Result<T, <T as TryFrom<U>>::Error>

Performs the conversion.

impl<T, U> TryInto<U> for T
where U: TryFrom<T>,

type Error = <U as TryFrom<T>>::Error

The type returned in the event of a conversion error.

fn try_into(self) -> Result<U, <U as TryFrom<T>>::Error>

Performs the conversion.