yade-users team mailing list archive

[Idan <question706221@xxxxxxxxxxxxxxxxxxxxx>]

New question #706221 on Yade:
https://answers.launchpad.net/yade/+question/706221

Hello to everyone,

I thought about trying ChatGPT 4 for writing new material laws and physics into YADE DEM, 

as I have little experience with the YADE code, it would be helpful if someone can give me a feedback.

Here is a code for temperature dependent cohesive burgers model written by ChatGPT 4:

HEADER:

#pragma once
#include<yade/pkg/dem/CohesiveFrictionalContactLaw.hpp>
#include<yade/pkg/dem/ScGeom.hpp>
#include<yade/pkg/dem/FrictPhys.hpp>

class TempDependentCohesiveBurgersViscoElasticModel: public CohFrictMat {
public:
    TempDependentCohesiveBurgersViscoElasticModel(Real young, Real poisson, Real frictionAngle, Real density, Real cohesion, Real viscousDamping, Real viscosity, Real tempRef);
    virtual ~TempDependentCohesiveBurgersViscoElasticModel();

    // Temperature reference for the model
    Real tempRef;

    // Viscosity for the viscoelastic component of the model
    Real viscosity;

    YADE_CLASS_BASE_DOC_ATTRS_CTOR(TempDependentCohesiveBurgersViscoElasticModel, CohFrictMat, "Temperature-dependent cohesive Burger's viscoelastic model inheriting from CohFrictMat",
        ((Real, viscousDamping, 0, ,"Viscous damping for the viscoelastic component of the model"))
        ((Real, tempRef, 298.15, ,"Reference temperature for the model"))
        ((Real, viscosity, 0, ,"Viscosity for the viscoelastic component of the model")),
        /* ctor */ TempDependentCohesiveBurgersViscoElasticModel(young, poisson, frictionAngle, density, cohesion, viscousDamping, viscosity, tempRef); );
};
REGISTER_SERIALIZABLE(TempDependentCohesiveBurgersViscoElasticModel);

class TempDependentCohesiveBurgersViscoElasticContactLaw: public Law2_ScGeom_FrictPhys_CundallStrack {
public:
    virtual void go(shared_ptr<IGeom>& ig, shared_ptr<IPhys>& ip, Interaction* contact);
    YADE_CLASS_BASE_DOC(TempDependentCohesiveBurgersViscoElasticContactLaw, Law2_ScGeom_FrictPhys_CundallStrack, "Contact law for the temperature-dependent cohesive Burger's viscoelastic model");
};
REGISTER_SERIALIZABLE(TempDependentCohesiveBurgersViscoElasticContactLaw);


CCP File:

#include "TempDependentCohesiveBurgersViscoElastic.hpp"

TempDependentCohesiveBurgersViscoElasticModel::TempDependentCohesiveBurgersViscoElasticModel(Real young, Real poisson, Real frictionAngle, Real density, Real cohesion, Real viscousDamping, Real viscosity, Real tempRef)
    : CohFrictMat(young, poisson, frictionAngle, density, cohesion)
    , viscousDamping(viscousDamping)
    , viscosity(viscosity)
    , tempRef(tempRef) {
}

TempDependentCohesiveBurgersViscoElasticModel::~TempDependentCohesiveBurgersViscoElasticModel() {
}

YADE_PLUGIN((TempDependentCohesiveBurgersViscoElasticModel));

void TempDependentCohesiveBurgersViscoElasticContactLaw::go(shared_ptr<IGeom>& ig, shared_ptr<IPhys>& ip, Interaction* contact) {
    // Cast the geometry and physical parameters to the correct types
    ScGeom* geom = static_cast<ScGeom*>(ig.get());
    FrictPhys* phys = static_cast<FrictPhys*>(ip.get());

    // Get the associated bodies and materials
    Body* b1 = Body::byId(contact->getId1(), scene);
    Body* b2 = Body::byId(contact->getId2(), scene);
    TempDependentCohesiveBurgersViscoElasticModel* mat1 = static_cast<TempDependentCohesiveBurgersViscoElasticModel*>(b1->material.get());
    TempDependentCohesiveBurgersViscoElasticModel* mat2 = static_cast<TempDependentCohesiveBurgersViscoElasticModel*>(b2->material.get());

    // Compute the normal and tangential forces based on the temperature-dependent cohesive Burger's viscoelastic model
    // Note: You will need to implement the specific calculations for the model
    Vector3r normalForce, tangentialForce;
    // ...

    // Apply the forces to the bodies
    scene->forces.addForce(contact->getId1(), normalForce + tangentialForce);
    scene->forces.addForce(contact->getId2(), -(normalForce + tangentialForce));

    // Update the interaction geometry
    // ...
}

YADE_PLUGIN((TempDependentCohesiveBurgersViscoElasticContactLaw));


It is not fully completed but does this look fine? what is the best IDE to be used debugging such model? thanks.



-- 
You received this question notification because your team yade-users is
an answer contact for Yade.