! --- Element and Material Definition --- ET,1,PLANE183 ! Element Type: 2D 8-Node Structural Solid ! Set KEYOPT to Plane Stress, which assumes a unit thickness of 1.0 by default KEYOPT,1,3,2 ! Option: Plane stress (thickness = 1.0)
! --- Geometry Creation --- ! Create the main plate outline by adding a rectangle and a circle RECTNG,0,150,0,100 ! Area 1: Rectangle CYL4,150,50,50 ! Area 2: Circle for the rounded end AADD,1,2 ! Add areas 1 and 2. Result is Area 3.
! Create the circular areas for the holes to be subtracted CYL4,150,50,25 ! Area 4: Large central hole (Radius = 50/2 = 25) CYL4,20,80,7.5 ! Area 5: Top-left small hole (Radius = 15/2 = 7.5) CYL4,20,20,7.5 ! Area 6: Bottom-left small hole (Radius = 15/2 = 7.5)
! Subtract the holes from the main plate ! The numbering is sequential. The plate is Area 3, holes are 4, 5, 6. ASBA,3,4 ! Subtract large hole (4) from plate (3) ASBA,3,5 ! Subtract top-left hole (5) from plate (3) ASBA,3,6 ! Subtract bottom-left hole (6) from plate (3)
! --- Meshing --- ESIZE,3 ! Set element edge length to 3 mm for a finer mesh MSHAPE,0,2D ! Use quadrilateral elements for the mesh MSHKEY,0 ! Use free meshing AMESH,ALL ! Mesh the final area (all)
FINISH
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! 2. SOLUTION !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ /SOLU ANTYPE,STATIC ! Set analysis type to Static
! --- Apply Constraints (Boundary Conditions) --- ! Fix the left edge of the plate (X=0) LSEL,S,LOC,X,0 ! Select the line at X=0 DL,ALL,,ALL,0 ! Constrain all DOF (UX & UY) on the selected line
! --- Apply Loads --- ! Apply a 500N downward force distributed on the bottom arc of the large hole ! Define a local cylindrical coordinate system at the hole center LOCAL,11,1,150,50,0 CSYS,11 ! Activate the local system
! Select nodes on a 60-degree arc at the bottom of the hole NSEL,S,LOC,X,25 ! Select nodes with a radius of 25 from the local origin NSEL,R,LOC,Y,240,300! Refine selection to an angle between 240 and 300 degrees
! Distribute the total force among the selected nodes *GET,N_COUNT,NODE,0,COUNT ! Get the number of selected nodes F_PER_NODE = -500/N_COUNT ! Calculate force per node
! Apply force in the global Y-direction to each selected node F,ALL,FY,F_PER_NODE
! --- Return to Global Settings & Solve --- CSYS,0 ! Return to global cartesian system ALLSEL,ALL ! Reselect all entities SOLVE ! Solve the system FINISH
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! 3. POST-PROCESSING !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ /POST1 SET,LAST ! Read the last result set from the results file
! --- Element and Material Definition --- ET,1,PLANE183 ! Element Type: 2D 8-Node Structural Solid KEYOPT,1,3,2 ! Option: Plane stress (assumes unit thickness of 1.0)
! --- Geometry Creation (CORRECTED LOGIC) --- ! Create the main plate outline by adding a rectangle and a circle RECTNG,0,150,0,100 ! Area 1: Rectangle CYL4,150,50,50 ! Area 2: Circle for the rounded end AADD,1,2 ! Add areas 1 and 2 to make the main plate *GET,main_plate,AREA,,NUM,MAX ! Get the area number of the combined plate
! Create the circular areas for the holes and store their IDs large_hole = CYL4(150,50,25) small_hole_1 = CYL4(20,80,7.5) small_hole_2 = CYL4(20,20,7.5)
! --- Subtract Holes Sequentially --- ! We must update the 'main_plate' variable after each subtraction because ! the ASBA command creates a new area with a new ID number.
! Subtract the large hole ASBA,main_plate,large_hole *GET,main_plate,AREA,,NUM,MAX ! IMPORTANT: Update the plate ID to the new area number
! Subtract the first small hole ASBA,main_plate,small_hole_1 *GET,main_plate,AREA,,NUM,MAX ! IMPORTANT: Update the plate ID again
! Subtract the second small hole ASBA,main_plate,small_hole_2
! --- Meshing --- ESIZE,3 ! Set element edge length to 3 mm for a finer mesh MSHAPE,0,2D ! Use quadrilateral elements for the mesh MSHKEY,0 ! Use free meshing AMESH,ALL ! Mesh the final area (all)
FINISH
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! 2. SOLUTION !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ /SOLU ANTYPE,STATIC ! Set analysis type to Static
! --- Apply Constraints (Boundary Conditions) --- ! Fix the left edge of the plate (X=0) LSEL,S,LOC,X,0 ! Select the line at X=0 DL,ALL,,ALL,0 ! Constrain all DOF (UX & UY) on the selected line
! --- Apply Loads --- ! Apply a 500N downward force distributed on the bottom arc of the large hole ! Define a local cylindrical coordinate system at the hole center LOCAL,11,1,150,50,0 CSYS,11 ! Activate the local system
! Select nodes on a 60-degree arc at the bottom of the hole NSEL,S,LOC,X,25 ! Select nodes with a radius of 25 from the local origin NSEL,R,LOC,Y,240,300! Refine selection to an angle between 240 and 300 degrees
! Distribute the total force among the selected nodes *GET,N_COUNT,NODE,0,COUNT ! Get the number of selected nodes F_PER_NODE = -500/N_COUNT ! Calculate force per node
! Apply force in the global Y-direction to each selected node F,ALL,FY,F_PER_NODE
! --- Return to Global Settings & Solve --- CSYS,0 ! Return to global cartesian system ALLSEL,ALL ! Reselect all entities SOLVE ! Solve the system FINISH
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! 3. POST-PROCESSING !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ /POST1 SET,LAST ! Read the last result set from the results file
/CLEAR
/FILNAME,Plate_Analysis_NoReal
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 1. PRE-PROCESSING
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
/PREP7
! --- Element and Material Definition ---
ET,1,PLANE183 ! Element Type: 2D 8-Node Structural Solid
! Set KEYOPT to Plane Stress, which assumes a unit thickness of 1.0 by default
KEYOPT,1,3,2 ! Option: Plane stress (thickness = 1.0)
! --- Material Properties (as specified in image) ---
MP,EX,1,2E5 ! Young's Modulus: 2e5 MPa (N/mm^2)
MP,PRXY,1,0.3 ! Poisson's Ratio: 0.3
! --- Geometry Creation ---
! Create the main plate outline by adding a rectangle and a circle
RECTNG,0,150,0,100 ! Area 1: Rectangle
CYL4,150,50,50 ! Area 2: Circle for the rounded end
AADD,1,2 ! Add areas 1 and 2. Result is Area 3.
! Create the circular areas for the holes to be subtracted
CYL4,150,50,25 ! Area 4: Large central hole (Radius = 50/2 = 25)
CYL4,20,80,7.5 ! Area 5: Top-left small hole (Radius = 15/2 = 7.5)
CYL4,20,20,7.5 ! Area 6: Bottom-left small hole (Radius = 15/2 = 7.5)
! Subtract the holes from the main plate
! The numbering is sequential. The plate is Area 3, holes are 4, 5, 6.
ASBA,3,4 ! Subtract large hole (4) from plate (3)
ASBA,3,5 ! Subtract top-left hole (5) from plate (3)
ASBA,3,6 ! Subtract bottom-left hole (6) from plate (3)
! --- Meshing ---
ESIZE,3 ! Set element edge length to 3 mm for a finer mesh
MSHAPE,0,2D ! Use quadrilateral elements for the mesh
MSHKEY,0 ! Use free meshing
AMESH,ALL ! Mesh the final area (all)
FINISH
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 2. SOLUTION
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
/SOLU
ANTYPE,STATIC ! Set analysis type to Static
! --- Apply Constraints (Boundary Conditions) ---
! Fix the left edge of the plate (X=0)
LSEL,S,LOC,X,0 ! Select the line at X=0
DL,ALL,,ALL,0 ! Constrain all DOF (UX & UY) on the selected line
! --- Apply Loads ---
! Apply a 500N downward force distributed on the bottom arc of the large hole
! Define a local cylindrical coordinate system at the hole center
LOCAL,11,1,150,50,0
CSYS,11 ! Activate the local system
! Select nodes on a 60-degree arc at the bottom of the hole
NSEL,S,LOC,X,25 ! Select nodes with a radius of 25 from the local origin
NSEL,R,LOC,Y,240,300! Refine selection to an angle between 240 and 300 degrees
! Distribute the total force among the selected nodes
*GET,N_COUNT,NODE,0,COUNT ! Get the number of selected nodes
F_PER_NODE = -500/N_COUNT ! Calculate force per node
! Apply force in the global Y-direction to each selected node
F,ALL,FY,F_PER_NODE
! --- Return to Global Settings & Solve ---
CSYS,0 ! Return to global cartesian system
ALLSEL,ALL ! Reselect all entities
SOLVE ! Solve the system
FINISH
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 3. POST-PROCESSING
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
/POST1
SET,LAST ! Read the last result set from the results file
! --- Plot Stress Results ---
PLNSOL,S,EQV,0,2.0
/TITLE, Von-Mises Stress (MPa)
!
! --- Plot Deformation Results ---
PLNSOL,U,SUM,0,2.0
/TITLE, Total Deformation (mm)
!
/CLEAR
/FILNAME,Plate_Analysis_Corrected
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 1. PRE-PROCESSING
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
/PREP7
! --- Element and Material Definition ---
ET,1,PLANE183 ! Element Type: 2D 8-Node Structural Solid
KEYOPT,1,3,2 ! Option: Plane stress (assumes unit thickness of 1.0)
! --- Material Properties (as specified in image) ---
MP,EX,1,2E5 ! Young's Modulus: 2e5 MPa (N/mm^2)
MP,PRXY,1,0.3 ! Poisson's Ratio: 0.3
! --- Geometry Creation (CORRECTED LOGIC) ---
! Create the main plate outline by adding a rectangle and a circle
RECTNG,0,150,0,100 ! Area 1: Rectangle
CYL4,150,50,50 ! Area 2: Circle for the rounded end
AADD,1,2 ! Add areas 1 and 2 to make the main plate
*GET,main_plate,AREA,,NUM,MAX ! Get the area number of the combined plate
! Create the circular areas for the holes and store their IDs
large_hole = CYL4(150,50,25)
small_hole_1 = CYL4(20,80,7.5)
small_hole_2 = CYL4(20,20,7.5)
! --- Subtract Holes Sequentially ---
! We must update the 'main_plate' variable after each subtraction because
! the ASBA command creates a new area with a new ID number.
! Subtract the large hole
ASBA,main_plate,large_hole
*GET,main_plate,AREA,,NUM,MAX ! IMPORTANT: Update the plate ID to the new area number
! Subtract the first small hole
ASBA,main_plate,small_hole_1
*GET,main_plate,AREA,,NUM,MAX ! IMPORTANT: Update the plate ID again
! Subtract the second small hole
ASBA,main_plate,small_hole_2
! --- Meshing ---
ESIZE,3 ! Set element edge length to 3 mm for a finer mesh
MSHAPE,0,2D ! Use quadrilateral elements for the mesh
MSHKEY,0 ! Use free meshing
AMESH,ALL ! Mesh the final area (all)
FINISH
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 2. SOLUTION
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
/SOLU
ANTYPE,STATIC ! Set analysis type to Static
! --- Apply Constraints (Boundary Conditions) ---
! Fix the left edge of the plate (X=0)
LSEL,S,LOC,X,0 ! Select the line at X=0
DL,ALL,,ALL,0 ! Constrain all DOF (UX & UY) on the selected line
! --- Apply Loads ---
! Apply a 500N downward force distributed on the bottom arc of the large hole
! Define a local cylindrical coordinate system at the hole center
LOCAL,11,1,150,50,0
CSYS,11 ! Activate the local system
! Select nodes on a 60-degree arc at the bottom of the hole
NSEL,S,LOC,X,25 ! Select nodes with a radius of 25 from the local origin
NSEL,R,LOC,Y,240,300! Refine selection to an angle between 240 and 300 degrees
! Distribute the total force among the selected nodes
*GET,N_COUNT,NODE,0,COUNT ! Get the number of selected nodes
F_PER_NODE = -500/N_COUNT ! Calculate force per node
! Apply force in the global Y-direction to each selected node
F,ALL,FY,F_PER_NODE
! --- Return to Global Settings & Solve ---
CSYS,0 ! Return to global cartesian system
ALLSEL,ALL ! Reselect all entities
SOLVE ! Solve the system
FINISH
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 3. POST-PROCESSING
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
/POST1
SET,LAST ! Read the last result set from the results file
! --- Plot Stress Results ---
PLNSOL,S,EQV,0,2.0
/TITLE, Von-Mises Stress (MPa)
! --- Plot Deformation Results ---
PLNSOL,U,SUM,0,2.0
/TITLE, Total Deformation (mm)