Merge pull request #213 from DongWuTUM/hotfix/thin_shell_model_

Xiangyu-Hu · web-flow · commit 5188552fa16e · 2023-02-28T14:51:20.000+01:00
Bug fix for thin shell.
diff --git a/SPHINXsys/src/shared/particle_dynamics/solid_dynamics/elastic_dynamics.cpp b/SPHINXsys/src/shared/particle_dynamics/solid_dynamics/elastic_dynamics.cpp
@@ -12,15 +12,15 @@ namespace SPH
 		AcousticTimeStepSize::AcousticTimeStepSize(SPHBody &sph_body, Real CFL)
 			: LocalDynamicsReduce<Real, ReduceMin>(sph_body, Real(MaxRealNumber)),
 			  ElasticSolidDataSimple(sph_body), CFL_(CFL),
-			  vel_(particles_->vel_), acc_(particles_->acc_),
+			  vel_(particles_->vel_), acc_(particles_->acc_), acc_prior_(particles_->acc_prior_),
 			  smoothing_length_(sph_body.sph_adaptation_->ReferenceSmoothingLength()),
 			  c0_(particles_->elastic_solid_.ReferenceSoundSpeed()) {}
 		//=================================================================================================//
 		Real AcousticTimeStepSize::reduce(size_t index_i, Real dt)
 		{
 			// since the particle does not change its configuration in pressure relaxation step
 			// I chose a time-step size according to Eulerian method
-			return CFL_ * SMIN(sqrt(smoothing_length_ / (acc_[index_i].norm() + TinyReal)),
+			return CFL_ * SMIN(sqrt(smoothing_length_ / ((acc_[index_i] + acc_prior_[index_i]).norm() + TinyReal)),
 							   smoothing_length_ / (c0_ + vel_[index_i].norm()));
 		}
 		//=================================================================================================//
diff --git a/SPHINXsys/src/shared/particle_dynamics/solid_dynamics/elastic_dynamics.h b/SPHINXsys/src/shared/particle_dynamics/solid_dynamics/elastic_dynamics.h
@@ -90,7 +90,7 @@ namespace SPH
 		{
 		protected:
 			Real CFL_;
-			StdLargeVec<Vecd> &vel_, &acc_;
+			StdLargeVec<Vecd> &vel_, &acc_, &acc_prior_;
 			Real smoothing_length_, c0_;
 
 		public:
diff --git a/SPHINXsys/src/shared/particle_dynamics/solid_dynamics/thin_structure_dynamics.cpp b/SPHINXsys/src/shared/particle_dynamics/solid_dynamics/thin_structure_dynamics.cpp
@@ -12,10 +12,11 @@ namespace SPH
 			  n0_(particles_->n0_), n_(particles_->n_), pseudo_n_(particles_->pseudo_n_),
 			  pos0_(particles_->pos0_), transformation_matrix_(particles_->transformation_matrix_) {}
 		//=================================================================================================//
-		ShellAcousticTimeStepSize::ShellAcousticTimeStepSize(SPHBody &sph_body, Real CFL)
+		ShellAcousticTimeStepSize::ShellAcousticTimeStepSize(SPHBody& sph_body, Real CFL)
 			: LocalDynamicsReduce<Real, ReduceMin>(sph_body, Real(MaxRealNumber)),
-			  ShellDataSimple(sph_body), CFL_(CFL), vel_(particles_->vel_), acc_(particles_->acc_),
+			ShellDataSimple(sph_body), CFL_(CFL), vel_(particles_->vel_), acc_(particles_->acc_),
 			  angular_vel_(particles_->angular_vel_), dangular_vel_dt_(particles_->dangular_vel_dt_),
+			  acc_prior_(particles_->acc_prior_),
 			  thickness_(particles_->thickness_),
 			  rho0_(particles_->elastic_solid_.ReferenceDensity()),
 			  E0_(particles_->elastic_solid_.YoungsModulus()),
@@ -27,14 +28,14 @@ namespace SPH
 		{
 			// Since the particle does not change its configuration in pressure relaxation step,
 			// I chose a time-step size according to Eulerian method.
-			Real time_setp_0 = CFL_ * SMIN(sqrt(smoothing_length_ / (acc_[index_i].norm() + TinyReal)),
+			Real time_setp_0 = SMIN(sqrt(smoothing_length_ / ((acc_[index_i] + acc_prior_[index_i]).norm() + TinyReal)),
 										   smoothing_length_ / (c0_ + vel_[index_i].norm()));
-			Real time_setp_1 = CFL_ * SMIN(sqrt(1.0 / (dangular_vel_dt_[index_i].norm() + TinyReal)),
+			Real time_setp_1 = SMIN(sqrt(1.0 / (dangular_vel_dt_[index_i].norm() + TinyReal)),
 										   1.0 / (angular_vel_[index_i].norm() + TinyReal));
 			Real time_setp_2 = smoothing_length_ * sqrt(rho0_ * (1.0 - nu_ * nu_) / E0_ /
 														(2.0 + (Pi * Pi / 12.0) * (1.0 - nu_) *
 																   (1.0 + 1.5 * powerN(smoothing_length_ / thickness_[index_i], 2))));
-			return SMIN(time_setp_0, time_setp_1, time_setp_2);
+			return CFL_ * SMIN(time_setp_0, time_setp_1, time_setp_2);
 		}
 		//=================================================================================================//
 		ShellCorrectConfiguration::
@@ -179,12 +180,17 @@ namespace SPH
 					transformation_matrix_[index_i] * current_transformation_matrix.transpose();
 
 				/** correct Almansi strain tensor according to plane stress problem. */
-				current_local_almansi_strain = getCorrectedAlmansiStrain(current_local_almansi_strain, nu_);
+				current_local_almansi_strain = getCorrectedAlmansiStrain(current_local_almansi_strain, nu_); 
+
+				/** correct out-plane numerical damping. */
+				Matd numerical_damping = current_transformation_matrix * transformation_matrix_[index_i].transpose() * F_gaussian_point *
+					elastic_solid_.NumericalDampingRightCauchy(F_gaussian_point, dF_gaussian_point_dt, smoothing_length_, index_i)
+					* F_gaussian_point.transpose() * transformation_matrix_[index_i] * current_transformation_matrix.transpose() 
+					/ F_gaussian_point.determinant();
+				numerical_damping.col(Dimensions - 1) *= thickness_[index_i] / smoothing_length_;
 				
-				Matd cauchy_stress = elastic_solid_.StressCauchy(current_local_almansi_strain, F_gaussian_point, index_i) + 
-					current_transformation_matrix * transformation_matrix_[index_i].transpose() * F_gaussian_point * 
-					elastic_solid_.NumericalDampingRightCauchy(F_gaussian_point, dF_gaussian_point_dt, smoothing_length_, index_i) 
-					* F_gaussian_point.transpose() * transformation_matrix_[index_i] * current_transformation_matrix.transpose() / F_gaussian_point.determinant();
+				Matd cauchy_stress = elastic_solid_.StressCauchy(current_local_almansi_strain, F_gaussian_point, index_i) + numerical_damping;
+					
 
 				/** Impose modeling assumptions. */
 				cauchy_stress.col(Dimensions - 1) *= shear_correction_factor_;
@@ -247,8 +253,8 @@ namespace SPH
 						pseudo_n_[index_j] - n0_[index_j],
 						transformation_matrix_[index_j].transpose() * F_bending_[index_j] * transformation_matrix_[index_j]);
 					Vecd rotation_jump = getRotationJump(pseudo_n_jump, transformation_matrix_[index_i]);
-					pseudo_normal_acceleration += hourglass_control_factor_ / 3.0 * weight * Dimensions * r_ij * G0_ * rotation_jump * 
-												  inner_neighborhood.dW_ijV_j_[n] * thickness_[index_i];
+					pseudo_normal_acceleration += hourglass_control_factor_ * weight * G0_ * rotation_jump * dim_inv_r_ij *
+												  inner_neighborhood.dW_ijV_j_[n] * powerN(thickness_[index_i], 3);
 				}
 
 				acceleration += (global_stress_i + global_stress_[index_j]) * inner_neighborhood.dW_ijV_j_[n] * inner_neighborhood.e_ij_[n];
diff --git a/SPHINXsys/src/shared/particle_dynamics/solid_dynamics/thin_structure_dynamics.h b/SPHINXsys/src/shared/particle_dynamics/solid_dynamics/thin_structure_dynamics.h
@@ -69,7 +69,7 @@ namespace SPH
 		{
 		protected:
 			Real CFL_;
-			StdLargeVec<Vecd> &vel_, &acc_, &angular_vel_, &dangular_vel_dt_;
+			StdLargeVec<Vecd> &vel_, &acc_, &angular_vel_, &dangular_vel_dt_, &acc_prior_;
 			StdLargeVec<Real> &thickness_;
 			Real rho0_, E0_, nu_, c0_;
 			Real smoothing_length_;
diff --git a/tests/2d_examples/test_2d_plate/2d_plate.cpp b/tests/2d_examples/test_2d_plate/2d_plate.cpp
@@ -83,6 +83,7 @@ int main()
 	//	Build up -- a SPHSystem
 	//----------------------------------------------------------------------
 	SPHSystem system(system_domain_bounds, resolution_ref);
+	system.generate_regression_data_ = false;
 	//----------------------------------------------------------------------
 	//	Creating body, materials and particles.
 	//----------------------------------------------------------------------
@@ -187,7 +188,13 @@ int main()
 	tt = t4 - t1 - interval;
 	std::cout << "Total wall time for computation: " << tt.seconds() << " seconds." << std::endl;
 
-	write_plate_max_displacement.newResultTest();
-
+	if (system.generate_regression_data_)
+	{
+		write_plate_max_displacement.generateDataBase(0.005);
+	}
+	else
+	{
+		write_plate_max_displacement.newResultTest();
+	}
 	return 0;
 }
diff --git a/tests/2d_examples/test_2d_plate/regression_test_tool/PlateObserver_Position_Run_0_result.xml b/tests/2d_examples/test_2d_plate/regression_test_tool/PlateObserver_Position_Run_0_result.xml
@@ -1,9 +1,9 @@
 <?xml version="1.0" encoding="UTF-8" ?>
 <result>
     <Snapshot_Element>
-        <Snapshot number_of_snapshot_for_local_result_="100" />
+        <Snapshot number_of_snapshot_for_local_result_="101" />
     </Snapshot_Element>
     <Result_Element>
-        <Particle_0 snapshot_0="~[5,0]" snapshot_1="~[5,0.0347921]" snapshot_2="~[5,0.165463]" snapshot_3="~[4.99999,0.481587]" snapshot_4="~[4.99994,0.928834]" snapshot_5="~[4.99995,1.40539]" snapshot_6="~[5,1.67986]" snapshot_7="~[4.99998,1.63683]" snapshot_8="~[5.00013,1.41895]" snapshot_9="~[4.99997,1.25354]" snapshot_10="~[5,1.25198]" snapshot_11="~[4.99997,1.38671]" snapshot_12="~[4.99997,1.53638]" snapshot_13="~[5,1.56036]" snapshot_14="~[5.00007,1.46893]" snapshot_15="~[4.99999,1.36053]" snapshot_16="~[5.00004,1.33213]" snapshot_17="~[4.99994,1.38557]" snapshot_18="~[4.99999,1.45912]" snapshot_19="~[5.00001,1.49221]" snapshot_20="~[4.99998,1.46565]" snapshot_21="~[5.00004,1.41428]" snapshot_22="~[4.99997,1.38418]" snapshot_23="~[5.00002,1.39711]" snapshot_24="~[4.99998,1.43441]" snapshot_25="~[4.99999,1.45928]" snapshot_26="~[5.00002,1.45457]" snapshot_27="~[4.99998,1.4302]" snapshot_28="~[5.00002,1.41028]" snapshot_29="~[4.99999,1.41064]" snapshot_30="~[4.99999,1.427]" snapshot_31="~[5,1.44195]" snapshot_32="~[5,1.44353]" snapshot_33="~[5.00001,1.43368]" snapshot_34="~[4.99999,1.42226]" snapshot_35="~[5.00001,1.41986]" snapshot_36="~[4.99999,1.42658]" snapshot_37="~[5,1.43515]" snapshot_38="~[5,1.43773]" snapshot_39="~[5,1.43343]" snapshot_40="~[5,1.42738]" snapshot_41="~[5,1.42484]" snapshot_42="~[5,1.42721]" snapshot_43="~[5,1.43147]" snapshot_44="~[5,1.43355]" snapshot_45="~[5,1.4323]" snapshot_46="~[5,1.42947]" snapshot_47="~[5,1.42772]" snapshot_48="~[5,1.42833]" snapshot_49="~[5,1.43033]" snapshot_50="~[5,1.4316]" snapshot_51="~[5,1.43135]" snapshot_52="~[5,1.43013]" snapshot_53="~[5,1.42918]" snapshot_54="~[5,1.42925]" snapshot_55="~[5,1.43004]" snapshot_56="~[5,1.43075]" snapshot_57="~[5,1.43082]" snapshot_58="~[5,1.43031]" snapshot_59="~[5,1.42979]" snapshot_60="~[5,1.42967]" snapshot_61="~[5,1.42997]" snapshot_62="~[5,1.43037]" snapshot_63="~[5,1.4305]" snapshot_64="~[5,1.43031]" snapshot_65="~[5,1.43]" snapshot_66="~[5,1.42986]" snapshot_67="~[5,1.42998]" snapshot_68="~[5,1.43021]" snapshot_69="~[5,1.43033]" snapshot_70="~[5,1.43027]" snapshot_71="~[5,1.43012]" snapshot_72="~[5,1.43002]" snapshot_73="~[5,1.43004]" snapshot_74="~[5,1.43015]" snapshot_75="~[5,1.43023]" snapshot_76="~[5,1.43023]" snapshot_77="~[5,1.43015]" snapshot_78="~[5,1.43009]" snapshot_79="~[5,1.43008]" snapshot_80="~[5,1.43013]" snapshot_81="~[5,1.43018]" snapshot_82="~[5,1.43019]" snapshot_83="~[5,1.43016]" snapshot_84="~[5,1.43012]" snapshot_85="~[5,1.43011]" snapshot_86="~[5,1.43013]" snapshot_87="~[5,1.43015]" snapshot_88="~[5,1.43016]" snapshot_89="~[5,1.43015]" snapshot_90="~[5,1.43014]" snapshot_91="~[5,1.43013]" snapshot_92="~[5,1.43013]" snapshot_93="~[5,1.43015]" snapshot_94="~[5,1.43015]" snapshot_95="~[5,1.43015]" snapshot_96="~[5,1.43014]" snapshot_97="~[5,1.43014]" snapshot_98="~[5,1.43014]" snapshot_99="~[5,1.43014]" />
+        <Particle_0 snapshot_0="~[5,0]" snapshot_1="~[5,0.0290544]" snapshot_2="~[5,0.158556]" snapshot_3="~[4.99998,0.410054]" snapshot_4="~[4.99994,0.792584]" snapshot_5="~[4.9999,1.2289]" snapshot_6="~[4.99996,1.5685]" snapshot_7="~[4.99998,1.67171]" snapshot_8="~[5.00009,1.53014]" snapshot_9="~[5.00005,1.34726]" snapshot_10="~[5,1.27724]" snapshot_11="~[4.99996,1.34091]" snapshot_12="~[5.00003,1.46365]" snapshot_13="~[4.99991,1.5287]" snapshot_14="~[5.00007,1.49678]" snapshot_15="~[5.00002,1.42217]" snapshot_16="~[4.99998,1.37387]" snapshot_17="~[5.00003,1.38203]" snapshot_18="~[4.99996,1.42619]" snapshot_19="~[5,1.46231]" snapshot_20="~[5,1.46138]" snapshot_21="~[5.00001,1.43603]" snapshot_22="~[5.00001,1.41335]" snapshot_23="~[4.99999,1.40991]" snapshot_24="~[5,1.4234]" snapshot_25="~[4.99999,1.4387]" snapshot_26="~[5,1.44271]" snapshot_27="~[5,1.43552]" snapshot_28="~[5,1.42639]" snapshot_29="~[5,1.42224]" snapshot_30="~[5,1.42513]" snapshot_31="~[5,1.43098]" snapshot_32="~[5,1.43421]" snapshot_33="~[5,1.43307]" snapshot_34="~[5,1.42953]" snapshot_35="~[5,1.42695]" snapshot_36="~[5,1.42729]" snapshot_37="~[5,1.42946]" snapshot_38="~[5,1.4314]" snapshot_39="~[5,1.43155]" snapshot_40="~[5,1.43023]" snapshot_41="~[5,1.42886]" snapshot_42="~[5,1.42853]" snapshot_43="~[5,1.42926]" snapshot_44="~[5,1.4302]" snapshot_45="~[5,1.43053]" snapshot_46="~[5,1.43016]" snapshot_47="~[5,1.42956]" snapshot_48="~[5,1.4293]" snapshot_49="~[5,1.42948]" snapshot_50="~[5,1.42983]" snapshot_51="~[5,1.43002]" snapshot_52="~[5,1.42995]" snapshot_53="~[5,1.42975]" snapshot_54="~[5,1.42962]" snapshot_55="~[5,1.42964]" snapshot_56="~[5,1.42975]" snapshot_57="~[5,1.42985]" snapshot_58="~[5,1.42986]" snapshot_59="~[5,1.42979]" snapshot_60="~[5,1.42972]" snapshot_61="~[5,1.4297]" snapshot_62="~[5,1.42974]" snapshot_63="~[5,1.42979]" snapshot_64="~[5,1.42981]" snapshot_65="~[5,1.42979]" snapshot_66="~[5,1.42976]" snapshot_67="~[5,1.42974]" snapshot_68="~[5,1.42975]" snapshot_69="~[5,1.42977]" snapshot_70="~[5,1.42978]" snapshot_71="~[5,1.42978]" snapshot_72="~[5,1.42976]" snapshot_73="~[5,1.42976]" snapshot_74="~[5,1.42976]" snapshot_75="~[5,1.42976]" snapshot_76="~[5,1.42977]" snapshot_77="~[5,1.42977]" snapshot_78="~[5,1.42977]" snapshot_79="~[5,1.42976]" snapshot_80="~[5,1.42976]" snapshot_81="~[5,1.42976]" snapshot_82="~[5,1.42977]" snapshot_83="~[5,1.42977]" snapshot_84="~[5,1.42977]" snapshot_85="~[5,1.42976]" snapshot_86="~[5,1.42976]" snapshot_87="~[5,1.42976]" snapshot_88="~[5,1.42977]" snapshot_89="~[5,1.42977]" snapshot_90="~[5,1.42977]" snapshot_91="~[5,1.42977]" snapshot_92="~[5,1.42976]" snapshot_93="~[5,1.42976]" snapshot_94="~[5,1.42977]" snapshot_95="~[5,1.42977]" snapshot_96="~[5,1.42977]" snapshot_97="~[5,1.42977]" snapshot_98="~[5,1.42977]" snapshot_99="~[5,1.42977]" snapshot_100="~[5,1.42977]" />
     </Result_Element>
 </result>
diff --git a/tests/2d_examples/test_2d_plate/regression_test_tool/PlateObserver_Position_Run_3_result.xml b/tests/2d_examples/test_2d_plate/regression_test_tool/PlateObserver_Position_Run_3_result.xml
@@ -1,9 +1,9 @@
 <?xml version="1.0" encoding="UTF-8" ?>
 <result>
     <Snapshot_Element>
-        <Snapshot number_of_snapshot_for_local_result_="100" />
+        <Snapshot number_of_snapshot_for_local_result_="101" />
     </Snapshot_Element>
     <Result_Element>
-        <Particle_0 snapshot_0="~[5,0]" snapshot_1="~[5,0.0347921]" snapshot_2="~[5,0.165463]" snapshot_3="~[4.99999,0.481587]" snapshot_4="~[4.99994,0.928834]" snapshot_5="~[4.99995,1.40539]" snapshot_6="~[5,1.67986]" snapshot_7="~[4.99998,1.63683]" snapshot_8="~[5.00013,1.41895]" snapshot_9="~[4.99997,1.25354]" snapshot_10="~[5,1.25198]" snapshot_11="~[4.99997,1.38671]" snapshot_12="~[4.99997,1.53638]" snapshot_13="~[5,1.56036]" snapshot_14="~[5.00007,1.46893]" snapshot_15="~[4.99999,1.36053]" snapshot_16="~[5.00004,1.33213]" snapshot_17="~[4.99994,1.38557]" snapshot_18="~[4.99999,1.45912]" snapshot_19="~[5.00001,1.49221]" snapshot_20="~[4.99998,1.46565]" snapshot_21="~[5.00004,1.41428]" snapshot_22="~[4.99997,1.38418]" snapshot_23="~[5.00002,1.39711]" snapshot_24="~[4.99998,1.43441]" snapshot_25="~[4.99999,1.45928]" snapshot_26="~[5.00002,1.45457]" snapshot_27="~[4.99998,1.4302]" snapshot_28="~[5.00002,1.41028]" snapshot_29="~[4.99999,1.41064]" snapshot_30="~[4.99999,1.427]" snapshot_31="~[5,1.44195]" snapshot_32="~[5,1.44353]" snapshot_33="~[5.00001,1.43368]" snapshot_34="~[4.99999,1.42226]" snapshot_35="~[5.00001,1.41986]" snapshot_36="~[4.99999,1.42658]" snapshot_37="~[5,1.43515]" snapshot_38="~[5,1.43773]" snapshot_39="~[5,1.43343]" snapshot_40="~[5,1.42738]" snapshot_41="~[5,1.42484]" snapshot_42="~[5,1.42721]" snapshot_43="~[5,1.43147]" snapshot_44="~[5,1.43355]" snapshot_45="~[5,1.4323]" snapshot_46="~[5,1.42947]" snapshot_47="~[5,1.42772]" snapshot_48="~[5,1.42833]" snapshot_49="~[5,1.43033]" snapshot_50="~[5,1.4316]" snapshot_51="~[5,1.43135]" snapshot_52="~[5,1.43013]" snapshot_53="~[5,1.42918]" snapshot_54="~[5,1.42925]" snapshot_55="~[5,1.43004]" snapshot_56="~[5,1.43075]" snapshot_57="~[5,1.43082]" snapshot_58="~[5,1.43031]" snapshot_59="~[5,1.42979]" snapshot_60="~[5,1.42967]" snapshot_61="~[5,1.42997]" snapshot_62="~[5,1.43037]" snapshot_63="~[5,1.4305]" snapshot_64="~[5,1.43031]" snapshot_65="~[5,1.43]" snapshot_66="~[5,1.42986]" snapshot_67="~[5,1.42998]" snapshot_68="~[5,1.43021]" snapshot_69="~[5,1.43033]" snapshot_70="~[5,1.43027]" snapshot_71="~[5,1.43012]" snapshot_72="~[5,1.43002]" snapshot_73="~[5,1.43004]" snapshot_74="~[5,1.43015]" snapshot_75="~[5,1.43023]" snapshot_76="~[5,1.43023]" snapshot_77="~[5,1.43015]" snapshot_78="~[5,1.43009]" snapshot_79="~[5,1.43008]" snapshot_80="~[5,1.43013]" snapshot_81="~[5,1.43018]" snapshot_82="~[5,1.43019]" snapshot_83="~[5,1.43016]" snapshot_84="~[5,1.43012]" snapshot_85="~[5,1.43011]" snapshot_86="~[5,1.43013]" snapshot_87="~[5,1.43015]" snapshot_88="~[5,1.43016]" snapshot_89="~[5,1.43015]" snapshot_90="~[5,1.43014]" snapshot_91="~[5,1.43013]" snapshot_92="~[5,1.43013]" snapshot_93="~[5,1.43015]" snapshot_94="~[5,1.43015]" snapshot_95="~[5,1.43015]" snapshot_96="~[5,1.43014]" snapshot_97="~[5,1.43014]" snapshot_98="~[5,1.43014]" snapshot_99="~[5,1.43014]" />
+        <Particle_0 snapshot_0="~[5,0]" snapshot_1="~[5,0.0290765]" snapshot_2="~[5,0.158627]" snapshot_3="~[4.99998,0.410192]" snapshot_4="~[4.99994,0.792871]" snapshot_5="~[4.9999,1.22918]" snapshot_6="~[4.99996,1.56871]" snapshot_7="~[4.99999,1.67175]" snapshot_8="~[5.00008,1.52999]" snapshot_9="~[5.00005,1.34705]" snapshot_10="~[5,1.2771]" snapshot_11="~[4.99996,1.34093]" snapshot_12="~[5.00003,1.4638]" snapshot_13="~[4.99991,1.52886]" snapshot_14="~[5.00007,1.49682]" snapshot_15="~[5.00002,1.42208]" snapshot_16="~[4.99998,1.37373]" snapshot_17="~[5.00003,1.38197]" snapshot_18="~[4.99996,1.42623]" snapshot_19="~[5.00001,1.46242]" snapshot_20="~[5,1.46145]" snapshot_21="~[5.00001,1.43602]" snapshot_22="~[5.00001,1.41327]" snapshot_23="~[4.99999,1.40984]" snapshot_24="~[5,1.42339]" snapshot_25="~[4.99999,1.43875]" snapshot_26="~[5,1.44276]" snapshot_27="~[5,1.43553]" snapshot_28="~[5,1.42636]" snapshot_29="~[5,1.4222]" snapshot_30="~[5,1.42512]" snapshot_31="~[5,1.43099]" snapshot_32="~[5,1.43424]" snapshot_33="~[5,1.43309]" snapshot_34="~[5,1.42952]" snapshot_35="~[5,1.42693]" snapshot_36="~[5,1.42727]" snapshot_37="~[5,1.42947]" snapshot_38="~[5,1.43141]" snapshot_39="~[5,1.43156]" snapshot_40="~[5,1.43023]" snapshot_41="~[5,1.42885]" snapshot_42="~[5,1.42852]" snapshot_43="~[5,1.42926]" snapshot_44="~[5,1.4302]" snapshot_45="~[5,1.43053]" snapshot_46="~[5,1.43016]" snapshot_47="~[5,1.42956]" snapshot_48="~[5,1.42929]" snapshot_49="~[5,1.42947]" snapshot_50="~[5,1.42983]" snapshot_51="~[5,1.43002]" snapshot_52="~[5,1.42995]" snapshot_53="~[5,1.42975]" snapshot_54="~[5,1.42962]" snapshot_55="~[5,1.42964]" snapshot_56="~[5,1.42975]" snapshot_57="~[5,1.42985]" snapshot_58="~[5,1.42986]" snapshot_59="~[5,1.42979]" snapshot_60="~[5,1.42972]" snapshot_61="~[5,1.4297]" snapshot_62="~[5,1.42974]" snapshot_63="~[5,1.42979]" snapshot_64="~[5,1.42981]" snapshot_65="~[5,1.42979]" snapshot_66="~[5,1.42975]" snapshot_67="~[5,1.42974]" snapshot_68="~[5,1.42975]" snapshot_69="~[5,1.42977]" snapshot_70="~[5,1.42978]" snapshot_71="~[5,1.42978]" snapshot_72="~[5,1.42976]" snapshot_73="~[5,1.42976]" snapshot_74="~[5,1.42976]" snapshot_75="~[5,1.42976]" snapshot_76="~[5,1.42977]" snapshot_77="~[5,1.42977]" snapshot_78="~[5,1.42977]" snapshot_79="~[5,1.42976]" snapshot_80="~[5,1.42976]" snapshot_81="~[5,1.42976]" snapshot_82="~[5,1.42977]" snapshot_83="~[5,1.42977]" snapshot_84="~[5,1.42977]" snapshot_85="~[5,1.42976]" snapshot_86="~[5,1.42976]" snapshot_87="~[5,1.42976]" snapshot_88="~[5,1.42977]" snapshot_89="~[5,1.42977]" snapshot_90="~[5,1.42977]" snapshot_91="~[5,1.42977]" snapshot_92="~[5,1.42976]" snapshot_93="~[5,1.42976]" snapshot_94="~[5,1.42977]" snapshot_95="~[5,1.42977]" snapshot_96="~[5,1.42977]" snapshot_97="~[5,1.42977]" snapshot_98="~[5,1.42977]" snapshot_99="~[5,1.42977]" snapshot_100="~[5,1.42977]" />
     </Result_Element>
 </result>
diff --git a/tests/2d_examples/test_2d_plate/regression_test_tool/PlateObserver_Position_Run_5_result.xml b/tests/2d_examples/test_2d_plate/regression_test_tool/PlateObserver_Position_Run_5_result.xml
diff --git a/tests/3d_examples/test_3d_roof/3d_roof.cpp b/tests/3d_examples/test_3d_roof/3d_roof.cpp
diff --git a/tests/3d_examples/test_3d_shell_stability_half_sphere/test_3d_shell_stability_half_sphere.cpp b/tests/3d_examples/test_3d_shell_stability_half_sphere/test_3d_shell_stability_half_sphere.cpp

Original file line number	Diff line number	Diff line change
`@@ -90,7 +90,7 @@ namespace SPH`
`90`	`90`	`{`
`91`	`91`	`protected:`
`92`	`92`	`Real CFL_;`
`93`		`- StdLargeVec<Vecd> &vel_, &acc_;`
	`93`	`+ StdLargeVec<Vecd> &vel_, &acc_, &acc_prior_;`
`94`	`94`	`Real smoothing_length_, c0_;`
`95`	`95`
`96`	`96`	`public:`
Original file line number	Diff line number	Diff line change
`@@ -69,7 +69,7 @@ namespace SPH`
`69`	`69`	`{`
`70`	`70`	`protected:`
`71`	`71`	`Real CFL_;`
`72`		`- StdLargeVec<Vecd> &vel_, &acc_, &angular_vel_, &dangular_vel_dt_;`
	`72`	`+ StdLargeVec<Vecd> &vel_, &acc_, &angular_vel_, &dangular_vel_dt_, &acc_prior_;`
`73`	`73`	`StdLargeVec<Real> &thickness_;`
`74`	`74`	`Real rho0_, E0_, nu_, c0_;`
`75`	`75`	`Real smoothing_length_;`