...did not start and the following error was observed:
“Permeability is specified on 61 elements but the initial void ratio is not defined. Use *initial conditions to the value. The elements are identified in element set ErrElemPermNoInitVoidRatio.
Unidentified or invalid output request...
Thank you for your responses.
I need field output.
The filed output generates the resultant shear force. Is it possible to define a field output that dived this resultant shear force by the corresponding element area?
Since I need resultant shear stress rather than resultant shear force.
Thank you
Hello
I am wondering if there is any way to visualise the resultant shear stress on a part?
I have two objects in contact and produce shear stress in different directions.
Like CSHEARF magnitude force, I wonder whether there is any way to calculate the resultant shear stress on the surface...
Hi
I have attached the CAE file.
You can download the ODB results from here.
Yes, it can be said that object A is a bone, object B is a residual limb and object C is a socket for that.
Object A has steel properties.
Objects B and C have hyperelastic properties...
Thanks for your response.
I have checked with displacement control, but there was no difference. Since the change in the stress on contact between objects B and C is absolute zero.
there should be an initial clearance, then after applying the force and the object B comes down, B and C should...
Hello
I am trying to simulate a force being applied on a cylinder inside an object encompassed by another semi-conical shape object.
Sorry, to make the explanation clearer, I named the parts (objects) inside the model and attached its image.
Similar to the attached image:
Object A is tied inside...
According to Matlab fmincon documentation, the 'trust-region-reflective' algorithm needs to specify the objective gradient.
I read the algorithm explanation but yet, I cannot understand the reason for that.
1- Why does it need gradient as input?
2- Does the objective function have to be...
...which I am sure it is way more than enough and ABAQUS solve it in 120 increment (0.5 time step, and no backward increment).
Let`s assume that touch occurred at the increment 40 (equal to time 20), in this case the pressure at the contact is: P= (500*40)/120 = 167 kPa
is it what you mean...
Hello
According to the attached picture, I have two surfaces (Master surface and slave surface) where there is a gap between them. The slave surface can be inflated with the pressure "P".
I aim to find an appropriate "P" magnitude to fill the gap and only TOUCH the master surface without...
Hello
In MATLAB, I am optimising a multiobjective function using the lsqnonlin function with the trust-region reflective algorithm. I was wondering how I can get the hessian matrix as an output? Matlab provides a matrix of Jacobian,is it possible to calculate Hessian matrix using the Jacobian...
Thank you so much, I will have a look into that.
I am wondering does it interpolate the result for the nodes that I have not defined?
For instance, my analytical solution returns the result for 60 nodes out of 500 nodes. Does Abaqus interpolate to visualise the rest of 440 nodes, or it will not...
Thank you for your response.
However, I am still not sure what I should do, or even I am doing it right.
I define my model (including mesh and loads etc), then in the property module, I define my map field. XYZ matrix coordinate and filed is my analytical displacement (gold standard).
Then I run...
Hello
I have 3D point clouds obtained from 3D image reconstruction, and I have also run a simulation with Abaqus and got the result.
I am trying to map the points obtained from the point cloud on the mesh and show them.
I followed the instruction in the documentation, but I do not know how can I...
Hello
I am trying to use python script to generate new surface in order to apply pressure on it using script.
I use the:
session.journalOptions.replayGeometry
session.journalOptions.setValues(replayGeometry=INDEX)
to achieve the index number, but the problem is that I cannot understand the...
...19.526271), ))
but I expect to get the following code work: (note that 0.785398 is the RAD of 45 DEG and I am assuming that it is (x,y,z)
Tsin=85*sin(0.785398)
Tcos=85*cos(0.785398)
s1 = a.instances['CylnderPart-1'].faces
side1Faces1 = s1.findAt(((85.0, 0, 21.0), ), ((-85.0, 0...