public void AddElementType(ElementType type) { if (_elementTypes.ContainsKey(type.eid)) return; _elementTypes[type.eid] = type; }
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); }