Mercurial > repos > public > sbplib

changeset 339:3135e13ae509 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Implemented free BC. Not tested yet.
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 24 Oct 2016 17:01:34 +0200
parents 3974dccff55b
children 25ae8333ce2b
files +scheme/Beam.m
diffstat 1 files changed, 14 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
diff -r 3974dccff55b -r 3135e13ae509 +scheme/Beam.m
--- a/+scheme/Beam.m	Mon Oct 24 17:01:05 2016 +0200
+++ b/+scheme/Beam.m	Mon Oct 24 17:01:34 2016 +0200
@@ -91,7 +91,7 @@
             delt = obj.delt;
 
             switch type
-                case {'dn'} % Dirichlet-neumann boundary condition
+                case {'dn', 'clamped'} % Dirichlet-neumann boundary condition
                     alpha = obj.alpha;
 
                     % tau1 < -alpha^2/gamma
@@ -111,6 +111,19 @@
 
                     penalty{1} = -obj.Hi*tau;
                     penalty{2} = -obj.Hi*sig;
+
+
+                case {'free'}
+                    a = -obj.alpha;
+
+                    tau =  s*a*d1;
+                    sig = -s*a*e;
+
+                    closure = obj.Hi*(tau*d2' + sig*d3');
+                    penalty{1} = -obj.Hi*tau;
+                    penalty{1} = -obj.Hi*sig;
+
+
                 otherwise % Unknown, boundary condition
                     error('No such boundary condition: type = %s',type);
             end