I haven't used this... but, there's a Prony series viscoelastic model in ABAQUS, if I am not mistaken. You should be able to specify the shear relaxation modulus in ABAQUS so that it follows your...
Type: Posts; User: matsgd
I haven't used this... but, there's a Prony series viscoelastic model in ABAQUS, if I am not mistaken. You should be able to specify the shear relaxation modulus in ABAQUS so that it follows your...
Print DSTRAN and DSTRES and make sure your routine does the right thing.
You'd have to use multi-point constraints. But the easiest way to do what you need to do is to write a piece of Matlab code in which you construct a regular, nice, mesh on the surface of the cube,...
Sounds like a sign error in the constratint equation. If you're using *EQUATION to accomplish this, be advised that all terms are on the lefthand side of the equation, meaning that a constraint A=B...
Ziguo,
The strain increment is used by ABAQUS to estimate a globally stable time step. Without going into detail (mostly because I really don't know these things well enough!), the stable time...
I *think* that F is also, in fact, in the local coordinates of the element. You should not use the 33 component of F, I think (not sure what it contains when passed to the VUMAT). Instead, F(3,3)...
If you want to, you could send me the UMAT. I wouldn't spread it. Write me a PM with your email address if you want.
Mats
Is this a a case of plane stress then?
What strain measure are you trying to calculate? If you are trying to calculate the Green-Lagrange strain, then do it this way:
E=0.5*(U*U - I)=0.5*(F^t*F...
The internal force vector of a solid element is given by
F=integral(BL^T*sigma)dV
where BL is the linear strain displacement matrix, and sigma is the stress tensor. In a UL formulation, BL is...
1) Not an expert on beams, but in general, stresses are given by the elastic strains and the elastic moduli. Take a look at Hooke's law for the 3-D case, and impose sig22=sig33=0. Then you should get...
No, they are the transverse normal strains, i.e. affecting the thickness of your beam. And, no, ABAQUS does not do this by itself, because it cannot. The transverse normal STRESSES are, by...
It must be that your argument list in the UMAT is incorrect, or something along those lines. I notice that there are spaces that separate certain variable names. Is that ok? I don't know. By the way,...
I agree! Consider the case of imposing a deformation such that I2=0. Does it seem alright to produce a negative strain energy..?
Mats
Be more specific. Give us a piece of your code so that we can take a look at it...
M
Here's a suggestion:
Print, to a text file, the stress at some time t+dt, and the six stress vectors corresponding to the perturbed states. Then post them here. Let's take a look at it all.
Also,...
Yes, that won't work. Here's why: The *PLASTIC assumes that plastic flow is incompressible, i.e. plastic volumetric strains remain zero, and the plastic behavior only depends on the deviatoric part...
No it is not possible. In a VUMAT, however, it is possible to deactivate elements (integration points).
Mats
This is steel? With a factor of 35 between tensile and compressive yield stresses? Do you have some more info on it?
Mats
You need to use one of the material models with pressure dependent yield surfaces / flow rules. Examples are the Mohr-Coulomb and Crushable foam models.
Mats
*depvar
1
... allocates memory for one state variable that your UMAT can use. The length of STATEV is 1, since NSTATV = 1 in this case (in the UMAT).
The code uses the utility routine DBFILE to read from the .fil-file. Each time the DBFILE is called, it reads one record, which is an array. This array contains values, both integers and double...
The dimensions of ARRAY and JRRAY are given so that they'll match the sizes that the various utility routines require. See the examples in the URDFIL documentation. The equivalence statement is used...
And one more thing before I start anything. The UEL documentation states, regarding the force vector RHS, that
"For most nonlinear analysis procedures, NRHS=1 and RHS should contain the residual...
... is it possible?
I realize that nodal dofs can be tied to dummy elements, provided that there IS a standard element with the same topology. But what about integration point values? Is there a...
I don't have time to run your problem, but I notice that you're assigning elements of ARRAY to integers (KD1, KD2, KD3). Instead, use names for the output starting with (A-H,O-Z), or simply replace
...