-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathinput
More file actions
185 lines (162 loc) · 4.92 KB
/
input
File metadata and controls
185 lines (162 loc) · 4.92 KB
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
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
#===============
#=====================================
# Materials types, properties, phases
#---------------------------------------------------
# Number of phases
N_phases 2
# System size
# CellDim 128 128 128
CellDim 256 256 1
# Data file of initial microstructure
#initial_ms /public/home/pzhao/work/VPSC/evpfft/msDataBase/AlSingleXtal_8x8x8.ms
#initial_ms /public/home/pzhao/work/VPSC/evpfft/msDataBase/poly_4N_Cu.ms
#initial_ms /public/home/pzhao/work/VPSC/evpfft/msDataBase/4N_Cu_many_32.ms
initial_ms Parent_Twin_Hori.txt
# Type of phase-1 and phase-2 (in accordance with the following order)
# 0(xtal) 1(gas)
Type_phases 0 0
# Elastic constants of phase 1: c11, c12, c44, nu
# Copper, in MPa 473K
ElastConst_PhaseI 162700. 119000. 71150. 0.35
# Copper, in MPa 723K
#ElastConst_PhaseI 152500. 114000. 64360. 0.35
# The above is not actually used because we use T-dependent elastic stiffness
# Data file containing info of phase 1
Slip_PhaseI Cu_p_hard.sx
# Elastic constants of phase 1: c11, c12, c44, nu
# Gas, ignored!
ElastConst_PhaseII 80000.0 55000.0 35150.0 0.35
# Data file containing info of phase 2
Slip_PhaseII Cu_p.sx
#====================================================
# Boundary conditions
#----------------------------------------------------ss
# Flags for velocity gradient (0:unknown, 1:known) Voigt notation
VelGrad_BC_Flag 1 1 1 1 1 1
# Velocity gradient
VelGrad_BC 0.0E-3 0.0E-3 0.0E-3 0.0 0.0 1.0E-3
# Flags for Cauchy stress
Stress_BC_Flag 0 0 0 0 0 0
# Applied Cauchy stress
Stress_BC 0. 0. 0. 0. 0. 0.
# elasticity BC: 0 -(e_homo=0), 1 -(relaxed bc)
ElastBC 1
#====================================================
# Simulation setup
#----------------------------------------------------
# Time increment
TimeStep 0.1
# Total simulation steps
N_steps 201
# Error tolerance
Err 2E-5
# Maximum iteration steps allowed to solve micromechanical problem
IterMax 1000
# write fields?(0:No, 1:Yes) write step interval
PrintControl 1 50
# Update texture flag
Update_Flag 1
# Hardening flag
Hard_Flag 1
# Texture output flag
Tex_Flag 1
#====================================================
# Dislocation density based constitutive model
#----------------------------------------------------
# A. Ma et al. Acta Mater. 54(2006) 2169
# lattice parameters for each phase [nm]
a0 0.362 0.362
# voxel length used in GND [mm]
#L0 0.03
L0 0.0045
# Used to scale L0 to saturate GND
#GND_ScaleRatio 1000.0
GND_ScaleRatio 200.0
# Activation energy for slip [1E-19J]
Q_slip 2.0 2.0
# Activation energy for climb [1E-19J]
Q_bulk 1.6 1.6
# Initial SSD density [1E12/m^2]
rho_SSD_initial 0.1
# Temperature [Ki
T__K 300
# Constant for passing stress
C_1 1.0 1.0
# Constant for jump width
C_2 1.0 1.0
# Constant for obstacle width
C_3 1.0 1.0
# Constant for lock forming rate [1E6/m]
C_4 600.0 600.0
## Constant for 6thermal annihilation rate
C_5 17.6 17.6
## Constant for thermal annihilation rate [1/nm]
C_6 0E0 0.0
## Constant for dipole forming rate [1E-37m^5*s^C_8/s]
C_7 0.0E-9 0.0
#C_7 0.0 0.0
# Constant for nonlinear climb of edge dislocations
C_8 0.24 0.24
# Self interaction
Selfinter0 0.2 0.2
# coplanar interaction
Coplanar0 0.2 0.2
# cross slip
CrossSlip0 0.2 0.2
#CrossSlip 0.8 0.0
# glissile junction
GlissileJunction0 0.3 0.3
#GlissileJunction 0.8 0
# Hirth loc
HirthLock0 0.2 0.2
#HirthLock 0.8 0.0
# Lomer-Cottrell Lock
LomerCottrellLock0 0.4 0.4
#LomerCottrellLock 0.8 0.0
#Diffusion coefficient for climb equation
D_0 1.5E-3
#nucleus_radius
nucleus_radius 4
#====================================================
# Phase field of recrystallization
#----------------------------------------------------
# System size for Phase-field (not need to be equal to FFDynamic recrystallization of copper polycrystals with different puritiesT size)
#CellDim_pf 64 64 64
CellDim_pf 128 128 128
# Seed f32 GSL random number generator
RandSeed 123
# GB energy [J/m^2]
E_gb 0.625
# GB mob2lity [m^4/(MJ*s)]
M_gb 7.4E-11
# Burgers vector [nm]
bb_len 0.256
# Scaler of stored energy [1]
Scale_Fdeform 0.5
# Characteristic nucleation strength (disl. density difference) for DRX [1E12/m^2]
k_c 320000.0
# Expon0nt in the statistical model for DRX nucleation [1]
alpha 8.5
# Norm3lized GB mobility in phase-field [in grid unit]
# inversely proportional to the # of phase-field steps
#M_bar 12.0
# Coefficient to be multiplied to k_c for re-nucleation
zeta_kc 1.0
# Aging time for DRX re-nucleation
age_drx 1E0
# Scaler used to be multiplied to SSDs for DRX grains after pahse-field
ssdScaler_DRX 0.74
gndScaler_DRX 0.66
sigScaler_DRX 0.66
# Static or dynamic implementation of nucleation? (1: static, 0: dynamic)
Nucl_Static_Flag 1
# Characteristic aging time for newly-formed DRX grains, normalized by TimeStep
tau_DRX 1000.0
# Additional barrier for slip for new DRX grains [J/m^2]
DeltaQ_DRX 0.0
# Additional passing stress for new DRX grains
DeltaTau_DRX 0.0
#pf_length_scale
pf_length_scale 2.5E-06
#pf_time_step
pf_time_step 0.004