The OpenSees Community

You need to set penalty a much larger number to reinforce contact law. for example, please change the following line: set k 10000; # penalty factor If you record the displacement, you can check by hand at what time, and at what acceleration that the sliding starts, theoretically. I cannot test your ...

The example pushes two blocks with applied axial load = 0.8, and lateral load=100, with frictional coefficient=0.6 between blocks. The static solution is non-existent since the system cannot stay in equilibrium. My suggestion: (1) reduce lateral load to a reasonable number (2) set k as a large penal...

I am curious about your remark. But, on what occasion do you really need hanging nodes? Impact analysis?

I donnot have whole lot of time to look into the code, but it seems to me the code contains some errors and a lot of misunderstanding. Just some answers to your questions: See the system below: **************** 7**** **************** |**** **************** |**** **************** |**** 4_________ 3 _...

So if use uniform excitation, the absolute acceleration = base input in uniform excitation + recorded acceleration, right? I want to clearify it is a plus not a minus.

I saw you already had a good answer in mind: use two 2D nodes for then zeroLength, then equalDOF the 1st and 2nd DOF of one of the 2D node to the 3D node of your pile.

Hello Perou, All your three questions can boil down to one answer: I have modified bbarQuad formulation to include (1) the axisymmetric case (2) self weight, and (3) a stress recorder to give stress state at Gaussian points, since the current version is incapable of doing all these. However, the mod...

#***************************************** # # # Cyclic Simple Shear Test for Clays # # # # MultiAxialCyclicPlasticity Model # # # # Gang Wang and Nicholas Siter # # # # CEE, University of California, Berkeley # # # # June 4, 2004 # # # #*****************************************# # #################...

as far as I checked the code today, the answer is no.
One has code a little bit to use shape function to project Gaussian point values to nodal values and modify function, eg. getResponse() etc., accordingly for that element. Not sure who is going to do that.

This command constructs a MultiaxialCyclicPlasticity material object, which can be used in 3D, 2D plane strain, and 2D axisymmetric analysis. The model is based on bounding surface concept. For the volumetric response, the stress-strain relationship is linearly elastic, as specified by a constant bu...

and add the command in the mannual. It will happen in a couple of days.
To be more correct, the model is developed for 3D, planestrain, and
axisymmetric cases, that is why it is called a multi-axial model.

Of course some materials have stress-strain information. The material
just takes trial strain from the element.
Why we have to use FEF to apply an element load?