Recent Posts

Pages: [1] 2 3 ... 10
1
General questions / Re: Pressure-Inflation Boundary Condition Query
« Last post by Prof. R.L. Taylor on May 30, 2025, 01:02:30 PM »
It is not clear what the shape of the patient specific meshes looks like.  Are they irregular or some smooth surface?

Are there equations for the 3-d boundary conditions to be imposed?  Can you describe exactly what your colleagues implementation is?  The boundary conditions is nonlinear and most likely in feap can be implemented as a user element using either Lagrange multiplier constraints or aa penalty approach.
2
General questions / Re: Pressure-Inflation Boundary Condition Query
« Last post by chrismi on May 29, 2025, 01:02:02 PM »
Dear  Prof. S. Govindjee,

Thank you for your response.

That is indeed something I've considered and will likely attempt as a last-case scenario. The problem is that the meshes/geometries that I'm using are patient-specific and have been derived from corneal topography images, so generating additional elements to represent the surrounding tissues beyond the limbus would mean adding a degree of uncertainty to an already well-established geometry. Also, the sclera is a complicated tissue in its own right, and if included, I'd have to properly account for how it integrates/anchors with the cornea at the limbus.

The free rotation boundary conditions I described have been successfully implemented in my colleague's purpose-built in-house FE software, and multiple studies have demonstrated its suitability in replicating the in vivo deformation of the cornea. So I would very much like to employ the same boundary conditions in FEAP, using the material model I've already implemented using the user element functionality (a previous model I developed that I've repurposed to describe the cornea), to form an effective comparison with my group's prior work.

With that in mind, would it be possible to implement these boundary conditions in FEAP? And if so, could you please point me in the right direction of how to do so? I would greatly appreciate any help you could provide me.
3
General questions / Re: Pressure-Inflation Boundary Condition Query
« Last post by Prof. S. Govindjee on May 28, 2025, 12:42:34 PM »
I would suggest that you construct a model of the limbus around your cornea model using elastic elements with the appropriate stiffness but without all the geometric detail, and then fix the outer part of your limbus model.  The support for the cornea seems like an elastic foundation and it is probably better to just model a bit of the foundation so as to not to over or under constrain your cornea model.  This will be a lot easier than try to enforce a complex boundary condition and it will probably be closer to the physical reality.
4
General questions / Re: Pressure-Inflation Boundary Condition Query
« Last post by chrismi on May 28, 2025, 09:18:04 AM »
Sorry, my initial post wasn't very clear. I'm trying to replicate the in vivo deformation of the cornea when subjected to intraocular pressure within the eye. The boundary conditions are supposed to reflect the action of the surrounding tissues (iris, sclera), which are relatively soft and allow for the rotational freedom of the cornea under physiological deformation.

Therefore, the purple point (which in 3D represents a line around the periphery) should be fixed in space, and orthogonality conditions should be preserved with respect to the deformed mid-section of the cornea. There are no other restrictions.

Thank you for response, and taking the time to answer my questions.
5
General questions / Re: Pressure-Inflation Boundary Condition Query
« Last post by FEAP_Admin on May 25, 2025, 03:31:53 PM »
What else is allowed to occur at the boundary edge?  It appears that there is free rotation.  Is the purple point fixed or can it translate?  Can the edge warp? etc.?  Also are you trying to model a particular experiment? if so can you describe the actual experiment?
6
General questions / Pressure-Inflation Boundary Condition Query
« Last post by chrismi on May 24, 2025, 10:22:01 AM »
Dear Professors,

I'm having difficulty determining how to apply the desired boundary conditions for the problem I'm trying to solve. I wasn't sure where exactly to post this query, so I chose "General Questions".

I am modelling the in vivo deformation of the cornea (eye tissue), which is essentially equivalent to a classical bulge-inflation test. I've provided a few figures for illustrative purposes. I'm applying the pressure using surface elements (highlighted red) with a follower load, applied to the interior surface.

However, the boundary conditions at the periphery are not quite as trivial as defining the loading. I want to enforce that the tangent to the central plane (red dotted line) where it meets the periphery (purple marker) is always orthogonal to the periphery (blue line), which would provide some rotational freedom. The provided diagram is for the 2D section, as it's difficult to visualise the issue properly in 3D.
 
I am using version 8.4 of FEAP, if that is of significance.

I look forward to hearing back from you!
7
Programming / Re: Ogden Material more with more than 3 coefficients
« Last post by duchoangnguyen on May 19, 2025, 11:21:25 PM »
For completeness:
d(31) seems to be occupied elsewhere. Accordingly, I have skipped this and continued with d(32). this must also be adapted in wder3f.
In the event that more than 10 inputs are given, I have read in another line via PINPUT (in the first block), as I have only seen reading in via TINPUT.
I do not know if that was nessesary but it seems to work.

Thank you very much!
8
Programming / Re: Ogden Material more with more than 3 coefficients
« Last post by Prof. R.L. Taylor on May 19, 2025, 09:08:17 AM »
I am not able to access the files for ver8.1 -- it is very old and compiler changes make it risky to use.  You should consider an update.  So I am not sure if the following applies to your version.

The changes are made to "inmate.f" and "wder3f.f".

In 'inmate.f'

1. search for where 'ogde' is input and input the additional values you want to include into new local values.

2. search for where the Ogden is output and modify the output of the parameters to the file.  In my version this is done with format 2024.

3.  below the output the parameters are stored in d(21) to d(28).  I would move the number of parameters to d(22) and the place all the pairs after this in d(23) to the max.  There is a remote risk you may overlap some other use of the parameter but let's hope not!

In "wder3f.f" search for where the d(22) etc. are used and adjust to match what you reset in "inmate.f"

Good luck.
9
Input File Issues / Re: Ogden materials
« Last post by Prof. R.L. Taylor on May 19, 2025, 08:44:07 AM »
The main thing you are missing is a PROP to ramp the load in several load steps, there are also some extra commands I have removed.  You should always look at the "O" file to see what the analysis did with the commands you provided.
10
Input File Issues / Ogden materials
« Last post by sara.grbcic on May 19, 2025, 06:33:31 AM »
Hello,

I am trying to model a hyperelastic membrane problem consisting of a rubber sheet in a simple tension test under large stretch by using the Ogden model integrated in FEAP. To be more precise, I am trying to simulate the numerical solution presented in the book of prof. Ibrahimbegovic - Nonlinear Solid Mechanics, pg. 335, Fig. 4.19. which I am providing in the attachment (DeformedConfigOgdenAIBook.pdf)

When modeling the problem in FEAP, I encounter several problems:
1.) Nodes which are and should be restrained in the definition of the problem, move, (left-hand edge, restrained)
2.) No matter of the choice of the applied loading, or prescribed displacement, on the free (right-hand edge), the configuration deforms only in the lateral direction, while the longitudinal direction remains unaffected (no stretch in the x direction)

I am attaching two versions of my input files and screenshots of the deformed configurations.

Could you please help me with providing an input file which is able to model the deformation presented in DeformedConfigOgdenAIBook.pdf?

I tried adding the mixed formulation, but without any success.

Thanks in advance.

Best regards.
Pages: [1] 2 3 ... 10