// برخی از مقادیر پارامتر هندسی را از Parameters بخوانید
Double voxel_size = model.param().evaluate(“vox_size”);
Double geom_tol = model.param().evaluate(“geom_tol”);
Double xmin = model.param().evaluate(“xmin”)+geom_tol;
دابل xmax = model.param().evaluate(“xmax”)+geom_tol;
Double ymin = model.param().evaluate(“ymin”)+geom_tol;
Double ymax = model.param().evaluate(“ymax”)+geom_tol;

// قسمت پاک کردن هندسه
model.geom(“part1”).feature().clear();

int i = 0; // شمارنده ردیف برای ماتریس

// متغیر برای محاسبه ضخامت الکترود
Double maxz = new Double(0.0);

// متغیرهایی برای محاسبه شعاع کره با میانگین جرم
Double Vsum = new Double(0.0);
Double rVsum = new Double(0.0);

// حلقه برای تمام سطرهای جدول
برای (String[] row : tableData) {
Double xminData = Double.parseDouble(row[2])*voxel_size;
Double yminData = Double.parseDouble(row[3])*voxel_size;
Double xmaxData = Double.parseDouble(row[5])*voxel_size;
Double ymaxData = Double.parseDouble(row[6])*voxel_size;

// بررسی کنید که آیا در کادر محدود است، و اگر شعاع بزرگتر از تلورانس است، کره اضافه کنید
if (((xminData < xmax && xminData > xmin) || (xmaxData < xmax && xmaxData > xmin)) && ((
yminData < ymax && yminData > ymin) || (ymaxData < ymax && ymaxData > ymin))
{ volume = Double.parseDouble(row[1])*Math.pow(voxel_size, 3.0);
شعاع دو برابر = Math.pow(0.75*volume/Math.PI، 0.333333)؛

if (radius > geom_tol) {
Double xpos = Double.parseDouble(row[8])*voxel_size;
Double ypos = Double.parseDouble(row[9])*voxel_size;
Double zpos = Double.parseDouble(row[10])*voxel_size;

// افزودن داده به ماتریس داده های کره ای
sphereData[i][0] = xpos;
sphereData[i][1] = ypos;
sphereData[i][2] = zpos;
sphereData[i][3] = شعاع;
i++;

maxz = maxz < (zpos+radius) ? zpos+radius : maxz;

// متغیرهای شعاع میانگین جرم را به‌روزرسانی کنید، اگر مرکز کره در کادر محدودکننده باشد
if ((xpos < xmax && xpos > xmin) && (ypos < ymax && ypos > ymin))) {
Vsum = Vsum+volume;
rVsum = rVsum+شعاع*حجم;
}

// افزودن کره به بخش هندسه
GeomFeature sph = model.geom(“part1”).create(“sph”+i, “Sphere”);
sph.set(“pos”, new String[]{Double.toString(xpos)+”[m]”, Double.toString(ypos)+”[m]”, Double.toString(zpos)+”[m] “})؛
sph.label(“Particle “+(int) Math.round(Double.parseDouble(row[0])));
sph.set(“r”, radius+”[m]”);

double hmin = model.param().evaluate(“hmin”); // حداقل اندازه عنصر مش
int cyl = 1;

// فواصل و شکاف‌های بین تمام ترکیب‌های دوتایی کره‌ها را
برای (int j = 0; j < i; j++) بررسی کنید {
for (int k = j+1; k <i; k++) {
Double x = sphereData[k][ 0] -sphereData[j][0];
Double y = sphereData[k][1]-sphereData[j][1];
z Double = sphereData[k][2]-sphereData[j][2];

فاصله دو برابر = Math.sqrt(x*x+y*y+z*z);

دو شکاف = dist-(sphereData[k][3]+sphereData[j][3]);
int index1 = sphereData[k][3] <sphereData[j][3] ? k : j;
int index2 = sphereData[k][3] > sphereData[j][3] ? k : j;
دابل r_sphere = sphereData[index1][3];

// اگر کره ها خیلی نزدیک یا همپوشانی دارند، استوانه اضافه کنید
if (gap <hmin && gap > -0.667*r_sphere) {

Double x_sphere = sphereData[index1][0];
Double y_sphere = sphereData[index1][1];
z_sphere دو برابر = sphereData[index1][2];

جهت دوگانه = sphereData[k][3] <sphereData[j][3] ? -1.0 : 1.0;

دابل r_help = r_sphere+gap-hmin;
Double r_cyl = Math.max(Math.sqrt(r_sphere*r_sphere-r_help*r_help), hmin*2);

// معادله ریشه به غیر خطی را برای شعاع استوانه حل کنید (به طوری که فاصله کره به کره = hmin در امتداد سطح استوانه).
دو برابر r1 = r_sphere;
Double r2 = sphereData[index2][3];
دو تا = 2*r2;
tb دو برابر = 2 * r1;
double tc = Math.sqrt(-((gap-hmin+tb+ta)*(gap-hmin+tb)*(gap-hmin+ta)*(gap-hmin)));
r_cyl = ((1/(شکاف-hmin+r1+r2))*tc)/2;

Double dist_corr = Math.sqrt(r_sphere*r_sphere-r_cyl*r_cyl)-2*hmin;
دابل posx = x_sphere+x/dist*dist_corr*direction;
دابل posy = y_sphere+y/dist*dist_corr*direction;
دابل posz = z_sphere+z/dist*dist_corr*direction;

// ایجاد سیلندر
GeomFeature cylinder = model.geom(“part1”).create(“cyl”+cyl++, “Cylinder”);
با (سیلندر)؛
set(“pos”، new String[]{Double.toString(posx)+”[m]”، Double.toString(posy)+”[m]، Double.toString(posz)+”[m]”} )
set(“axitype”، “cartesian”);
set(“ax3″، new String[]{Double.toString(x*direction)+”[m]”, Double.toString(y*direction)+”[m]”، Double. toString(z*direction)+”[m]”});
set(“r”، r_cyl+”[m]”);
set(“h”، hmin*5+”[m]”);
پایان یافتن با()؛
}
}
}

// تنظیم طول الکترود در مدل
model.param().set(“L_elec”, maxz+geom_tol+”[m]”);

// ساخت هندسه در model.geom
(“part1”).run();
model.component(“comp1”).geom(“geom1”).run();

// اندازه گیری و به روز رسانی پارامتر حجم
model.component(“comp1”).geom(“geom1”).measure().selection().init(3);
model.component(“comp1”).geom(“geom1”).measure().selection().named(“uni1”);
model.param().set(“Vp”, model.geom(“geom1″).measure().getVolume()+”[m^3]”);

//
مدل پارامتر ناحیه را اندازه گیری و به روز کنید.
model.component(“comp1”).geom(“geom1”).measure().selection().named(“intsel1”);
model.param().set(“Ap”, model.geom(“geom1″).measure().getVolume()+”[m^2]”);

// به روز رسانی پارامتر شعاع میانگین جرمی
model.param().set(“rp_avg_spheres”, rVsum/Vsum+”[m]”);