-
Notifications
You must be signed in to change notification settings - Fork 0
/
EffectiveModel.hpp
127 lines (86 loc) · 2.15 KB
/
EffectiveModel.hpp
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
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
//
// EffectiveModel.hpp
// CompleteModelCollection
//
// Created by Will Zhang on 9/26/17.
// Copyright © 2017 Will Zhang. All rights reserved.
//
#ifndef EffectiveModel_hpp
#define EffectiveModel_hpp
#include <stdio.h>
#include <vector>
#include <cmath>
#endif /* EffectiveModel_hpp */
class smpfl_model
{
private:
int ex_i;
int ex_j;
int ex_k;
int deg[3];
double vec_M[2];
double vec_S[2];
double gamma_1;
double gamma_2;
double gamma_3;
double energy;
public:
void degrees(int degrees);
void calc_vec_M(double phi);
void calc_vec_S(double phi);
void calc_lamdbda_M2();
void calc_lamdbda_S2();
void calc_phi();
void calc_gammas();
void stress(double parm[], double strain[], int n);
double partial_i(double w, double gamma1, double i);
smpfl_model(double phi, double rank, double degrees);
smpfl_model(): ex_i{2}, ex_j{2}, ex_k{2}, deg{10,10,5}, vec_M{1,0}, vec_S{0,1} {};
};
void smpfl_model::degrees(int degrees)
{
deg[0] = degrees;
deg[1] = degrees;
deg[2] = degrees/2;
return;
}
smpfl_model::smpfl_model(double theta, int rank, int exponent)
{
degrees(rank);
ex_i = exponent/deg[0]/deg[1];
ex_j = (exponent % (deg[0]/deg[1])) / deg[2];
ex_k = (exponent % (deg[0]/deg[1])) % deg[2];
calc_vec_M(theta);
calc_vec_S(theta);
return;
}
void smpfl_model::calc_vec_M(double theta)
{
vec_M[0] = cos(theta);
vec_M[1] = sin(theta);
return;
}
void smpfl_model::calc_vec_S(double theta)
{
vec_S[0] = -sin(theta);
vec_S[1] = cos(theta);
return;
}
double mat_multiply(double tsr[4], double vec_M[2], double vec_S[2])
{
double result = 0.0;
for (int i = 0; i<2; i++) {
for (int j = 0; j<2; j++) {
result += vec_M[i]*tsr[2*i + j]*vec_S[j];
}
}
return result;
}
double strainenergy(double gamma1, double gamma2, double gamma3, int i, int j, int k)
{
return pow(gamma1, i) * pow(gamma2, j) * pow(gamma3, k);
}
double smpfl_model::partial_i(double w, double gamma1, double i)
{
return (i - 1) * w / gamma1;
}