sciforge-core 0.0.4

Shared engineering primitives: materials, fasteners, fluids, propulsion taxonomies, std components.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67

use crate::moleculars::Material;
use sciforge_hub::prelude::constants::physics::solid_mechanics::friction;
use sciforge_hub::prelude::constants::physics::solid_mechanics::wear;
use sciforge_hub::prelude::physics::solid_mechanics::contact as sf_contact;

#[derive(Debug, Clone, Copy)]
pub struct Tribology {
    pub material_a: Material,
    pub material_b: Material,
}

impl Tribology {
    pub fn new(material_a: Material, material_b: Material) -> Self {
        Self { material_a, material_b }
    }

    pub fn mu_static(&self, condition: &str) -> Option<f64> {
        friction::by_pair(self.material_a.formula, self.material_b.formula, condition)
            .map(|f| f.mu_static)
    }

    pub fn mu_kinetic(&self, condition: &str) -> Option<f64> {
        friction::by_pair(self.material_a.formula, self.material_b.formula, condition)
            .map(|f| f.mu_kinetic)
    }

    pub fn friction_force_n(
        &self,
        normal_force_n: f64,
        condition: &str,
        moving: bool,
    ) -> Option<f64> {
        let mu = if moving {
            self.mu_kinetic(condition)?
        } else {
            self.mu_static(condition)?
        };
        Some(mu * normal_force_n)
    }

    pub fn archard_wear_constant(&self, condition: &str) -> Option<f64> {
        wear::by_pair(self.material_a.formula, self.material_b.formula, condition)
            .map(|w| w.archard_k)
    }

    pub fn archard_wear_volume_m3(
        &self,
        normal_force_n: f64,
        sliding_distance_m: f64,
        condition: &str,
    ) -> Option<f64> {
        let k = self.archard_wear_constant(condition)?;
        let hv_softer = self
            .material_a
            .vickers_hv()
            .into_iter()
            .chain(self.material_b.vickers_hv())
            .reduce(f64::min)?;
        let hardness_pa = hv_softer * 9.81e6;
        Some(sf_contact::archard_wear_volume(
            k,
            normal_force_n,
            sliding_distance_m,
            hardness_pa,
        ))
    }
}