Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Nonlinear bending of an isotropic plate

Please login with a confirmed email address before reporting spam

Hello,

I m trying to simulate a very known benchmark problem which also apears on Reddy and Arciniega's article, ''Tensor-based finite element formulation for geometrically nonlinear analysis of shell structures''.

The problem consists of a plate with width=12 and height=1 (as defined in Comsol), E (module Young) = 1.2*10^6,
poisson ratio is 0 and the thickness of the plate is 0.1. The boundary conditions are fixed constraint at one edge and moment at the opposite edge. (the moment is defined in the Local Edge System and not in the Global Coordinate System).

While in the article Reddy and Arciniega are able to to reach fully circular rotation, using comsol i only can reach much lower rotations when the applied moment is 42. applying moments higher then 42 makes the solution false.

When applying moments lower then 42 , i get rotations that are accurate when comparing to analytical results, but i still wonder why i cant fully replicate what Reddy and Rrciniega did.

thx a lot
Iftah Nudel

5 Replies Last Post Feb 23, 2012, 7:39 a.m. EST
Nagi Elabbasi Facebook Reality Labs

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Feb 19, 2012, 5:13 p.m. EST
If you post your model maybe someone on the Forum can help. Also, COMSOL support can help with such problems.
If you post your model maybe someone on the Forum can help. Also, COMSOL support can help with such problems.

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Feb 20, 2012, 5:24 a.m. EST
Hi,

I m attaching the model.

thx a lot
Hi, I m attaching the model. thx a lot


Ivar KJELBERG COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Feb 20, 2012, 11:44 a.m. EST
Hi

first of all define the material, or remove it and use the "User Definition fileds" (but all of them)

Then turn on the plot while solving (solver node) and use a parameter and a continuation parameter to ramp up your load from typically Param=0..1 * -45 (define Parameter Param = 0, select the stationary solver node, go the the Extensions tab and select "continuation" add (+) the Param and set a range such as " range(0,0.1,1)"

Still I'm not sure the shell physics accept such large deformations, one would need to go down and dig into the doc first


by the way check if you cannot update to latest version, there are a few improvement in the latest patch

you could also try with a 2D solid first

--
Good luck
Ivar
Hi first of all define the material, or remove it and use the "User Definition fileds" (but all of them) Then turn on the plot while solving (solver node) and use a parameter and a continuation parameter to ramp up your load from typically Param=0..1 * -45 (define Parameter Param = 0, select the stationary solver node, go the the Extensions tab and select "continuation" add (+) the Param and set a range such as " range(0,0.1,1)" Still I'm not sure the shell physics accept such large deformations, one would need to go down and dig into the doc first by the way check if you cannot update to latest version, there are a few improvement in the latest patch you could also try with a 2D solid first -- Good luck Ivar

Nagi Elabbasi Facebook Reality Labs

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Feb 22, 2012, 10:49 a.m. EST
Use the parametric feature that Ivar mentioned and also put a symmetry condition on Line 3. Right now after some deformation the plate switches to another deformation mode. It’s a valid mode but that’s not the one you want. Adding the symmetry condition gives you the benchmark solution you’re expecting.

Nagi Elabbasi
Veryst Engineering
Use the parametric feature that Ivar mentioned and also put a symmetry condition on Line 3. Right now after some deformation the plate switches to another deformation mode. It’s a valid mode but that’s not the one you want. Adding the symmetry condition gives you the benchmark solution you’re expecting. Nagi Elabbasi Veryst Engineering

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Feb 23, 2012, 7:39 a.m. EST
Thx a lot! it's working.
Thx a lot! it's working.

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.