Stuck? 
support@mecway.com
support@mecway.com
Mecway
FEA
Non-linear analysis of a masonry arch bridge
Hello,
I'm trying to analyze a masonry arch bridge by using a plane model. As of now, I have modelled all the components with linear elastic materials but I'd like to use a compression only model at least for the masonry arch.
Is this possible in Mecway?
Secondly, is any constitutive model for concrete available in Mecway?
Thanks,
p.s. I attach for reference the file I'm working on.
I'm trying to analyze a masonry arch bridge by using a plane model. As of now, I have modelled all the components with linear elastic materials but I'd like to use a compression only model at least for the masonry arch.
Is this possible in Mecway?
Secondly, is any constitutive model for concrete available in Mecway?
Thanks,
p.s. I attach for reference the file I'm working on.
Howdy, Stranger!
It looks like you're new here. If you want to get involved, click one of these buttons!
Comments
Mmartin has a nice example with only compression and very similar to yours in the forum.
Clay arch simulation
https://mecway.com/forum/discussion/1089/clay-arch-simulation
Be careful , I have seen some nodes disconnected on the lower corners. You can use the tool Open cracks in the View menu to find them.
I saw the clay arch simulation file but it is a 3d model. I was wondering whether it is possible to model compression only materials also in 2d models.
Regards,
I am experienced on simulating 3d archs using compression only material.
Maybe soon, I expect to get a new mfront compilation linked to MW with mazars material model ready to use. Then we will be able to reproduce cracks on brittle materials.
Manuel
Line load doesn't work .
Find example with random material properties for a compression only arch.
You need to erase the material on the main tree to avoid conflict.
Check it carefully. I never try this.
I've tried to modify the model as suggested but it doesn't run the analysis. I'm still testing the free edition and the model has less than 1000 nodes. However, when I try to save the file, I get an error message due to the number of nodes.
Please find here attached my modified file. Doesn't it run just because I'm using the free edition or for some other reason?
The error "element set ARCH has not yet been defined." is because Mecway didn't generate an ELSET for ARCH which is because it doesn't generate duplicate ELSETs and already has Element_Selection which is the same. So modify or remove the Element_Selection named selection.
It's also lacking Z displacement constraints but that's a subsequent issue.
*ERROR: too many cutbacks best solution and residuals are in the frd file
In addition, after running the analysis, I can't save the file anymore and I get this message:
Cannot save.
Free edition: Cannot solve, save or run sripts on models with more than 1000 nodes.
The last version of the file is here attached. Thanks for your support.
Andrea
By your description seems your analysis is running but not converging. That's a different problem.
Check the time settings in Nonlinear.
Regarding the number of nodes, I would try symmetry if you are just playing. Loads may not be symmetric but... just for testing could be enough.
Line load do not work. Try pressure or traction. Expand the geometry to see the thickness of the shell and you can select the area.
File attached.
On the other hand, if I try to set a new file with the same features you provided in the one you sent, the new analysis doesn't converge.
I still get the error message when I try to save and I don't understand why since the model has less than 1000 nodes, but this is a second order issue at the moment.
Since I'm novice with Mecway, is there any reference to help me through the Custom Model Definition?
To be able to use those new capabilities that Calculix offer, you need to learn how to set up an analysis for Calculix (or at least have some general idea)
COMPRESSION_ONLY set up can be found in the ccx Manual v 2.19 page 245
https://www.dhondt.de/ccx_2.19.pdf
Regarding the number of nodes, in the log file produced after the analysis I can see 56747 nodes, which corresponds to what you were saying. However, I think I will work out this issue with the complete software.
I'd like to have another clarification on the Mecway functioning that I can't find in the manual.
In the arch analysis there are mainly two load conditions: self weight of the model and variable load due to traffic. I'd like to run an analysis by applying the self weight first and then the traffic load. This because the self weight stabilizes the arch, while the traffic loads tend to produce tensile stresses.
As far as I understand, the CCX solver allows to increase all loads at the same time.
Isn't it possibile to run two different load steps, the first by applying the selt weight only and the second by applying the traffic loads after the self weight is already applied?
Thanks,
When Nonlinear Quasistatic analysis is active, loads can be entered as a function of time by means of a function or a table. Wait some seconds before activating the pressure load. You have a lot of possible combinations. Load functions are also available like heaviside, sin, cos, ...if you want to introduce them in a more complex way.
Once you familiarize with the software you can also introduce a load that cross the bridge. It is explained in the forum. Mpg video Attached.
https://mecway.com/forum/discussion/comment/5741/
If I have understood correctly, in defining the Compression_only material, just the Young modulus and the tensile strength are specified. Doesn't the Poisson coefficient influence the analysis too?
For questions more related with the theory involved you better ask at the Calculix forum directly. There are real experts there who can give you an answer or address you to some book or paper related with the subject.
https://calculix.discourse.group/
By calculating the mass of the arch on the basis of total vertical reactions, I get a value of 42628 kg, about 30% smaller than the real one.
I don't understand why I get this discrepancy. Since the weight of the arch is much smaller than the weight of the backfill, the global error in terms of reaction is low (about 5%), but I can't understand why there should be any difference.
Looking at the value of the error (30%), it might be due to something related to the shell thickness (3 m). Any hint on this?
It will be easier to follow.
The correction term is in fact identical to the term used to cut off tensile stresses for penalty contact in Equation(203) and Figure (130). Replacing “overclosure” and “pressure” by “principal strain” and “principal stress” in that figure yields the function f.
I guess “absolute value of the maximum allowed pressure” parameter should be read as : “absolute value of the maximum allowed principal stress”
It agrees with the result. If you just show the Arch in the Solution tree the values adjust for the arch.
2- ¿How do you get the reaction forces from shell elements?
2) I was looking at the sum of external forces in the y direction over the named selection BASE-SX.
You can create a New table in solution.
Check in the manual.
You have not defined the element selection called ARCH.
Your custom model ELSET is pointing to nowhere.
Maybe Victor can take a look at this as it should fail and it is not.
Take into consideration these are kind of tricky things are out of normal Mecway procedures.
Did you manage to obtain different results?
My model still provides a wrong result of the arch weight, that I measure by using the function Sum applied to External forces in the y direction over the named selection BASE-SX.
In the same way, the maximum principal tensile stress at midspan (and in the same way the tangential stress, which should be equivalent in that point) at step 10 (t=1 s) is equal to 0.032 MPa, larger than the defined 0.02 MPa.
i.e. The tensile stress strain behavior is not accurate for the compression only material, but one that allows for a solution. The method I use is to run as a normal material, and accept the load at which the modulus of rupture is achieved as the limit, or run as a compression only material with the tensile strength approximately equal to that provided by any reinforcing.
BASE-SX is only one side's support. Not sure if you doubled that? I get 801kN as the sum of external force Y on all the Y constraint nodes (displacement_nodes(2)) at t=0.5 s. Is that the erroneous value?
I'm pretty sure density is always mass per unit volume, not area, so I wouldn't expect thickness to be a problem.
I started again from Andrea file and one additional refinement makes the big difference.
@Andrea.
I have set up the whole model as Compression only and everything keeps below 0.02 MPa.
If I set only the Arch, it exceeds punctually to 0.024MPa but it is in a corner where the node is shared with a regular element (Node value is averaged so it can exceed 0.02Mpa.
However, my main concern on this point was that both the compression-only analysis and the elastic analysis provide the same tensile stress of 0.032 MPa.
This sounded strange to me but probably it is just a matter of refinement.
The problem on mass is of major importance to me and I need to work it out.
@Victor If I measure on BASE-SX the sum of external force Y at t = 0.5 s, I get 1610 kN. BASE-SX does correspond to one side's support.
On the other hand, if I measure on Displacement_nodes(2) the sum of external force Y at t = 0.5 s, I get 3204 kN, which is approximately the double of the value I obtained on BASE-SX.
Now, if I run an elastic analysis, the sum of external force Y over BASE-SX is 1702 kN, which should be the correct value from hand calculations.
What do you think?
Thanks,