cmp

Author: claudioperez

content/en/examples/Example8/index.md

@@ -22,9 +22,9 @@ model -ndm 3 -ndf 3
 {{% /tab %}}
 {{% tab name="Python (RT)" %}}
 ```python
-import opensees.openseespy as ops
+import xara
 
-model = ops.Model(ndm=3, ndf=3)
+model = xara.Model(ndm=3, ndf=3)
 ```
 {{% /tab %}}
 {{< /tabs >}}
@@ -43,7 +43,7 @@ The final argument of the `block3D` command is a block specifying the location o
 {{% tab name="Tcl" %}}
 ```tcl
 # mesh generation
-block3D $nx $ny $nz   1 1  $element  $eleArgs {
+block3D $nx $ny $nz   1 1  $element  $args {
     1   -1     -1      0
     2    1     -1      0
     3    1      1      0
@@ -57,7 +57,7 @@ block3D $nx $ny $nz   1 1  $element  $eleArgs {
 {{% /tab %}}
 {{% tab name="Python (RT)" %}}
 ```python
-model.block3D(nx, ny, nz, 1, 1, Brick, 1, {
+model.block3D((nx, ny, nz), (1, 1), "Brick", 1, {
               1: [-1.0, -1.0,  0.0],
               2: [ 1.0, -1.0,  0.0],
               3: [ 1.0,  1.0,  0.0],

content/en/examples/FiniteTension/index.md

@@ -15,23 +15,29 @@ $$
 The deformation gradient is, in coordinates:
 
 $$
-\boldsymbol{F}=\left[\begin{array}{lll}
+\boldsymbol{F}
+=\begin{bmatrix}
 \Lambda & 0 & 0 \\
 0 & 1 & 0 \\
 0 & 0 & 1
-\end{array}\right]
+\end{bmatrix}
 $$
-
+And the Cauchy-Green tensor is
 $$
-\boldsymbol{C}=\boldsymbol{F}^T \boldsymbol{F}=\left[\begin{array}{ccc}
+\boldsymbol{C}=\boldsymbol{F}^T \boldsymbol{F}
+= \Lambda^2 \mathbf{e}\otimes\mathbf{e} + \mathbf{e}_{\alpha} \otimes \mathbf{e}_{\alpha}
+=\left[\begin{array}{ccc}
 \Lambda^2 & 0 & 0 \\
 0 & 1 & 0 \\
 0 & 0 & 1
 \end{array}\right]
 $$
-furnishing the Green-Lagrange strain \(\boldsymbol{E}\)
+furnishing the Green-Lagrange strain $\boldsymbol{E}$
 $$
-\boldsymbol{E}=\frac{1}{2}(\boldsymbol{C}-\boldsymbol{I})=\left[\begin{array}{ccc}
+\boldsymbol{E}
+=\frac{1}{2}(\boldsymbol{C}-\boldsymbol{1})
+= \frac{1}{2}\left(\Lambda^2-1\right) \mathbf{e}\otimes\mathbf{e}
+=\left[\begin{array}{ccc}
 \frac{1}{2}\left(\Lambda^2-1\right) & 0 & 0 \\
 0 & 0 & 0 \\
 0 & 0 & 0
@@ -40,7 +46,7 @@ $$
 
 For a SVK material the 2PK stress is
 $$
-\boldsymbol{S}=\lambda \operatorname{tr}(\boldsymbol{E}) \boldsymbol{I}+2 \mu \boldsymbol{E}=\left[\begin{array}{ccc}
+\boldsymbol{S}=\lambda \operatorname{tr}(\boldsymbol{E}) \boldsymbol{1}+2 \mu \boldsymbol{E}=\left[\begin{array}{ccc}
 \left(\Lambda^2-1\right)\left(\frac{\lambda}{2}+\mu\right) & 0 & 0 \\
 0 & \frac{\lambda}{2}\left(\Lambda^2-1\right) & 0 \\
 0 & 0 & \frac{\lambda}{2}\left(\Lambda^2-1\right)

content/en/examples/FrameVecXZ/index.md

@@ -1,5 +1,6 @@
 ---
 title: Orienting frames in 3D
+weight: 10
 description: This example demonstrates how to orient beam elements in 3D. 
 thumbnail: img/examples/vecxz.png
 downloads:
@@ -82,4 +83,4 @@ $$
 & \frac{2}{L}    EI_g 
 & \frac{4}{L}    EI_g + \frac{4}{H} EI_2 \\
 \end{array}\right]
-$$
+$$

content/en/examples/PlaneFlexure/index.md

@@ -27,9 +27,9 @@ $$
 
 with $R = \tfrac{1}{12} Ed^3$.
 
-<figure>
-<img src="img/cross.png" width="40%" >
-<img src="img/beam.png" width="40%" >
+<figure style="text-align: center">
+<img src="img/cross.png" width="45%" >
+<img src="img/beam.png" width="45%" >
 <figcaption>Deformed cross section from 3D solution of pure flexure.</figcaption>
 </figure>
 

content/en/examples/PlaneTaper/plane_ledge.py

@@ -0,0 +1,65 @@
+
+import xara
+from xara.helpers import find_node, find_nodes
+from veux.stress import node_average
+
+def create_ledge(mesh,
+                thickness=1,
+                element: str = "LagrangeQuad"):
+
+    nx, ny = mesh
+
+    # create model in two dimensions with 2 DOFs per node
+    model = ops.Model(ndm=2, ndf=2)
+
+    # Define the material
+    # -------------------
+    #                                 tag  E      nu   rho
+    model.material("ElasticIsotropic", 1, 1000.0, 0.25, 0, "-plane-strain")
+
+    # Define geometry
+    # ---------------
+
+    load = 1000.0
+    angle = 0
+    b = 30
+    L = 100.0
+    r = 15/L
+
+    # Create the nodes and elements using the surface command
+    # {"quad", "enhancedQuad", "tri31", "LagrangeQuad"}:
+    args = (thickness, "PlaneStrain", 1)
+
+    surface = model.surface((nx, ny),
+                  element=element, args=args,
+                  points={
+                    1: [  0.0,   0.0],
+                    5: [  L/2,   10.],
+                    2: [   L,    L*r],
+                    3: [   L,     b ],
+                    4: [  0.0,    b ]
+                  })
+
+    # Single-point constraints
+    #            x   (u1 u2)
+    for node in find_nodes(model, x=0.0):
+        print("Fixing node ", node)
+        model.fix(node, (1, 1))
+
+    # Define gravity loads
+    # create a Plain load pattern with a linear time series
+    model.pattern("Plain", 1, "Linear")
+
+#   # Find the node at the tip center
+#   place = find_node(model, x=L, y=15.0)
+#   print(f"Placing load at node {place}")
+#   force = (1.0, 0.0)
+#   model.load(place, force, pattern=1)
+
+    tip = list(find_nodes(model, x=L))
+    for node in tip:
+        model.load(node, (load/len(tip), 0.0), pattern=1)
+
+    return model
+
+

content/en/examples/Wrench/index.md

@@ -19,7 +19,7 @@ For the triangular mesh, the [`Tri31`](https://xara.so/user/manual/model/element
 {{% tab name="Quadrilateral" %}}
 ```python
 def create_quads():
-    model = ops.Model(ndm=2, ndf=2)
+    model = xara.Model(ndm=2, ndf=2)
     model.nDMaterial("ElasticIsotropic", 1, 200e3, 0.25)
 
     for num,block in blocks.items():
@@ -34,7 +34,7 @@ def create_quads():
 {{% tab name="Triangle" %}}
 ```python
 def create_tris():
-    model = ops.Model(ndm=2, ndf=2)
+    model = xara.Model(ndm=2, ndf=2)
     model.nDMaterial("ElasticIsotropic", 1, 200e3, 0.25)
 
     elem = 1