forked from calculix/Cubit-CalculiX
-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathbolted_connection.jou
148 lines (124 loc) · 5.8 KB
/
bolted_connection.jou
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
#!cubit
reset
create brick x 40 y 20 z 5
volume 1 copy move x 20 z -5.001
create cylinder height 11 radius 4
move volume 3 x 10 z -2.5 include_merged
create cylinder height 4 radius 6
move volume 4 x 10 z 4.501 include_merged
create cylinder height 4 radius 6
move volume 5 x 10 z -9.502 include_merged
subtract volume 4 from volume 3 keep_tool
subtract volume 5 from volume 3 keep_tool
create cylinder height 20 radius 4.25
move volume 6 x 10 z -2.5 include_merged
subtract volume 6 from volume 1 2
imprint vol all
merge vol all
#vol all size 1
mesh vol all
block 1 add vol 1
block 2 add vol 2
block 3 add vol 3 to 5
#block all element type hex20
#nodeset 1 add surface 4
#nodeset 1 name "plate_1"
#nodeset 2 add surface 12
#nodeset 2 name "plate_2"
sideset 1 add surface 30
sideset 1 name "plate_1_top"
sideset 2 add surface 31
sideset 2 name "plate_1_bottom"
sideset 3 add surface 33
sideset 3 name "plate_2_top"
sideset 4 add surface 34
sideset 4 name "plate_2_bottom"
sideset 5 add surface 36
sideset 5 name "bolt_top"
sideset 6 add surface 38
sideset 6 name "bolt_bottom"
sideset 7 add surface 25
sideset 7 name "bolt_shaft"
sideset 8 add surface 29 32
sideset 8 name "plates_holes"
create material "steel" property_group "CalculiX-FEA"
modify material "steel" scalar_properties "CCX_ELASTIC_USE_CARD" 1
modify material "steel" scalar_properties "CCX_ELASTIC_ISO_USE_CARD" 1
modify material "steel" matrix_property "CCX_ELASTIC_ISO_MODULUS_VS_POISSON_VS_TEMPERATURE" 2.1e+5 0.3 0
modify material "steel" scalar_properties "CCX_PLASTIC_ISO_USE_CARD" 1
modify material "steel" scalar_properties "CCX_DENSITY_USE_CARD" 1
modify material "steel" matrix_property "CCX_DENSITY_DENSITY" 7.85e-06 0
modify material "steel" scalar_properties "CCX_EXPANSION_USE_CARD" 1
modify material "steel" scalar_properties "CCX_EXPANSION_ISO_USE_CARD" 1
modify material "steel" matrix_property "CCX_EXPANSION_ISO_A_TEMPERATURE" 1.2e-05 0
modify material "steel" scalar_properties "CCX_PLASTIC_USE_CARD" 1
modify material "steel" matrix_property "CCX_PLASTIC_ISO_YIELD_STRESS_VS_STRAIN_VS_TEMPERATURE" 235 0 0
ccx create section solid block all material 1
ccx create surfaceinteraction name "interaction" linear slopeK 1e+7 sigmaINF 1 c0 1e-3
ccx modify surfaceinteraction 1 friction mu 0.1 lambda 500
ccx create contactpair surfaceinteraction 1 surfacetosurface master 1 slave 5
ccx create contactpair surfaceinteraction 1 surfacetosurface master 2 slave 3
ccx create contactpair surfaceinteraction 1 surfacetosurface master 4 slave 6
ccx create contactpair surfaceinteraction 1 surfacetosurface master 7 slave 8
create temperature on volume all value 0
ccx create initialcondition temperature
ccx modify initialcondition 1 temperature bc_id 1
#preload
create temperature on volume 3 value -50
ccx modify temperature 2 keyword_type temp
ccx modify temperature 1 keyword_type temp
# plate bc
create displacement name "plate_1_fix" on surface 4 dof 1 dof 2 dof 3 fix 0
create displacement name "plate_2_fix" on surface 12 dof 1 dof 2 dof 3 fix 0
create displacement name "plate_1_pull" on surface 4 dof 1 fix -0.35
create displacement name "plate_2_pull" on surface 12 dof 1 fix 0.35
ccx create historyoutput name "ho_1" element
ccx modify historyoutput 1 element block 1 totals_yes
ccx modify historyoutput 1 element key_on S E
ccx modify historyoutput 1 element key_off SVF ME PEEQ CEEQ ENER SDV HFL HFLF COORD ELSE ELKE EVOL EMAS EBHE CENT
ccx create historyoutput name "ho_2" element
ccx modify historyoutput 2 element block 2
ccx modify historyoutput 2 element key_on S E
ccx modify historyoutput 2 element key_off SVF ME PEEQ CEEQ ENER SDV HFL HFLF COORD ELSE ELKE EVOL EMAS EBHE CENT
ccx create historyoutput name "ho_3" element
ccx modify historyoutput 3 element block 3
ccx modify historyoutput 3 element key_on S E
ccx modify historyoutput 3 element key_off SVF ME PEEQ CEEQ ENER SDV HFL HFLF COORD ELSE ELKE EVOL EMAS EBHE CENT
ccx create historyoutput name "ho_4" contact
ccx modify historyoutput 4 contact contactpair 1 totals_yes
ccx modify historyoutput 4 contact key_on CDIS CSTR CELS CNUM CF CFN CFS
ccx create historyoutput name "ho_5" node
ccx modify historyoutput 5 node nodeset 1 totals_yes
ccx modify historyoutput 5 node key_on U RF
ccx modify historyoutput 5 node key_off NT TSF TTF PN PSF PTF MACH CP VF DEPF TURB MF RFL
ccx create fieldoutput name "fo_1" node
ccx create fieldoutput name "fo_2" element
ccx create fieldoutput name "fo_3" contact
ccx modify fieldoutput 1 node
ccx modify fieldoutput 1 node key_on RF U NT
ccx modify fieldoutput 1 node key_off CP DEPF DEPT DTF HCRI KEQ MACH MAXU MF PNT POT PRF PS PSF PT PTF PU RFL SEN TS TSF TT TTF TURB V VF
ccx modify fieldoutput 2 element
ccx modify fieldoutput 2 element key_on E S
ccx modify fieldoutput 2 element key_off CEEQ ECD EMFB EMFE ENER ERR HER HFL HFLF MAXE MAXS ME PEEQ PHS SF SMID SNEG SPOS SVF SDV THE ZZS
ccx modify fieldoutput 3 contact
ccx modify fieldoutput 3 contact key_on CDIS CSTR CELS
ccx modify fieldoutput 3 contact key_off PCON
ccx create step name "preload" static
ccx modify step 1 parameter nlgeom_yes
ccx modify step 1 static totaltimeatstart 0 initialtimeincrement 0.5 timeperiodofstep 1 minimumtimeincrement 1e-5 maximumtimeincrement 0.5
ccx step 1 add fieldoutput 1 2 3
ccx step 1 add historyoutput 1 2 3 4 5
ccx step 1 add bc temperature 1 2
ccx step 1 add bc displacement 1 2
ccx create step name "load" static
ccx modify step 2 parameter nlgeom_yes
ccx modify step 2 static totaltimeatstart 1 initialtimeincrement 0.01 timeperiodofstep 1 minimumtimeincrement 1e-5 maximumtimeincrement 0.02
ccx step 2 add fieldoutput 1 2 3
ccx step 1 add historyoutput 1 2 3 4 5
ccx step 2 add bc displacement 3 4
ccx create job name "bolted_connection"
#ccx run job 1 no_conversion
#ccx wait job 1
#ccx result load job 1
#ccx use logfile core on
#ccx result convert job 1 partial block 1 2 3