cppSimplifiedExpStrucModel/ModelPSstruc.hpp at master · willwiz/cppSimplifiedExpStrucModel · GitHub

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
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
//
//  ModelPSstruc.hpp
//  fullstructuralmodel
//
//  Created by Will Zhang on 6/12/18.
//  Copyright © 2018 Will Zhang. All rights reserved.
//

#ifndef ModelPSstruc_hpp
#define ModelPSstruc_hpp

#include <stdio.h>
#include "ModelMatrixWZ.hpp"
#include "ModelCollagenEXLstruc.hpp"


class model_ps_struc
{
protected:

    double rate_PS;

    double time_elapsed;

    double kM;
    double ratio;
    double alpham1;
    double betam1;

    double strainhistory[4];

    double kC;
    double ODFmean;
    double ODFstdev;
    double Dxmean;
    double Dxstdev;
    double Dxlb;
    double Dxub;
    double Dx_aniso;
    double kI;

    double F11_pushforward;
    double F12_pushforward;
    double F21_pushforward;
    double F22_pushforward;

public:

    double fraction_old;
    double fraction_new;

    model_matrix_yeoh_wz matrixmodel;

    // This is the updated standard structural model.
//    model_col_tissue_struc_PF collagenmodel;

    // This give the full exogenously crosslinked constitutive model
    model_col_int_tissue_struc_PF collagenmodel;

    void calc_new_old_fraction(double rate, double time);

    void update_reference(double arg[4]);

    void update_history(double arg[4]);

    void updata_parameters(double arg[13]);

    model_ps_struc():matrixmodel(), collagenmodel(),
    kM{1.0}, ratio{1.0}, alpham1{1.0}, betam1{1.0},
    kC{100000.0}, ODFmean{0.0}, ODFstdev{0.3},
    Dxmean{1.2}, Dxstdev{0.015}, Dxlb{1.0}, Dxub{1.24}, Dx_aniso{1.0},
    kI{0.0},
    F11_pushforward{1.0}, F12_pushforward{0.0},F21_pushforward{0.0},F22_pushforward{1.0},
    strainhistory{1.0, 0.0, 0.0, 1.0}
    {calc_new_old_fraction(0.01, 0.0);};

    model_ps_struc(double arg[]);

    int stress(double vF[4], double res[4]);

};


#endif /* ModelPSstruc_hpp */