Victor

About

Username
Victor
Joined
Visits
4,719
Last Active
Roles
Administrator

Comments

  • You can enter formulas in the direction vector components. Use y and -x for the X and Y components respectively to constraint in the circumferential direction about Z. Use a separate displacement constraint for the axial displacement.
  • @ahjiuhdotcom There must be 3D solid elements in the solution. So stress linearization won't work for axisymmetric with the internal solver but will with CCX where the output is the expanded thin solid wedge. Does that help? Otherwise, if the line …
  • Stress appears for me. You do have to solve it again if you added stress after it solved - unlike the internal solver. Are you using the Mystran that's installed with Mecway, not one you obtained separately?
  • Mystran's 8-node hex is just the classic compatible mode formulation so it's terrible in bending. Hex20 should be OK though. Stress should be available for both solids and shells. Perhaps post your .liml file and I'll see what's wrong.
  • The CCX manual identifies it as "Keyword type: model definition, material" which means it goes in a material data block, which is the cards immediately following *MATERIAL. Do that in Mecway using the CCX -> Custom material type by copying all th…
  • No, sorry. Some ideas, none of which are easy: - Fit it to a curve (somehow) and enter it as a formula in Mecway. - Convert the pressure data directly to CCX *DLOAD cards and put them in CCX -> Custom step contents in Mecway. The biggest obstac…
  • Not in Mecway's post-processor but CCX can output stress at integration points to the human-readable .dat file by putting the following in the *STEP section (CCX -> custom step contents in Mecway): *EL PRINT,ELSET=EALL S and having an ELSET …
  • For the beam model, just add a thermal stress item and unsuppress the temperatures.
  • Hello Marcos 1. Don't split the bolt. Make it a single part. Mecway finds the faces for the split and CalculiX splits it. 2. I'm not sure what you mean but the split happens roughly in the middle of the prism (cylinder in this case). The ends of t…
  • I made 4 changes: - Set the analysis type to Static 3D. - Unsuppressed the mechanical constraints. - Added a thermal stress item. I chose 0°C as the Reference temperature so you might have to change that. - Added displacement and reaction force to t…
  • I should also add that items displayed in red text in the outline tree have errors. To see what it is, right click it and pick What's wrong? Items with orange text are warnings which are also worth looking at if something's not working right.
  • Yes, shells can too. Could you post the .liml file and I can identify all the problems. I notice there's no thermal stress item in Loads & Constraints. You need to include that when you do the static solve.
  • Yes, beam elements work with thermal stress. Your screenshot shows analysis type is thermal, so change it to static which has displacement DOFs. That's also why many of the boundary conditions are shown in red, which means error (displacement const…
  • From your pictures, it looks like Spooles doesn't converge in either version so there's no inconsistency there, and Pastix is the one that gives different results. Is that right? Pastix can give random results with mixed precision, which is the def…
  • File not attached. That is a worry if they both converge with different solutions. I'm looking forward to seeing your example. You can use PARDISO by building from the source included with Mecway, or possibly download a compiled binary from dhondt…
  • @Fatmac C3D8I performs much better than C3D8, which is practically useless for anything with bending. C3D8R has some problems according to the manual - under-stiff in bending and constant strain due to 1 integration point. C3D20 is pretty good alrea…
    in Mecway 33 Comment by Victor March 22
  • Hello windys53. The forum can be flakey sometimes and might need a reload or cache clearing. To get radial displacement in the solution, use Solution -> Formula and enter a formula like ( x*u.x + y*u.y ) / sqrt(x^2+y^2)
  • I've seen this before. It seems to be common for elastic contact on curved surfaces with incompatible meshes, even when there's no interference. I wonder if it might be a bug in CCX.
    in CXX Contact Comment by Victor March 13
  • Hello Youk 1. It can't import OpenFoam pressure, and having a different pressure on each element/node is difficult in Mecway because each one appears as a separate pressure item in the outline tree. You could write a script to generate the .liml fi…
  • There's something strange about your picture showing time = 1000000000 s. I didn't wait for a solution that way though and got a more reasonable looking one by doing this: * Add No extra contact step to the CCX branch in the outline tree. * Ch…
  • I see. It looks like the AI is trying to editing the mesh instead of applying a displacement constraint. You can do a radial prescribed displacement without a script, using formulas in Displacement like this: If you need to adjust the distance in…
  • You can do that with Mesh tools -> Move/copy, Set Coordinate system to Cylindrical (Z axis). But if you do need a script, could you please format the script to "Code" in the forum otherwise indenting is lost and other things get changed. Is it …
  • Shell offset doesn't currently work with OpenRadioss. I'm not sure if OpenRadioss even has that capability? I see that Mecway allows you to define shell offset and doesn't give any error message. This is wrong and I'll make it clear that it won't wo…
  • I'm not sure. It looks like it has read equivalent plastic strain from the CCX output because it's showing zero instead of nothing. Perhaps have a look at the .frd file to see if there should be non-zero values on those nodes.
  • This looks like a bug Disla found https://github.com/Dhondtguido/CalculiX/issues/129 for heat flow which has apparently been fixed in CCX 2.23.
  • Happy to hear that Jachym! You can also come here if you have any questions - there are some very knowledgeable and helpful people around.
  • 40/(0.9*1.565) is half the length of the bridge (unit is meter) so that it starts from the -x end instead of the origin. unittriangle(w,x) is this with area = 1. You add a constant to x to shift it in space, and add a multiple of t to make it shift…
  • I get the same displacement at the center deck node at t=20 s. What is the difference you're seeing?
  • You may have to defeature the geometry (in CAD) by removing the tiny radius on the trailing edge. Looks like that's what Fatmac did too?
  • I'll change the traction unit from MPa to Pa because there's a unit conversion factor in there that makes it more confusing to me: -(180/4.3) * unittriangle(0.5, (0.9*1.565)*t-40/(0.9*1.565)-x ) the unittriangle() function defines a force per unit…
Avatar

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!