Hi all,
I am doing a contact simulation and would like to obtain the total area in contact. I am able to get an output using:
*************************Output, history
*CONTACT OUTPUT, MASTER=masterSurf, SLAVE=slaveSurf
CAREA
************************
However, the output is a single integer...
Eureka! I just finished implemented the changes to the Mathematica code and checking against Abaqus. Using the 'sigma-transverse == 0' equations, the system is harder to solve, so I have to work out a new method for fitting the parameters. I can deal with that however.
@IceBreakerSours:
I...
I'm using cylindrical because the final mesh will be round (see images above). Using cylindrical coordinates, I don't have to define a local material orientation for each element; I simply assign the coordinate system and move forward.
I'm doing it in the validation just to make sure it works...
Excellent! I never thought to use sigma_2 == 0 when solving for sigma_1 as a function of lambda_1! Yes, Mathematica should have no problem with this.
I can't wait to try this. I owe you both, you have been incredibly patient and helpful!
...is anisotropic, my intuition suggests this to be incorrect. The constitutive equation for incompressibilitiy should still hold however:
labmda_1*lambda_2*lambda_3 == 1
OR
labmda_t*lambda_r*lambda_z == 1
The problem is, how do I relate the lambda's? In order to use the starting SEDF, I need...
@sdebock,
Thank you for trying the input file! By lambda_t, do you mean the theta direction? I stretch the 10x10mm element 1.5mm in each direction (one per file), which should yield a lambda_t of 1.15. Can you clarify why you get a lambda_t of 1.3?
How did you get lambda_x from the output? Do...
Hey again,
I added the constraints for the rotational degrees of freedom where they were missing as IceBreakerSours suggested. The 'kappa' I am using is consistent with the example I found in the documentation. The mathematical model also seems to work okay with this value.
I'm still getting...
Thank you both! Yes, I meant to ask what 'S' in Abaqus CAE represents. Thank you for confirming sdebock.
I'll try constraining the rotational degrees of freedom as suggested by IceBreakerSours. I didn't think that a uniaxial displacement would have any influence on the node rotations.
Not yet. I'm reading Holzapfel's book over Christmas break so I understand the fundamentals better.
I found other sources saying that S is Cauchy, but I need to check the documentation.
I figured out which elements may be used with incompressibility by trying a C3D8 element and getting this error:
Looks like shell elements are okay. The final geometry we will be using is a 3D surface with uniform thickness, so shell elements are the most convenient. That's why I've been...
I checked both the .msg and .dat files and couldn't find the suggested warning. Also, in the Abaqus documentation example, they suggest that D=0 is acceptable:
Anyways, small changes to the compressibility do yield large changes in material response. Have you seen this paper? Its very recent...
Ah yes, E_alpha is positive for all lambda except for extreme cases of gamma, so that should be okay.
In that case I must be doing something else wrong; I'm verifying against Abaqus using a single 10x10mm shell element in tension. I displace one side of the element by 1.5mm, yielding a stretch...
...function seems a little weird to me.
If I understand correctly, the Cauchy stress 'sigma' in the 11 direction would be:
sigma_11 = lambda_1 * partial(U)/partial(lambda_1)
where 'U' is the strain energy density function as defined in the Abaqus documentation (neglecting the compressibility...
Edit 2: I'm looking in the 6.12 documentation for HGO now to see how the fiber orientations are specified.
Here is the desired behavior of the material model. You can see that both curves are highly non-linear, with the longitudinal direction about twice as distensible as the circumferential...
@IceBreakerSours
By "bifurcations", I mean that the vein branches. I should actually call them anastomoses I guess (joining of two smaller veins to form a larger vein). I've been defining a local cylindrical coordinate system for each branch up to now, which has worked out okay. I'll see how it...
@ IceBreakerSours:
The reason I moved away from the HGO model is that I had trouble getting two material orientations for the fibers using CAE.
When you say "allow two directions of anisotropy", do you mean specifying "Number of local directions" to be two? Would I then specify two sets of...
...material could be modelled by fitting two uniaxial stress-strain curves to Fung-orthotropic material constants. I guess I will need to substitute *reasonable* parameters for the out-of-plane behavior in order to define the 'transverse shear stiffness'. I'm having trouble finding any...
I am working on modelling the response of an orthotropic hyperelastic material (vein) to a given load. Because the local thickness of the material is unknown, shell elements with a constant thickness are used. I am able to get the simulation to run. However, I had some questions about the...