public static void GenerateAllParts(MatSetting matSetting, Model model)
{
//lining element, shell143
Part partShell = new Part(1, 1, 1, 1);
ElementType typeShell = new ElementType(1, EType.shell143);
Mat matShell = new Mat(1, matSetting.shell_coe, matSetting.shell_pr, matSetting.shell_dens);
Shell143 sectionShell = new Shell143(1);
sectionShell.thickness = matSetting.shell_thickness;
model.AddPart(partShell);
model.AddElementType(typeShell);
model.AddMat(matShell);
model.AddSection(sectionShell);
//radial ground spring, link180
Part partRadialGround = new Part(2, 2, 2, 2);
ElementType typeRadialGround = new ElementType(2, EType.link180);
Mat matRadialGround = new Mat(2, matSetting.radial_spring_coe, matSetting.radial_spring_pr, 0);
Link180 linkRadialGround = new Link180(2);
linkRadialGround.area = matSetting.radial_spring_area;
linkRadialGround.keyOption1 = 3;
linkRadialGround.keyOption2 = 2;
model.AddPart(partRadialGround);
model.AddElementType(typeRadialGround);
model.AddMat(matRadialGround);
model.AddSection(linkRadialGround);
//tangential ground spring, link180
Part partTangentialGround = new Part(3, 3, 3, 3);
ElementType typeTangentialGround = new ElementType(3, EType.link180);
Mat matTangentialGround = new Mat(3, matSetting.tangential_spring_coe, matSetting.tangential_spring_pr, 0);
Link180 linkTangentialGround = new Link180(3);
linkTangentialGround.area = matSetting.tangential_spring_area;
model.AddPart(partTangentialGround);
model.AddElementType(typeTangentialGround);
model.AddMat(matTangentialGround);
model.AddSection(linkTangentialGround);
//circumferential joint
Part partCircumferentialJoint = new Part(4, 4, 0, 4);
ElementType typeCircumferentialJoint = new ElementType(4, EType.combin39);
//Mat matCircumferentialJoint = new Mat(4, 0, 0, 0);
Combin39 combinCircumferentialJoint = new Combin39(4);
combinCircumferentialJoint.keyOption1 = 3;
combinCircumferentialJoint.keyOption2 = 1;
combinCircumferentialJoint.stiffness = GeneratePairs(
-0.001 * matSetting.circumferential_joint_compression, -0.001,
0, 0,
0.001 * matSetting.circumferential_joint_tensile, 0.001);
model.AddPart(partCircumferentialJoint);
model.AddElementType(typeCircumferentialJoint);
//model.AddMat(matCircumferentialJoint);
model.AddSection(combinCircumferentialJoint);
//longitudinal joint rotation X
Part partLongitJointZ = new Part(5, 5, 0, 5);
ElementType typeLongitJointZ = new ElementType(5, EType.combin39);
//Mat matLongitJointZ = new Mat(5, 0, 0, 0);
Combin39 combinLongitJointZ = new Combin39(5);
combinLongitJointZ.keyOption1 = 3;
combinLongitJointZ.keyOption2 = 4;
combinLongitJointZ.stiffness = GeneratePairs(GenerateCoordinate(
matSetting.longitudinal_joint_rotation_k_1, matSetting.longitudinal_joint_rotation_strain_1,
matSetting.longitudinal_joint_rotation_k_2, matSetting.longitudinal_joint_rotation_strain_2));
model.AddPart(partLongitJointZ);
model.AddElementType(typeLongitJointZ);
//model.AddMat(matLongitJointZ);
model.AddSection(combinLongitJointZ);
//circumferential joint Y
Part partLongitJointX = new Part(6, 6, 0, 6);
ElementType typeLongitJointX = new ElementType(6, EType.combin39);
//Mat matLongitJointX = new Mat(6);
Combin39 combinLongitJointX = new Combin39(6);
combinLongitJointX.keyOption1 = 3;
combinLongitJointX.keyOption2 = 2;
combinLongitJointX.stiffness = GeneratePairs(GenerateCoordinate(
matSetting.longitudinal_joint_rotation_k_1 / 2, matSetting.longitudinal_joint_rotation_strain_1,
matSetting.longitudinal_joint_rotation_k_2 / 2, matSetting.longitudinal_joint_rotation_strain_2));
model.AddPart(partLongitJointX);
model.AddElementType(typeLongitJointX);
//model.AddMat(matLongitJointX);
model.AddSection(combinLongitJointX);
//circumferential joint Z
Part partLongitJointY = new Part(7, 7, 0, 7);
ElementType typeLongitJointY = new ElementType(7, EType.combin39);
//Mat matLongitJointY = new Mat(7);
Combin39 combinLongitJointY = new Combin39(7);
combinLongitJointY.keyOption1 = 3;
combinLongitJointY.keyOption2 = 3;
combinLongitJointY.stiffness = GeneratePairs(GenerateCoordinate(
matSetting.longitudinal_joint_rotation_k_1 / 2, matSetting.longitudinal_joint_rotation_strain_1,
matSetting.longitudinal_joint_rotation_k_2 / 2, matSetting.longitudinal_joint_rotation_strain_2));
model.AddPart(partLongitJointY);
model.AddElementType(typeLongitJointY);
//model.AddMat(matLongitJointY);
model.AddSection(combinLongitJointY);
//pre-load element, in axial direction, link180
Part partAxial = new Part(8, 8, 8, 8);
ElementType typeAxial = new ElementType(8, EType.link180);
Mat matAxial = new Mat(8, matSetting.radial_spring_coe, matSetting.radial_spring_pr, 0);
Link180 linkAxial = new Link180(8);
linkRadialGround.area = matSetting.radial_spring_area;
linkRadialGround.keyOption1 = 3;
linkRadialGround.keyOption2 = 2;
model.AddPart(partRadialGround);
model.AddElementType(typeRadialGround);
model.AddMat(matRadialGround);
model.AddSection(linkRadialGround);
}
private void SettingInitial()
{
//Model Setting
modelSetting = new ModelSetting();
modelSetting.outerRadius = 3.1;
modelSetting.thickness = 0.35;
modelSetting.width = 1.2;
modelSetting.num_ring = 2;
modelSetting.num_longit = 1;
modelSetting.num_circum = 90;
modelSetting.pos_joint.Add(8);
modelSetting.pos_joint.Add(73);
modelSetting.pos_joint.Add(138);
modelSetting.pos_joint.Add(222);
modelSetting.pos_joint.Add(287);
modelSetting.pos_joint.Add(352);
modelSetting.axialForce = 15000;
modelSetting.circumferential_joint_length = 0.4;
modelSetting.circumferential_joint_diameter = 0.03;
modelSetting.circumferential_joint_Es = 200000000;
modelSetting.circumferential_joint_Ec = 35000000;
modelSetting.circumferential_ring_friction_factor = 0.6;
modelSetting.isCurve = true;
modelSetting.axis = TunnelTools.GetTunnelFootprintAxis(121800);
modelSetting.Initial();
//Mat Setting
matSetting = new MatSetting();
matSetting.shell_coe = 35000000;
matSetting.shell_pr = 0.25;
matSetting.shell_dens = 2.5;
matSetting.shell_thickness = modelSetting.thickness;
matSetting.radial_spring_coe = 5000;
matSetting.radial_spring_pr = 0.25;
matSetting.radial_spring_area = 2 * PI * modelSetting.modelRadius
/ modelSetting.num_circum * modelSetting.width / modelSetting.num_longit;
matSetting.tangential_spring_coe = 100;
matSetting.tangential_spring_pr = 0.25;
matSetting.tangential_spring_area = matSetting.radial_spring_area / 3;
matSetting.circumferential_joint_tensile = modelSetting.circumferential_joint_Es
* modelSetting.circumferential_joint_area / modelSetting.circumferential_joint_length;
matSetting.circumferential_joint_compression = matSetting.circumferential_joint_tensile * 1000;
matSetting.longitudinal_joint_rotation_k_1 = 200000 / modelSetting.num_longit;
matSetting.longitudinal_joint_rotation_strain_1 = 0.0015;
matSetting.longitudinal_joint_rotation_k_2 = 20000 / modelSetting.num_longit;
matSetting.longitudinal_joint_rotation_strain_2 = 0.002;
matSetting.longitudinal_joint_shear_force = modelSetting.axialForce * modelSetting.circumferential_ring_friction_factor / modelSetting.num_circum;
}
