v0.14.0
Loading...
Searching...
No Matches
pinched_cylinder_convetgence_test.sh
Go to the documentation of this file.
1#!/bin/sh
2
3NB_PROC=4
4
5BIN_PATH="/Users/likask/MyBuild/um_debug/solid_shell_prism_element"
6SRC_PATH="/Users/likask/MyBuild/um_debug/users_modules/solid_shell_prism_element"
7
8MPIRUN="/opt/petsc/arch-darwin-c-debug/bin/mpirun -np $NB_PROC"
9
10set -x
11cd $BIN_PATH
12
13# MESH 9
14
15# echo LOG: Mesh 9 - Convergence in plane - locking
16#
17# for ODISP in `seq 1 8`;
18# do
19# $MPIRUN $BIN_PATH/solid_shell \
20# -my_file $SRC_PATH/meshes/pinched_cylinder_9.h5m \
21# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
22# -my_order $ODISP -my_order_thickness 0 -my_order_3d 0 \
23# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
24# 2>&1 | tee tmp_log
25# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 0 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
26# done
27#
28# echo LOG: Mesh 9 - Convergence in plane - thinckes 2 - plane locking
29#
30# for ODISP in `seq 1 8`;
31# do
32# $MPIRUN $BIN_PATH/solid_shell \
33# -my_file $SRC_PATH/meshes/pinched_cylinder_9.h5m \
34# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
35# -my_order $ODISP -my_order_thickness 2 -my_order_3d 0 \
36# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
37# 2>&1 | tee tmp_log
38# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 2 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
39# done
40#
41# echo LOG: Mesh 9 - Convergence in plane - thinckes 2 - plane 3
42#
43# for ODISP in `seq 1 8`;
44# do
45# $MPIRUN $BIN_PATH/solid_shell \
46# -my_file $SRC_PATH/meshes/pinched_cylinder_9.h5m \
47# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
48# -my_order $ODISP -my_order_thickness 3 -my_order_3d 2 \
49# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
50# 2>&1 | tee tmp_log
51# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 2 2 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
52# done
53#
54# echo LOG: Mesh 6 - Convergence in plane - locking
55#
56for ODISP in `seq 5 5`;
57do
58 FILE=$SRC_PATH/meshes/pinched_cylinder_6.h5m
59 ./solid_shell_partition_mesh -my_file $FILE -my_out_file tmp.h5m -my_nb_parts $NB_PROC
60 $MPIRUN $BIN_PATH/solid_shell \
61 -my_file tmp.h5m -my_is_partitioned 1 \
62 -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
63 -my_order $ODISP -my_order_thickness 3 -my_order_3d 2 \
64 -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
65 2>&1 | tee tmp_log
66 echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 0 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
67done
68#
69# echo LOG: Mesh 6 - Convergence in plane - thinckes 2 - plane locking
70#
71# for ODISP in `seq 1 8`;
72# do
73# $MPIRUN $BIN_PATH/solid_shell \
74# -my_file $SRC_PATH/meshes/pinched_cylinder_6.h5m \
75# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
76# -my_order $ODISP -my_order_thickness 2 -my_order_3d 0 \
77# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
78# 2>&1 | tee tmp_log
79# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 2 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
80# done
81#
82# echo LOG: Mesh 6 - Convergence in plane - thinckes 3 - plane 2
83#
84# for ODISP in `seq 1 8`;
85# do
86# $MPIRUN $BIN_PATH/solid_shell \
87# -my_file $SRC_PATH/meshes/pinched_cylinder_6.h5m \
88# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
89# -my_order $ODISP -my_order_thickness 2 -my_order_3d 3 \
90# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
91# 2>&1 | tee tmp_log
92# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 3 2 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
93# done
94
95# Mesh 4
96
97# echo LOG: Mesh 4 - Convergence in plane - locking
98#
99# for ODISP in `seq 1 8`;
100# do
101# $MPIRUN $BIN_PATH/solid_shell \
102# -my_file $SRC_PATH/meshes/pinched_cylinder_4.h5m \
103# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
104# -my_order $ODISP -my_order_thickness 0 -my_order_3d 0 \
105# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
106# 2>&1 | tee tmp_log
107# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 0 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
108# done
109#
110# echo LOG: Mesh 4 - Convergence in plane - thinckes 2 - plane locking
111#
112# for ODISP in `seq 1 8`;
113# do
114# $MPIRUN $BIN_PATH/solid_shell \
115# -my_file $SRC_PATH/meshes/pinched_cylinder_4.h5m \
116# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
117# -my_order $ODISP -my_order_thickness 2 -my_order_3d 0 \
118# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
119# 2>&1 | tee tmp_log
120# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 2 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
121# done
122#
123# echo LOG: Mesh 4 - Convergence in plane - thinckes 3 - plane 2
124#
125# for ODISP in `seq 1 8`;
126# do
127# $MPIRUN $BIN_PATH/solid_shell \
128# -my_file $SRC_PATH/meshes/pinched_cylinder_4.h5m \
129# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
130# -my_order $ODISP -my_order_thickness 3 -my_order_3d 2 \
131# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
132# 2>&1 | tee tmp_log
133# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 3 2 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
134# done
135#
136# # Mesh 2
137#
138# echo LOG: Mesh 2 - Convergence in plane - locking
139#
140# for ODISP in `seq 1 6`;
141# do
142# $MPIRUN $BIN_PATH/solid_shell \
143# -my_file $SRC_PATH/meshes/pinched_cylinder_2.h5m \
144# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
145# -my_order $ODISP -my_order_thickness 0 -my_order_3d 0 \
146# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
147# 2>&1 | tee tmp_log
148# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 0 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
149# done
150#
151# echo LOG: Mesh 2 - Convergence in plane - thinckes 2 - plane locking
152#
153# for ODISP in `seq 1 6`;
154# do
155# $MPIRUN $BIN_PATH/solid_shell \
156# -my_file $SRC_PATH/meshes/pinched_cylinder_2.h5m \
157# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
158# -my_order $ODISP -my_order_thickness 2 -my_order_3d 0 \
159# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
160# 2>&1 | tee tmp_log
161# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 2 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
162# done
163#
164# echo LOG: Mesh 2 - Convergence in plane - thinckes 3 - plane 2
165#
166# for ODISP in `seq 1 6`;
167# do
168# $MPIRUN $BIN_PATH/solid_shell \
169# -my_file $SRC_PATH/meshes/pinched_cylinder_2.h5m \
170# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
171# -my_order $ODISP -my_order_thickness 3 -my_order_3d 2 \
172# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
173# 2>&1 | tee tmp_log
174# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 3 2 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
175# done
176
177# rm -f tmp_log
178
179#/opt/local/bin/valgrind --dsymutil=yes --track-origins=yes --trace-children=yes
180
181# $MPIRUN $BIN_PATH/solid_shell \
182# -my_file $SRC_PATH/meshes/pinched_cylinder_4.h5m \
183# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
184# -my_order 5 -my_order_thickness 2 -my_order_3d 0 \
185# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3e6 \
186# 2>&1 | tee tmp_log
187
188# $MPIRUN $BIN_PATH/solid_shell \
189# -my_file $SRC_PATH/meshes/simple_plate.h5m \
190# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
191# -my_order 2 -my_order_thickness 0 -my_order_3d 0 \
192# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3e6 \
193# 2>&1 | tee tmp_log
194
195# nonlinear analysis
196
197# Tension test
198# $MPIRUN $BIN_PATH/solid_shell_nonlinear \
199# -my_file $SRC_PATH/meshes/pinched_cylinder_4_long_tension.h5m \
200# -snes_monitor -snes_max_it 24 -snes_atol 1e-6 -snes_rtol 1e-5 -snes_stol 0 \
201# -snes_linesearch_type l2 -snes_linesearch_monitor -snes_linesearch_minlambda 0.01 \
202# -snes_converged_reason \
203# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
204# -ksp_converged_reason \
205# -my_step_size 1 -my_arc_beta 0 \
206# -my_numerical_damping 0 \
207# -my_order 6 -my_order_thickness 2 -my_order_3d 0 \
208# -my_shell_thickness 0.094 -my_poisson_ratio 0.3125 -my_young_modulus 10.5 \
209# -ts_monitor -ts_dt 0.001 -ts_final_time 100000 \
210# -ts_max_snes_failures -1 \
211# 2>&1 | tee tmp_log
212
213# Compresion test
214# $MPIRUN $BIN_PATH/solid_shell_nonlinear \
215# -my_file $SRC_PATH/meshes/pinched_cylinder_6_long_compression.h5m \
216# -snes_monitor -snes_max_it 24 -snes_atol 1e-6 -snes_rtol 1e-5 -snes_stol 0 \
217# -snes_linesearch_type l2 -snes_linesearch_monitor -snes_linesearch_minlambda 0.01 \
218# -snes_converged_reason \
219# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
220# -ksp_converged_reason \
221# -my_step_size 1 -my_arc_beta 0 \
222# -my_numerical_damping 0 \
223# -my_order 5 -my_order_thickness 2 -my_order_3d 0 \
224# -my_shell_thickness 1 -my_poisson_ratio 0.3 -my_young_modulus 30 \
225# -ts_monitor -ts_dt 0.001 -ts_final_time 100000 \
226# -ts_max_snes_failures -1 \
227# 2>&1 | tee tmp_log