Skip to content
This repository was archived by the owner on Nov 24, 2021. It is now read-only.

Commit 8949e0e

Browse files
author
jorgepz
committed
add corotational frame element diagram
1 parent fcb5c4c commit 8949e0e

8 files changed

+1263
-8
lines changed

docs/make.jl

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -36,6 +36,7 @@ elseif outputFormat == "html"
3636
"Cantilever Beam" => "cantileverBeam.md"],
3737
"User guide" => Any["Installation" => "howtouse/install.md",
3838
"Creating Models" => "howtouse/creatingModels.md",
39+
"Corotational frame" => "corotationalFrameElement.md",
3940
"References" => "theory/references.md"]
4041
#"Developer guide" => Any["ONSAS_solve" => "ONSAS_solve.md"]
4142
],

docs/src/ONSAS_solve.md

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -27,7 +27,7 @@ while finalTimeReachedBoolean == false
2727
2828
% check if final time was reached
2929
finalTimeReachedBoolean = ( modelNextSol.currTime - modelProperties.analysisSettings.finalTime ) ...
30-
>= ( -(modelProperties.analysisSettings.finalTime) * 1e-8 ) ;
30+
>= ( -(modelProperties.analysisSettings.finalTime) * 1e-8 )
3131
3232
% store results and update structs
3333
modelCurrSol = modelNextSol ;
@@ -36,7 +36,7 @@ while finalTimeReachedBoolean == false
3636
3737
% generate vtk file for the new state
3838
if strcmp( modelProperties.plotsFormat, 'vtk' )
39-
vtkMainWriter( modelCurrSol, modelProperties )
39+
vtkMainWriter( modelCurrSol, modelProperties );
4040
end % if vtk output format
4141
4242
end %while time

docs/src/cantileverBeam.md

Lines changed: 4 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -103,7 +103,7 @@ end
103103
```
104104
analysisSettings.methodName = 'newtonRaphson' ;
105105
analysisSettings.deltaT = 0.1 ;
106-
analysisSettings.finalTime = 1.0 ;
106+
analysisSettings.finalTime = 0.2 ;
107107
analysisSettings.stopTolDeltau = 1e-6 ;
108108
analysisSettings.stopTolForces = 1e-6 ;
109109
analysisSettings.stopTolIts = 10 ;
@@ -113,13 +113,14 @@ analysisSettings.stopTolIts = 10 ;
113113
```
114114
otherParams.problemName = 'uniformCurvatureCantilever';
115115
otherParams.controlDofs = [ numElements+1 4 ] ;
116+
otherParams.plotsFormat = 'vtk' ;
116117
```
117118
## Analysis case 1: NR with Rotated Eng Strain
118119
In the first case ONSAS is run and the solution at the dof (angle of node B) of interest is stored:
119120
```
120121
[matUs, loadFactorsMat] = ONSAS( materials, elements, boundaryConds, initialConds, mesh, analysisSettings, otherParams ) ;
121-
122122
```
123+
123124
the control dof to verificate the solution is the node angle B, this corresponds to the following dof number:
124125
```
125126
angleControlDof = (numElements+1)*6 - 2;
@@ -136,7 +137,7 @@ analyticLoadFactorsNREngRot = @(w) E * Iy * w / l ;
136137
```
137138
verifBoolean = norm( analyticLoadFactorsNREngRot( controlDispsNREngRot) ...
138139
- loadFactorsNREngRot' ) ...
139-
< ( norm( analyticLoadFactorsNREngRot( controlDispsNREngRot) ) * 1e-4 )
140+
< ( norm( analyticLoadFactorsNREngRot( controlDispsNREngRot) ) * 1e-4 ) ;
140141
```
141142

142143

docs/src/corotationalFrameElement.md

Lines changed: 5 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,5 @@
1+
# Co-rotational frame system
2+
3+
```@raw html
4+
<img src="diagramCorotationalFrameSystem.svg" alt="corotational diagram" width="800"/>
5+
```

0 commit comments

Comments
 (0)