New point generator functions

Project.toml

@@ -1,7 +1,7 @@
 name = "Peridynamics"
 uuid = "4dc47793-80f3-4232-b30e-ca78ca9d621b"
 authors = ["Kai Partmann"]
-version = "0.3.3"
+version = "0.4.0-DEV"
 
 [deps]
 AbaqusReader = "bc6b9049-e460-56d6-94b4-a597b2c0390d"

docs/src/api_reference.md

@@ -37,6 +37,9 @@ precrack!
 contact!
 uniform_box
 uniform_sphere
+uniform_cylinder
+round_sphere
+round_cylinder
 n_points
 ```
 

docs/src/literate/tutorial_logo.jl

@@ -39,9 +39,9 @@ material!(plate; horizon=3.015ΔX₀ₚ, E=30e9, rho=2000, Gc=10)
 cz = Ø / 2 + lz / 2 + 1.1 * ΔX₀ₛ
 r_logo = Ø / 2 + 0.2 * Ø
 sxy, cxy = r_logo * sin(30π / 180), r_logo * cos(30π / 180)
-posₛ₁, volₛ₁ = uniform_sphere(Ø, ΔX₀ₛ; center_y=r_logo, center_z=cz)
-posₛ₂, volₛ₂ = uniform_sphere(Ø, ΔX₀ₛ; center_x=cxy, center_y=-sxy, center_z=cz)
-posₛ₃, volₛ₃ = uniform_sphere(Ø, ΔX₀ₛ; center_x=-cxy, center_y=-sxy, center_z=cz)
+posₛ₁, volₛ₁ = uniform_sphere(Ø, ΔX₀ₛ; center=(0, r_logo, cz))
+posₛ₂, volₛ₂ = uniform_sphere(Ø, ΔX₀ₛ; center=(cxy, -sxy, cz))
+posₛ₃, volₛ₃ = uniform_sphere(Ø, ΔX₀ₛ; center=(-cxy, -sxy, cz))
 sphere₁ = Body(BBMaterial(), posₛ₁, volₛ₁)
 sphere₂ = Body(BBMaterial(), posₛ₂, volₛ₂)
 sphere₃ = Body(BBMaterial(), posₛ₃, volₛ₃)

src/Peridynamics.jl

@@ -17,7 +17,7 @@ export NoCorrection, EnergySurfaceCorrection
 # Discretization
 export Body, point_set!, point_sets, failure_permit!, material!, velocity_bc!, velocity_ic!,
        forcedensity_bc!, precrack!, MultibodySetup, contact!, uniform_box, uniform_sphere,
-       n_points
+       uniform_cylinder, round_sphere, round_cylinder, n_points
 
 # Running simulations
 export VelocityVerlet, DynamicRelaxation, Job, submit

src/auxiliary/precompile_workload.jl

@@ -6,8 +6,8 @@
         error(msg)
     end
     root = joinpath(@__DIR__, "temp_precompilation")
-    pos1, vol1 = uniform_box(1, 1, 1, 0.5; center_x=0.5, center_y=0.5, center_z=0.5)
-    pos2, vol2 = uniform_box(1, 1, 1, 0.5; center_x=-0.5, center_y=0.5, center_z=0.5)
+    pos1, vol1 = uniform_box(1, 1, 1, 0.5; center=(0.5, 0.5, 0.5))
+    pos2, vol2 = uniform_box(1, 1, 1, 0.5; center=(-0.5, 0.5, 0.5))
     path_bb = joinpath(root, "BB")
     path_osb = joinpath(root, "OSB")
     path_nosb = joinpath(root, "NOSB")