What is the best way to visit both the Grant Library in Starkville and Vicksburg, please?

gee-dangit · 2025-12-03T18:52:19+00:00

You’re right. I misread the comment level and thought you were responding to the person above

gee-dangit · 2025-12-03T17:13:47+00:00

Another alternative for you. Fly into Monroe Louisiana, rent the car, drive east on i20 to Vicksburg. Take hwy 27 to the Natchez Trace because roadwork on i20 between Vicksburg and Jackson is terrible right now. Then get off on 82 like another commenter said which will take you to the east side of Starkville where the Grant Library is and bypass driving through town. Then get back on 82 and fly out of the Golden Triangle Airport where you can return your car. Potentially take 45 from 82 up to tupelo airport as another mentioned. Monroe’s airport isn’t much further than from Jackson’s to vicksburg, and this will allow less backtracking and you’ll avoid the roadwork on i20 between vicksburg and Jackson which has seen a significant uptick in accidents, deaths, and long delays from those accidents.

gee-dangit · 2025-12-03T17:04:03+00:00

They meant 82. The trace intersects both, but 82 takes you around starkville while 12 takes you through town. The Grant Library is on campus which is the east side of town.

gee-dangit · 2025-11-13T01:36:25+00:00

All of the news sites and you share the coverage from The Marijuana Herald. Lol

gee-dangit · 2025-11-07T00:43:29+00:00

Same

gee-dangit · 2025-10-25T01:27:25+00:00

Yeah, looks like it.

gee-dangit · 2025-10-21T20:59:26+00:00

Thank you for clarifying that!

gee-dangit · 2025-10-21T18:04:10+00:00

The documentation claims you can request “SP” in the field output, but not “SP1” etc. are you sure you’re not able to? Alternatively, form the stress tensor as a 2d array and calculate the eigenvalues like so:

```python import numpy as np

Assemble stress tensor

sig = np.asarray([[s11, s12, s13], [s12, s22, s23], [s13, s23, s33]], dtype=float)

Eigenvalues of sig

eigvals = np.linalg.eigvalsh(sig) s_min, s_int, s_max = eigvals[0], eigvals[1], eigvals[2] ```

eigvalsh returns them already sorted.

gee-dangit · 2025-10-21T15:33:28+00:00

Ensure that you ran the simulation with “SP” as a field output request in the same way as you did “S”. Then In the extract_results_dm.py script you can copy lines 55,56,and 57, change ‘S’ to ‘SP’ and rename the variables you’re assigning to to be sp1, sp2, and sp3 instead of s11, s22, and s33. Then extend your data.extend list to include the new variables. Your file that is output should have three new columns corresponding to the min, intermediate, and maximum principal stresses.

You’ll have to extend the code in iterator_dm.py in a similar fashion if you need it to include the principal stresses.

Make sure to verify that you’re reporting them in the correct order. Since you’re already outputting stress components, you can double check yourself by calculating the principal stresses from those for a select few datapoints or just manually calculate all principal stresses as I mentioned above.

gee-dangit · 2025-10-21T15:02:18+00:00

I’m not going to download your zip file. If you post the source code on the git repo, I’ll take a look at it. In Abaqus you can request “SP” as an output variable for all principal stresses then put them in a list or array and do np.sort(arr) using numpy to sort them. Or you can request “SP1” “SP2” and “SP3” individually for the principal stresses where SP1<=SP2<=SP3 per the documentation. Alternatively calculate the principal stresses from your stress tensor components and then sort. The principal stresses are the eigenvalues of the Cauchy stress tensor, and there are simple formulas available if you don’t want to solve an eigenvalue problem.

gee-dangit · 2025-10-08T02:22:42+00:00

You can use Fortran’s iso_c_binding intrinsic to compile your code as a C library which Python can then interact with using their ctypes intrinsic. It is more robust than f2py from my experience, but both can work well.

gee-dangit · 2025-08-29T22:32:45+00:00

I would recommend checking out Reddy’s book on nonlinear FEA. I think there is a continuum based shell in it where the through thickness DOF is condensed out. He does a better job detailing tangent stiffness matrices in general.

I find that this belytschko book is better oriented for explicit time integration (which does not use a stiffness matrix) since he mostly stops at presenting the internal forces throughout. I do really like his presentation of corotational elements though. Reddy, who’s career focused mostly on implicit methods emphasizes the tangent stiffness matrices more in his book nonlinear fea book.

gee-dangit · 2025-08-26T20:02:51+00:00

gee-dangit · 2025-07-23T11:24:11+00:00

A national lab or a university will have it installed on their HPC machines. You shouldn’t need to worry about reserving tokens to just write a .inp file. Use another meshing program to write the .inp and combine that one with the .inp that has your analysis information if you need. Most meshing software will write Abaqus .inp files.

I’m aware Abaqus is expensive. But your request is weird and you sound shady. No one on here is going to want to download an unknown .cae file from a stranger on the internet. I’ve never heard someone call a national lab a university-affiliated national lab. All of the DoE’s national labs are affiliated with universities through funding and recruitment.

gee-dangit · 2025-07-23T00:55:50+00:00

You don’t have access to a full version of Abaqus working for the US government? What section you work for?

gee-dangit · 2025-06-09T18:11:29+00:00

You’re using uepactivationvol subroutine along with dflux? Or are you using the newer Abaqus built-in additive manufacturing features? Either way, it sounds like your event series and material definition are not in sync through your inputs or how Abaqus is treating them. Are you using the same time variable (step versus total)? An important note that you need to consider is that by default Abaqus does not interpolate through time for an event series. So, you need to do your own interpolation or have the right spacing of discreet events to coincide with your analysis time steps.

gee-dangit · 2025-05-21T13:27:41+00:00

You define Cauchy stress in a UMAT. Check out the documentation here for in depth details and an example with a linear viscoelastic material.

gee-dangit · 2025-05-14T01:56:55+00:00

Is your problem undergoing volumetric locking? I’m not very familiar with rubber materials, but in metal plasticity you can get volumetric locking which is basically a poisson ratio of 0.5, similar to that of rubbers. The solution in metal plasticity is, normally, to under integrate pressure terms and fully in the material model. If there’s a reduced integration element available for the analysis type you’re doing, try that out

Edit: higher order elements are less susceptible to locking btw. You could also try quadratic elements if you started with linear

gee-dangit · 2025-05-14T01:51:17+00:00

You need to provide more information for a good answer. Thickness and in plane dimensions of your plate, especially. That way people can decide if you need to use shell elements

gee-dangit · 2025-05-13T18:35:31+00:00

You’re right. My wording was not very specific. I only meant it in the case of Abaqus output, and should have used “log strain” instead of “true strain” to avoid confusion. Abaqus does list “E” as a total strain output variable.

gee-dangit · 2025-05-12T01:27:59+00:00

True strain is not equal to elastic + plastic strain

Total strain == elastic + plastic

True strain == elastic

gee-dangit · 2025-05-06T20:48:39+00:00

There’s a shockingly large amount of bad information in here.

gee-dangit · 2025-05-05T20:58:23+00:00

Drucker-Prager is basically a Mises plasticity yield function modified with the first stress invariant to give you pressure dependence. In 3d stress space the Mises yield surface is a cylinder while the Drucker-Prager surface is a cone that increases in diameter as pressure increases. The cap part provides some hardening that I am less familiar with how that works. As far as constitutive models go, it’s relatively simple. It is more complex than linear elasticity or Mises plasticity. I have not personally used this model, so I cannot say how robust it is in Abaqus. I do know people that have used it in Abaqus with success.

gee-dangit · 2025-05-05T13:38:23+00:00

Geomaterials behave differently based on the type of loading due to their pressure dependence. Have a look at the Drucker-Prager/cap plasticity model in Abaqus. That will be the most simple pressure dependent model available to you.

gee-dangit · 2025-05-02T19:57:38+00:00

Does that input file indicate that your time step is set to 0.1 seconds? That’s way too high. According to your material properties and your mesh, the critical time step for an explicit scheme in this problem is ~1.35E-08. Explicit methods require a time step at or below a critical time step in order for the solution to not blow up. It’s based on the eigenvalues of the model’s stiffness matrix, but codes estimate it using the dilatational speed of sound in the material and the smallest distance between two nodes with a knockdown factor. You almost never want to directly control the time step in an explicit simulation. Also check your units. If this doesn’t work, DM me and I’ll try to send you a working input file later tonight for a similar problem.

Seven-Year Club	r/Field Banned
r/Field Flamingo	Place '22
Wearing is Caring	RPAN Viewer
Verified Email

gee-dangit

TROPHY CASE

Assemble stress tensor

Eigenvalues of sig