VM51

VM51
Electrostatic Forces Between Charged Spheres

Overview

Reference: Any General Physics Textbook
Analysis Type(s): Electrostatic Analysis (ANTYPE = 0)
Element Type(s):
3D 20-Node Electrostatic Solid Element (SOLID122)
3D 10-Node Tetrahedral Electrostatic Solid Element (SOLID123)
3D Infinite Solid Element (INFIN111)
2D 8-Node Electrostatic Solid (PLANE121)
Meshing Facet (MESH200)
Input Listing:

vm51.dat

VM51 requires a supplemental .cdb input file which is too long to include full input listings. This file must be downloaded and placed in your working directory for the test case to run properly. Additionally, the geometry and mesh should be regenerated.

Download link: MAPDL Test Case Files for 2024 R2

vm51_case1.cdb

vm51_case2.cdb

Test Case

Two spheres with radii = 1 m, separated by a distance of 3 m, are subjected to a surface charge. Find the resultant electrostatic force between the spheres.

Figure 68: Charged Spheres Problem Sketch

Charged Spheres Problem Sketch

Material PropertiesGeometric PropertiesLoading
ε = 8.854E-12 F/m
r1 = 1.0 m
r2 = 3.0 m
r3 = 6.0 m
r4 = 1.25 m
Surface charge = 8.854 E-12 C/m2

Analysis Assumptions and Modeling Notes

The finite element mesh and the infinite element mesh are generated from the solid model. A planar section of the model is meshed with PLANE121 elements and rotated through 30 degrees. The 2D mesh creates a 3D mesh of SOLID122 elements. SHELL281 elements are generated over the outer surface of the 3D mesh and extruded in the radial direction to complete the finite element domain. The process is repeated to extrude the INFIN111 mesh in the radial direction.

It can be assumed that a symmetry plane exists at Y = 0, at which a zero voltage constraint is imposed. Infinite flags are set for the outer surface of the INFIN111 elements.

The same problem is then solved using SOLID123 elements

The theoretical solution is:

where ε = 8.854E-12, q1 = q2 = 4* π* ε. This charge corresponds to a surface force of ε on the sphere.

surface charge = q/area = (4* π* ε )/4* π*(r1**2) = ε

Results Comparison

SOLID122
TargetMechanical APDLRatio
YFORCE-0.1236E-10-0.1236E-101.000
SOLID123
TargetMechanical APDLRatio
YFORCE-0.1236E-10-0.1235E-100.999

Figure 69: Solid Model

Solid Model

Figure 70: FEA Model with SOLID122 and INFIN111 Elements

FEA Model with SOLID122 and INFIN111 Elements

Figure 71: Electric Field Plot with SOLID122 and INFIN111 Elements

Electric Field Plot with SOLID122 and INFIN111 Elements

Figure 72: Voltage Plot with SOLID122 and INFIN111 Elements

Voltage Plot with SOLID122 and INFIN111 Elements

Figure 73: FEA Model with SOLID123 and INFIN111 Elements

FEA Model with SOLID123 and INFIN111 Elements

Figure 74: Electric Field Plot with SOLID123 and INFIN111 elements

Electric Field Plot with SOLID123 and INFIN111 elements

Figure 75: Voltage Plot with SOLID123 and INFIN111 Elements

Voltage Plot with SOLID123 and INFIN111 Elements